Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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 559 records · Page 31

Implementation of extrinsic cohesive zone model (ECZM) in 2D finite-discrete element method (FDEM) using node binding scheme

The combined finite-discrete element method (FDEM) has been widely used for rock fracturing simulations. Conventionally, FDEM is realized using the intrinsic cohesive zone model (ICZM); however, it has the drawback of artificial compliance and high computational expense. As a complement, the extrinsic cohesive zone model (ECZM) is seen to be realized in FDEM recently, whereas the node splitting scheme utilized is cumbersome. Here, within the framework of ICZM-based FDEM, we propose a node binding scheme to efficiently bind the pre-discretized finite elements and thus guarantee the continuum behavior of materials in the elastic stage. The yield surfaces, controlled by ECZM, are dynamically embedded by invoking the pre-inserted cohesive elements. The effectiveness and efficiency of the proposed approach are validated and tested by performing a suite of numerical experiments. Compared with ICZM-based FDEM, the proposed approach can correctly capture material deformation and reduce the computation cost. In contrast to the existing ECZM-based FDEM, the proposed approach can overcome the frequent and complex element topology updating. Finally, this work provides a novel perspective that fully inherits the advantages of both ICZM and ECZM, but circumvents their shortcomings, which guarantees a more efficient and effective simulation of brittle material evolution from continuum to discontinuum.

58 GEOSCIENCES↗

Water-dispersible nanocolloids and higher temperatures promote the release of carbon from riparian soil

Increasing temperatures in alpine regions accompanied by glacial retreat is occurring rapidly due to climate change. This may affect riparian soils by increasing weathering rates, resulting in greater organic carbon (OC) release to rivers via movement of iron-containing colloids and nanonminerals. Increased concentrations of iron- or silcate-nanominerals would result in higher surface area for OC adsorption. To test the influence of temperature on OC leaching, we examined mineral weathering and nanocolloid facilitated release of OC through a series of controlled laboratory batch and column experiments using sediment from the banks of the Nisqually River, Mount Rainier in Washington State (USA). Additional experiments were conduscted using the same sediments, but with an illite amendment added to test the influence of additional surface area and nanominerals that many sediments along the Nisqually River contain. These higher- and lower-surface-area sediments (i.e., sediments with and with the illite amendment) were incubated for three months at 4 or 20 °C, followed by batch and column OC leaching tests. Results show that OC leaching rates for 20 °C were two to three times greater than for 4 °C. Further, our results suggest that nanocolloids are responsible for moving this increased OC load from these sediments. When hydrologically connected, OC is released from back sediments to rivers faster than presently anticipated in fluvial environments experiencing climate change-induced glacial retreat. Further, a one-dimensional, finite-element computational model developed for this study estimates that a 1 °C increase in temperature over a 3-mo summer runoff period increases the OC release rate from sediments by 79%.

54 ENVIRONMENTAL SCIENCES↗

Advanced multimaterial shape optimization methods as applied to advanced manufacturing of wind turbine generators

Abstract Currently, many utility‐scale wind turbine generator original equipment manufacturers are dependent on imported rare earth permanent magnets, which are susceptible to market risks from cost instability. To lower the production costs of these generators and stay competitive in the market, several small wind manufacturers are pursuing continuous improvements to both generator design and manufacturing. However, traditional design and manufacturing methods have yielded marginal improvements in wind power performance. This work presents novel methods to redesign a baseline 15‐kW wind turbine generator with reduced rare‐earth permanent magnets by leveraging cutting‐edge three‐dimensional (3D) printed polymer‐bonded permanent magnets and steel. Symmetric, asymmetric, and multimaterial‐magnet parametrization methods are introduced for shape optimization. We extend the symmetric and asymmetric methods to the back iron in the stator to further investigate the impact and opportunities for performance improvements with lesser active materials. We employ a design‐of‐experiments approach with parametric computer‐aided design for shape generation and evaluate different designs by magneto‐thermal modeling and finite‐element analysis. We use adaptive sampling technique to identify better performing designs with lesser magnet mass, higher efficiency, and lower cogging torque when compared with the baseline generator. Asymmetric pole designs resulted in a magnet mass in the range of 4.77–5.37 kg, which was 27%–35% lighter than the baseline generator, suggesting that a new design freedom exists that can be enabled by advanced manufacturing, such as 3D printing. Shaping the back iron in the stator resulted in material savings in electrical steel of up to 14.62 kg, which was 20% lighter than the baseline stator. We conducted a structural analysis to evaluate an optimized asymmetric rotor design from the point of view of mechanical integrity and air‐gap stiffness. The magnetically optimal shape profile was shown as having a positive impact on the radial stiffness, and an optimal solution was discovered to reduce the structural mass by nearly 30 kg, which was 29% lighter than the baseline.

17 WIND ENERGY↗

A combined experimental and numerical approach that eliminates the non-uniqueness associated with the Johnson-Cook parameters obtained using inverse methods

Abstract Johnson-Cook constitutive model is a commonly used material model for machining simulations. The model includes five parameters that capture the initial yield stress, strain-hardening, strain-rate hardening, and thermal softening behavior of the material. These parameters are difficult to determine using experiments since the conditions observed during machining (such as high strain-rates of the order of $$10^5$$ 10 5 /sec - $$10^6$$ 10 6 /sec) are challenging to recreate in the laboratory. To address this problem, several researchers have recently proposed inverse approaches where a combination of experiments and analytical models are used to predict the Johnson-Cook parameters. The errors between the measured cutting forces, chip thicknesses and temperatures and those predicted by analytical models are minimized and the parameters are determined. In this work, it is shown that only two of the five Johnson-Cook parameters can be determined uniquely using inverse approaches. Two different algorithms, namely, Adaptive Memory Programming for Global Optimization (AMPGO) and Particle Swarm Optimization (PSO), are used for this purpose. The extended Oxley’s model is used as the analytical tool for optimization. For determining a parameter’s value, a large range for each parameter is provided as an input to the algorithms. The algorithms converge to several different sets of values for the five Johnson-Cook parameters when all the five parameters are considered as unknown in the optimization algorithm. All of these sets, however, yield the same chip shape and cutting forces in FEM simulations. Further analyses show that only the strain-rate and thermal softening parameters can be determined uniquely and the three parameters present in the strain-hardening term of the Johnson-Cook model cannot be determined uniquely using the inverse method. A combined experimental and numerical approach is proposed to eliminate this determine all parameters uniquely.

42 ENGINEERING↗

A comparative analysis of residual stresses from friction stir processing of aluminum cast 380 and wrought 7075 alloy sheets: experimental characterization and modeling

Residual stresses are often overlooked in friction stir processing (FSP), but their significant impact on fatigue performance necessitates their consideration in optimizing processing parameters. The first step in this effort is understanding how process conditions influence residual stress distributions, especially across different alloys. This study focuses on determining and explaining the through-thickness residual stress variations and the effect of process temperature on the residual stress magnitude in wrought AA7075 and cast AA380.0 alloys. Additionally, for AA380.0, the impact of a second FSP pass was investigated. To achieve this, hole-drilling electronic speckle pattern interferometry (ESPI) and the thermal pseudo-mechanical (TPM) model within finite element analysis were employed to study the 3D distributions of in-plane residual stresses in processed samples under various conditions. A key finding was the varying impact of process temperatures on residual stress magnitudes. Higher process temperatures reduced stresses in AA380.0 but increased them in AA7075. Additionally, the through-thickness stress distributions differed between the two alloys. Further analysis revealed that yield stresses are crucial in explaining these phenomena and the effects of additional FSP passes. Further, this fundamental understanding will be vital in guiding the efforts to mitigate residual stresses and assess their impact on the performance of FSP aluminum alloys.

36 MATERIALS SCIENCE↗

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗

Impacts of Mode Mixity on Controlled Spalling of (100)-Oriented Germanium

Controlled spalling is a technology to prepare single-crystal thin films of semiconductors by fracture with a subsurface crack propagating nearly parallel to the substrate surface. Practical applications require uniform thickness and a smooth surface across the whole film. Both wafer-scale and patterned-stressor-defined small-area spalling of germanium substrates are conducted experimentally and numerically. River line features are observed on spalled surfaces close to lateral edges of the spall, regardless of the spall direction and the size of the spalled area. Three-dimensional finite element method modeling shows the river lines are caused by mixed mode I?+?III loading near the lateral edges of spall and predicts a spall depth variation near the lateral edges of spall due to mixed mode I?+?II loading. The absolute range of river lines increases with lateral size of spall, while the relative range of river lines decreases, consistent with variations in mode mixity.

36 MATERIALS SCIENCE↗

MURR LEU structural and thermal hydraulics analyses: Part I – Preliminary irradiation thermo-mechanical behavior

The University of Missouri Research Reactor (MURR) is expected to be converted from highly enriched uranium (HEU, ≥ 20 wt% U-235) U-Al x dispersion fuel to low-enriched uranium (LEU, < 20 wt% U-235) with U-10Mo monolithic fuel. This work introduces high-fidelity irradiation thermo–mechanical (T-M) analysis of the MURR LEU focusing on changes in coolant channel gap thickness. Three-dimensional (3D) finite element (FE) models were developed to simulate the irradiation T-M behavior of the MURR LEU element with all 23 curved fuel plates, the two side plates, and the combs. It was shown that channel gap thickness changes were influenced not only by plate thickness variations due to fuel swelling and creep but also by the radial displacement of consecutive MURR LEU plates. Modeling the fuel element assembly captured side plate displacements, which were shown to reduce radial fuel plate displacements towards the convex side. The maximum local radial displacement in the element was predicted at the end of life (EOL) as 23.7 mil (602.0 µm) on the lateral centerline of plate 23 towards the convex side. The maximum stripe-averaged reduction in channel gap thickness, particularly relevant for thermal hydraulics (TH) safety analysis, was calculated as 15.9 mil (403.9 µm) in single-side heated channel 24 (the outermost channel). These results account for the thermal resistance from the oxide build-up on cladding surfaces which was shown to be up to 0.82 mil (20.8 µm) thick. It was demonstrated that accounting for oxide layer thermal resistance led to a 10 °C higher peak fuel temperature and a 4.4 mil (111.8 µm) greater maximum local radial displacement. The impact of the calculated channel gap thickness changes on the MURR LEU TH safety analysis is evaluated in Part II.

Irradiation thermo-mechanical analysis↗

Anisotropic thermal diffusivity and conductivity in SiC/SiC tubes studied by infrared imaging and X-ray computed tomography

Traditional flash diffusivity evaluation of thermal diffusivity/conductivity of composite tubes require machining of specimens. For a thin-wall tube, this method can only be used to obtain through-thickness transport property. A novel method to evaluate anisotropic thermal diffusivity in a composite tube has been developed. Braided SiC/SiC composite tubes were subjected to a xenon flash heating pulse. A high-speed, high-sensitivity infrared camera was used to measure surface temperature changes as a function of time and nondestructively detect subsurface defects/damages, such as macroscopic pores. Standard reference material (Pyroceram 9606) and curved SiC/SiC composite tube specimens were used to validate thermal diffusivity obtained from infrared imaging. Unlike the traditional method, there is no need prepare special specimens, and thermal diffusivity values in three orientations are obtained after a single flash. A finite element analysis model based on x-ray computed tomography scans was developed to simulate the heat transfer. In conclusion, this technique is significant in assessing thermal conductivity and inspecting the health of ceramic tubes during and after service.

36 MATERIALS SCIENCE↗

Influence of copper content on the high temperature tensile and low cycle fatigue behavior of cast Al-Cu-Mn-Zr alloys

The relationship between low cycle fatigue (LCF) and monotonic tensile fracture strain was investigated for cast Al-Cu-Mn-Zr alloys containing 6 and 9 wt% Cu at 250 °C. The 9% Cu alloy consisted of larger size and volume fraction of brittle intermetallic grain boundary particles that fractured during tensile tests to reduce the fracture strain by 50% compared to the 6% Cu alloy. We found that LCF life was similar between the two alloys and weakly influenced by the particles. LCF life and fracture strain were poorly correlated due to the dominant influence of particles on monotonic tensile crack but not on fatigue crack.

36 MATERIALS SCIENCE↗

Directed energy deposition of radially functionally graded nickel based superalloys IN625 and IN738: Microstructure and mechanical properties

This study investigates the microstructure and mechanical properties of radially functionally graded (RFG) materials fabricated from nickel-based superalloys IN625 and IN738 using directed energy deposition (DED). RFG cylindrical deposits with radial gradients of IN625-IN738 and IN738-IN625 were successfully fabricated and confirmed to be crack-free through tailored composition grading, scan strategy and process parameters that minimized thermal gradients and solidification velocities. Microstructural analysis revealed distinct variations across the gradient, characterized by differences in dendrite arm spacing, microsegregation, and the formation of distinct γ' precipitates and carbide/boride phases (MC, M 6 C, M 2 C, and MB 2 ) in as-deposited (AD) conditions, which evolved into M 23 C 6 and M 2 B upon subsequent heat treatment (HT). EBSD showed predominant columnar grains with a [100] texture, while microhardness profiles demonstrated consistent variation across the gradients. Mechanical testing highlighted the synergy of the graded structure, combining IN738's strength with IN625's ductility, resulting in enhanced properties such as a yield strength of 1131 MPa, ultimate tensile strength of 1382 MPa, and elongation of 12 % for HT IN625-IN738 grading. Numerical simulations employing CALPHAD thermodynamic calculations and finite element thermal modeling are used to establish phase evolution and microstructural feature correlation with thermal profiles, aiding in the optimization of process-structure-property relationship.

Additive manufacturing↗

Achieving clean production with nanostructured coated milling tools dry machining low carbon steel

The advancement of clean production using nanostructured materials in subtractive manufacturing processes has focused on replacing the use of liquid lubrication with solid lubricants coated to conventional cutting tools machining low carbon steel (~ 0.2 wt. % carbon). However, little is known about the wear mechanisms that dominate such tools when coated with functional graded nanostructured coatings. The present study follows an international standard for characterizing such wear (ISO 8688) but uses advanced measurement techniques such as x-ray fluorescence technology to measure the diffusion of chemical species from tool to workpiece and vice versa, that contributes to tool wear with the view to enhance the principles of cleaner production. The results show that when commercially-available tools coated with the appropriate functionally graded industrial coating, diffusion does not take place and creates an effective chemical barrier during machining. The experimental methods used in the study include machining workpieces using a CNC milling machine, tachometer for measuring spindle speed, dynamometer for measuring cutting forces, infra-red camera for measuring cutting temperatures, chemical species diffusion wear using x-ray florescence detector, and an optical microscope for measuring tool wear. Coupled with the physical measurement of wear, it is concluded that mechanical wear dominates the milling of a low carbon steel and that thermal properties of the coatings do not correlate with such wear. It is also concluded that flank wear as a function of volume of workpiece removed and machining duration (machining index) is a sound method for assessing the transition of the stages of wear during clean machining operations. The present study not only has implication for the design of better cutting tools, but also advances dry machining processes that eliminate the use of liquid lubricants making subtractive production processes cleaner. The main conclusions drawn from this study show that coated tools have a longer life under dry conditions compared to uncoated tools and that there is no correlation between coating properties and tool wear. In conclusion, diffusion of chemical species into the cut chips did not occur because of the lack of thermally induced chemical wear of the cutting tool and that wear is caused by the gradual erosion of the flank face due to the chip abrading the cutting tool.

36 MATERIALS SCIENCE↗

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗

Strength-ductility synergy through microstructural and compositional heterogeneity in directed energy deposition additive manufacturing of face-centered cubic materials

Directed energy deposition (DED) is an additive manufacturing (AM) process based on welding technology and offers the advantages of large build volume, high deposition rate, and ability to fabricate multi-material parts. Epitaxial continuous columnar grain growth is a characteristic microstructural feature of DED processed alloys. In this study, a bamboo-like microstructure (periodic alternation of equiaxed and columnar structure) was produced by adopting an intermittent deposition strategy in 316L stainless steel and Inconel 625. The formation of a bamboo-like alternating microstructure was confirmed through electron backscattered diffraction (EBSD) analysis. Hardness mapping showed that the columnar to equiaxed transition (CET) occurred at the region right below the fusion line. A finite element (FE) model was used to investigate the relationship between the temperature gradient (G) and the solidification rate (R). The FE model showed a low G/R ratio at the region right below the interface promoting the CET. The grain size and material-dependent deformation behaviors are analyzed using digital image correlation (DIC). The lower deformation on the fine-grain regions observed in DIC analysis is attributed to a higher strain hardening rate, which is confirmed through dislocation density analysis on a tensile-interrupted specimen. The periodically alternating grain size coupled with the microstructural changes caused by intermittent deposition strategy result in a better strength-ductility synergy in both single-material and bimetallic specimens.

36 MATERIALS SCIENCE↗

Mechanical behavior analysis using small-size graphite disc compression testing with digital image correlation methods

The ASTM International Standard Test Method for Tensile Stress-Strain of Carbon and Graphite (ASTM C749) requires dog-bone shaped specimens. Unfortunately, specimen sizes required in ASTM C749 are much too large for irradiation within a test reactor and are incompatible with irradiation capsule volumes, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed as an alternative that provides a means for testing tensile properties using smaller specimens when geometric constraints are present. However, the important stress–strain relation throughout the test loading history that could be recorded using ASTM C749 is lost by using ASTM D8289 because of the difficulties involved with using traditional strain gauges or extensometers on such a small specimen. In this work, a digital image correlation (DIC) method, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 splitting tensile test on eight small samples of graphite grade Mersen 2114 to track their deformation during compressive loading. Initial results showed that the DIC technique can measure the surface displacement/strain on these small (Ø6 × 3 mm) graphite specimens with high accuracy and good repeatability, indicating that the DIC technique, combined with the small-disc splitting tensile test in ASTM D8289, can capture the strain evolution in the test. In conclusion, further comparison of measured strains with analytical and numerical simulation results provides additional understanding of the mechanical behaviour of this nuclear graphite material in complex stress states.

36 MATERIALS SCIENCE↗

A numerical and experimental approach to compare the effect of sample thickness in small in-situ SEM and large ex-situ tensile testing in Alloy 709

ASTM standards for tensile tests define specific sample size requirements regardless of grain size of the material. However, sample size requirements for testing should be considered in conjunction with the number of grains within its reduced cross-sectional area. This is particularly important for in-situ Scanning Electron Microscope (SEM) tensile tests, as they have to be conducted on smaller samples. In this study, a comprehensive experimental and numerical evaluation of the effect of specimen thickness (and the resulting number of grains within the cross-section) were conducted using in-situ SEM tensile test (on submillimeter thick samples) and ex-situ tensile tests (on samples of 0.68–5.9 mm thickness) and the results are compared with FEM simulations outcome. All tests were conducted at room temperature and the results are correlated to the number of grains within the thickness of specimens. The tensile test results indicated that even though the 0.2% proof stress and the tensile strength do not vary for the current range of thicknesses (with number of grains in the cross-section varying from 13 to 118), a difference in necking mechanism exists. Post tensile strength, thinner samples undergo more shear failure and diagonal localized necking whereas thicker samples experience more diffuse necking indicated by decreased area of shear failure at the edges and increased area of dimpled tensile failure at the center of the sample. Further, FEM results complement the experimental findings by showing the formation of conjugated localized shear bands on the upper and lower surface in the 0.68 mm thick sample and shear bands combining to form singular shear bands in thicker samples. These results also confirm the validity of the in-situ SEM tensile tests conducted on thinner samples as long as the required minimum number of grains exist within the cross-section.

36 MATERIALS SCIENCE↗

Irradiation performance of a U-7Mo in Al-Si matrix dispersion full-size fuel plate assembly

The Korea Atomic Energy Research Institute (KAERI) is leading the Ki-Jang Research Reactor (KJRR) project with the intent to develop a new reactor for medical isotope production and other nuclear research purposes. The KJRR core is designed to use high density fuel system where uranium alloyed with 7 wt% molybdenum (U-7Mo) particles are dispersed in a matrix of aluminum alloyed with 5 wt% silicon (Al-5Si) and clad in aluminum alloy 6061 (Al-6061) to form fuel plates. KAERI developed a fabrication facility to construct KJRR fuel assemblies and partnered with the Idaho National Laboratory (INL) to irradiate a full-size fuel assembly, with 21 total fuel plates, in the Advanced Test Reactor (ATR). Irradiation testing and subsequent Post Irradiation Exam (PIE) campaigns were performed successfully over a multi-year project. Monte Carlo neutronic calculations, coupling with a depletion code, were performed based on ATR’s as-run power history which showed that the highest power plate (plate 20) reached 83.1 % end-of-life (EOL) local burnup based on initial 235 U content. Additionally, finite element thermal modeling was performed based as-run power history which showed a beginning-of-life (BOL) peak local heat flux of 184 W/cm 2 . No anomalous fuel performance was observed during the irradiation and target test conditions were achieved. PIE showed favorable performance of the fuel assembly regarding all important phenomena. This paper describes the KJRR fuel assembly irradiation conditions and PIE data to support the conclusion that it performed well, without evidence of unexpected or problematic fuel performance, within an irradiation test designed to bound the KJRR design environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗