Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Methodology for Sensitivity Analysis of Homogenized Cross-Sections to Instantaneous and Historical Lattice Conditions with Application to AP1000® PWR Lattice

In the two-step method for nuclear reactor simulation, lattice physics calculations are performed to compute homogenized cross-sections for a variety of burnups and lattice configurations. A nodal code is then used to perform full-core analysis using the pre-calculated homogenized cross-sections. One source of uncertainty introduced in this method is that the lattice configuration or depletion conditions typically do not match a pre-calculated one from the lattice physics simulations. Therefore, some interpolation model must be used to estimate the homogenized cross-sections in the nodal code. This current study provides a methodology for sensitivity analysis to quantify the impact of state variables on the homogenized cross-sections. This methodology also allows for analyses of the historical effect that the state variables have on homogenized cross-sections. An application of this methodology on a lattice for the Westinghouse AP1000® reactor is presented where coolant density, fuel temperature, soluble boron concentration, and control rod insertion are the state variables of interest. The effects of considering the instantaneous values of the state variables, historical values of the state variables, and burnup-averaged values of the state variables are analyzed. Using these methods, it was found that a linear model that only considers the instantaneous and burnup-averaged values of state variables can fail to capture some variations in the homogenized cross-sections.

Price, Dean (ORCID:0000000309990111)↗
NASA NTRS2011

A Groundmass Composition for EET 79001A Using a Novel Microprobe Technique for Estimating Bulk Compositions. Lithology A as an Impact Melt?

Petrologic investigation of the shergottites has been hampered by the fact that most of these meteorites are partial cumulates. Two lines of inquiry have been used to evaluate the compositions of parental liquids: (i) perform melting experiments at different pressures and temperatures until the compositions of cumulate crystal cores are reproduced [e.g., 1]; and (ii) use point-counting techniques to reconstruct the compositions of intercumulus liquids [e.g., 2]. The second of these methods is hampered by the approximate nature of the technique. In effect, element maps are used to construct mineral modes; and average mineral compositions are then converted into bulk compositions. This method works well when the mineral phases are homogeneous [3]. However, when minerals are zoned, with narrow rims contributing disproportionately to the mineral volume, this method becomes problematic. Decisions need to be made about the average composition of the various zones within crystals. And, further, the proportions of those zones also need to be defined. We have developed a new microprobe technique to see whether the point-count method of determining intercumulus liquid composition is realistic. In our technique, the approximating decisions of earlier methods are unnecessary because each pixel of our x-ray maps is turned into a complete eleven-element quantitative analysis. The success or failure of our technique can then be determined by experimentation. As discussed earlier, experiments on our point-count composition can then be used to see whether experimental liquidus phases successfully reproduce natural mineral compositions. Regardless of our ultimate outcome in retrieving shergottite parent liquids, we believe our pixel-bypixel analysis technique represents a giant step forward in documenting thin-section modes and compositions. For a third time, we have analyzed the groundmass composition of EET 79001, 68 [Eg]. The first estimate of Eg was made by [4] and later modified by [5], to take phase diagram considerations into account. The Eg composition of [4] was too olivine normative to be the true Eg composition, because the ,68 groundmass contains no forsteritic olivine. A later mapping by [2] basically reconfirmed the modifications of [5]. However, even the modified composition of [5] has olivine on the liquidus for ~50 C before low-Ca pyroxene appears [6].

Jones, John H.↗
NASA NTRS2009

Cast Aluminum Alloys for High Temperature Applications Using Nanoparticles Al2O3 and Al3-X Compounds (X = Ti, V, Zr)

In this paper, the effect of nanoparticles Al2O3 and Al3-X compounds (X = Ti, V, Zr) on the improvement of mechanical properties of aluminum alloys for elevated temperature applications is presented. These nanoparticles were selected based on their low cost, chemical stability and low diffusions rates in aluminum at high temperatures. The strengthening mechanism at high temperature for aluminum alloy is based on the mechanical blocking of dislocation movements by these nanoparticles. For Al2O3 nanoparticles, the test samples were prepared from special Al2O3 preforms, which were produced using ceramic injection molding process and then pressure infiltrated by molten aluminum. In another method, Al2O3 nanoparticles can also be homogeneously mixed with fine aluminum powder and consolidated into test samples through hot pressing and sintering. With the Al3-X nanoparticles, the test samples are produced as precipitates from in-situ reactions with molten aluminum using conventional permanent mold or die casting techniques. It is found that cast aluminum alloy using nanoparticles Al3-X is the most cost effective method to produce high strength aluminum alloys for high temperature applications in comparison to nanoparticles Al2O3. Furthermore, significant mechanical properties retention in high temperature environment could be achieved with Al3-X nanoparticles, resulting in tensile strength of nearly 3 times higher than most 300- series conventional cast aluminum alloys tested at 600 F.

Lee, Jonathan A.↗
DOE OSTI2022

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The basic objective of this research is to extend the capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. In the efforts related to LES, we were primarily involved with assessing the performance of the various modern methods based on the Probability Density Function (PDF) methods for providing closures for treating the subgrid fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we concentrated on understanding some of the relevant physics of compressible reacting flows by means of statistical analysis of the data generated by DNS of such flows. In the research conducted in the second year of this program, our efforts focused on the modeling of homogeneous compressible turbulent flows by PDF methods, and on DNS of non-equilibrium reacting high speed mixing layers. Some preliminary work is also in progress on PDF modeling of shear flows, and also on LES of such flows.

Givi, P.↗
NASA NTRS1992

LAI inversion from optical reflectance using a neural network trained with a multiple scattering model

The inversion of the leaf area index (LAI) canopy parameter from optical spectral reflectance measurements is obtained using a backpropagation artificial neural network trained using input-output pairs generated by a multiple scattering reflectance model. The problem of LAI estimation over sparse canopies (LAI < 1.0) with varying soil reflectance backgrounds is particularly difficult. Standard multiple regression methods applied to canopies within a single homogeneous soil type yield good results but perform unacceptably when applied across soil boundaries, resulting in absolute percentage errors of >1000 percent for low LAI. Minimization methods applied to merit functions constructed from differences between measured reflectances and predicted reflectances using multiple-scattering models are unacceptably sensitive to a good initial guess for the desired parameter. In contrast, the neural network reported generally yields absolute percentage errors of <30 percent when weighting coefficients trained on one soil type were applied to predicted canopy reflectance at a different soil background.

Smith, James A.↗
NASA NTRSDate not supplied

Active Thermography Based on Tensor Rank Decomposition

Principal Component Thermography applies Singular Value Decomposition (SVD) to post-process data that are derived from active thermographic inspections. SVD provides useful compression of the data and allows for better understanding of substructure and indications of potential damage. In the standard approach, SVD is applied to a certain reshaping of a three-dimensional data stack into a two-dimensional array. This work applies the CANDECOMP-PARAFAC (CP) tensor rank decomposition directly to the three-dimensional data to avoid the initial reshaping step in order to begin to develop an inspection method that can more accurately detect defects in non-homogeneous and anisotropic materials. Tests against simulated data that compare the CP decomposition method with traditional Principal Component Thermography based on SVD are described. Finally, the method of Proper Generalized Decomposition (PGD) is used to derive the CP decomposition, and its performance against other algorithms is also discussed.

Thermography↗
NASA NTRS1974

Buckling of cylindrical shells with random imperfections

Buckling of long cylindrical shells with random imperfections subjected to axial loads is treated using the method of truncated hierarchy. A system of homogeneous variational equations is set up to examine the existence of bifurcation in the neighborhood of the equilibrium state. These equations are linear and involve stationary random coefficients. The truncated-hierarchy method is applied, and characteristic equations are obtained. Various exponential cosine correlation functions associated with asymmetric imperfections are examined numerically. Qualitatively, the results obtained are as anticipated.

Fersht, R. S.↗
DOE OSTI2021

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Lignin-based polymers with enhanced melt extrusion ability

A solid polymer blend material comprising: (i) a lignin-acrylonitrile component containing a homogeneous blend of a lignin component and an acrylonitrile-containing rubber component; and (ii) a styrene-containing thermoplastic component that is non-elastomeric; wherein components (i) and (ii) are homogeneously dispersed in the polymer blend material. Methods for producing the blend material are also described. Methods for producing objects made of the blend material by melt extrusion are also described, comprising: (a) melt blending components (i) and (ii) to form a polymer blend in which components (i) and (ii) are homogeneously blended, wherein the polymer blend exhibits a melt viscosity of no more than 2000 Pa·s at a shear rate of 100-1000 s −1 and when heated to a temperature of no more than 240° C.; and (b) forming an object made of said polymer blend material.

Naskar, Amit K.↗
DOE OSTI2023

Lignin-based polymers with enhanced melt extrusion ability

A solid polymer blend material comprising: (i) a lignin-acrylonitrile component containing a homogeneous blend of a lignin component and an acrylonitrile-containing rubber component; and (ii) a styrene-containing thermoplastic component that is non-elastomeric; wherein components (i) and (ii) are homogeneously dispersed in the polymer blend material. Methods for producing the blend material are also described. Methods for producing objects made of the blend material by melt extrusion are also described, comprising: (a) melt blending components (i) and (ii) to form a polymer blend in which components (i) and (ii) are homogeneously blended, wherein the polymer blend exhibits a melt viscosity of no more than 2000 Pa·s at a shear rate of 100-1000 s −1 and when heated to a temperature of no more than 240° C.; and (b) forming an object made of said polymer blend material.

Naskar, Amit K.↗
NASA NTRS1992

Improved inhomogeneous finite elements for fabric reinforced composite mechanics analysis

There is a need to do routine stress/failure analysis of fabric reinforced composite microstructures to provide additional confidence in critical applications and guide materials development. Conventional methods of 3-D stress analysis are time consuming to set up, run and interpret. A need exists for simpler methods of modeling these structures and analyzing the models. The principal difficulty is the discrete element mesh generation problem. Inhomogeneous finite elements are worth investigating for application to these problems because they eliminate the mesh generation problem. However, there are penalties associated with these elements. Their convergence rates can be slow compared to homogeneous elements. Also, there is no accepted method for obtaining detailed stresses in the constituent materials of each element. This paper shows that the convergence rate can be significantly improved by a simple device which substitutes homogeneous elements for the inhomogeneous ones. The device is shown to work well in simple one and two dimensional problems. However, demonstration of the application to more complex two and three dimensional problems remains to be done. Work is also progressing toward more realistic fabric microstructural geometries.

Foye, R. L.↗
DOE OSTI2023

Computational design of metamaterials with self contact

Inverse homogenization in combination with contact modeling, topology optimization and shape optimization is used to design metamaterials with optimized macroscopic response. The homogenization assumes length scale separation which allows the non-linear macroscopic behavior to be obtained by analyzing a single unit cell in a lattice structure. Self contact in the unit cell, which is modeled using a third medium contact method, is leveraged to obtain a complex homogenized response. The inverse homogenization problem is initially formulated as a topology optimization problem, where the macroscopic stress–strain behavior is tuned to our liking. However, it is well known that boundary phenomena are difficult to model in topology optimization and that interface modeling is crucial to accurately analyze contact. For that reason, the boundary representation of the topology optimized design is extracted and used as initial design in a subsequent shape optimization. The behaviors of our designs are verified by performing rigorous post-processing analyzes using conforming meshes and conventional contact formulations.

42 ENGINEERING↗
DOE OSTI2024

Mechanical and Electrochemical Implications of Drying Temperature on Lithium‐Ion Battery Electrodes

Abstract Lithium‐ion battery (LIB) electrodes are typically produced with n‐methyl‐2‐pyrrolidone, a toxic solvent that is a known carcinogen and reproductive hazard. Accordingly, aqueous processing has been an expanding area of research interest in the field of LIB manufacturing. Although aqueous processing has been widely successful in anode processing, serious challenges remain in processing the cathode. In this work, the drying mechanics of cathode processed with both solvents is investigated though implementation of a chemical‐engineering‐based model to better understand the utilization of heat provided by experimentally determining the heat and mass transfer coefficients. Electrochemical performance is also evaluated to determine the impact of drying temperature on cycling performance. Binder distribution is determined via various methods to confirm differences in binder homogeneity as a function of both solvent and drying temperature. Identified is the large difference in the efficiency in which the heat is used as well as an ideal drying temperature for both aqueous and non‐aqueous processed cathodes. Also identified is the increased sensitivity to processing temperature for aqueous processed electrodes compared to non‐aqueous processed counterparts, pointing to the possibility of tuned drying regimes which would capitalize on the potential cost savings of aqueous processing for cathodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗