Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalization”

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 973 records · Page 54

Exploring structural transitions at grain boundaries in Nb using a generalized embedded atom interatomic potential

The advancement in experimental techniques, like the atom probe tomography and high resolution electron microscopy, is fueling interest in studying structural transformations of grain boundaries in metal and alloys to uncover correlations between mechanical properties and solute or impurity segregation to grain boundaries. Atomistic modeling is an important tool that can pinpoint the intricate dynamics of grain boundary phase transitions, but the lack of accurate interatomic potentials needed to simulate the complex dynamics of grain boundary structural transitions and identify different metastable phases has been the bottleneck. To this end, we use niobium as a model body centered cubic (BCC) metal and develop an interatomic potential to study grain boundary phase transitions. The potential for Nb is based on a generalization of the embedded atomic method potential and has sufficient flexibility to learn complex energy landscapes using a small set of training structures. We systematically test and validate the using data from ab initio density functional theory calculations and experiments. Using this potential, we calculate energies of multiple symmetric-tilt grain boundaries spanning a wide range of misorientation angles. Additionally, we explore different metastable structures of the Σ 27(552) [$1\overline{1}0$] grain boundary and use molecular dynamic simulations to study the coexistence of metastable phases and grain boundary transitions at finite temperature.

36 MATERIALS SCIENCE↗

Analysis of phase stability and chemical segregation in the Mo-V alloys using a generalized embedded atom method potential

A new interatomic potential for the Mo-V system is introduced to facilitate the study of phase stability and mechanical properties at lower temperatures. This potential is based on a generalization of the embedded atom method and includes contributions from embedding energy, explicit two- and three-body interactions and nonlocal many-body interaction terms. The parameters of the potential are optimized by using data from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, such as elastic constants, equation of states, phonon dispersion curves, point defect properties and melting temperatures for different compositions. Even though our potential is trained on a small dataset, its accuracy is comparable to available machine learning potentials for Mo and V. Furthermore, our results show that an ordered B2 phase is stable at low temperatures in alloys containing 50% V, but the solid solution phase is stable above 800 K. However, such long-range ordering is not observed in V-rich or Mo-rich alloys. In addition, our results show that V segregates to dislocation cores and grain boundaries.

36 MATERIALS SCIENCE↗

A general-purpose method for Pareto optimal placement of flow rate and concentration sensors in networked systems – With application to wastewater treatment plants

The advent of affordable computing, low-cost sensor hardware, and high-speed and reliable communications have spurred installation of ubiquitous sensors in complex engineered systems. However, ensuring reliable data quality remains a challenge. Exploitation of redundancy among sensor signals can help improving the precision of measured variables, detecting the presence of gross errors, and identifying faulty sensors. The cost of sensor ownership, maintenance efforts in particular, can still be cost-prohibitive however. Maximizing the ability to assess and control data quality while minimizing the cost of ownership thus requires a careful sensor placement. To solve this challenge, in this work we develop a generally applicable method to solve the multi-objective sensor placement problem in systems governed by linear and bilinear balance equations. Importantly, the method computes all Pareto-optimal sensor layouts with conventional computational resources and requires no information about the expected sensor quality.

42 ENGINEERING↗

LAVA 1.0: A general-purpose python toolkit for calculation of material properties with LAMMPS and VASP

Here, we introduce LAVA, a general-purpose python toolkit to provide user-friendly and high throughput calculations for material properties using both LAMMPS and VASP. It contains a set of classes as well as pre-processing and post-processing functions to prepare, execute and extract information from LAMMPS/VASP simulations. An overview of the program structure is provided. The current version contains modules to calculate elastic and mechanical properties such as lattice constant, cohesive energy, cold curve, elastic constants, bulk and shear modulus, volume conserving/non-conserving deformation path, vacancy/interstitial formation energy, surface energy, stacking fault energy, melting point, radial distribution function, and thermal expansion. These functionalities are demonstrated for different interatomic potentials and DFT calculations in Al.

36 MATERIALS SCIENCE↗

Kernel methods for evolution of generalized parton distributions

Generalized parton distributions (GPDs) characterize the 3-dimensional structure of hadrons, combining information about their internal quark and gluon longitudinal momentum distributions and transverse position within the hadron. The dependence of GPDs on the factorization scale Q 2 allows one to connect hard exclusive processes involving GPDs at disparate energy and momentum scales, which is needed in global analyses of experimental data. Here, in this work, we explore how finite element methods can be used to construct fast and differentiable Q 2 evolution codes for GPDs in momentum space, which can be used in a machine learning framework. We show numerical benchmarks of the methods' accuracy, including a comparison to an existing evolution code from PARTONS/APFEL++, and provide a repository where the code can be accessed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

General solution to the Kohn–Luttinger nonconvergence problem

A simple, but general solution is proposed for the Kohn–Luttinger problem, i.e., the nonconvergence of the finite-temperature many-body perturbation theory with its zero-temperature counterpart as temperature is lowered to zero under some circumstances. How this nonconvergence can be avoided by altering the reference wave function is illustrated numerically by using up to the fifth order of the perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗

A general predictive methodology for fuel-mixture properties up to supercritical conditions

A predictive thermodynamic model is utilized for the calculation of fuel properties of oxymethylene dimethyl ethers (OME 3–4 ), surrogates for gasoline, diesel and aviation fuel, as well as alcohol blends with gasoline and diesel. The alcohols used for these blends are methanol, ethanol, propanol, butanol and pentanol; their mixing ratio ranges from 10 to 50% by volume. The model is based on the Perturbed-Chain Statistical Association Fluid Theory (PC-SAFT) equation of state (EoS) and Vapor Liquid Equilibrium (VLE) calculations at constant temperature, density and composition. The model includes the association term, with the assumption of two association sites (2B scheme), to enable the modeling of alcohols. The pure-component parameters are estimated based on the Group Contribution (GC) method of various sources, as well as a parametrization model specifically designed for the case of OME 3–4 . The results of the computational model for the density, vapor pressure and distillation curves at various conditions, including high-pressure, high-temperature (HPHT), are compared to experimental and computational data available in the literature. In the cases where no measurements are available for the surrogates, experimental data for the corresponding target fuel are used, taking into consideration the inherent deviation in properties between real and surrogate fuel. Overall, the results are in good agreement with the data from the literature, with the average deviation not exceeding 12% for temperature (Kelvin) on the distillation curves, 10% for density and 46% for vapor pressure and the general trend being captured successfully. The use of different pure component parameter estimation techniques can further improve the prediction quality in the cases of OME3–4 and the aviation fuel surrogate, especially for the vapor pressure, leading to an average deviation lower than 18%. These results demonstrate the predictive capabilities of the model, which extend to a wide range of fuel types and pressure/temperature conditions. Through this investigation, the present work aims to establish the limits of applicability of this thermodynamic property prediction methodology.

33 ADVANCED PROPULSION SYSTEMS↗

A generalized analytical energy balance model for evaluating agglomeration from a binary collision of wet particles

Agglomeration of wet particles, i.e., particles coated with a thin liquid layer, is a common phenomenon in many processes like fluidized bed combustion of low rank fuels. The availability of an agglomeration model that can evaluate the outcome of a binary collision between wet particles differing in solid particle properties, liquid layer thicknesses, and initial collision (impact) speeds is essential for obtaining a comprehensive understanding on the existing processes experiencing wet particle agglomeration or for a successful development of new processes with high chances of wet particle agglomeration. This study presents a generalized agglomeration model on the basis of energy conservation before and after collision when colliding wet particles may differ in solid particle properties, liquid layer thicknesses, and impact speeds. The model was established based on the approximate values of energy losses that may happen during the collision. It incorporates body forces, solid-solid contacting, liquid capillary, and viscous contributions, as well as the liquid bridge volume effect. Predictions of the new model for collision outcomes of identical wet particles were like those from an analytical energy balance model developed recently by the group for identical wet particles. We also validated the new model by experimental data from literature. The results of a collision direction analysis indicated that the direction often has a minimal effect on the collision outcome in many practical scenarios. The results of Monte Carlo uncertainty analyses with the new model revealed that proper estimations of impact speed, under capillary limiting conditions, and thickness of coating layers and asperity heights, under viscous limiting conditions, are critical for the realistic prediction of collision outcomes at impact speeds close to critical impact speed, i.e., the minimum particle speed required for the particles to rebound.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A time-varying vulnerability index for COVID-19 in New Mexico, USA using generalized propensity scores

The coronavirus disease (COVID-19) pandemic has highlighted systemic inequities in the United States and resulted in a larger burden of negative social outcomes for marginalized communities. New Mexico, a state in the southwestern US, has a unique population with a large racial minority population and a high rate of poverty that may make communities more vulnerable to negative social outcomes from COVID-19. To identify which communities may be at the highest relative risk, we created a county-level vulnerability index. After the first COVID-19 case was reported in New Mexico on March 11, 2020, we fit a generalized propensity score model that incorporates sociodemographic factors to predict county-level viral exposure and thus, the generic risk to negative social outcomes such as unemployment or mental health impacts. We used four static sociodemographic covariates important for the state of New Mexico—population, poverty, household size, and minority population—and weekly cumulative case counts to iteratively run our model each week and normalize the exposure score to create a time-varying vulnerability index. We found the relative vulnerability between counties varied in the first eight weeks from the initial COVID-19 case before stabilizing. This framework for creating a location-specific vulnerability index in response to an ongoing disaster may be used as a quick, deployable metric to inform health policy decisions such as allocating state resources to the county level.

59 BASIC BIOLOGICAL SCIENCES↗

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗

A cavitation and dynamic void growth model for a general class of strain-softening amorphous materials

In this work, we derive a homogenized framework for cavitation and dynamic void growth in a general class of strain-softening materials with particular emphasis on amorphous materials. Analytic solutions are derived by considering an assemblage of spherical shells subject to dynamic hydrostatic tensile pressure. The isochoric response of the matrix material is captured by a free volume theory based amorphous viscoplasticity model. The framework accounts for the coupled rate-dependent, strain-softening behavior typically exhibited in amorphous materials. Validation of the homogenized theory is carried out against numerous direct numerical simulations with excellent agreement. Furthermore, we find that a reduction in cavitation strength and significant acceleration of dynamic void growth rates is induced as a consequence of strain-softening in the amorphous matrix material. Lastly, the model is utilized to interpret some non-intuitive experimental observations of amorphous materials subject to shock compression and subsequent spall failure.

36 MATERIALS SCIENCE↗

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗

Steady-state fuel performance analyses for the preliminary fuel concept of general atomics fast modular reactor

Here this manuscript presents the fuel performance analysis results of the General Atomics Fast Modular Reactor (FMR) based on an axi-symmetric (2D-RZ) geometry. Three fuel performance model sets that fit the FMR fuel specifications best, i.e., a BISON baseline model set, a BISON diffusion enhancement model set, and a BISON-FASTGRASS model set, were identified and evaluated against a series of relevant experimental cases featuring high burnup and low irradiation temperature conditions. The three BISON-based model sets were then utilized to conduct a comprehensive fuel performance analysis of the FMR fuel under normal operation including the shutdown/restarting periods for refueling. The evaluation of the fuel performance parameters, represented by temperature, internal pressure, stress, and strain, shows that the FMR fuel maintains its thermal and mechanical integrity during normal operation. Technology gaps and limitations are also discussed to guide future efforts for extending the performance analysis to transient scenarios as well as improving the fuel performance evaluation through experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantitative measurement of pebble fuel anisotropy using generalized ellipsometry

The renewed development of high-temperature gas-cooled reactors has created a need to modernize pebble fuel quality control methods. One important specification in these methods is the coefficient of thermal expansion (CTE) anisotropy, which reflects graphite grain alignment and affects fuel performance. This study investigates the feasibility of using two-modulator generalized ellipsometry microscopy (2-MGEM) as a potential surrogate metric for CTE anisotropy. Graphite matrix compacts with controlled levels of anisotropy were fabricated by varying pre-pressures during forming and were subsequently characterized using both optical dilatometry and 2-MGEM. Measured responses were compared to an isotropic reference pebble produced using an industrial manufacturing process to provide a practical baseline. The resulting measurements revealed a strong quantitative correlation between the CTE anisotropy ratio and the corresponding diattenuation ratio. However, propagated measurement uncertainties in the optical method were found to limit definitive pass/fail discrimination, highlighting the need for improved precision before qualification as a quality control tool. Still, the results establish a quantitative relationship between optical and CTE anisotropy and support further development of 2-MGEM as a potential surrogate method for evaluating graphite matrix anisotropy in pebble fuel manufacturing.

Raftery, Alicia [ORNL] (ORCID:0000000259076427)↗

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗

Suzuki type estimates for exponentiated sums and generalized Lie-Trotter formulas in JB-algebras

Lie-Trotter-Suzuki product formulas are ubiquitous in quantum mechanics, computing, and simulations. Approximating exponentiated sums with such formulas are investigated in the JB-algebraic setting. We show that the Suzuki type approximation for exponentiated sums holds in JB-algebras, we give explicit estimation formulas, and we deduce three generalizations of Lie-Trotter formulas for arbitrary number elements in such algebras. In conclusion, we also extended the Lie-Trotter formulas in a Jordan Banach algebra from three elements to an arbitrary number of elements.

97 MATHEMATICS AND COMPUTING↗

Directional amorphization of covalently-bonded solids: A generalized deformation mechanism in extreme loading

Shock compression subjects materials to a unique regime of high quasi-hydrostatic pressure and coupled shear stresses for durations on the order of 1–10 nanoseconds for laser-driven loading of samples. There is, additionally, an attendant temperature increase due to the shock and the mechanisms of plastic deformation in metals whereby dislocations, twins, and phase transitions nucleate and propagate at velocities near the sound speed. Covalently bonded materials have, by virtue of the directionality of their bonds, great difficulty in responding by conventional plastic deformation to this extreme regime of shock compression. Here, we propose that the shear from shock compression induces amorphization, as observed in Si, Ge, B4C, SiC, and olivine ((Mg, Fe) 2 SO 4 ) and that this is a general deformation mechanism in a broad class of covalently bonded materials. The crystalline structure transforms to amorphous along regions of maximum shear stress, forming nanoscale bands, and thereby relaxing the shear component of the imposed shock stress. This process is usually preceded by the emission and propagation of a critical concentration of dislocations.

36 MATERIALS SCIENCE↗

Deep Kronecker neural networks: A general framework for neural networks with adaptive activation functions

Here we propose a new type of neural networks, Kronecker neural networks (KNNs), that form a general framework for neural networks with adaptive activation functions. KNNs employ the Kronecker product, which provides an efficient way of constructing a very wide network while keeping the number of parameters low. Our theoretical analysis reveals that under suitable conditions, KNNs induce a faster decay of the loss than that by the feed-forward networks. This is also empirically verified through a set of computational examples. Furthermore, under certain technical assumptions, we establish global convergence of gradient descent for KNNs. As a specific case, we propose the Rowdy activation function that is designed to get rid of any saturation region by injecting sinusoidal fluctuations, which include trainable parameters. The proposed Rowdy activation function can be employed in any neural network architecture like feed-forward neural networks, Recurrent neural networks, Convolutional neural networks etc. The effectiveness of KNNs with Rowdy activation is demonstrated through various computational experiments including function approximation using feed-forward neural networks, solution inference of partial differential equations using the physics-informed neural networks, and standard deep learning benchmark problems using convolutional and fully-connected neural networks.

97 MATHEMATICS AND COMPUTING↗