Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods”

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 1,189 records · Page 66

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗

CAMERA: A method for cost-aware, adaptive, multifidelity, efficient reliability analysis

Estimating probability of failure in aerospace systems is a critical requirement for flight certification and qualification. Failure probability estimation involves resolving tails of probability distributions, and Monte Carlo sampling methods are intractable when expensive high-fidelity simulations have to be queried. Here, we propose a method to use models of multiple fidelities that trade accuracy for computational efficiency. Specifically, we propose the use of multifidelity Gaussian process models to efficiently fuse models at multiple fidelity, thereby offering a cheap surrogate model that emulates the original model at all fidelities. Furthermore, we propose a novel sequential acquisition function based experiment design framework that can automatically select samples from appropriate fidelity models to make predictions about quantities of interest at the highest fidelity. We use our proposed approach in an importance sampling setting and demonstrate our method on the failure level set and probability estimation on synthetic test functions and two real-world applications, namely, the reliability analysis of a gas turbine engine blade using a finite element method and a transonic aerodynamic wing test case using Reynolds-averaged Navier-Stokes equations. We show that our method predicts the failure boundary and probability more accurately and at a fraction of the computational cost compared with using just a single expensive high-fidelity model. Finally, we show that our sequential approach is guaranteed to asymptotically converge to the true failure boundary with high probability.

97 MATHEMATICS AND COMPUTING↗

Sparse-grid discontinuous Galerkin methods for the Vlasov–Poisson–Lenard–Bernstein model

Sparse-grid methods have recently gained interest in reducing the computational cost of solving high-dimensional kinetic equations. In this paper, we construct adaptive and hybrid sparse-grid methods for the Vlasov–Poisson–Lenard–Bernstein (VPLB) model. This model has applications to plasma physics and is simulated in two reduced geometries: a 0x3v space homogeneous geometry and a 1x3v slab geometry. Here we use the discontinuous Galerkin (DG) method as a base discretization due to its high-order accuracy and ability to preserve important structural properties of partial differential equations. We utilize a multiwavelet basis expansion to determine the sparse-grid basis and the adaptive mesh criteria. We analyze the proposed sparse-grid methods on a suite of three test problems by computing the savings afforded by sparse-grids in comparison to standard solutions of the DG method. The results are obtained using the adaptive sparse-grid discretization library ASGarD.

97 MATHEMATICS AND COMPUTING↗

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗

Performance evaluation of automated data-driven feature extraction and selection methods for practical and scalable building energy consumption prediction models

Here, this study quantifies the impact of automated feature engineering methods (feature extraction and selection) on the quality and accuracy of machine learning models that predict building energy consumption. The case study compares model performance for three main scenarios: baseline (no feature extraction and selection), feature extraction only, and feature extraction combined with feature selection (filter and/or wrapper methods) for fully trained machine learning models for 200 metered/sub-metered energy measurements across 118 real buildings. For consistency, the same machine learning model architecture (a black box deep learning neural network with probabilistic forecast output) was used for all scenarios. Based on results, all feature engineering methods provided noticeable prediction accuracy improvements (e.g., 29%-68% median prediction improvement) compared to baseline scenarios. However, in this application, feature selection methods provide little practical value due to their limited performance gains and high computational cost. Smarter algorithm development supported by better computational environments will be needed before feature selection methods can reliably and efficiently improve predictive model performance.

97 MATHEMATICS AND COMPUTING↗

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗

Accurately measuring the ionic conductivity of membranes via the direct contact method

Accurately measuring the ionic conductivity of membranes is important both for assessing the performance of new membranes as well as for advancing the fundamental understanding of electric field-driven ion transport in membranes. Yet, despite a considerable number of reports on the topic, there is no standardized method for ionic conductivity measurements. One of the main challenges with ionic conductivity measurements is isolating the membrane resistance from the other resistances in the electrochemical cell. In this study, we present a detailed discussion of factors that affect ionic conductivity measurements with the commonly used direct contact method, in which the membrane contacts the electrodes. Such factors include external resistances due to imperfect contact between the membrane and electrodes, residual impedances due to the sample holder and wiring, and fringe effects. Further, we also present and validate simple methods for isolating the true membrane ionic conductivity in both the through-plane and in-plane orientations. These methods can be implemented with virtually any membrane to obtain accurate and reproducible ionic conductivity values. Finally, we outline recommendations for performing ionic conductivity measurements with the direct contact method and a set of parameters to verify when reporting these results.

36 MATERIALS SCIENCE↗

Novel method for accurately estimating membrane transport properties and mass transfer coefficients in reverse osmosis

Here, we present a simple and robust method to simultaneously characterize the water and salt permeability (A, B) of reverse osmosis (RO) membranes and mass transfer coefficient (k) in membrane modules. The proposed methodology comprises a set of RO experiments performed at different operating pressures or stages. The measured water and salt fluxes in each stage are simultaneously fitted to the RO transport equations by performing a non-linear regression, using A, B, and k as regression parameters. We first perform a systematic accuracy analysis of the proposed method across the full operational range of RO. The assessment shows that the method accuracy is substantially higher than current methods and increases with number of experimental stages and driving forces. This assessment is used to inform the design of an experimental protocol that minimizes errors in estimated A, B, and k. We then evaluate two commercial RO membranes following the new protocol. For both membranes, A and B parameters decrease by 17% and 15% from the dilute solution to seawater concentrations, whereas the k parameter remains constant. Our study demonstrates that the proposed method, informed by data-driven experimental designs, provides a new approach for accurately characterizing transport phenomenon in membrane processes with feeds of less than 100 g/L total dissolved solids.

36 MATERIALS SCIENCE↗

Evaluation of a high-throughput method for processing sponge-stick samples to detect viable, non-spore-forming biothreat agents

After a bioterrorism incident, surface sampling is often used to determine the extent of contamination and exposure, guiding decontamination efforts and decisions for re-occupancy of affected sites. The sponge-stick (SS) is a preferred and commonly used device for sample collection to detect both spore-forming and non-spore-forming biothreat agents from non-porous surfaces. Here, in this study, a recently developed high-throughput method (HTM) for processing SS samples to detect viable Bacillus anthracis spores was adapted for detection of non-spore-forming biothreat agents, Yersinia pestis and Francisella tularensis. The scalable HTM was used to process up to 20 SS samples simultaneously, compared to the current stomacher-based method which processes one SS at a time. Comparisons of the HTM and the stomacher-based method were statistically indistinguishable for most experiments (P > 0.05) with HTM recoveries of 37–60 % for Y. pestis inoculated at 102–103 cells/SS and held 48 h at 4 °C to mimic sample transport/storage. The HTM was integrated with Rapid Viability-Polymerase Chain Reaction (RV-PCR) analysis to detect viable Y. pestis in the presence of particulate contamination (Arizona Test Dust, ATD). This approach detected Y. pestis inoculated at 20 cells/SS and ATD did not impact detection (P > 0.05). F. tularensis showed significantly lower recoveries between no-hold time and 48-h hold time (4 °C, P < 0.05) using the HTM, which further testing showed could be due to toxicity of the neutralizing buffer used for SS pre-wetting. With modifications, this method could enhance throughput capacity while maintaining similar recovery efficiencies to current methods for other non-spore-forming bacterial pathogens.

Biological and medical sciences↗

New method for level-lifetime measurements with thick scintillators

Level lifetimes provide key insight into the structure of atomic nuclei and serve as stringent tests of theoretical descriptions. Though several methods for determining level lifetimes exist for both reaction measurements and decay studies, here the focus is on techniques involving the direct measurement of time differences between population and subsequent depopulation of excited states. The techniques presented herein are broadly applicable across multiple timing ranges, but the approach is specifically described for the βγ timing method. A multi-step, amplitude-dependent time walk correction was employed to address the data analysis complications that arise from using thick scintillators for electron detection. Further, additional corrections for the depth of the interaction must also be performed when the parent isotope has been implanted into the detector. Techniques are presented for performing these time walk and depth of interaction corrections. Subsequently, a new Monte Carlo method utilizing measured detector responses obtained from the data, coupled with chi-square minimization, is presented for extracting excited state lifetimes ≳ 100 picoseconds. The framework of this Monte Carlo method is developed for the decay of a state in 68 Zn with a known 1.6 ps half life, which is considered prompt given the detection sensitivity of the technique, and then benchmarked using two other excited states in neutron-rich Ni isotopes with 120(34) ps and 1.05(3) ns half lives. Using this new method which takes into account the thick scintillator used, these same half lives were measured to be 135(10) ps and 1.04(24) ns, respectively. The overall good agreement demonstrates the validity of the technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Diagnosing nuclear power plant pipe wall thinning due to flow accelerated corrosion using a passive, thermal non-destructive evaluation method: Feasibility assessment via numerical experiments

Flow accelerated corrosion (FAC) in nuclear power plant pipes is one of the leading causes of accidents, fatalities, damage and outages. Current FAC identification methods employ expensive sensing technology and are “active” methods, where the response of the piping system to an externally-generated thermal, mechanical or optical excitation must be measured. As a result, these techniques require a disruptive and time-consuming setup. Here we propose a method that utilizes pipe surface temperature measurements to passively monitor for FAC-induced pipe wall thinning without the need for expensive equipment or post-installation setup time. This diagnostic method utilizes a simulation data-driven diagnostic model to estimate the amount of thickness reduction in a pipe based on changes in measured steady-state pipe temperatures. In order to reduce the computational burden of generating large, simulation-based datasets, the behavior of the insulation of the pipe was modeled using a suitably calibrated heat transfer boundary parameter. Additionally, global sensitivity analysis was performed to determine system parameter(s), such as the temperature of water flowing inside the pipe, which significantly affect the steady state pipe wall temperature and could cause errors in diagnosis. Two diagnostic models, one using only the change in steady-state temperature as an indicator for FAC-induced pipe wall thinning and the other using water temperature as an additional diagnostic model input were evaluated for their ability to estimate thickness reductions in a pipe using simulated pipe wall temperature data. For the numerical experiments conducted in this work, both models estimated wall thickness with errors within 0.5 mm, indicating that the proposed technique can potentially be used as a low-cost, first-pass method for FAC monitoring.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Use of high throughput screening methods to study dual-functional crosslinkable latexes

Herein, structure-formulation-property relationships for crosslinkable core-shell latexes were studied initially using high throughput screening methods and then using conventional experimental methods to determine optimum coating formulations. Two approaches that are used within coatings technology to improve sustainability are reduction in volatile organic compounds (VOC) and reduction in energy needed to cure coatings. Hybrid crosslinkable core-shell latexes have shown potential in meeting these challenges where film properties provided by particle coalescence are not sufficient and chemical crosslinking is needed. However, with numerous parameters involved such as latex particle architecture, nature and amount of crosslinker, and curing conditions, the screening of formulations via conventional coating testing methods becomes laborious and time-consuming. The results from high-throughput dye extraction and puncture testing were compared to data from conventional testing methods such as chemical resistance (MEK double rubs), pendulum hardness, and tensile testing to identify correlations. The dye extraction results have shown a strong correlation to solvent resistance and pendulum hardness data whereas toughness from puncture testing data correlated to tensile testing. Both high throughput dye extraction and puncture testing methods can be used for the accelerated screening of crosslinkable coating formulations to optimize curing parameters to achieve the best coating performance.

36 MATERIALS SCIENCE↗

Unpacking the drivers of diurnal dynamics of sun-induced chlorophyll fluorescence (SIF): Canopy structure, plant physiology, instrument configuration and retrieval methods

Sun-induced chlorophyll fluorescence (SIF) from spaceborne sensors is a promising tool for global carbon cycle monitoring, but its application is constrained by insufficient understanding of the drivers underlying diurnal SIF dynamics. SIF measurements from ground-based towers can reveal diurnal SIF dynamics across biomes and environmental conditions; however, meaningful interpretation of diurnal variations requires disentangling impacts from canopy structure, plant physiology, instrument configuration and retrieval methods, which often interact with and confound each other. This study aims to unpack these drivers using 1) concurrent ground and airborne canopy-scale and leaf-scale measurements at a corn field, 2) a mechanistic SIF model that explicitly considers the dynamics of photochemistry (via the fraction of open photosystem II reaction centers, qL) and photoprotection (via nonphotochemical quenching, NPQ) as well as their interactive dependence on the sub-canopy light environment, and 3) cross-comparison of SIF instrument configurations and retrieval methods. We found that crop row orientations and sun angles can introduce a distinctive midday dip in SIF in absence of stress, due to a midday drop of absorbed photosynthetically active radiation (APAR) when crop rows are north-south oriented. Canopy structure caused distinctive responses in both qL and NPQ at different positions within the vertical canopy that collectively influenced fluorescence quantum yield (Φ F ) at the leaf scale. Once integrated at the canopy scale, diurnal dynamics of both APAR and canopy escape probability (ε) are critical for accurately shaping diurnal SIF variations. While leaf-level qL and NPQ exhibited strong diurnal dynamics, their influence was attenuated at the canopy scale due to opposing effects on SIF at different canopy layers. Furthermore, different system configurations (i.e., bi-hemispherical vs. hemispherical-conical) and retrieval methods can bias the SIF magnitude and distort its diurnal shapes, therefore confounding the interpretation of inherent strength and dynamics of SIF emission. Our findings demonstrate the importance of crop row structures, interactive variations in canopy structure and plant physiology, instrument configuration, and retrieval method in shaping the measured dynamics of diurnal SIF. This study highlights the necessity to account for these factors to accurately interpret satellite SIF, and informs future synthesis work with different SIF instrumentation and retrieval methods across sites.

59 BASIC BIOLOGICAL SCIENCES↗

A polarimetry-based field-deployable non-interruptive mirror soiling detection method

The soiling level of heliostat mirrors in Concentrated Solar Power (CSP) fields is one of the key factors that significantly influences optical efficiency. State-of-the-art methods of monitoring heliostats soiling levels still face various challenges, including slow speed, labor-intensive operations, resolution and accuracy constraints or interruptions to solar field operations. Here, we present a rapid, cost-effective, and non-intrusive method for mirror soiling detection based on polarimetric imaging, referred to as Polarimetric Imaging-based Mirror Soiling (PIMS). The compact PIMS device is designed for integration with unmanned aerial vehicles (UAVs), enabling rapid, large-area assessments of heliostat mirrors for efficient soiling detection. Our method utilizes the correlation between the Degree of Linear Polarization (DoLP) and surface soiling level based on Mie scattering theory and Monte Carlo simulations. Field deployment of the PIMS method requires minimal device installation, and its UAV-based operation allows for soiling detection without interrupting plant activities. The PIMS method holds the potential for mirror soiling detection across various concentrated solar power (CSP) plants and can be further adapted for other types of solar fields, such as parabolic trough systems.

CSP Field↗

A Quantitative Method for In Situ Pump-Beam Metrology in 4D Ultrafast Electron Microscopy

In this work, we report a method for measuring spot size and focusing conditions of the femtosecond (fs) excitation laser in situ at the specimen location in 4D ultrafast electron microscopy (UEM). The method makes use of threshold laser ablation behaviors of thin amorphous carbon membranes. For Gaussian beam profiles and for ablation threshold fluence values, we analytically derive expressions describing the relationship between ablated hole size and the actual laser spot size. Using these expressions, we developed experimental procedures for characterizing the shape and spot size of the pump beam at the specimen. We demonstrate the viability of the approach for incident excitation wavelengths of 343 nm and 515 nm, thus illustrating the applicability of the method to a range of optical wavelengths without modification. Further, we show that by measuring ablated hole size as a function of focusing condition, a full metrological characterization of the Gaussian beam propagation properties can be performed. Finally, we find good agreement for spot sizes determined with this method and with those determined by extrapolation from measurements taken outside the microscope column. Overall, this method is a simple, cost-effective means for accurate and precise determination of key pump-beam parameters in situ at the specimen location in UEM experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development and experimental verification of an active noise cancellation (ANC) method for non-invasive sensing

Ultrasonic wave propagation may be used in non-invasive sensing to identify the properties of a fluid in a container. An ultrasonic excitation can be coupled into the liquid or the gas through the walls of the container and detected by a receiver on the other side. However, in addition to this fluid path, another signal will also propagate through the structure and is referred as noise in this study. For most materials, this structural noise will be of higher amplitude compared to the fluid signal due to the large acoustic impedance difference between the solid container and the fluid. In this paper, an active noise cancellation (ANC) method is developed and demonstrated on an experimental setup. It was shown that this method cancels more than 90% of the noise signal in a given time window of interest. A sequential version of this cancellation method is also proposed, and it was shown that its performance improves as the number of canceling excitations is increased. The performance was found to be independent of the location of canceling transducers and position of the silence window but dependent on the width of the silence window. As an example, the ANC method was applied to noise cancellation of a signal travelling through water in a stainless-steel container. The signal-to-noise ratio (SNR) of the water signal in the window of interest was increased by 148% from 13.98 dB to 34.68 dB by using the proposed ANC method.

42 ENGINEERING↗

Quantifying Phospholipids in Organic Samples Using a Hydrophilic Interaction Liquid Chromatography–Inductively Coupled Plasma High-Resolution Mass Spectrometry (HILIC-ICP-HRMS) Method

Here, in this study, a novel method using hydrophilic interaction liquid chromatography (HILIC) coupled with inductively coupled plasma high-resolution mass spectrometry (ICP-HRMS) was introduced for the quantification of phospholipids in oil samples. The method employed a bridged ethyl hybrid (BEH) stationary phase HILIC column with a tetrahydrofuran (THF)/water mobile phase, enhancing the solubility and detection of phospholipids. During the study, a gradient/matrix effect on ICP-HRMS sensitivity was observed and successfully compensated for experimentally, ensuring reliable quantification results. This approach has proven effective for a wide range of different oil samples including vegetable oils, animal fats, and phospholipid supplements. Notably, this method allowed the direct quantification of phospholipids in oil samples, bypassing the need for prior sample preparation methods, such as solid phase extraction (SPE), thereby streamlining the analytical process. The precision, accuracy, and reduced need for extensive sample preparation offered by this method mark a significant advancement in lipids analysis. Its robustness and broad applicability have substantial implications for industries such as food and renewable energy production, where both efficient and accurate lipid identification and quantification are crucial.

09 BIOMASS FUELS↗