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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Minimizing thickness variation in monolithic U-10Mo fuel foil and Zr interlayer during hot rolling: A microstructure-based finite element method analysis

Low-enriched uranium alloyed with 10 wt. % molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium fuel for the United States’ high performance research reactors. The monolithic U-10Mo fuel plate consists of a metallic U-10Mo fuel foil with a 25 µm Zr interlayer and a relatively thick cladding of aluminum alloy 6061. The Zr interlayer is typically applied during the hot co-rolling process, and this process dictates the uniformity of the Zr interlayer. Thickness variation observed in the U-10Mo and Zr interlayer has been attributed to several sources: the initial grain size of the U-10Mo castings, can materials, rolling temperature, inhomogeneous molybdenum content, and porosity in the cast U-10Mo. This thickness variation limits the ability to meet the dimensional specification; thus, a better understanding of the factors causing the nonuniform thickness is needed. In this work, we used a novel, microstructure-based finite element method to model the hot rolling process to address these concerns. Grain microstructures in U-10Mo were tessellated and explicitly considered in the finite element model. Each grain was assigned a random material property to mimic the grain strength variations induced by different grain orientations. Simulations were performed using six steel can thicknesses, four grain sizes, and with or without a Zr interlayer to investigate the influences of those variables on the thickness nonuniformity. The simulation results showed that a thinner steel can and finer U-10Mo grain size reduce thickness variations in both the U-10Mo fuel foil and Zr interlayer. The direct findings from the simulations and analysis can be used to optimize the hot rolling schedule, reduce fabrication defects, and meet the dimensional specifications. The proposed microstructure-based finite element model can be also coupled with experimental microstructure characterization data, images, and models to simulate multi-pass hot rolling.

36 MATERIALS SCIENCE↗
NASA NTRS1979

Interfacial geometry and D-variation effects in two-phase systems

Numerical solutions of the governing diffusion equation for two-phase concentration dependent diffusion coefficients are examined. Solutions were also calculated for planar, cylindrical, and spherical geometries to compare the effect of interface geometries with those caused by concentration-dependent diffusion coefficients, and two methods of averaging D were considered to determine the best averaging method for different types of D-variations. The effects of interface-location criteria on mass conservation and convergence of interface location, diffusion coefficient variation in the alpha and beta-phases of a two-phase binary alloy system, effect of D(alpha) variation in a cylindrical couple on beta-phase thickness, and geometry and D-variation effects on the degree of homogenization were determined. It is concluded that typical D(alpha)-variations can have a greater influence on the kinetics of interdiffusion than the geometry.

Tenney, D. R.↗
NASA NTRS1977

Spatially resolved excitation temperature measurements in a hypersonic flow using the hook method

Following an analysis of the hook equations for variation of refractivity along a spectrograph slit a series of experiments is reported, in which CO2 was used as the test gas, and the flow variations were spatially resolved. The hooked spectra were used to measure the ground state number densities of the impurity atom of CrI present in the flow. The excitation temperature of the CrI 1-eV level was obtained from the separation of the hooks of the CrI 425.4-nm transition resonance line and simultaneously the 434.4-nm CrI eV lower level line, from the shock to the body of the hypersonic shock tunnel flow. This excitation temperature is compared with theoretical calculations of the CO2 gas temperature. Some further applications of the hook type methods to gas diagnostics are discussed.

Sandeman, R. J.↗
DOE OSTI2024

The collisional particle-in-cell method for the Vlasov–Maxwell–Landau equations

We introduce an extension of the particle-in-cell method that captures the Landau collisional effects in the Vlasov–Maxwell–Landau equations. The method arises from a regularisation of the variational formulation of the Landau equation, leading to a discretisation of the collision operator that conserves mass, charge, momentum and energy, while increasing the (regularised) entropy. The collisional effects appear as a fully deterministic effective force, thus the method does not require any transport–collision splitting. The scheme can be used in arbitrary dimension, and for a general interaction, including the Coulomb case. We validate the scheme on scenarios such as the Landau damping, the two-stream instability and the Weibel instability, demonstrating its effectiveness in the numerical simulation of plasma.

Bailo, Rafael (ORCID:0000000180183799)↗
DOE OSTI2024

Uncertainty quantification of graph convolution neural network models of evolving processes

The application of neural network models to scientific machine learning tasks has proliferated in recent years. In particular, neural networks have proved to be adept at modeling processes with spatial–temporal complexity. Nevertheless, these highly parameterized models have garnered skepticism in their ability to produce outputs with quantified error bounds over the regimes of interest. Hence there is a need to find uncertainty quantification methods that are suitable for neural networks. In this work we present comparisons of the parametric uncertainty quantification of neural networks modeling complex spatial–temporal processes with Hamiltonian Monte Carlo and Stein variational gradient descent and its projected variant. Specifically we apply these methods to graph convolutional neural network models of evolving systems modeled with recurrent neural network and neural ordinary differential equations architectures. We show that Stein variational inference is a viable alternative to Monte Carlo methods with some clear advantages for complex neural network models. For our exemplars, Stein variational interference gave similar pushed forward uncertainty profiles through time compared to Hamiltonian Monte Carlo, albeit with generally more generous variance. As a result, projected Stein variational gradient descent also produced similar uncertainty profiles to the non-projected counterpart, but large reductions in the active weight space were confounded by the stability of the neural network predictions and the convoluted likelihood landscape.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Ice Accretion Roughness Variations on a Hybrid CRM65-Midspan Wing Model

Ice accretion roughness measurements were performed in the Icing Research Tunnel (IRT) at NASA Glenn Research Center for the Hybrid CRM65-Midspan model in a range of icing conditions. The Hybrid CRM65-Midspan model was chosen for this investigation because 1) the model exhibits high sweep relative to models previously explored in the roughness investigations, 2) the model has leading edge characteristics similar to wing shapes currently used in mid-size commercial airliners, and 3) the sweep and thickness ratios relate better to hybrid lifting body designs for N+2 and N+3 vehicles than other models available. The investigation consisted of multiple sets of tests which focused on 1) 0-angle of attack cases replicating the conditions employed by Anderson et al. (1998) using both Appendix C and SLD cloud conditions, 2) cases based on the “Max Scallop” case by Broeren et al. (2016) and a “High Temperature” case with cloud properties similar to the “Max Scallop” case. Additional tests were performed 1) based on the “Max Scallop” case with variations in freestream static temperature and 2) using test section speeds near 10,000-hold flight conditions. The point clouds were characterized using the approach of McClain and Kreeger (2013) for the ice roughness variations and using the approach of McClain (2016) for the mean ice thickness variations. The resulting roughness and mean thickness variations generally follow the temporal scaling previously identified using on airfoil models without sweep, but the collapse of the time progression profiles is not as tight as found for past measurements on models without sweep. LEWICE and modified panel-method predictions were used explore spatial roughness variations and to compare to the roughness correlations developed by McClain et al. (2021) for the “Max Scallop” cases.

Icing↗
NASA NTRS1951

An Instrument Employing a Coronal Discharge for the Determination of Droplet-Size Distribution in Clouds

A flight instrument that uses electric means for measuring the droplet-size distribution in above-freezing clouds has been devised and given preliminary evaluation in flight. An electric charge is placed on the droplets and they are separated aerodynamically according to their mass. Because the charge placed on the droplets is a. function of the droplet size, the size spectrum can 'be determined by measurement of the charge deposited on cylinders of several different sizes placed to intercept the charged droplets. An expression for the rate of charge acquisition by a water droplet in a field of coronal discharge is derived. The results obtained in flight with an instrument based on the method described indicate that continuous records of droplet-size spectrum variations in clouds can be obtained. The experimental instrument was used to evaluate the method and was not refined to the extent necessary for obtaining conclusive meteorological data. The desirable features of an instrument based on the method described are (i) The instrument can be used in clouds with temperatures above freezing; (2) the size and the shape of the cylinders do not change during the exposure time; (3) the readings are instantaneous and continuous; (4) the available sensitivity permits the study of variations in cloud structures of less than 200 feet in extent.

Brun, Rinaldo J.↗
DOE OSTI2025

Joint Modeling of Wind Speed and Wind Direction Through a Conditional Approach

Atmospheric near surface wind speed and wind direction play an important role in many applications, ranging from air quality modeling, building design, wind turbine placement to climate change research. It is therefore crucial to accurately estimate the joint probability distribution of wind speed and direction. In this work, we develop a conditional approach to model these two variables, where the joint distribution is decomposed into the product of the marginal distribution of wind direction and the conditional distribution of wind speed given wind direction. To accommodate the circular nature of wind direction, a von Mises mixture model is used; the conditional wind speed distribution is modeled as a directional dependent Weibull distribution via a two-stage estimation procedure, consisting of a directional binned Weibull parameter estimation, followed by a harmonic regression to estimate the dependence of the Weibull parameters on wind direction. A Monte Carlo simulation study indicates that our method outperforms two other approaches in estimation efficiency: one that utilizes periodic spline quantile regression and another that generates data from the commonly used Abe-Ley distribution for cylindrical data. We illustrate our method by using the output from a regional climate model to investigate how the joint distribution of wind speed and direction may change under some future climate scenarios. Our method indicates significant changes in the variation of wind speed with respect to some directions.

17 WIND ENERGY↗
DOE OSTI2021

Correlation of Early-Stage Growth Process Conditions with Dislocation Evolution in MOCVD-Based GaP/Si Heteroepitaxy

To identify the complex relationships between early-stage growth processes and the resultant defect microstructure in GaP/Si heteroepitaxy, a holistic study of several key metal-organic chemical vapor deposition (MOCVD) parameters was conducted, focusing on Si surface preparation and GaP atomic layer epitaxy (ALE) based nucleation processes. Here, crystalline defects related to the lattice mismatch and/or interfacial heterovalency, namely misfit dislocations (MD), threading dislocations (TD), and stacking fault pyramids (SFP), were quantitatively characterized via electron channeling contrast imaging (ECCI) and correlated against the different process variations. Choice of Si surface preparation method between the two examined (dilute SiH4 annealing versus Si2H6 based homoepitaxy) had little impact on resultant GaP film morphology and defect content, whereas differing GaP ALE nucleation conditions produced much more substantial changes. In particular, the initial precursor species (tert-butylphosphine versus triethygallium) and ALE cycle purge times both yielded significant influence over threading dislocation densities (TDD) in thin (100 nm), post-critical thickness GaP/Si films, with TDD spanning two orders of magnitude, from 6.7×10 7 cm -2 to 7.1×10 5 cm -2 , depending on the specific process conditions employed. SFP densities were also found to follow a similar trend, ranging from 2.0×10 7 cm -2 to 1.8×10 5 cm -2 , but with no apparent causal relationship between SFP density and TDD. To help explain the dramatic differences observed, detailed, large-area MD network characterization was used to provide statistically-relevant quantitative analyses of the critical dislocation dynamics (introduction rates and glide velocities) associated with the different process variants. These extracted values are then correlated against the ALE process variants to provide insight into the potential mechanistic roles of the different growth processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory

The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density-matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies, relative to reference values from configuration-interaction-based calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Dehydrogenated Polyethylene from Discarded Plastics as a Synthon for Functional Polyolefins

Many valuable specialty chemicals and drug candidates rely on synthesizing reactive intermediates to generate the final product(s). An approach that leverages C-H activation chemistry to achieve chemically active reactive handles for polymer diversification is desirable. Further, this report describes the synthesis and functionalization of dehydrogenated high-density polyethylene (HDPE) from discarded post-consumer plastic via acceptorless Ir-pincer-catalyzed non-oxidative dehydrogenation followed by acid-catalyzed hydroamination and hydroalkoxylation. The final products show significant material property changes with tunability and scalability through variation of the olefin content. This method highlights the potential value of dehydrogenated intermediates as platforms to enable carbon circularity and polyolefin end-of-life management strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A nanosecond pulsed discharge circuit model for engine applications

Non-equilibrium plasma discharges in spark gaps have been an increasingly studied method for alleviating cycle to cycle variation in lean and dilute combustion environments. However, ignition models that account for streamer propagation, cathode fall, and transmission line amplification over nanosecond time scales have so far not been developed. The present study develops such a model, with emphasis on the energy delivered from circuit to cylinder. Here, key pieces of the relevant physics and chemistry are summarized, simplified, and systematically coupled to one another. The set of parameters is limited to a handful of key observables and modeled using Modelica. Results show non-trivial behavior in the energy delivery characteristics of such discharges with important implications for ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Three-dimensional modelling of a self-sustained atmospheric pressure glow discharge

The atmospheric pressure glow discharge (APGD) is a relatively simple and versatile plasma source used in diverse applications. Stable APGD operation at high currents, generally a challenge due to instabilities leading to glow-to-arc transition, has been demonstrated using actively-controlled cathodic cooling. This article presents the computational modelling and simulation of a self-sustained direct-current APGD in helium within a 10 mm pin-to-plate inter-electrode gap for currents ranging from 4 to 40 mA. The APGD model is comprised of the conservation equations for total mass, chemical species, momentum, thermal energy of heavy-species and of free electrons, and electric charge. The model equations are discretized using a nonlinear variational multi-scale finite element method that has demonstrated superior accuracy in other plasma flow problems, on a temporal and three-dimensional computational domain suitable to unveil the potential occurrence of instabilities. Modelling results show good agreement with experimental measurements of voltage drop and the same trend but higher values of temperature. The higher temperatures obtained by the simulations appear to be due to the absence of a near-cathode heat dissipation model. Here, the results also reveal that the distribution of electron density and of the ratio of atomic helium ions to total ions transitions from monotonically increasing away from the cathode to presenting a minimum near the centre of the gap with increasing current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantum error mitigation by hidden inverses protocol in superconducting quantum devices

We present a method to improve the convergence of variational algorithms based on hidden inverses (HIs) to mitigate coherent errors. In the context of error mitigation, this means replacing the hardware implementation of certain Hermitian gates with their inverses. Doing so results in noise cancellation and a more resilient quantum circuit. This approach improves performance in a variety of two-qubit error models where the noise operator also inverts with the gate inversion. We apply the mitigation scheme on superconducting quantum processors running the variational quantum eigensolver (VQE) algorithm to find the H 2 ground-state energy. When implemented on superconducting hardware we find that the mitigation scheme effectively reduces the energy fluctuations in the parameter learning path in VQE, reducing the number of iterations for a converged value. We also provide a detailed numerical simulation of VQE performance under different noise models and explore how HIs and randomized compiling affect the underlying loss landscape of the learning problem. These simulations help explain our experimental hardware outcomes, helping to connect lower-level gate performance to application-specific behavior in contrast to metrics like fidelity which often do not provide an intuitive insight into observed high level performance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Human papillomavirus–associated anal squamous cell carcinoma: sociodemographic, geographic, and county-level economic trends in incidence rates—United States, 2001-2019

Abstract Background Incidence of anal squamous cell carcinoma is increasing, but vaccination against human papillomavirus (HPV) and removal of precancerous anal lesions could prevent new cases. The overall HPV-associated cancer incidence is reported to be higher in rural populations and in counties with lower economic status. We assessed these differences specifically for HPV-associated anal squamous cell carcinoma and described the geographic, county-level economic, and sociodemographic variations in incidence rates and trends. Methods We analyzed data from the US Cancer Statistics to assess age-standardized incidence rates of HPV-associated squamous cell carcinomas among adults aged 18 years and older from 2001 to 2019. We calculated rate ratios and 95% confidence intervals to examine differences in incidence rates. We also quantified changes in incidence rates over time using joinpoint regression. Results From 2001 to 2019, 72 421 new cases of HPV-associated anal squamous cell carcinoma were diagnosed among women (2.8 per 100 000) and 37 147 among men (1.7 per 100 000). Age-standardized incidence rates were higher in the South compared with other census regions and in counties ranked in the bottom 25% and 25%-75% economically than in the top 25%. The overall incidence rate increased in women but remained stable in men during 2009-2019. Incidence rates increased in adults aged 50 years and older but decreased among those aged 40-44 years from 2001 to 2019 in women and from 2007 to 2019 in men. Conclusions There were inequities in HPV-associated anal squamous cell carcinoma incidence by geographic and county-level economic characteristics. Failure to improve vaccine and treatment equity may widen existing disparities.

Oncology↗
DOE OSTI2021

Variational study of the Kitaev-Heisenberg-Gamma model

Here we compute the low-energy excitation spectrum and the dynamical spin structure factor of the Kitaev-Heisenberg-Gamma model through a variational approach based on the exact fractionalized excitations of the pure Kitaev honeycomb model. This novel approach reveals the physical reason for the asymmetric stability of the Kitaev spin liquid phases around the ferromagnetic and antiferromagnetic Kitaev limits. Moreover, we demonstrate that the fractionalized excitations form bound states in specific regions of each Kitaev spin liquid phase and that certain phase transitions induced by Heisenberg and Gamma interactions are driven by the condensation of such a bound state. Remarkably, this bound state appears as a sharp mode in the dynamical spin structure factor, while its condensation patterns at the appropriate phase transitions provide a simple explanation for the magnetically ordered phases surrounding each Kitaev spin liquid phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fixed-phase diffusion Monte Carlo study of activation gap and skyrmion excitations of a ν = 1 system in the presence of charged impurities

The discrepancy between the theoretically calculated and experimentally measured activation gaps in quantum Hall effect has long been a puzzle. We revisit this issue in the context of the v = 1 quantum Hall state, while also incorporating the skyrmion physics. We find that the finite width and the Landau level mixing effects are not sufficient to explain the observed activation gap. We further show that the presence of charged impurities located adjacent to the quantum well can cause a significant reduction in the activation gap, while also causing a suppression of the skyrmion size.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Phase transitions of the ferroelectric (ND 4 ) 2 FeCl 5 ·D 2 O under magnetic field

Due to the strong coupling between magnetism and ferroelectricity, (ND 4 ) 2 FeCl 5 ·D 2 O exhibits several intriguing magnetic and electric phases. In this Letter, we include higher-order onsite anisotropic spin interactions to explain the ferroelectric phase transitions of (ND 4 ) 2 FeCl 5 ·D 2 O in a magnetic field and to produce the large weights of high-order harmonic components in the cycloid structure that are observed from neutron diffraction experiments. Moreover, we predict a new ferroelectric phase sandwiched between the FE II and FE III phases in a magnetic field. Furthermore, by emphasizing the importance of the higher-order spin anisotropic interactions, our work provides a framework to understand multiferroic materials with rich phase diagrams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗