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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗
DOE OSTI2024

A conservative discontinuous-Galerkin-in-time (DGiT) multirate time integration framework for interface-coupled problems with applications to solid–solid interaction and air–sea models

In this paper we extend the DGiT multirate framework, developed in Connors and Sockwell (2022) for scalar transmission problems, to a solid–solid interaction (SSI) problem involving two coupled elastic solids and a coupled air–sea model with the rotating, thermal shallow water equations. In so doing we aim to demonstrate the broad applicability of the mathematical theory and governing principles established in Connors and Sockwell (2022) to coupled problems characterized by subproblems evolving at different temporal scales. Further, multirate time integration algorithms employing different time steps, optimized for the dynamics of each subproblem, can significantly improve simulation efficiency for such coupled problems. However, development of multirate algorithms is a highly non-trivial task due to the coupling, which can impact accuracy, stability or other desired properties such as preservation of system invariants. DGiT provides a general template for multirate time integration that can achieve these properties. To elucidate the manner in which DGiT accomplishes this task, we fully detail each step in the application of the framework to the SSI and air–sea coupled problems. Numerical examples illustrate key properties of the resulting multirate schemes for both problems.

42 ENGINEERING↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Soil carbon change in intensive agriculture after 25 years of conservation management

Changes in soil organic carbon (SOC) and nitrogen (SON) are strongly affected by land management but few long-term comparative studies have surveyed changes throughout the whole soil profile. We quantified 25-year SOC and SON changes to 1 m in 10 replicate ecosystems at an Upper Midwest, USA site. We compared four annual cropping systems in maize (Zea mays)-soybean (Glycine max)-winter wheat (Triticum aestivum) rotations, each managed differently (Conventional, No-till, Reduced input, and Biologically based); in three managed perennial systems (hybrid Poplar (Populus × euramericana), Alfalfa (Medicago sativa), and Conifer (Pinus spp.); and in three successional systems (Early, Mid- and Late succession undergoing a gradual change in species composition and structure over time). Both Reduced input and Biologically based systems included winter cover crops. Neither SOC nor SON changed significantly in the Conventional or Late successional systems over 25 years. All other systems gained SOC and SON to different degrees. SOC accrual was fastest in the Early successional system (0.8 ± 0.1 Mg C ha –1 y –1 ) followed by Alfalfa and Conifer (avg. 0.7 ± 0.1 Mg C ha –1 y –1 ), Poplar, Reduced input, and Biologically based systems (avg. 0.4 ± 0.1 Mg C ha –1 y –1 ), and Mid-successional and No-till systems (0.3 and 0.2 Mg C ha –1 y –1 , respectively). Over the most recent 12 years, rates of SOC accrual slowed in all systems except Reduced input and Mid-successional. There was no evidence of SOC loss at depth in any system, including No-till. Rates of SON accrual ranged from 64.7 to 0.8 kg N ha –1 y –1 in the order Alfalfa ≥ Early successional > Reduced input and Biologically based ≥ Poplar > No-till and Conifer > Mid-successional systems. Pyrogenic C levels in the Conventional, Early, and Late successional systems were similar despite 17 years of annual burning in the Early successional system (~ 15 % of SOC to 50 cm, on average, and ~40 % of SOC from 50 to 100 cm). Results underscore the importance of cover crops, perennial crops, and no-till options for sequestering whole profile C in intensively managed croplands.

60 APPLIED LIFE SCIENCES↗