Search NASA⌕ Search

SEARCH · Search NASA

Results for “RVE”

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 19 records

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗

An adaptive granular representative volume element model with an evolutionary periodic boundary for hierarchical multiscale analysis

The hierarchical multiscale analysis normally utilizes a microscopic representative volume element (RVE) model to capture path/history-dependent macroscopic responses instead of using phenomenological constitutive models. However, for problems involving large deformation, the current RVE model used in geomechanics may lose representative properties due to the progressive distortion of the RVE box, unless a proper reinitialization is applied. This work develops an adaptive RVE model in conjunction with an evolutionary periodic boundary (EPB) algorithm for hierarchical multiscale analysis of granular materials undergoing large deformation based on a recent RVE model proposed for coupling molecular dynamics and the material point method. The proposed adaptive RVE model avoids the reinitialization of the RVE box that even undergoes extremely large shear deformation; meanwhile, it accounts for the deformation history of the RVE model and treats the interaction between boundary particles and other image particles in a more efficient way. Numerical examples with extremely large deformation are used to illustrate the adaptive granular RVE model enhanced by the proposed EPB algorithm. Furthermore, some key features of this new methodology are further discussed for clarification.

42 ENGINEERING↗

Effects of Resistive Vibration Exercise Combined with Whey Protein and KHCO3 on Bone Tturnover Markers in Head-down Tilt Bed Rest (MTBR-MNX Study)

High protein intake further increases bone resorption markers in head-down tilt bed rest (HDBR), most likely induced by low-grade metabolic acidosis. Adding an alkaline salt to a diet with high protein content prevents this additional rise of bone resorption markers in HDBR. In addition, high protein intake, specifically whey protein, increases muscle protein synthesis and improves glucose tolerance, which both are affected by HDBR. Resistive vibration exercise (RVE) training counteracts the inactivity-induced bone resorption during HDBR. To test the hypothesis that WP plus alkaline salt (KHCO3) together with RVE during HDBR will improve bone turnover markers, we conducted a randomized, three-campaign crossover design study with 12 healthy, moderately fit male subjects (age 34+/-8 y, body mass [BM] 70 +/- 8 kg). All study campaigns consisted of a 7-d ambulatory period, 21days of -6 deg. head-down tilt bed rest (HDBR), and a 6-d recovery period. Diet was standardized and identical across phases. In the control (CON) campaign, subjects received no supplement or RVE. In the intervention campaigns, subjects received either RVE alone or combined with WP and KHCO3 (NEX). WP was applied in 3 doses per day of 0.6 g WP/kg BM together with 6 doses of 15 mmol KHCO3 per day. Eleven subjects completed the RVE and CON campaign, 8 subjects completed all three campaigns. On day 21 of HDBR excretion of the bone resorption marker C-telopeptide (CTX) was 80+/-28% (p<0.001) higher than baseline, serum calcium concentrations increased by 12 +/- 29% (p<0.001) and serum osteocalcin concentrations decreased by 6+/-12% (p=0.001). Urinary CTX excretion was 11+/- 25% (p=0.02) lower on day 21 of HDBR in the RVE- and tended to decrease by 3+/- 22% (p=0.06) in the NEX campaign compared to CON. Urinary calcium excretion was higher on day 21 in HDBR in the RVE and NEX (24+/- 43% p=0.01; 25+/- 37% p=0.03) compared to the CON campaign. We conclude that combination of RVE with WP/KHCO3 was not superior to RVE alone in any of these results.

Graf, Sonja↗

On the Representativity of Electrode Microstructure Parameters and Their Electrochemical Response for Lithium Ion Batteries

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the microstructure analysis. A microstructure scale electrochemical model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. To evaluate the representativeness, and thus relevance, of these model predictions, a coupled representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response. Electrode microstructure parameters representativeness has been first quantified using the representative volume element (RVE) methodology. The RVE major flaw is that ultimately it can only conclude if a FOV contains representative subvolumes of the FOV, but not if the FOV itself is representative of the electrode volume. Analysis can conclude negatively ('FOV is not representative'), but not positively ('FOV is representative'). One major contribution of this work was to quantify the convergence of the RVE size with the FOV, to actually investigate the FOV representativeness and thus partly remedy this intrinsic limitation. The analysis determined that performing a standard RVE calculation, without exploring its FOV convergence, is likely to strongly underestimate the actual RVE size. The new RVE methodology has been automated in the NREL open-source Microstructure Analysis Toolbox (MATBOX) and is available to the battery community. Representativeness of microstructure parameters is however only an intermediate step, as the end-results of an electrochemical model are performances predictions. Indeed, what is the practical consequence of a given deviation for a microstructure parameter? The microstructure parameter deviation propagations to the 3D microstructure scale electrochemical response have been then quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Such deviation propagation analysis is analogous to error propagation analysis and is necessary to determine the relevance of microstructure scale model predictions for macroscale predictions. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Therefore, we introduced the novel concept of electrochemical RVE (eRVE) that is a function of the operating conditions (thus defined as a dynamic RVE), with an increasing dependence with the C-rate. Representativity analysis of the investigated cell determined a FOV of 144.4 x 54.4 m2 is large enough to establish a convergence on the representative section areas for low to intermediate C-rate (=2.5C), but not large enough to conclude for higher rates. This work aims to emphasize the importance of representativity analysis for LIB electrode microstructures, as it is required to estimate the error, and thus the relevance, of microstructure parameters intended to be used in macroscale models. The methodology and results can help researchers to select the relevant imaging and associated FOV required to provide accurate enough microstructure parameters.

ADVANCED PROPULSION SYSTEMS↗

Measuring representative volume elements from high-resolution grain-scale strain fields

Most crystalline materials present a highly heterogeneous response at the microscale, which can be affected by both internal factors (such as microstructural parameters) and external factors (such as loading). Relating microscale inhomogeneities to the macroscale response of a material requires the use of homogenisation techniques, usually based on the concept of a representative volume element (RVE)—the smallest volume of material that represents the global average response. Here, in this work, we present a new and robust experimental method of measuring the size of a strain-based RVE from high-resolution grain-scale strain fields obtained using digital image correlation (DIC). The proposed method is based on the statistical (stereological) nature of the RVE, which has been widely adopted in numerical studies, and involves dividing a strain field into randomly selected regions of varying sizes and statistically analysing the distributions of average strains within them. To validate the new method, we generate a large number of synthetic strain fields from a fractional Gaussian noise algorithm. The proposed stereological method is shown to be capable of producing reliable RVE measurements from a very large range of possible microscale strain fields while at the same time being robust in that it can produce RVE measurement results even in cases where other existing methods may be unable to do so. The proposed method has a low field-of-view requirement, only needing a field-of-view about 1.2 times as large as the RVE to produce reliable measurements. In addition, the stereological method offers significant flexibility since its statistical nature allows for control over how strict the RVE measurement should be in each case.

36 MATERIALS SCIENCE↗

Discrete dislocation dynamics for crystal RVEs. Part 1: Periodic network kinematics

A novel implementation of the dislocation flux boundary condition in discrete dislocation dynamics is presented. The continuity of the individual dislocation loops in a periodic representative crystal volume (RVE) is enforced across the boundary of the RVE with the help of a dual topological description for representing dislocation line kinematics in two equivalent spaces representing the deforming crystal, the RVE and the unbounded crystal spaces. The former describes the motion of the dislocations in the simulated crystal RVE whereas the latter represents the motion of dislocations in an infinite space containing all replicas of the RVE. A mapping between the two spaces forms the basis of the implementation of flux boundary condition. The implementation details are discussed in the context of statistical homogeneity of bulk crystals undergoing macroscopically homogeneous plastic deformation. In this case, the boundary nodes associated with dislocation segments bear no relevance in the motion of the dislocations. Finally, some test cases are presented and discussed to establish the proposed approach.

42 ENGINEERING↗

Experimental investigation of multi-scale strain-field heterogeneity in rocks

Rocks have inherently heterogeneous micro-structures. The mechanical response of these heterogeneous micro-structures primarily governs the overall macroscopic behavior of intact rocks, as these features cause heterogeneity in the distribution of stress and strain within a rock volume. These attributes have generally been studied through numerical tools, but in this study, the stress-induced multi-scale spatial heterogeneity in the strain-field of three types of rocks was analyzed experimentally. The full-field non-contact strain measurement approach of 2D Digital Image Correlation (DIC) was employed for real-time evaluation of the progression of strain-field heterogeneity across the surface of the rock specimens when subjected to uniaxial compression. The results show that the strain-field heterogeneity is a characteristic of the inherent heterogeneity of the rock specimens, and that it amplifies with the magnitude of stress applied on the rock specimens. The changes in the strain-field heterogeneity with loading were consistent with established microcrack evolution processes in the rock specimens. The statistical fluctuation in the strain-field heterogeneity of the three rocks at varying gauge lengths (GL) was then analyzed for the determination of the Representative Volume Element (RVE) length-scale at increasing load levels. The evolution of the RVE length-scale with loading facilitates an objective way of determination of the stress at which the strain-field heterogeneity increases and begins to influence the overall global response of the material. The RVE length-scale variation with loading showed that at the RVE length-scale reaches the specimen-scale at load levels well below the peak strength.

2D-digital image correlation↗

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗

Evaluation of Constitutive Relationship for Silicon Carbide Cladding Material

A curved representative volume element (RVE) was employed to analyze SiC/SiC tubular composites. By imposing periodicity on the tubular model through appropriate boundary conditions, the elastic constants of the composite was evaluated. It was found that as the radius of the tubular composite increases, the more the flat RVE model and curved RVE model disagreed.

Thandaga Nagaraju, Hemanth↗

Stress Analysis and Failure Behavior of SiC/SiC Textile Composite Tubes

A novel method for analyzing tubular textile composites was developed by obtaining a curved representative volume element (RVE) from the transformation of a flat RVE. Periodic boundary conditions were imposed on the curved RVE. The stresses in a transeversely isotropic yarn were compared between the two RVEs. It was found that there was a discrepancy between the two models.

Thandaga Nagaraju, Hemanth↗

Statistical Homogenization of Elastic and Fracture Properties of a Sample Selective Laser Melting Material

Selective laser melting (SLM) is an additive manufacturing technique commonly used in the rapid prototyping of components. The complexity of the SLM microstructure poses a unique challenge to deriving effective mechanical properties at different length scales. Representative volume elements (RVEs) are often used to homogenize the material properties of composites. Instead of RVEs, we use statistical volume elements (SVEs) to homogenize the elastic and fracture properties of the material. This relates the inherent variation of a material’s microstructure to the variation in its mechanical properties at different observation scales. The convergence to the RVE limit is examined from two perspectives: the stability of the mean value as the SVE size increases for the mean-based approach, and the tendency of the normalized variation in homogenized properties to zero as the SVE size increases for the variation-based approach. Fracture properties tend to make the RVE limit slower than do elastic properties from both perspectives. There are also differences between vertical (normal to printing plane) and horizontal (in-plane) properties. While the elastic properties tend to make the RVE limit faster for the horizontal direction, i.e., having a smaller variation and more stable mean value, the fracture properties exhibit the opposite effect. We attributed these differences to the geometry of the melt pools.

36 MATERIALS SCIENCE↗

Meso-scale modeling and damage analysis of carbon/epoxy woven fabric composite under in-plane tension and compression loadings

The mechanical properties and damage behaviors of carbon/epoxy woven fabric composite under in-plane tension and compression are studied at the meso-scale level through experiment and simulation. An efficient representative volume element (RVE) modeling method with consistent mesh, high yarn volume fraction and realistic geometry is proposed. The material constitutive laws with plasticity, tension-compression asymmetry and damage evolution are established for the three components - yarn, matrix and interface, respectively. Significantly different mechanical properties and damage evolutions are observed depending on loading conditions and initial geometry characteristics. It shows a non-linear stress-strain curve with clear transition region and intensive damage in tension, while a quasi-linear behavior up to facture is observed in compression with little damage prior to final fracture. Furthermore, compared to the constant Poisson's ratio with straining in compression, a dramatic increase in Poisson's ratio appears in tension. Simulation shows damage mechanisms including transverse damage, matrix damage and delamination, which all play critical roles in the property evolution. In particular, the rapid damage accumulation after elastic deformation destroys the strong bonds and causes the easy deformation of transverse yarns which results in the transition region and large Poisson's ratio in tension. All the mechanical behaviors and damage evolutions are well captured and explained with the current RVE model.

42 ENGINEERING↗

The SEEL motif and members of the MYB-related REVEILLE transcription factor family are important for the expression of LORELEI in the synergid cells of the Arabidopsis female gametophyte

Synergid cells in the micropylar end of the female gametophyte are required for critical cell–cell signaling interactions between the pollen tube and the ovule that precede double fertilization and seed formation in flowering plants. LORELEI (LRE) encodes a putative GPI-anchored protein that is expressed primarily in the synergid cells, and together with FERONIA, a receptor-like kinase, it controls pollen tube reception by the receptive synergid cell. Still, how LRE expression is controlled in synergid cells remains poorly characterized. We identified candidate cis-regulatory elements enriched in LRE and other synergid cell-expressed genes. One of the candidate motifs (‘TAATATCT’) in the LRE promoter was an uncharacterized variant of the Evening Element motif that we named as the Short Evening Element-like (SEEL) motif. Deletion or point mutations in the SEEL motif of the LRE promoter resulted in decreased reporter expression in synergid cells, demonstrating that the SEEL motif is important for expression of LRE in synergid cells. Additionally, we found that LRE expression is decreased in the loss of function mutants of REVEILLE (RVE) transcription factors, which are clock genes known to bind the SEEL and other closely related motifs. Here, we propose that RVE transcription factors regulate LRE expression in synergid cells by binding to the SEEL motif in the LRE promoter. Identification of cis-regulatory elements and transcription factors involved in the expression of LRE will serve as a foundation to characterize the gene regulatory networks in synergid cells.

59 BASIC BIOLOGICAL SCIENCES↗

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗

An integrated computational materials engineering framework to analyze the failure behaviors of carbon fiber reinforced polymer composites for lightweight vehicle applications

A bottom-up multi-scale modeling approach is used to develop an Integrated Computational Materials Engineering (ICME) framework for carbon fiber reinforced polymer (CFRP) composites, which has the potential to reduce development to deployment lead time for structural applications in lightweight vehicles. In this work, we develop and integrate computational models comprising of four size scales to fully describe and characterize three types of CFRP composites. In detail, the properties of the interphase region are determined by an analytical gradient model and molecular dynamics analysis at the nano-scale, which is then incorporated into micro-scale unidirectional (UD) representative volume element (RVE) models to characterize the failure strengths and envelopes of UD CFRP composites. Then, the results are leveraged to propose an elasto-plastic-damage constitutive law for UD composites to study the fiber tows of woven composites as well as the chips of sheet molding compound (SMC) composites. Subsequently, the failure mechanisms and failure strengths of woven and SMC composites are predicted by the meso-scale RVE models. Finally, building upon the models and results from lower scales, we show that a homogenized macro-scale model can capture the mechanical performance of a hat-section-shaped part under four-point bending. Along with the model integration, we will also demonstrate that the computational results are in good agreement with experiments conducted at different scales. The present work illustrates the potential and significance of integrated multi-scale computational modeling tools that can virtually evaluate the performance of CFRP composites and provide design guidance for CFRP composites used in structural applications.

36 MATERIALS SCIENCE↗

Computational micromechanics model based failure criteria for chopped carbon fiber sheet molding compound composites

Chopped carbon fiber sheet molding compound has a great potential in lightweight automotive, marine, and aerospace applications. One of the most challenging tasks is to predict the failure strength of the material due to its anisotropy and heterogeneity, as well as complex stress states in real-world working conditions. In this work, a novel constitutive model of carbon fiber chip is proposed to capture the pre- and post-failure behaviors under different loading modes. On this basis, we propose a new computational micromechanics model, which is calibrated and validated by uniaxial tensile, compressive, and in-plane shear experiments. Furthermore, a set of microstructures representative volume element (RVE) models under complex loading conditions are reconstructed to understand the relationship between the microstructure characteristics and the failure envelopes. Finally, several modified versions of classical failure criteria are proposed for anisotropic materials with consideration of the fiber orientation tensor. The modified Tsai-Wu failure criterion, which shows the best accuracy among all failure criteria, is highlighted in the comparative study.

36 MATERIALS SCIENCE↗

Shade signals alter the expression of circadian clock genes in newly-formed bioenergy sorghum internodes

Stem internodes of bioenergy sorghum inbred R.07020 are longer at high plant density (shade) than at low plant density (control). Initially, the youngest newly-formed subapical stem internodes of shade-treated and control plants are comparable in length. However, full-length internodes of shade-treated plants are three times longer than the internodes of the control plants. To identify the early molecular events associated with internode elongation in response to shade, we analyzed the transcriptome of the newly-formed internodes of shade-treated and control plants sampled between 4 and 6 hr after the start of the light period (14 hr light/10 hr dark). Sorghum genes homologous to the Arabidopsis shade marker genes ATHB2 and PIL1 were not differentially expressed. The results indicate that shade signals promote internode elongation indirectly because sorghum internodes are not illuminated and grow while enclosed with leaf sheaths. Sorghum genes homologous to the Arabidopsis morning-phased circadian clock genes LHY, RVE, and LNK were downregulated and evening-phased genes such as TOC1, PRR5, and GI were upregulated in young internodes in response to shade. We hypothesize that a change in the function or patterns of expression of the circadian clock genes is the earliest molecular event associated with internode elongation in response to shade in bioenergy sorghum. Increased expression of CycD1, which promotes cell division, and decreased expression of cell wall-loosening and MBF1-like genes, which promote cell expansion, suggest that shade signals promote internode elongation in bioenergy sorghum in part through increasing cell number by delaying transition from cell division to cell expansion.

54 ENVIRONMENTAL SCIENCES↗

An FFT-based approach for Bloch wave analysis: application to polycrystals

A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.

97 MATHEMATICS AND COMPUTING↗