Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 775 records · Page 43

Exponential Decay in the Sensitivity Analysis of Nonlinear Dynamic Programming

In this report, we study the sensitivity of discrete-time dynamic programs with nonlinear dynamics and objective to perturbations in the initial conditions and reference parameters. Under uniform controllability and boundedness assumptions for the problem data, we prove that the directional derivative of the optimal state and control at time $k$, $x^*_k$, and $u^*_k$, with respect to the reference signal at time $i$, $d_i$, will have exponential decay in terms of $|k-i|$ with a decay rate $\rho$ independent of the temporal horizon length. The key technical step is to prove that a version of the convexification approach proposed by Verschueren et al. can be applied to the KKT conditions and results in a convex quadratic program with uniformly bounded data. In turn, Riccati techniques can be further employed to obtain the sensitivity result, borne from the observation that the directional derivatives are solutions of quadratic programs with structure similar to the KKT conditions themselves. We validate our findings with numerical experiments on a small nonlinear, nonconvex, dynamic program.

97 MATHEMATICS AND COMPUTING↗

Leveraging Application Data Constraints to Optimize Database-Backed Web Applications

Exploiting the relationships among data is a classical query optimization technique. As persistent data is increasingly being created and maintained programmatically, prior work that infers data relationships from data statistics misses an important opportunity. We present Coco, the first tool that identifies data relationships by analyzing database-backed applications. Once identified, Coco leverages the constraints to optimize the application's physical design and query execution. Instead of developing a fixed set of predefined rewriting rules, Coco employs an enumerate-test-verify technique to automatically exploit the discovered data constraints to improve query execution. Each resulting rewrite is provably equivalent to the original query. Using 14 real-world web applications, our experiments show that Coco can discover numerous data constraints from code analysis and improve real-world application performance significantly.

Computer Science↗

Measurement of Close-in Ground Motion from an Underground Chemical Explosion

Understanding the geophysical response near an underground explosion is crucial for generating insights into the source and emplacement conditions that produce distinct observations in monitoring scenarios occurring at greater distances. Recently, Shot A of the Low Yield Nuclear Monitoring (LYNM) Physics Experiment 1 (PE1) series was conducted at the Nevada National Security Site to provide ground truth for subsurface explosion signal models. This experiment resulted in measuring near-source ground motion at distances ranging from 70 to 1000 m/kt with a 99% success rate, yielding high-fidelity knowledge of the near-field response that can serve as benchmarks for future numerical modeling and experiment planning. However, technical challenges exist in observing near-source phenomena while safeguarding sensitive data acquisition components from the detrimental effects of ground motion in the subsurface. This report outlines tools and techniques to address challenges associated with observing near-source accelerations and within the tunnel drift of the PE1 test bed. Additionally, we describe key systems designed with both modern advancements and legacy guidance to maximize the collection of high-quality ground motion data, which may be applied to constitutive and computational models, leading to new or improved understanding of the near- and far-field signals produced by underground explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

To have value, comparisons of high-throughput phenotyping methods need statistical tests of bias and variance

The gap between genomics and phenomics is narrowing. The rate at which it is narrowing, however, is being slowed by improper statistical comparison of methods. Quantification using Pearson’s correlation coefficient ( r ) is commonly used to assess method quality, but it is an often misleading statistic for this purpose as it is unable to provide information about the relative quality of two methods. Using r can both erroneously discount methods that are inherently more precise and validate methods that are less accurate. These errors occur because of logical flaws inherent in the use of r when comparing methods, not as a problem of limited sample size or the unavoidable possibility of a type I error. A popular alternative to using r is to measure the limits of agreement (LOA). However both r and LOA fail to identify which instrument is more or less variable than the other and can lead to incorrect conclusions about method quality. An alternative approach, comparing variances of methods, requires repeated measurements of the same subject, but avoids incorrect conclusions. Variance comparison is arguably the most important component of method validation and, thus, when repeated measurements are possible, variance comparison provides considerable value to these studies. Statistical tests to compare variances presented here are well established, easy to interpret and ubiquitously available. The widespread use of r has potentially led to numerous incorrect conclusions about method quality, hampering development, and the approach described here would be useful to advance high throughput phenotyping methods but can also extend into any branch of science. The adoption of the statistical techniques outlined in this paper will help speed the adoption of new high throughput phenotyping techniques by indicating when one should reject a new method, outright replace an old method or conditionally use a new method.

59 BASIC BIOLOGICAL SCIENCES↗

PyDA: A Hands-On Introduction to Dynamical Data Assimilation with Python

Dynamic data assimilation offers a suite of algorithms that merge measurement data with numerical simulations to predict accurate state trajectories. Meteorological centers rely heavily on data assimilation to achieve trustworthy weather forecast. With the advance in measurement systems, as well as the reduction in sensor prices, data assimilation (DA) techniques are applicable to various fields, other than meteorology. However, beginners usually face hardships digesting the core ideas from the available sophisticated resources requiring a steep learning curve. In this tutorial, we lay out the mathematical principles behind DA with easy-to-follow Python module implementations so that this group of newcomers can quickly feel the essence of DA algorithms. We explore a series of common variational, and sequential techniques, and highlight major differences and potential extensions. We demonstrate the presented approaches using an array of fluid flow applications with varying levels of complexity.

97 MATHEMATICS AND COMPUTING↗

Weak-lensing Detection of Intercluster Filaments in Three Nearby Cluster Systems

Abstract Direct detection of intercluster filaments is challenging due to their low surface density, resulting in a weak deflection field. We present weak-lensing detections of intercluster filaments using wide-field Dark Energy Camera observations from the Local Volume Complete Cluster Survey. A matched-filter method was applied to identify filamentary structures in three nearby ( z < 0.1) systems centered on A401, A2029, and A3558. We discover two filaments (>3 σ ) in each system, with the strongest detections (5.2 σ –5.8 σ ) around A401 and A2029. In particular, we report the first robust weak-lensing detections (≳5 σ ) of the intercluster bridges connecting the cluster pairs A401/399, A2029/2033, A2029/SIG, and A3558/3556. Adopting a filament convergence model motivated by numerical simulations, we infer the maximum convergence ( κ 0 ) and characteristic width ( h c ) for all six filaments, yielding κ 0 ∼ 0.016–0.040 and h c ∼ 0.23–0.43 Mpc. The performance of the matched-filter technique is validated using mock shear catalogs and further tested on a null field around A2351. We explore the potential of using the B-mode lensing signal of filaments to suppress cluster-induced shear contamination. We also quantify the biasing effect of closely separated terminal clusters to the filament signal. These results demonstrate the feasibility of directly mapping dark matter filaments with current and future wide-field weak-lensing datasets.

Shinde, Rahul [Brown University] (ORCID:0000000273↗

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

Hakobyan, Hayk [Flatiron Institute, New York, NY (↗

The impact of stratospheric aerosol intervention on the North Atlantic and Quasi-Biennial Oscillations in the Geoengineering Model Intercomparison Project (GeoMIP) G6sulfur experiment

As part of the Geoengineering Model Intercomparison Project a numerical experiment known as G6sulfur has been designed in which temperatures under a high-forcing future scenario (SSP5-8.5) are reduced to those under a medium-forcing scenario (SSP2-4.5) using the proposed geoengineering technique of stratospheric aerosol intervention (SAI). G6sulfur involves introducing sulfuric acid aerosol into the tropical stratosphere where it reflects incoming sunlight back to space, thus cooling the planet. Here, we compare the results from six Earth-system models that have performed the G6sulfur experiment and examine how SAI affects two important modes of natural variability, the northern wintertime North Atlantic Oscillation (NAO) and the Quasi-Biennial Oscillation (QBO). Although all models show that SAI is successful in reducing global mean temperature as designed, they are also consistent in showing that it forces an increasingly positive phase of the NAO as the injection rate increases over the course of the 21st century, exacerbating precipitation reductions over parts of southern Europe compared with SSP5-8.5. In contrast to the robust result for the NAO, there is less consistency for the impact on the QBO, but the results nevertheless indicate a risk that equatorial SAI could cause the QBO to stall and become locked in a phase with permanent westerly winds in the lower stratosphere.

54 ENVIRONMENTAL SCIENCES↗

Laser-Induced Breakdown Spectroscopy (LIBS) Sensing for Environmental and Subsurface Monitoring

Groundwater monitoring is essential to timely and accurately reflect the current situations and the trends of water quality. However, the long-term stability and survival of any potential in-situ monitoring method is threatened by the harsh conditions and the groundwater monitoring in downhole environments poses numerous challenges to the sensor community. Laser induced breakdown spectroscopy (LIBS) has been demonstrated as a promising technology for chemical monitoring in high-temperature high-pressure (HTHP) environments and hard to reach places. The technique demonstrates several advantages, such as rapid, real-time, in-situ, and simultaneous detection of multiple elements with simple or no sample preparation. This chapter includes a brief review of the field-portable LIBS systems, and development of a compact, robust, and simple LIBS instrument for downhole HTHP water quality monitoring. The fieldable prototype sensor is tested in an onsite monitoring well where trace elements’ concentrations are tracked over an extended period. The testing has verified that the fiber coupled design performs as desired. The system shows good calibration linearity for tested elements and collection times, and Limits of Detection (LODs) that are comparable to those of tabletop LIBS instruments. In addition to groundwater quality monitoring, the fabricated LIBS-based sensor could have widespread sub-surface detection applications.

Jain, Jinesh [NETL Site Support Contractor, Nation↗

Long noncoding RNA MEG3 expressed in human dental pulp regulates LPS-Induced inflammation and odontogenic differentiation in pulpitis

Highlights: • LncRNA MEG3, upregulated in inflamed dental pulp, is involved in the inflammation and regeneration in pulpitis. • The expression of lncRNA MEG3 was semiquantitatively detected in the pulp tissue by the innovative RNAscope® technique. • MEG3 downregulation can inhibit inflammatory cytokines secretion and also promote the odontogenic differentiation of hDPCs. • LncRNA MEG3 can be a potential therapeutic target in the inflammation and regeneration of dentin-pulp complex in pulpitis. Pulpitis refers to inflammation of the inner pulp by invading microbes, and tissue repair occurs due to odontogenic differentiation of human dental pulp cells (hDPCs) with multidifferentiation potential. Long noncoding RNAs (lncRNAs) can modulate numerous pathological and biological processes; however, the role of lncRNAs in the inflammation and regeneration of the dentin-pulp complex in pulpitis is unclear. Here, we performed high-throughput sequencing to identify differentially expressed lncRNAs between human normal and inflamed pulp and concluded that lncMEG3 (lncRNA maternally expressed gene 3, MEG3) was signicantly upregulated in both inflamed pulp and LPS-treated hDPCs. MEG3 expression in the pulp tissue was detected using the RNAscope® technique. RNA pulldown assays identified the MEG3-interacting proteins and the potential mechanisms. With MEG3 knockdown, we investigated the role of MEG3 in the secretion of inflammatory cytokines in LPS-treated hDPCs and odontogenic differentiation of hDPCs. MEG3 downregulation inhibited the secretion of TNF-α, IL-1β and IL-6 in LPS-treated hDPCs, and the p38/MAPK signaling pathway may be related to this effect. MEG3 knockdown promoted odontogenic differentiation of hDPCs by regulating the Wnt/β-catenin signaling pathway. Our study suggested that MEG3 has a negative effect on inflammation and regeneration of the dentin-pulp complex in pulpitis.

60 APPLIED LIFE SCIENCES↗

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗

Dynamic Performance Enhancement for Nonlinear Stochastic Systems Using RBF Driven Nonlinear Compensation with Extended Kalman Filter

In this paper, a novel hybrid control method is proposed to enhance the control performance of the PI based control system for a class of nonlinear and non-Gaussian stochastic dynamic processes with unmeasurable states. Firstly, to enhance the tracking performance of the PI controller where the PI parameters are fixed in many actual control systems, the compensative signal is formed using the extended Kalman filter(EKF) based state estimator and driven by the radial basis function (RBF) neural network based compensator. In addition, the weights of RBF is trained to minimize the entropy criterion of tracking error as the process is subjected to non-Gaussian disturbances. Meanwhile, since the precise statistical property of noises is hard to obtain for many industrial processes, the kernel density estimation (KDE) technique is employed in this paper to estimate the entropy. The convergence of RBF network is discussed and the stability of the resulted closed-loop hybrid control system is analyzed in mean square sense. Finally, a numerical example and a practical system testing are given to illustrate the effectiveness of the proposed control method.

stochastic systems, PI Controller, Entropy, Extend↗

Reduced Order Modeling conditioned on monitored features for response and error bounds estimation in engineered systems

Reduced Order Models (ROMs) form essential tools across engineering domains by virtue of their function as surrogates for computationally intensive digital twinning simulators. Although purely data-driven methods are available for ROM construction, schemes that allow to retain a portion of the physics tend to enhance the interpretability and generalization of ROMs. However, physics-based techniques can adversely scale when dealing with nonlinear systems that feature parametric dependencies. This study introduces a generative physics-based ROM that is suited for nonlinear systems with parametric dependencies and is additionally able to provide numerical error bounds associated with the respective estimates. A main contribution of this work is the conditioning of these parametric ROMs to features that can be derived from monitoring measurements, feasibly in an online fashion. This is contrary to most existing ROM schemes, which remain restricted to the prescription of the physics-based, and usually a priori unknown, system parameters. Our work utilizes conditional Variational Autoencoders to continuously map the required reduction bases to a feature vector extracted from limited output measurements, while additionally allowing for a probabilistic assessment of the ROM-estimated Quantities of Interest. An auxiliary task using a neural network-based parametrization of suitable probability distributions is introduced to re-establish the link with physical model parameters. We verify the proposed scheme on a series of simulated case studies incorporating effects of geometric and material nonlinearity under parametric dependencies related to system properties and input load characteristics.

Conditional VAEs↗

Confocal Raman Microscopy as an In Situ Probe of Volume Change in Hydration-Sensitive Polymer Membranes

An approach is described for measuring hydration-induced volume change within ion conductive polymer membranes and thin films by adapting the technique of confocal Raman microscopy. With careful consideration of factors that affect excitation and scattering within the confocal probe volume region, material swelling and deswelling were estimated from the intensities of polymer matrix spectral features. A high numerical aperture oil-immersion objective was used to achieve tight focusing within samples and the efficient collection of Raman scattered light. The approach is demonstrated for a fluorinated cation-exchange ionomer (Nafion) and the hydrocarbon anion-exchange ionomer Sustainion. Samples were monitored while under a nitrogen atmosphere that was cycled between dry and humid (∼50% relative humidity, RH) states. The volume changes estimated from the Raman spectra were in close agreement with those derived from conventional measurements. The reported work advances understanding needed to adapt confocal Raman microscopy for quantitative in situ and operando studies of ionomers within electrochemical devices, such as polymer electrolyte membrane fuel cells and electrolyzers, and lays a foundation for broader applications in the study of hydration-induced polymer swelling and the associated molecular level water and polymer framework structural changes.

Absorption↗

Hyper‐Local Temperature Prediction Using Detailed Urban Climate Informatics

The accurate modeling of urban microclimate is a challenging task given the high surface heterogeneity of urban land cover and the vertical structure of street morphology. Recent years have witnessed significant efforts in numerical modeling and data collection of the urban environment. Nonetheless, it is difficult for the physical‐based models to fully utilize the high‐resolution data under the constraints of computing resources. The advancement in machine learning (ML) techniques offers the computational strength to handle the massive volume of data. In this study, we proposed a modeling framework that uses ML approach to estimate point‐scale street‐level air temperature from the urban‐resolving meso‐scale climate model and a suite of hyper‐resolution urban geospatial data sets, including three‐dimensional urban morphology, parcel‐level land use inventory, and weather observations from a sensor network. We implemented this approach in the City of Chicago as a case study to demonstrate the capability of the framework. The proposed approach vastly improves the resolution of temperature predictions in cities, which will help the city with walkability, drivability, and heat‐related behavioral studies. Moreover, we tested the model's reliability on out‐of‐sample locations to investigate the modeling uncertainties and the application potentials to the other areas. This study aims to gain insights into next‐gen urban climate modeling and guide the observation efforts in cities to build the strength for the holistic understanding of urban microclimate dynamics.

54 ENVIRONMENTAL SCIENCES↗

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗

Embedding hard physical constraints in neural network coarse-graining of three-dimensional turbulence

In recent years, deep learning approaches have shown much promise in modeling complex systems in the physical sciences. A major challenge in deep learning of partial differential equations is enforcing physical constraints and boundary conditions. In this work, we propose a general framework to directly embed the notion of an incompressible fluid into convolutional neural networks, and apply this to coarse-graining of turbulent flow. These physics-embedded neural networks leverage interpretable strategies from numerical methods and computational fluid dynamics to enforce physical laws and boundary conditions by taking advantage the mathematical properties of the underlying equations. Here, we demonstrate results on three-dimensional fully developed turbulence, showing that this technique drastically improves local conservation of mass, without sacrificing performance according to several other metrics characterizing the fluid flow.

97 MATHEMATICS AND COMPUTING↗

Learning Conservation Laws in Unknown Quantum Dynamics

We present a learning algorithm for discovering conservation laws given as sums of geometrically local observables in quantum dynamics. This includes conserved quantities that arise from local and global symmetries in closed and open quantum many-body systems. The algorithm combines the classical shadow formalism for estimating expectation values of observable and data analysis techniques based on singular value decompositions and robust polynomial interpolation to discover all such conservation laws in unknown quantum dynamics with rigorous performance guarantees. Our method can be directly realized in quantum experiments, which we illustrate with numerical simulations, using closed and open quantum system dynamics in a ℤ 2 gauge theory and in many-body localized spin chains.

Machine learning↗