Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational 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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

A Framework Based on 2-D Taylor Expansion for Quantifying the Impacts of Subpixel Reflectance Variance and Covariance on Cloud Optical Thickness and Effective Radius Retrievals Based on the Bispectral Method

The bispectral method retrieves cloud optical thickness (τ) and cloud droplet effective radius (re) simultaneously from a pair of cloud reflectance observations, one in a visible or near-infrared (VIS/NIR) band and the other in a shortwave infrared (SWIR) band. A cloudy pixel is usually assumed to be horizontally homogeneous in the retrieval. Ignoring subpixel variations of cloud reflectances can lead to a significant bias in the retrieved τ and re. In the literature, the retrievals of τ and re are often assumed to be independent and considered separately when investigating the impact of subpixel cloud reflectance variations on the bispectral method. As a result, the impact on τ is contributed only by the subpixel variation of VIS/NIR band reflectance and the impact on re only by the subpixel variation of SWIR band reflectance. In our new framework, we use the Taylor expansion of a two-variable function to understand and quantify the impacts of subpixel variances of VIS/NIR and SWIR cloud reflectances and their covariance on the τ and re retrievals. This framework takes into account the fact that the retrievals are determined by both VIS/NIR and SWIR band observations in a mutually dependent way. In comparison with previous studies, it provides a more comprehensive understanding of how subpixel cloud reflectance variations impact the τ and re retrievals based on the bispectral method. In particular, our framework provides a mathematical explanation of how the subpixel variation in VIS/NIR band influences the re retrieval and why it can sometimes outweigh the influence of variations in the SWIR band and dominate the error in re retrievals, leading to a potential contribution of positive bias to the re retrieval. We test our framework using synthetic cloud fields from a large-eddy simulation and real observations from Moderate Resolution Imaging Spectroradiometer. The predicted results based on our framework agree very well with the numerical simulations. Our framework can be used to estimate the retrieval uncertainty from subpixel reflectance variations in operational satellite cloud products and to help understand the differences in τ and re retrievals between two instruments.

reflectance↗
NASA NTRS1978

Computation of three-dimensional nozzle-exhaust flow fields with the GIM code

A methodology is introduced for constructing numerical analogs of the partial differential equations of continuum mechanics. A general formulation is provided which permits classical finite element and many of the finite difference methods to be derived directly. The approach, termed the General Interpolants Method (GIM), can combined the best features of finite element and finite difference methods. A quasi-variational procedure is used to formulate the element equations, to introduce boundary conditions into the method and to provide a natural assembly sequence. A derivation is given in terms of general interpolation functions from this procedure. Example computations for transonic and supersonic flows in two and three dimensions are given to illustrate the utility of GIM. A three-dimensional nozzle-exhaust flow field is solved including interaction with the freestream and a coupled treatment of the shear layer. Potential applications of the GIM code to a variety of computational fluid dynamics problems is then discussed in terms of existing capability or by extension of the methodology.

Spradley, L. W.↗
DOE OSTI2020

Integrating Predictions for Improving Defect Classification Accuracy in NDT-based Assessment of Concrete - 20229

There is an increasing need to create predictive models for defect classification in concrete using the output of non-destructive testing (NDT) techniques. Recent advancement of machine-learning algorithms has offered several techniques for developing classification models for different types of data sets. However, the performance of these algorithms is very uncertain, mainly when applied to small and noisy datasets. For example, when human access is limited (e.g., nuclear facility), robot-based NDT is preferred. But compared to manual tests with humans present on site, the data sets are small, and more noise can exist. Therefore, it is imperative to develop new approaches to ensure a consistently high classification accuracy for inadequate data sets. This study explores the classification performance on NDT dataset using classifiers from different machine-learning algorithms, namely k-Nearest Neighbor (kNN), Decision Tree, Naive Bayes, Logistic Regression, and Support Vector Machine (SVM). The authors further integrated the predictions from these classifiers using proposed methods. The integration strategy combines the output of the classifiers based on two different measures, accuracy, and performance (ACC and PERF), using equations such as sum, average, and square-root-of-sums-of-squares (SRSS). Our results reveal varying classification accuracies across individual classifiers with different misclassifications across the test data set. The integration strategy provided significant improvement in the classification accuracy compared to the individual classifiers. Furthermore, the results indicate minimal variation across the integration methods as compared to the variation across the individual classifiers. To conclude, prediction integration offers a unique approach for combining the output of multiple classifiers to create redundancies with the potential of achieving high classification performance and improved reliability in predictive models for defect detection in concrete. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Capturing the Relative‐Humidity‐Sensitive Gas–Particle Partitioning of Organic Aerosols in a 2D Volatility Basis Set

Abstract Aerosol water affects the physicochemical properties and mass concentration of organic aerosols (OA), but it is typically omitted by air quality, weather, and climate models. We compare two classes of simplified models to estimate the OA water uptake and gas–particle partitioning of organic compounds. One class uses a single‐hygroscopicity‐parameter (κ) approach while the other is based on the reduced‐complexity Binary Activity Thermodynamics (BAT) model. We show that a BAT‐based two‐dimensional volatility basis set (VBS) model always predicts a higher OA mass concentration at elevated relative humidity (RH), for example, ∼16% at 80% RH, than any variation of theκ‐based method considered—even when BAT‐VBS predicts a lower water uptake. The main reason being that the BAT‐VBS model captures variations in effective saturation mass concentration of organics (C*) with RH, a feature that other VBS methods lack. The BAT‐VBS framework offers an efficient, RH‐sensitive treatment for reduced‐complexity OA modeling.

Geology↗
DOE OSTI2024

Residual Stress Measurements of Cr-coated Cladding

Accident-tolerant fuel concepts for light-water reactor applications have been developed and tested in diverse research programs around the world. Industry teams have developed coated cladding concepts using different coating application methods, such as physical vapor deposition and coldspraying. Recently, variations in thickness and application method were investigated, and their effects on the microstructure and mechanical properties were reported. The results suggest an increase in performance in an unirradiated state, however, it remains unclear what phenomenon drove this performance. One explanation to this performance is associated with the coating process imparting a residual stress. Residual stress analysis was performed on coatings produced with various process parameters to understand how a thin layer of Cr impacts the mechanical properties after coating application and after simulated operating conditions. The gathered information is intended to address safety concerns or support margin identification related to accident-tolerant fuels and to support subsequent modeling and simulation efforts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1988

Some notes on applying the one-step multiconic method of trajectory propagation

The one-step multiconic method (also called the pseudostate method) has been known for years to offer up to an order of magnitude of improvement over the patched-conic method for propagating trajectories under the gravitational influence of two bodies. In these notes, the method and some of its variations are described. Then, a crucial parameter of this method, the sweepback duration, is examined, and it is shown that the traditional wisdom concerning this parameter should be replaced by an algorithm which relates it to the trajectory characteristics. Remarkably, this algorithm is independent of the masses of the primary and secondary bodies.

Sweetser, Theodore H., III↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Wind regime peculiarities in the lower thermosphere in the winter of 1983/84

Temporal variations of prevailing winds at 90 to 100 km obtained from measurements carried out in winter 1983 to 1984 at three sites in the USSR and two sites in East Germany are reported. These variations are compared with those of the thermal stratospheric regime. Measurements were carried out using the drifts D2 method (meteor wind radar) and the D1 method (ionospheric drifts). Temporal variations of zonal and meridional prevailing wind components for all the sites are given. Also presented are zonal wind data obtained using the partial reflection wind radar. Wind velocity values were obtained by averaging data recorded at between 105 and 91 km altitude. Wind velocity data averaged in such a way can be related to about the same height interval to which the data obtained by the meteor radar and ionospheric methods at other sites, i.e., the mean height of the meteor zone (about 95 km). The results presented show that there are significant fluctuations about the seasonal course of both zonal and meridional prevailing winds.

Lysenko, I. A.↗
DOE OSTI2026

Reweighting Underlying Event and Colour Reconnection parameter variations in Sherpa

We propose and validate a new method to trace the impact of parameter variations in the simulation of multi-parton interactions and colour reconnections in the Sherpa event generator. They are reflected, at an event-by-event basis, through relative weights with respect to the central production parameters that give rise to the generated events and distributions. Our method facilitates the tuning of the Monte Carlo event generator at a dramatically reduced computational cost, alleviates parameter sensitivity studies, and enables robust quantification of parametric uncertainties on-the-fly, one of the missing ingredients for future simulations of high-energy particle collisions. The method can easily be adapted to and implemented in other event generators. To illustrate its potential, we here consider combined tunes of the multi-parton-interaction and colour-reconnection models in Sherpa using LHC proton-proton collision data at $\sqrt{s}=7\,\text{TeV}$. We furthermore calibrate the energy-scaling behaviour of dimensionful model parameters based on $\sqrt{s}=13\,\text{TeV}$ LHC data and Tevatron data taken at $\sqrt{s}=1.96\,\text{TeV}$.

Pabst, Moritz [Gottingen U.]↗
NASA NTRS1997

Evaluation of heart rate changes: electrocardiographic versus photoplethysmographic methods

The heart rate (HR) variation to forced deep breathing (HRDB) and to the Valsalva maneuver (Valsalva ratio; VR) are the two most widely used tests of cardiovagal function in human subjects. The HR is derived from a continuously running electrocardiographic (ECG) recording. Recently, HR derived from the arterial waveform became available on the Finapres device (FinapHR), but its ability to detect rapid changes in HR remains uncertain. We therefore evaluated HRDB and VR derived from FinapHR using ECG-derived HR (ECGHR) recordings as the standard. We also compared the averaged HR on Finapres (Finapav) with beat-to-beat Finapres (FinapBB) values. Studies were undertaken in 12 subjects with large HR variations: age, 34.5 +/- 9.3 (SD) years; six males and six females. FinapBB values were superimposable upon ECGHR for both HRDB and VR. In contrast, Finapav failed to follow ECGHR for HRDB and followed HRECG with a lag for the VR. To evaluate statistically how closely FinapHR approximated ECGHR, we undertook regression analysis, using mean values for each subject. To compare the two methods, we evaluated the significance of the difference between test and standard values. For HRDB, FinapBB reproducibly recorded HR (R2 = 0.998), and was significantly (p = 0.001) better than Finapav (R2 = 0.616; p < 0.001). For VR, HRBB generated a VR that was not significantly different from the correct values, while HRav generated a value that was slightly but consistently lower than the correct values (p < 0.001). We conclude that FinapHR reliably records HR variations in the beat-to-beat mode for cardiovascular HR tests.

Non-NASA Center↗
DOE OSTI2023

Simulating scalar field theories on quantum computers with limited resources

Here we present a quantum algorithm for implementing φ 4 lattice scalar field theory on qubit computers. The field is represented in the discretized field amplitude basis. The number of qubits and elementary gates required by the implementation of the evolution operator is proportional to the lattice size. The algorithm allows efficient φ 4 state preparation for a large range of input parameters in both the normal and broken-symmetry phases. The states are prepared using a combination of variational and adiabatic evolution methods. First, the ground state of a local Hamiltonian, which includes the φ 4 self-interaction, is prepared using short variational circuits. Next, this state is evolved by switching on the coupling between the lattice sites adiabatically. The parameters defining the local Hamiltonian are adjustable and constitute the input of our algorithm. We present a method to optimize these parameters in order to reduce the adiabatic time required for state preparation. For preparing broken-symmetry states, the adiabatic evolution problems caused by crossing the phase transition critical line and by the degeneracy of the broken-symmetry ground state can be addressed using an auxiliary external field which gradually turns off during the adiabatic process. We show that the time dependence of the external field during the adiabatic evolution is important for addressing the broken-symmetry ground state degeneracy. The adiabatic time dependence on the inverse error tolerance can be reduced from quadratic to linear by using a field strength that decreases exponentially in time relative to one that decreases linearly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Synchronous and Concurrent Multidomain Computing Method for Cloud Computing Platforms

We present a numerical method for synchronous and concurrent solution of transient elastodynamics problem where the computational domain is divided into subdomains that may reside on separate computational platforms. Here, this work employs the variational multiscale discontinuous Galerkin (VMDG) method to develop interdomain transmission conditions for transient problems. The fine-scale modeling concept leads to variationally consistent coupling terms at the common interfaces. The method admits a large class of time discretization schemes, and decoupling of the solution for each subdomain is achieved by selecting any explicit algorithm. Numerical tests with a manufactured solution problem show optimal convergence rates. The energy history in a free vibration problem is in agreement with that of the solution from a monolithic computational domain.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

A New Method to Retrieve the Diurnal Variability of Planetary Boundary Layer Height from Lidar Under Different Thermodynamic Stability Conditions

The planetary boundary layer height (PBLH) is an important parameter for understanding the accumulation of pollutants and the dynamics of the lower atmosphere. Lidar has been used for tracking the evolution of PBLH by using aerosol backscatter as a tracer, assuming aerosol is generally well-mixed in the PBL; however, the validity of this assumption actually varies with atmospheric stability. This is demonstrated here for stable boundary layers (SBL), neutral boundary layers (NBL), and convective boundary layers (CBL) using an 8-year dataset of micropulse lidar (MPL) and radiosonde (RS) measurements at the ARM Southern Great Plains, and MPL at the GSFC site. Due to weak thermal convection and complex aerosol stratification, traditional gradient and wavelet methods can have difficulty capturing the diurnal PBLH variations in the morning and forenoon, as well as under stable conditions generally. A new method is developed that combines lidar-measured aerosol backscatter with a stability dependent model of PBLH temporal variation (DTDS). The latter helps “recalibrate” the PBLH in the presence of a residual aerosol layer that does not change in harmony with PBL diurnal variation. The hybrid method offers significantly improved PBLH detection, with better correlation and smaller biases, under most thermodynamic conditions, especially for SBL and CBL. Relying on the physical process of PBL diurnal development, different schemes are developed for growing, maintenance, and decaying periods. Comprehensive evaluation of this new method shows much better tracking of diurnal PBLH variation and significantly smaller biases under various pollution levels.

Tianning Su↗
NASA NTRS1978

Random motion analysis of flexible satellite structures

A singular perturbation formulation is used to study the responses of a flexible satellite when random measurement errors can occur. The random variables, at different instants of time, are assumed to be uncorrelated. Procedures for obtaining maxima and minima are described, and a variation of the linear method is developed for the formal solution of the two-point boundary-value problems represented by the variational equations. Random and deterministic solutions for the structural position coordinates are studied, and an analytic algorithm for treating the force equation of motion is developed. Since the random system indicated by the variational equation will always be asymptotically unstable, any analysis of stability must be based on the deterministic system.

Huang, T. C.↗
NASA NTRSDate not supplied

Delta Method Application for the Correlation of IR Detector Thermal Parasitic Loads with Statistically Accurate Results

The Linear Etalon Imaging Spectral Array (LEISA) is a cryogenic IR detector that is passively cooled below 110K. LEISA is part of the L’Ralph instrument, which is going on the Lucy mission, with planned flyby’s of the Trojan asteroids around 5.5 AU. One of the primary thermal challenges is to successfully quantify each significant parasitic heat flow value with statistical certainty.Simulation of flight-like thermal environments and correlation of thermal models is an important part of this challenge, which can be critical for thermally sensitive systems such as cryogenic detectors that are passively cooled. It is common practice during TVAC to achieve several thermal balance points and then to correlate several design variables in order to achieve the closest match to all balance points tested. This approach lends itself to high levels of uncertainty in the correlated values due to the cumulative effect of absolute temperature uncertainty and differences between the sensors used for the measurement. Accuracy can be improved with the use of more accurate sensors, by calibrating all sensors with respect to each other, and by using Zero-Q methodologies, all of which add cost and complexity to the test.The Delta Method uses the difference in temperature caused by the changes in the boundary conditions to calculate heat flows and/or thermal resistances across temperature differentials. This eliminates the problem of high uncertainties with absolute temperature measurements, and yields statistically accurate results down to the sensor resolution of 1mK at < 100K, by taking the sensor’s calibration errors out of the equation. During testing, this method relies on parametric variations around an initial balance point, by varying the boundary conditions of the different parasitic heat sources. The Delta Method, as well as its application on the LEISA thermal test, will be discussed in detail.

Daniel G. Bae↗