Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical differentiation”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Differential Organic Carbon Mineralization Responses to Soil Moisture in Three Different Soil Orders Under Mixed Forested System: Supporting Data

This data contains data from 90-day long incubation study which aimed to look at the soil moisture-texture relationship on soil organic carbon (SOC) cycling. Soils were collected from three distinct soil textures from mixed forests in 2017: sandy (Georgia, 2017-05-01), loamy (Missouri, 2017-06-14) and clayey (Texas, December 2017) were incubated at different soil moisture levels (air-dried, 25% water holding capacity (WHC), 50% WHC, 100% WHC and 175% WHC) at room temperature for a period of 90 days. Files contain microbial respiration, active and slow SOC pools, and their respective mineralization rates, extractable organic carbon (C), and C-acquiring extracellular enzymes. Findings from these data were used in Singh et al. (2021). This study aimed to examine the interactive effect of soil moisture and texture on SOC mineralization. Soil samples of three distinct textures (sandy, loamy, and clayey) were collected from mixed forests of Georgia, Missouri, and Texas, respectively. Soil cores of 5 cm diameter were collected from numerous random locations at each site from 0-15 cm depth after scraping the litter layer and mixed thoroughly to obtain a composite sample per site. Three additional soil cores were collected to determine the WHC using pressure plate extractors. Soil samples were composited, and triplicate soil samples were incubated in mason jars for a period of 90 days at room temperature under different moisture regimes: air dried, 25% WHC, 50% WHC, at WHC and 100% saturation. Soil respiration was measured weekly, and destructive sampling was conducted at 1, 15, 60, and 90 days to determine extractable organic C, C acquiring enzyme activity, and active and slow SOC pools with their respective mineralization rates. The C acquiring enzyme activity was the total activity of α-glucosidase, β-glucosidase, cellobiohydrolase, and β-xylosidase enzymes. Gas samples for microbial respiration measurements were collected from headspace of incubation jars through the sampling ports on the lids and then analyzed using a Shimadzu Gas Chromatograph (GC-2014). Prior to sampling, the vials were evacuated. Blank correction was also done by collecting gas samples from empty incubation jars. Double pool exponential decay model was used in SigmaPlot to determine the active and slow SOC pools and their mineralization rates (Farrar et al., 2012; Jagadamma et al., 2014). The C-acquiring extracellular enzymes were measured using the microplate method by German et al., (2011). Microbial community structure was determined using the phospholipid fatty acid (PLFA) and neutral lipid fatty acid (NLFA) analyses (Buyer and Sasser, 2012). This dataset has seven data files provided in comma-separate (*.csv) format. Additional metadata are provided: seven data dictionaries and a file-level metadata file in comma separate (*.csv) format and a user guide in PDF (*.pdf) format.

Carbon acquiring enzyme activity↗
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

VALVE TIMING AND AREA PROFILE SELECTION FOR HYDRAULIC PUMPS AND MOTORS

A significant amount of research has been conducted to select valve timing and area profiles that create efficient and quiet hydraulic pumps and motors. Numerous active valve architectures have been modeled and optimized, but the rationale for the final solution is often unclear. The solution is usually highly dependent on the modeled valve geometry constraints and the duty cycle of the pump or motor for which the valve was optimized. This paper first presents an ideal valve area profile and what compromises are made to design for a working pressure and speed range. The ideal valve area profile is resolved using a piston-cylinder model with a pressure dependent bulk modulus and variable valve openings. The piston-cylinder is first modeled as a pump with passive (check) valves and a specified pressure drop across the valves to determine the ideal area profile. A generalization of the ideal timing is demonstrated for positive pressure differential and positive rotation direction, also known as the first quadrant. The effects on valve timing when varying pressure are shown for the first quadrant, but the trends can be extrapolated for all quadrants of operation. The ideal valve profile is implemented as fixed valve timing, meaning the valve area is only a function of the timing angle, in a pump-motor. The pump-motor is simulated in one rotation direction through a pressure range. Performance is high in pumping operation, but when the pressure differential is reversed, cylinder pressure spikes ensue. Two strategies to modify an ideal valve area profile are presented: timing grooves and a pressure shifted valve timing. Timing grooves reduce pressure spikes and cavitation in the cylinder but generally increase throttling losses. A pressure shifted valve timing has lower throttling energy losses, making it the favored solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics

We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

ECLEIRS: Exact conservation law embedded identification of reduced states for parameterized nonlinear conservation laws from sparse and noisy data

Multi-query applications such as parameter estimation, uncertainty quantification and design optimization for parameterized partial differential equation (PDE) systems are expensive. While reduced/latent state dynamics approaches for parameterized PDEs offer a viable alternative, these approaches rely on high-quality data and struggle with highly sparse spatiotemporal noisy measurements typically obtained from experiments. Furthermore, there is no guarantee that these models satisfy governing physical conservation laws. In this article, we propose a reduced state dynamics approach, referred to as ECLEIRS, that embeds exact conservation in the solution and flux representation by utilizing a space-time divergence-free neural network formulation. We compare ECLEIRS with other reduced state dynamics approaches, those that do not enforce any physical constraints and those with physics-informed loss functions, for three shock-propagation problems: 1-D advection, 1-D Burgers and 2-D Euler equations. In conclusion, the numerical experiments conducted in this study demonstrate that ECLEIRS provides the most accurate prediction of dynamics for unseen parameters even in the presence of highly sparse and noisy data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

On fast simulation of dynamical system with neural vector enhanced numerical solver

The large-scale simulation of dynamical systems is critical in numerous scientific and engineering disciplines. However, traditional numerical solvers are limited by the choice of step sizes when estimating integration, resulting in a trade-off between accuracy and computational efficiency. To address this challenge, we introduce a deep learning-based corrector called Neural Vector (NeurVec), which can compensate for integration errors and enable larger time step sizes in simulations. Our extensive experiments on a variety of complex dynamical system benchmarks demonstrate that NeurVec exhibits remarkable generalization capability on a continuous phase space, even when trained using limited and discrete data. NeurVec significantly accelerates traditional solvers, achieving speeds tens to hundreds of times faster while maintaining high levels of accuracy and stability. Moreover, NeurVec’s simple-yet-effective design, combined with its ease of implementation, has the potential to establish a new paradigm for fast-solving differential equations based on deep learning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Can lidars assess wind plant blockage in simple terrain? A WRF-LES study

Wind plant blockage reduces wind velocity upstream of wind plants, reducing the power generated by turbines adjacent to the inflow, and potentially throughout the plant as well. The nature of the mechanism that amplifies blockage as well as the velocity reductions in both the induction zone and potentially deeper into the array are not well understood. Field observations can provide valuable insight into the characteristics of the induction zone and the mechanisms that amplify it. However, the relatively small velocity reductions that have been measured experimentally pose a challenge in quantifying blockage, especially in onshore environments with flow heterogeneities that may be of the same scale as the blockage effect itself. We simulate the flow around the King Plains wind plant in the relatively simple terrain of Oklahoma, the location of the American WAKE experimeNt, to evaluate wind plant blockage in this environment. Using numerical simulations, we find the largest velocity deceleration (0.64 m s−1; 8%) immediately upstream of the wind plant, and 1% velocity deficits 24 rotor diameters upstream of the first turbine row. We also use virtual measurements upstream of the wind plant to analyze the uncertainties and difficulties in measuring blockage using a scanning lidar on shore. Based on our virtual lidar study, the induction zone of land-based wind plants can be incorrectly estimated using observations if the effects of nonuniform terrain on the flow are not carefully considered. Changes in terrain elevation produce local variations in wind speed (as measured by a scanning lidar) that exceed in magnitude the deceleration within the induction zone. We refer to these local changes in wind speed as terrain effects. A methodology to differentiate between terrain effects and blockage in experimental settings is proposed and evaluated herein, highlighting the difficulties and uncertainties associated with measurement and simulation of blockage in even relatively simple onshore environments.

17 WIND ENERGY↗
DOE OSTI2021

Multilevel Hierarchical Decomposition of Finite Element White Noise with Application to Multilevel Markov Chain Monte Carlo

In this work we develop a new hierarchical multilevel approach to generate Gaussian random field realizations in an algorithmically scalable manner that is well suited to incorporating into multilevel Markov chain Monte Carlo (MCMC) algorithms. This approach builds off of other partial differential equation (PDE) approaches for generating Gaussian random field realizations; in particular, a single field realization may be formed by solving a reaction-diffusion PDE with a spatial white noise source function as the right-hand side. While these approaches have been explored to accelerate forward uncertainty quantification tasks, e.g., multilevel Monte Carlo, the previous constructions are not directly applicable to multilevel MCMC frameworks which build fine-scale random fields in a hierarchical fashion from coarse-scale random fields. Our new hierarchical multilevel method relies on a hierarchical decomposition of the white noise source function in $L^2$ which allows us to form Gaussian random field realizations across multiple levels of discretization in a way that fits into multilevel MCMC algorithmic frameworks. After presenting our main theoretical results and numerical scaling results to showcase the utility of this new hierarchical PDE method for generating Gaussian random field realizations, this method is tested on a four-level MCMC algorithm to explore its feasibility.

algebraic multigrid↗
DOE OSTI2021

Implicit Multirate GARK Methods

This article considers multirate generalized-structure additively partitioned Runge–Kutta methods for solving stiff systems of ordinary differential equations with multiple time scales. These methods treat different partitions of the system with different timesteps for a more targeted and efficient solution compared to monolithic single rate approaches. With implicit methods used across all partitions, methods must find a balance between stability and the cost of solving nonlinear equations for the stages. In order to characterize this important trade-off, we explore multirate coupling strategies, problems for assessing linear stability, and techniques to efficiently implement Newton iterations for stage equations. Unlike much of the existing multirate stability analysis which is limited in scope to particular methods, we present general statements on stability and describe fundamental limitations for certain types of multirate schemes. New implicit multirate methods up to fourth order are derived, and their accuracy and efficiency properties are verified with numerical tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Small-𝑥 Factorization in the Target Fragmentation Region

We consider the differential cross-section for single-inclusive jet production with transverse momentum 𝑃 ⊥ in deep inelastic scattering at small Bjorken 𝑥 Bj , mediated by a virtual photon with virtuality 𝑄 2 . Unlike most studies at small 𝑥, which focus on particle production in the current fragmentation region, we investigate the kinematic regime where the jet is produced in the target fragmentation hemisphere of the Breit frame, and with 𝑃 ⊥ ≪ 𝑄. For a longitudinally polarized photon, we demonstrate that this cross section is not power suppressed in 𝑃 ⊥ /𝑄 and we derive a factorized expression in terms of extended quark and gluon jet fracture functions. Our formula, valid at next-to-leading order in 𝛼 𝑠 at small 𝑥, is akin to the Altarelli-Martinelli identity for the longitudinal DIS structure function. Numerical estimates show that the extended quark jet fracture function is the most sensitive to saturation effects in large nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The generation of a multiphase medium in ‘Splash’ bridge systems: towards an understanding of star formation suppression in turbulent galaxy systems

ABSTRACT Cloud–cloud collisions in splash bridges produced in gas-rich disc galaxy collisions offer a brief but interesting environment to study the effects of shocks and turbulence on star formation rates in the diffuse intergalactic medium, far from the significant feedback effects of massive star formation and active galactic nucleus. Expanding on our earlier work, we describe simulated collisions between counter-rotating disc galaxies of relatively similar mass, focusing on the thermal and kinematic effects of relative inclination and disc offset at the closest approach. This includes essential heating and cooling signatures, which go some way towards explaining the luminous power in H$_2$ and [C ii] emission in the Taffy bridge, as well as providing a partial explanation of the turbulent nature of the recently observed compact CO-emitting clouds observed in Taffy by the Atacama Large Millimeter Array (ALMA). The models show counter-rotating disc collisions result in swirling, shearing kinematics for the gas in much of the post-collision bridge. Gas with little specific angular momentum due to collisions between counter-rotating streams accumulates near the centre of mass. The disturbances and mixing in the bridge drive continuing cloud collisions, differential shock heating, and cooling throughout. A wide range of relative gas phases and line-of-sight velocity distributions are found in the bridges, depending sensitively on initial disc orientations, and the resulting variety of cloud collision histories. Most cloud collisions can occur promptly or persist for quite a long duration. Cold and hot phases can largely overlap throughout the bridge or can be separated into different parts of the bridge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Long–short-term memory encoder–decoder with regularized hidden dynamics for fault detection in industrial processes

The ability of recurrent neural networks (RNN) to model nonlinear dynamics of high dimensional process data has enabled data-driven RNN-based fault detection algorithms. Previous studies have focused on detecting faults by identifying the discrepancies in data distribution between the faulty and normal data, as reflected in prediction errors generated by RNN models. However, in industrial processes, variations in data distribution can also result from changes in normal control setpoints and compensatory control adjustments in response to disturbances, making it hard to differentiate between normal and faulty conditions. This paper proposes a fault detection method utilizing a long short-term memory (LSTM) encoder–decoder structure with regularized hidden dynamics and reversible instance normalization (RevIN) to compactly represent high-dimensional measurements for effective monitoring. During training, the hidden states of the model are regularized to form a low-dimensional latent space representation of the original multivariate time series data. As a result, the prediction errors of the latent states can be used to monitor the abnormal dynamic variations, while the reconstruction errors of the measured variables are used to monitor the abnormal static variations. Furthermore, the proposed indices can reflect operating conditions, even when the distribution of test data changes, which helps distinguish faults from normal adjustments and disturbances that controllers can settle. Here, data from numerical simulation and the Tennessee Eastman process are used to illustrate the effectiveness of the proposed fault detection method.

42 ENGINEERING↗
DOE OSTI2022

Calabi-Yau CFTs and random matrices

Using numerical methods for finding Ricci-flat metrics, we explore the spectrum of local operators in two-dimensional conformal field theories defined by sigma models on Calabi-Yau targets at large volume. Focusing on the examples of K3 and the quintic, we show that the spectrum, averaged over a region in complex structure moduli space, possesses the same statistical properties as the Gaussian orthogonal ensemble of random matrix theory.

superstring vacua↗
DOE OSTI2025

Discovery of Proteoforms Associated With Alzheimer's Disease Through Quantitative Top-Down Proteomics

The complex nature of Alzheimer's disease (AD) and its heterogenous clinical presentation has prompted numerous large-scale - omic analyses aimed at providing a global understanding of the pathophysiological processes involved. AD involves isoforms, proteolytic products, and posttranslationally modified proteins such as amyloid beta (Aβ) and microtubule-associated protein tau. Top-down proteomics directly measures these species and thus, offers a comprehensive view of pathologically relevant proteoforms that are difficult to analyze using traditional proteomic techniques. Here, we broadly explored associations between proteoforms and clinicopathological traits of AD by deploying a quantitative top-down proteomics approach across frontal cortex of 103 subjects selected from the ROS and MAP cohorts. The approach identified 1213 proteins and 11,782 proteoforms, of which 154 proteoforms had at least one significant association with a clinicopathological phenotype. One important finding included identifying Aβ C-terminal truncation state as the key property for differential association between amyloid plaques and cerebral amyloid angiopathy. Furthermore, various N-terminally truncated forms of Aβ had noticeably stronger association with amyloid plaques and global cognitive function. Additionally, we discovered six VGF neuropeptides that were positively associated with cognitive function independent of pathological burden. The database of brain cortex proteoforms provides a valuable context for functional characterization of the proteins involved in AD and other late-onset brain pathologies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

On the effects of quadrupolar relaxation in Earth’s field NMR spectra

There is growing interest in using low-field magnetic resonance experiments for routine chemical characterization. Earth’s field NMR is one such technique that can garner structural information and enable sample differentiation with low cost and highly portable designs. The resulting NMR spectra are primarily influenced by $J$-couplings, resulting in so-called $J$-coupled spectra (JCS). Many small molecules include atoms with NMR-active nuclei that are quadrupolar either at natural abundance or are often isotopically enriched (e.g., 2 H, 6 Li, 11 B, 14 N, 17 O, etc.) where the effects of quadrupolar J-couplings and relaxation on JCS of strongly- and weakly-coupled spin systems have not been explored to date. Herein, using a set of seven fluoropyridine samples with unique substitution and $J$-couplings, we demonstrate that the 14 N relaxation rates can induce drastic line-broadening in the JCS. This includes a previously unexplored unique line broadening mechanism enabled by strongly coupled spins at low-field. In conclusion, numerical simulations are used to model and refine the magnitudes and signs of $J$-couplings, as well as indirectly determine the 14 N relaxation rates in a single 1D experiment that has a higher fidelity than observed in high-field NMR experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗