Search NASA⌕ Search

SEARCH · Search NASA

Results for “PRESERVATION”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
DOE OSTI2020

On differentiable local bounds preserving stabilization for Euler equations

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

42 ENGINEERING↗
DOE OSTI2015

A freestream-preserving fourth-order finite-volume method in mapped coordinates with adaptive-mesh refinement

We present a fourth-order accurate finite-volume method for solving time-dependent hyperbolic systems of conservation laws on mapped grids that are adaptively refined in space and time. Some novel considerations for formulating the semi-discrete system of equations in computational space are combined with detailed mechanisms for accommodating the adapting grids. Furthermore, these considerations ensure that conservation is maintained and that the divergence of a constant vector field is always zero (freestream-preservation property). The solution in time is advanced with a fourth-order Runge-Kutta method. A series of tests verifies that the expected accuracy is achieved in smooth flows and the solution of a Mach reflection problem demonstrates the effectiveness of the algorithm in resolving strong discontinuities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Selective preservation of pyrogenic carbon across soil organic matter fractions and its influence on calculations of carbon mean residence times

The long-standing perspective that recalcitrance of soil organic carbon (SOC) controls its stability and persistence has shifted to one in which physical inaccessibility of SOC to microorganisms plays a predominant role. This paradigm shift has been facilitated by analytical techniques that isolate SOC into physical fractions protected from decomposers by different mechanisms. The correlation between these fractions and SOC age has reinforced the emphasis of SOC inaccessibility. Pyrogenic C (PyC; also called charcoal or black carbon), which has been thermally altered by fire, is known to contain highly recalcitrant components that decompose very slowly and could represent an exception to this paradigm shift. We employed hydrogen pyrolysis to quantify the contribution of PyC to total SOC across soil fractions from three long-term agricultural experiments with land use conversions that caused reductions in SOC. We show that all soil fractions contain PyC and up to one-fifth of SOC in soil fractions considered to have low accessibility is comprised of PyC. Regardless of the soil fraction in which it was located, PyC was relatively unaffected by land use conversion compared to biogenic C (organic C not altered by fire), which suggests that selective preservation, rather than physical protection, is the dominant mechanism limiting PyC decomposition in these sites. We accounted for PyC in calculations of C mean residence times (MRTs) using differences in stable C isotope ratios between PyC and SOC. Though results varied by site and soil fraction, MRTs for biogenic C were generally shorter than for total SOC. Based on these results, PyC decomposition is controlled by a different mechanism than biogenic C, and this should be considered in studies of soil C dynamics. In addition, methods based on physical fractionation alone may place too great an emphasis on the role of inaccessibility for long-term SOC persistence.

58 GEOSCIENCES↗
DOE OSTI2015

Total Laryngectomy Versus Larynx Preservation for T4a Larynx Cancer: Patterns of Care and Survival Outcomes

Purpose: To examine practice patterns and compare survival outcomes between total laryngectomy (TL) and larynx preservation chemoradiation (LP-CRT) in the setting of T4a larynx cancer, using a large national cancer registry. Methods and Materials: Using the National Cancer Database, we identified 969 patients from 2003 to 2006 with T4a squamous cell larynx cancer receiving definitive treatment with either initial TL plus adjuvant therapy or LP-CRT. Univariate and multivariable logistic regression were used to assess predictors of undergoing surgery. Survival outcomes were compared using Kaplan-Meier and propensity score–adjusted and inverse probability of treatment–weighted Cox proportional hazards methods. Sensitivity analyses were performed to account for unmeasured confounders. Results: A total of 616 patients (64%) received LP-CRT, and 353 (36%) received TL. On multivariable logistic regression, patients with advanced nodal disease were less likely to receive TL (N2 vs N0, 26.6% vs 43.4%, odds ratio [OR] 0.52, 95% confidence interval [CI] 0.37-0.73; N3 vs N0, 19.1% vs 43.4%, OR 0.23, 95% CI 0.07-0.77), whereas patients treated in high case-volume facilities were more likely to receive TL (46.1% vs 31.5%, OR 1.78, 95% CI 1.27-2.48). Median survival for TL versus LP was 61 versus 39 months (P<.001). After controlling for potential confounders, LP-CRT had inferior overall survival compared with TL (hazard ratio 1.31, 95% CI 1.10-1.57), and with the inverse probability of treatment–weighted model (hazard ratio 1.25, 95% CI 1.05-1.49). This survival difference was shown to be robust on additional sensitivity analyses. Conclusions: Most patients with T4a larynx cancer receive LP-CRT, despite guidelines suggesting TL as the preferred initial approach. Patients receiving LP-CRT had more advanced nodal disease and worse overall survival. Previous studies of (non-T4a) locally advanced larynx cancer showing no difference in survival between LP-CRT and TL may not apply to T4a disease, and patients should be counseled accordingly.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Early Radiosurgery Improves Hearing Preservation in Vestibular Schwannoma Patients With Normal Hearing at the Time of Diagnosis

Purpose: Vestibular schwannomas (VS) are increasingly diagnosed in patients with normal hearing because of advances in magnetic resonance imaging. We sought to evaluate whether stereotactic radiosurgery (SRS) performed earlier after diagnosis improved long-term hearing preservation in this population. Methods and Materials: We queried our quality assessment registry and found the records of 1134 acoustic neuroma patients who underwent SRS during a 15-year period (1997-2011). We identified 88 patients who had VS but normal hearing with no subjective hearing loss at the time of diagnosis. All patients were Gardner-Robertson (GR) class I at the time of SRS. Fifty-seven patients underwent early (≤2 years from diagnosis) SRS and 31 patients underwent late (>2 years after diagnosis) SRS. At a median follow-up time of 75 months, we evaluated patient outcomes. Results: Tumor control rates (decreased or stable in size) were similar in the early (95%) and late (90%) treatment groups (P=.73). Patients in the early treatment group retained serviceable (GR class I/II) hearing and normal (GR class I) hearing longer than did patients in the late treatment group (serviceable hearing, P=.006; normal hearing, P<.0001, respectively). At 5 years after SRS, an estimated 88% of the early treatment group retained serviceable hearing and 77% retained normal hearing, compared with 55% with serviceable hearing and 33% with normal hearing in the late treatment group. Conclusions: SRS within 2 years after diagnosis of VS in normal hearing patients resulted in improved retention of all hearing measures compared with later SRS.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

A diffuse interface model with immiscibility preservation

A new, simple, and computationally efficient interface capturing scheme based on a diffuse interface approach is presented for simulation of compressible multiphase flows. Multi-fluid interfaces are represented using field variables (interface functions) with associated transport equations that are augmented, with respect to an established formulation, to enforce a selected interface thickness. The resulting interface region can be set just thick enough to be resolved by the underlying mesh and numerical method, yet thin enough to provide an efficient model for dynamics of well-resolved scales. A key advance in the present method is that the interface regularization is asymptotically compatible with the thermodynamic mixture laws of the mixture model upon which it is constructed. It incorporates first-order pressure and velocity non-equilibrium effects while preserving interface conditions for equilibrium flows, even within the thin diffused mixture region. We first quantify the improved convergence of this formulation in some widely used one-dimensional configurations, then show that it enables fundamentally better simulations of bubble dynamics. Demonstrations include both a spherical-bubble collapse, which is shown to maintain excellent symmetry despite the Cartesian mesh, and a jetting bubble collapse adjacent a wall. Comparisons show that without the new formulation the jet is suppressed by numerical diffusion leading to qualitatively incorrect results.

42 ENGINEERING↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Here we present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Fast nonlinear iterative solver for an implicit, energy-conserving, asymptotic-preserving charged-particle orbit integrator

Here recently, an asymptotic-preserving (AP) particle orbit integrator has been proposed with remarkable properties including exact energy conservation, the ability to capture of all first-order drifts (including the ∇B-drift), the ability to capture trapped-passing boundaries with parallel velocity extremely close to the critical velocity, and the ability to transition from strongly to weakly magnetized spatial regions. The new AP orbit integrator is implicit, employing a Crank-Nicolson (CN) temporal discretization to ensure exact energy conservation. This, in turn, requires a local nonlinear iteration involving particle velocities and positions, and the local electromagnetic fields, to obtain the new-time solution. Ref. [1] did not attempt to provide an efficient solver for this system, and employed a brute-force GMRES-driven Jacobian-free Newton-Krylov (JFNK) solver to invert the particle orbit equations at every timestep for expediency. While JFNK is robust and reliable, it is also expensive and very intrusive for practical implementations of the method (it requires having the JFNK machinery available and solving a 6 x 6 Jacobian system iteratively once per iteration per particle).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A method for bounding high-order finite element functions: Applications to mesh validity and bounds-preserving limiters

We introduce a novel method for bounding high-order multi-dimensional polynomials in finite element approximations. The method involves precomputing optimal piecewise-linear bounding boxes for polynomial basis functions, which can then be used to locally bound any combination of these basis functions. This approach can be applied to any element/basis type at any approximation order, can provide local (i.e., subcell) extremum bounds to a desired level of accuracy, and can be evaluated efficiently on-the-fly in simulations. Furthermore, we show that this approach generally yields more accurate bounds in comparison to traditional methods based on convex hull properties (e.g., Bernstein polynomials). Furthermore, the efficacy of this technique is shown in applications such as mesh validity checks and optimization for high-order curved meshes, where positivity of the element Jacobian determinant can be ensured throughout the entire element, and continuously bounds-preserving limiters for hyperbolic systems, which can enforce maximum principle bounds across the entire solution polynomial.

Bounding box↗
DOE OSTI2017

Data preservation at the Fermilab Tevatron

The Fermilab Tevatron collider's data-taking run ended in September 2011, yielding a dataset with rich scientific potential. The CDF and D0 experiments each have approximately 9 PB of collider and simulated data stored on tape. A large computing infrastructure consisting of tape storage, disk cache, and distributed grid computing for physics analysis with the Tevatron data is present at Fermilab. The Fermilab Run II data preservation project intends to keep this analysis capability sustained through the year 2020 and beyond. To achieve this goal, we have implemented a system that utilizes virtualization, automated validation, and migration to new standards in both software and data storage technology and leverages resources available from currently-running experiments at Fermilab. Lastly, these efforts have also provided useful lessons in ensuring long-term data access for numerous experiments, and enable high-quality scientific output for years to come.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Tip radius preservation for high resolution imaging in amplitude modulation atomic force microscopy

The acquisition of high resolution images in atomic force microscopy (AFM) is correlated to the cantilever's tip shape, size, and imaging conditions. In this work, relative tip wear is quantified based on the evolution of a direct experimental observable in amplitude modulation atomic force microscopy, i.e., the critical amplitude. We further show that the scanning parameters required to guarantee a maximum compressive stress that is lower than the yield/fracture stress of the tip can be estimated via experimental observables. In both counts, the optimized parameters to acquire AFM images while preserving the tip are discussed. The results are validated experimentally by employing IgG antibodies as a model system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Data preservation at the Fermilab Tevatron

The Fermilab Tevatron collider's data-taking run ended in September 2011, yielding a dataset with rich scientific potential. The CDF and DO experiments each have nearly 9 PB of collider and simulated data stored on tape. A large computing infrastructure consisting of tape storage, disk cache, and distributed grid computing for physics analysis with the Tevatron data is present at Fermilab. The Fermilab Run II data preservation project intends to keep this analysis capability sustained through the year 2020 or beyond. To achieve this, we are implementing a system that utilizes virtualization, automated validation, and migration to new standards in both software and data storage technology as well as leveraging resources available from currently-running experiments at Fermilab. Furthermore, these efforts will provide useful lessons in ensuring long-term data access for numerous experiments throughout high-energy physics, and provide a roadmap for high-quality scientific output for years to come.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2014

Production and Early Preservation of Lipid Biomarkers in Iron Hot Springs

The bicarbonate-buffered anoxic vent waters at Chocolate Pots hot springs in Yellowstone National Park are 51–54°C, pH 5.5–6.0, and are very high in dissolved Fe(II) at 5.8–5.9 mg/L. The aqueous Fe(II) is oxidized by a combination of biotic and abiotic mechanisms and precipitated as primary siliceous nanophase iron oxyhydroxides (ferrihydrite). Four distinct prokaryotic photosynthetic microbial mat types grow on top of these iron deposits. Lipids were used to characterize the community composition of the microbial mats, link source organisms to geologically significant biomarkers, and investigate how iron mineralization degrades the lipid signature of the community. The phospholipid and glycolipid fatty acid profiles of the highest-temperature mats indicate that they are dominated by cyanobacteria and green nonsulfur filamentous anoxygenic phototrophs (FAPs). Diagnostic lipid biomarkers of the cyanobacteria include midchain branched mono- and dimethylalkanes and, most notably, 2-methylbacteriohopanepolyol. Diagnostic lipid biomarkers of the FAPs (Chloroflexus and Roseiflexus spp.) include wax esters and a long-chain tri-unsaturated alkene. Surprisingly, the lipid biomarkers resisted the earliest stages of microbial degradation and diagenesis to survive in the iron oxides beneath the mats. Understanding the potential of particular sedimentary environments to capture and preserve fossil biosignatures is of vital importance in the selection of the best landing sites for future astrobiological missions to Mars. Finally, this study explores the nature of organic degradation processes in moderately thermal Fe(II)-rich groundwater springs—environmental conditions that have been previously identified as highly relevant for Mars exploration.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2019

Beam emittance preservation using Gaussian density ramps in a beam-driven plasma wakefield accelerator

A current challenge that is facing the plasma wakefield accelerator (PWFA) community is transverse beam emittance preservation. This can be achieved by balancing the natural divergence of the beam against the strong focusing force provided by the PWFA plasma source in a scheme referred to as beam matching. One method to accomplish beam matching is through the gradual focusing of a beam with a plasma density ramp leading into the bulk plasma. Here, the beam dynamics in a Gaussian plasma density ramp are considered, and an empirical formula is identified that gives the ramp length and beam vacuum waist location needed to achieve near-perfect matching. The method uses only the beam vacuum waist beta function as an input. Numerical studies show that the Gaussian ramp focusing formula is robust for beta function demagnification factors spanning more than an order of magnitude with experimentally favourable tolerances for future PWFA research facilities. This article is part of the Theo Murphy meeting issue ‘Directions in particle beam-driven plasma wakefield acceleration’.

43 PARTICLE ACCELERATORS↗