Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational homogenization”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Porting hypre to heterogeneous computer architectures: Strategies and experiences

We report that linear systems are occurring in many applications, and solving them can take a large amount of the total simulation time. The high performance library hypre provides a variety of interfaces and linear solvers, including various multigrid methods, that have achieved good scalability on a variety of homogeneous parallel computer architectures. Heterogeneous architectures with nodes that have both CPUs and accelerators provide new challenges, since they require more fine-grained parallelism and reduced data movement between different memories on a single node as well as across nodes. We will discuss our experiences and strategies to port hypre to heterogeneous computers with accelerators, including the design of a new memory model, the use of abstractions, the BoxLoop macros in the structured and semi-structured interfaces, and the restructuring of algebraic multigrid (AMG) into modular components. We present numerical experiments comparing CPU and GPU performance for several test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Analyzing photovoltaic module mechanics using composite plate theories and finite element solutions

Deflection and stress calculated from an experimentally validated, high-fidelity finite element model (FEM) of a photovoltaic module experiencing mechanical load was compared to results from a simplified FEM treating the module laminate as a homogenized composite using a rule of mixtures approach, and further compared to analytical calculations treating the module as a Kirchoff-Love flat plate. The goal of this study was to determine the error incurred by analyzing module mechanics with varying levels of simplification, since resolving the aspect ratios of a module is computationally expensive. Homogenized FEMs were found to underpredict peak deflection under a 1.0 kPa load by between 13 and 19% for lower and upper bound application of the rule of mixtures. However, module shape was captured, implying that a useful replication of a resolved model could be achieved with a reduced, calibrated material stiffness. Homogenized stress results captured glass layer tensile stress components to within 46 to 52% at a sample location of interest, though agreement was poor through the remainder of the laminate due to the lack of material resolution. For plate theory, deflection was overpredicted by 45 to 67% for upper and lower bound homogenizations, and frame-adjacent module shapes were not adequately replicated. Stress results mirrored FEM trends but magnitudes were not well correlated to resolved model values. These results support the use of homogenized laminate models for module shape derivation, though resolved models remain necessary for stress analyses. The accuracy of plate theory was found to be inadequate for most applications.

14 SOLAR ENERGY↗
NASA NTRS1977

Improved methods for large-eddy simulations of turbulence

By using Fourier transforms for evaluating spatial derivatives, the accuracy of the large-eddy simulation of homogeneous isotropic turbulence is improved. Especially, the treatment of certain terms that arise in filtering the equations is considerably improved in both speed and accuracy. Use of vorticity as the principal variable is shown to be a viable and potentially useful alternative to the primitive variables. A method of deriving conservation properties of numerical schemes is given which is much simpler than previous methods and is widely applicable. The methods are applied to the computation of homogeneous isotropic turbulence, and it is found that the subgrid scale model is improved by using finite differences in place of 'exact' derivatives.

Mansour, N. N.↗
NASA NTRS2002

Graph Partitioning for Parallel Applications in Heterogeneous Grid Environments

The problem of partitioning irregular graphs and meshes for parallel computations on homogeneous systems has been extensively studied. However, these partitioning schemes fail when the target system architecture exhibits heterogeneity in resource characteristics. With the emergence of technologies such as the Grid, it is imperative to study the partitioning problem taking into consideration the differing capabilities of such distributed heterogeneous systems. In our model, the heterogeneous system consists of processors with varying processing power and an underlying non-uniform communication network. We present in this paper a novel multilevel partitioning scheme for irregular graphs and meshes, that takes into account issues pertinent to Grid computing environments. Our partitioning algorithm, called MiniMax, generates and maps partitions onto a heterogeneous system with the objective of minimizing the maximum execution time of the parallel distributed application. For experimental performance study, we have considered both a realistic mesh problem from NASA as well as synthetic workloads. Simulation results demonstrate that MiniMax generates high quality partitions for various classes of applications targeted for parallel execution in a distributed heterogeneous environment.

Bisws, Rupak↗
NASA NTRS1993

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An Alternative Casting Technique to Improve the Creep Resistance of Cast INCONEL Alloy 740H

The increasing performance requirements of power plant designs, such as advanced-ultra supercritical (A-USC), require the use of Ni-based superalloys to replace high-strength, ferritic-martensitic steels for components subjected to temperatures above 898 K (625 °C) and for austenitic stainless steels components at temperatures above 973 K (700 °C). To date, commercial Ni-based superalloy INCONEL 740H has been shown to be appropriate for use in A-USC power plants as boiler components in a wrought product. However, large complex components in boilers as well as other casings in the turbine and valve chest require castings of a thick-wall nature. Using the alloy in its cast form would be significantly valuable in terms of range of component size, geometry and complexity. Previous investigations revealed short creep lives from cast INCONEL alloy 740H. Here in this investigation, an alternative casting route that utilized a melt procedure resulting in a fine-grain casting, and in conjunction with a computationally optimized homogenization heat treatment, not only controlled the grain size and grain boundary structure but minimized chemistry variability and segregation. A primarily equiaxed and homogenous grain size distribution was obtained from this approach with better repartition of M 23 C 6 carbides along the grain boundaries. Furthermore, better than 38 pct increase was obtained for this material in comparison to the creep life obtained from the best performing conventionally cast material. More importantly, the fine-grain homogenized (FGH) casting route resulted in the Larson–Miller plot for this material that coincided with that of wrought alloy 740. At low creep stresses (with a test still in progress), the FGH casting is resulting in higher values of the LMP than the wrought alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Long-Term Creep Behavior of a CoCrFeNi Medium-Entropy Alloy

We report that current research shows that high-entropy alloys (HEAs) can offer superior oxidation resistance compared to more traditional austenitic stainless steels with equivalent tensile mechanical performance. While many reports have been published on creep of CoCrFeNiMn, few focus on the CoCrFeNi family of low to medium entropy alloys. In this manuscript, uniaxial tensile creep test results for more than 46,000 accumulated hours are presented on the equiatomic CoCrFeNi alloy. The alloy was manufactured at NETL using a vacuum induction melt solidification approach with a computationally designed homogenization cycle. After completing the homogenization cycle, the ingot was first hot forged and then hot rolled into 10 mm thick plate, a product form that is relevant to many industries. Creep tests were performed in the temperature range of 600-650 °C using applied stresses from 120 to 189 MPa. This study builds upon previous results on the CoCrFeNiMn family of high-entropy alloys and offers a modicum of clarification about performance compared to small button melts by using ASTM standard test procedures and mechanical test specimens. The apparent creep activation energy and creep stress exponent were similar to those for CoCrFe medium-entropy and CoCrFeNiMn high-entropy alloys. The post-test creep microstructure showed some Cr segregation, which has also been noted for other CoCrFe alloys. The CoCrFeNi alloy also presents a double minimum or two minimum creep rate phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Variation of γ′ Formers and Refractory Elements for Enhanced Creep resistance and Phase Stability of HAYNES® 282® Alloy

Modifications to the chemistry of alloy 282 were performed to improve the alloy’s resistance to long term creep deformation as well as phase stability. Alloys were prepared using vacuum induction melting, computationally optimized homogenization heat treatment and hot working. Phase stability studies were carried out for up to 5,000 hours at 800°C and 900°C while creep testing was performed at conditions leading to lives past 7,000 hours for temperatures ranging from 740°C to 900°C. The formation of σ and μ phases were reduced in the modified alloy. Atom probe tomography (APT) was performed on the specimens from the phase stability study to investigate changes in the elemental partitioning to the γ and γ′ phases. Post deformation microstructures were analyzed using EBSD and TEM to study the influence of the detrimental phases on damage accumulation during creep.

Detrois, Martin↗
DOE OSTI2025

The final WaZP galaxy cluster catalog of the Dark Energy Survey and comparison with SZE data

In this work, we present and characterize the galaxy cluster catalog detected by the WaZP cluster finder, which is not based on red-sequence identification, on the full six years of observations of the Dark Energy Survey (DES-Y6). The full catalog contains over 400k detected clusters with richnesses, Ngals, above 5 and that reach redshifts up to 1.3. We also provide a version of the catalog where the observation depth and richness computation are homogenized to be used for cosmology, containing 33k rich (Ngals >25) clusters. We compare our results with the previous WaZP catalog obtained from the DES first-year data release (DES-Y1). We find that essentially all clusters within the common footprint and depth limit are recovered. The deeper observations on DES-Y6 and the more complete available spectroscopic redshift sample lead to improvements in the redshifts of the clusters, resulting in an average scatter of 1.4% and offset of 0.2%. The optical clusters are also cross-matched with Sunyaev Zel'dovich Effect (SZE) cluster samples detected by the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT). We find that essentially all SZE clusters with reasonable overlapping footprint have a corresponding WaZP cluster. Conversely, 90% of the optical detections with richness greater than 150 have a counterpart in the deeper regions of the SZE surveys. Based on cross-match with the SZE catalogs, we also find that 15-20% of the SZE matched systems have more than one possible WaZP counterpart at the same redshift and within the SZE R500c, indicating possible interacting or unrelaxed systems. Finally, given the optical and SZE beams, WaZP and SZE centerings are found to be consistent. A more detailed study of the SZE-WaZP mass-richness relation will be presented in a separate paper.

Benoist, C. [OCA, Nice, Lab. Lagrange; LIneA, Rio ↗
NASA NTRS1985

The preparation and archiving of machine-readable astronomical data

The preparation of machine-readable astronomical catalogs for archiving and distribution from the ADC is described. Examples of specific format and data structure changes to improve compatibility with other computers, to homogenize the presentation of astronomical data among catalogs, and to increase storage efficiency are given. Various checks which are often made to insure internal consistency and for comparison with published versions are identified and the importance of permanently archiving astronomical data in the best possible condition is emphasized.

Warren, W. H., Jr.↗
NASA NTRS2001

Micromechanics-Based Inelastic Finite Element Analysis Accomplished Via Seamless Integration of MAC/GMC

A critical issue in the micromechanics-based analysis of composite structures becomes the availability of a computationally efficient homogenization technique: one that is 1) Capable of handling the sophisticated, physically based, viscoelastoplastic constitutive and life models for each constituent; 2) Able to generate accurate displacement and stress fields at both the macro and the micro levels; 3) Compatible with the finite element method. The Generalized Method of Cells (GMC) developed by Paley and Aboudi is one such micromechanical model that has been shown to predict accurately the overall macro behavior of various types of composites given the required constituent properties. Specifically, the method provides "closed-form" expressions for the macroscopic composite response in terms of the properties, size, shape, distribution, and response of the individual constituents or phases that make up the material. Furthermore, expressions relating the internal stress and strain fields in the individual constituents in terms of the macroscopically applied stresses and strains are available through strain or stress concentration matrices. These expressions make possible the investigation of failure processes at the microscopic level at each step of an applied load history.

Arnold, Steven M.↗
NASA NTRS2008

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗
NASA NTRS2019

Micromechanical Modeling of High-Strain Thin-Ply Composites

This paper presents a micromechanical analysis for the elastic and viscoelastic behavior of high-strain thin-ply composites. The modeling approach is based on unit cell homogenization. The geometry of the internal weave architecture and ply configuration is characterized via micrographic analysis and explicitly modeled in the unit cell. The composites are modeled as Kirchhoff plates and the homogenization analysis computes the effective relaxation ABD matrix represented by Prony series using the elastic and viscoelastic properties of the constituent fiber and matrix. The formulation of the micromechanical model and numerical implementation d. Composite laminates with 3-ply and 4-ply configurations are studied.

Hamillage, Milinda Yapa↗
NASA NTRS1992

Critical assessment of Reynolds stress turbulence models using homogeneous flows

In modeling the rapid part of the pressure correlation term in the Reynolds stress transport equations, extensive use has been made of its exact properties which were first suggested by Rotta. These, for example, have been employed in obtaining the widely used Launder, Reece and Rodi (LRR) model. Some recent proposals have dropped one of these properties to obtain new models. We demonstrate, by computing some simple homogeneous flows, that doing so does not lead to any significant improvements over the LRR model and it is not the right direction in improving the performance of existing models. The reason for this, in our opinion, is that violation of one of the exact properties can not bring in any new physics into the model. We compute thirteen homogeneous flows using LRR (with a recalibrated rapid term constant), IP and SSG models. The flows computed include the flow through axisymmetric contraction; axisymmetric expansion; distortion by plane strain; and homogeneous shear flows with and without rotation. Results show that for most general representation for a model linear in the anisotropic tensor, performs either better or as good as the other two models of the same level.

Shabbir, Aamir↗
NASA NTRS1975

Nonstationary homogeneous nucleation

A computer calculation was made of the relaxation time to steady-state nucleation in a supersaturated vapor. A binary-collision model of nucleation based on classical liquid-drop thermodynamics was used. The nondimensional nucleation relaxation time is plotted as a function of the critical drop size for various values of the surface-tension parameter. The curves permit a rough estimation of the limit of the validity of steady-state nucleation theory.

Harstad, K. G.↗
NASA NTRS1990

Modelling a nickel cadmium battery as a homogeneous device

A computer model of the nickel-cadmium (Ni-Cd) battery cell has been developed. This one-dimensional macrohomogeneous model predicts performance for Ni-Cd battery cells based on electrochemical phenomena. Reaction rates, concentrations, current densities, porosities, and potentials are predicted over a range of conditions. A description of the model, some initial results, and plans for continued development are presented.

Timmerman, Paul J.↗