Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Casting Modeling Simulation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Modeling and simulation framework for dynamic strain localization in elasto-viscoplastic metallic materials subject to large deformations

This study describes a theoretical and computational framework for the treatment of adiabatic shear band formation in rate-sensitive polycrystalline metallic materials. From a computational perspective, accurate representation of strain localization behavior has been a long-standing challenge. In addition, the underlying physical mechanisms leading to the localization of plastic deformation are still not fully understood. The proposed framework is built around an enhanced-strain finite element formulation, designed to alleviate numerical pathologies known to arise in localization problems, by allowing a localization band of given finite width (weak discontinuity) to be embedded within individual elements. The mechanical threshold strength (MTS) model is used to represent the temperature and strain rate-dependent viscoplastic response of the material. This classical flow stress model employs an internal state variable to quantify the effect of dislocation structure evolution (work hardening and recovery). In light of growing evidence suggesting that the softening effect of dynamic recrystallization may play a significant role, alongside thermal softening, in the process of shear band formation and growth, a simple dynamic recrystallization model is proposed and cast within the context of the MTS model with the aid of the aforementioned internal state variable. An initiation criterion for shear localization in rate and temperature-sensitive materials is introduced and used in the present context of high-rate loading, where material rate-dependence is pronounced and substantial temperature increases are achieved due to the dissipative nature of viscoplastic processes. In addition, explicit time integration is adopted to facilitate treatment of the dynamic problems under consideration, where strain rates in excess of 10 4 s –1 are typically attained. Two series of experiments are conducted on AISI 316L stainless steel, employing the commonly used top-hat sample geometry and the Split-Hopkinson Pressure Bar dynamic test system. Axi-symmetric finite element simulation results are compared to cross-sectional micrographs of recovered samples and experimental load–displacement results, in order to examine the performance of the proposed framework and demonstrate its effectiveness in treating the initiation and growth of adiabatic shear banding in dynamically loaded metallic materials. These comparisons demonstrate that thermal softening alone is insufficient to induce shear localization behaviors observed in some materials, such as stainless steel, and support the hypothesis that dynamic recrystallization and/or other softening mechanisms play an essential role in this process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Effects on image quality of a 2D antiscatter grid in x-ray digital breast tomosynthesis: Initial experience using the dual modality (x-ray and molecular) breast tomosynthesis scanner

Purpose: Radiation scattered from the breast in digital breast tomosynthesis (DBT) causes image degradation, including loss of contrast between cancerous and background tissue. Unlike in 2-dimensional (2D) mammography, an antiscatter grid cannot readily be used in DBT because changing alignment between the tube and detector during the scan would result in unacceptable loss of primary radiation. However, in the dual modality breast tomosynthesis (DMT) scanner, which combines DBT and molecular breast tomosynthesis, the tube and detector rotate around a common axis, thereby maintaining a fixed tube-detector alignment. This C-arm geometry raises the possibility of using a 2D (cellular) focused antiscatter grid. The purpose of this study is to assess change in image quality when using an antiscatter grid in the DBT portion of a DMT scan under conditions of fixed radiation dose. Methods: Two 2D focused prototype grids with 80 cm focal length were tested, one stack-laminated from copper (Cu) and one cast from a tungsten-polymer (W-poly). They were reciprocated using a motion scheme designed to maximize transmission of primary x-ray photons. Grid-in and grid-out scatter-to-primary ratios (SPRs) were measured for rectangular blocks of material simulating 30%, 50%, and 70% glandular tissue compositions. For assessment of changes in image quality through the addition of a grid, the Computerized Imaging Reference Systems, Inc., phantom Model 011A containing a set of 1 cm thick blocks simulating a range of glandular/adipose ratios from 0/100 to 100/0 was used. To simulate 6.5 and 8.5 cm thick compressed breasts, 1 cm thick slices of PMMA were added to the Model 011A phantom. DBT images were obtained with and without the grid, with exposure parameters fixed for a given compressed thickness. Signal-difference-to-noise ratios (SDNRs), contrast, and voxel value-based attenuation coefficients (μ) were measured for all blocks from reconstructed phantom images. Results: For 4, 6, and 8 cm tissue-equivalent block phantom thicknesses, the inclusion of the W-poly grid reduced the SPR by factors of 5, 6, and 5.8, respectively. For the same thicknesses, the copper grid reduced the SPR by factors of 3.9, 4.5, and 4.9. For the 011A phantom, the W-poly grid raised the SDNR of the 70/30 block from 0.8, −0.32, and −0.72 to 0.9, 0.76, and 0.062 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. It raised the SDNR of the 100/0 block from 3.78, 1.95, and 1.0 to 3.79, 3.67, and 3.25 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. Inclusion of the W-poly grid improved the accuracy of image-based μ values for all block compositions. However, smearing of attenuation across slices due to limited angular sampling decreases the sensitivity of voxel values to changing composition compared to theoretical μ values. Conclusions: Under conditions of fixed radiation dose to the breast, use of a 2D focused grid increased contrast, SDNR, and accuracy of estimated attenuation for mass-simulating block compositions in all phantom thicknesses tested, with the degree of improvement depending upon material composition. A 2D antiscatter grid can be usefully incorporated in DBT systems that employ fully isocentric tube-detector rotation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Robust metal-ceramic coaxial cable sensors for distributed temperature monitoring in fossil energy power systems

The project aims to develop a new type of low cost, robust metal-ceramic coaxial cable (MCCC) Fabry-Pérot interferometer (FPI) sensor and demonstrate the capability of cascading a series of FPI sensing points in a single MCCC for real-time distributed monitoring of temperature up to 1000°C. The technical accomplishments of this project include: (i) development of metal and ceramic materials for fabricating the MCCC suitable for use in high temperature (up to 1000oC) fossil fuel, (ii) design, construction and test of single point MCCC-FPI sensors for high temperature measurement at deviation within ±2.7%, (iv) demonstration of the MCCC-FPI sensors for real-time distributed temperature measurement with desired sensitivity, spatial resolution, stability, and response speed that are important to practical applications, and (v) development of a novel MCCC-FPI sensor platform for measuring high frequency dielectric constants for solids and liquids as function of temperature. Single point MCCC-FPI sensors are fabricated based on the FPI parameters, including the relationship between reflector width and reflectance as well as the inter-reflector distance for a FPI, estimated by model calculations. MCCC-FPI sensors have been made using insulators of packed ceramic powders, ceramic tubes, and air or inert gases; outer conductors of various metal tubes; inner conductors of wires of the same metal materials as the tubes used; and reflectors of dense ceramic discs, casted ceramic discs and air gaps. The insulator, conductor and reflector materials are investigated for their thermal, chemical, and structural stabilities under temperatures up to 1000°C in simulated coal-derived syngas. Based on the experimental findings, alumina tube and discs, and stainless steel tube and wire are selected for reflector, insulator and outer and inner conductors, respectively. Air or inert gases are found to be useful for reflectors (in the form of narrow groves) or insulators (in the form of empty chamber cavity between ceramic discs) in construction of MCCC-FPI sensors. Single point MCCC-FPI sensors have been fabricated using stainless steel as conductors, dense α-alumina tubes as dielectric insulators, and two air gaps in the insulator as reflectors. The MCCC-FPI sensor with air gap reflectors has been demonstrated for measuring temperature between 200 and 500°C using microwave as sensing signal in a frequency range of 2 – 8 GHz. The temperature measurement is achieved by monitoring the frequency shift (Δƒ) of a selected resonant peak of the microwave interferogram reflected from the pair of air gap reflectors. The MCCC-FPI sensor exhibited excellent linear temperature-dependence for Δƒ, good accuracy (deviation ±2.7%), and reasonable response speed (<180 s for stabilization). Thermal drifting of the sensing signal was also rather minimal for the MCCC-FPI sensors after operating for more than one week at around 500oC. MCCC-FPI sensors with each containing 2 or 3 pairs of reflectors (i.e. 2 or 3 FPIs in a single cable) have been fabricated and demonstrated for distributed temperature measurement between 250 and 550°C in a microwave frequency range of 3-6 GHz. The measurement is achieved by monitoring the frequency shift of a selected resonant peak of the reflected microwave interferogram, which is obtained by the join-time-domain signal processing method. The results demonstrate that the FPI sensing points in both the three-point and two-point MCCC sensors can measure distributed high temperature with good reproducibility and consistency in multiple runs of measurement tests. The sensor signal displays excellent linear dependence between temperature and the corresponding frequency shift. The 2-point and 3-point MCCC-FPI sensors also exhibit fast response to environmental temperature changes. A three-point MCCC-FPI sensor has been further examined for temperature measurement up to 1000°C for multiple heating-cooling cycles and the results prove that the sensor is suitable for distributed temperature measurement up to 1000°C. MCCC-FPI sensors with more than 4 FPIs point in a single line are made in two structures: the first uses alumina tube as insulation and air gaps as reflectors; and the second uses air chamber as insulation and complete circular alumina discs (2 mm-thick) as reflectors. MCCC-FPIs in neither of the two structures are able to exceed four FPIs in a cable because the reflections from the FPIs are too strong for the microwave to transmit through >4 FPIs over a long distance. Specially shaped alumina discs made by removing parts of the circular discs are found to decrease the microwave reflectance because of reduced reflection areas. However, the smallest reflectance (in terms of percentage of the incident signal) obtained in this work is as high as 7.5%, which is still too large to construct more than 10 FPIs in a single MCCC. Nevertheless, the feasibility of the proposed multipoint MCCC-FPI has been demonstrated using air insulation and 2 mm-thick low-ε r Teflon discs as reflectors on a 2 m-long cable containing 10 FPIs (with reflectance <5% on each FPI). These results indicate that the contrast of dielectric constant between the insulation and reflector has a more significant impact on the reflectance. Thus, future research efforts should also be directed to developing ceramic materials (dense or porous) with low dielectric constant values which, when used as insulation or reflector, can realize large numbers of FPIs in a single cable over a long distance. To assist the development of high temperature dielectric materials, a novel MCCC-FPI sensor platform is developed and successfully demonstrated for measuring high frequency dielectric constants for fluids, packed dielectric solids, and ceramic tube materials and determining ε r dependencies on temperature and chemical composition, etc..

36 MATERIALS SCIENCE↗
DOE OSTI2025

Macro-micro multiscale modeling to assist the design of HPDC Al castings microstructure and alloys for EV super-large body structures (Phase 1)

Implementation of High Pressure Die Casting (HPDC) Aluminum (Al) body structures for high volume electrified vehicles (EV) to improve electric efficiency remains a key strategy within many original equipment manufacturer (OEM)s. In addition to high strength for safety requirements, superior Self-Piercing Riveting (SPR) performance is demanded for HPDC Al alloys to be compatible with high volume SPR joining. In this work, it is proposed to extend and validate an existing Contractor finite element multiscale macro-micro modeling approach to quantify the influence of the microstructure of HPDC alloys on the fracture strain/displacement under 3-point bend and clinch testing. The success of this work will allow to replace solution treatment stage with low energy consumption heat treatment (HT) processes, or to design new non heat treatable (NHT) HPDC Al alloys to eliminate HT requirements. Ultimately, this project will facilitate the application of HPDC Al alloys for super-large vehicle structures to significantly reduce vehicle weight, and thus improving energy efficiency. The purpose of this project is to extend and validate an existing finite element code, which is based on the Contractor developed macro-micro multi-scale modeling approach, to numerically simulate the three-point bending and clinch test and study the influences of material microstructural characteristics and phase properties on the rivetability. The macro-micro modeling approach begins with a sample scale model and identify the location, which is mostly prone to failure, the deformation history of the boundaries of that location calculated will be used to drive a microstructure-based sub-models where the material microstructure and microscale properties are considered. Using this approach, the wrap-bending failure for two Al alloys are correctly predicted for the first time. This will start with phase I effort of building a framework of macro-micro three point bending test and clinch test of Al10SiMgMn HPDC alloy in the as-cast and T7 heat treated conditions. Those results will then be validated with experimental test results. The phase II effort will involve the utilization of the knowledge learned in phase I to establish the quantitative correlation between the microstructure characteristics and the riveting performance, which will be further used to guide the optimization of HPDC Al alloy microstructure using heat treatment process to achieve sufficient rivetability to join large thin-wall HPDC alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Greenhouse Gas Source Attribution: Measurements Modeling and Uncertainty Quantification

In this project we have developed atmospheric measurement capabilities and a suite of atmospheric modeling and analysis tools that are well suited for verifying emissions of green- house gases (GHGs) on an urban-through-regional scale. We have for the first time applied the Community Multiscale Air Quality (CMAQ) model to simulate atmospheric CO 2 . This will allow for the examination of regional-scale transport and distribution of CO 2 along with air pollutants traditionally studied using CMAQ at relatively high spatial and temporal resolution with the goal of leveraging emissions verification efforts for both air quality and climate. We have developed a bias-enhanced Bayesian inference approach that can remedy the well-known problem of transport model errors in atmospheric CO 2 inversions. We have tested the approach using data and model outputs from the TransCom3 global CO 2 inversion comparison project. We have also performed two prototyping studies on inversion approaches in the generalized convection-diffusion context. One of these studies employed Polynomial Chaos Expansion to accelerate the evaluation of a regional transport model and enable efficient Markov Chain Monte Carlo sampling of the posterior for Bayesian inference. The other approach uses de- terministic inversion of a convection-diffusion-reaction system in the presence of uncertainty. These approaches should, in principle, be applicable to realistic atmospheric problems with moderate adaptation. We outline a regional greenhouse gas source inference system that integrates (1) two ap- proaches of atmospheric dispersion simulation and (2) a class of Bayesian inference and un- certainty quantification algorithms. We use two different and complementary approaches to simulate atmospheric dispersion. Specifically, we use a Eulerian chemical transport model CMAQ and a Lagrangian Particle Dispersion Model - FLEXPART-WRF. These two models share the same WRF assimilated meteorology fields, making it possible to perform a hybrid simulation, in which the Eulerian model (CMAQ) can be used to compute the initial condi- tion needed by the Lagrangian model, while the source-receptor relationships for a large state vector can be efficiently computed using the Lagrangian model in its backward mode. In ad- dition, CMAQ has a complete treatment of atmospheric chemistry of a suite of traditional air pollutants, many of which could help attribute GHGs from different sources. The inference of emissions sources using atmospheric observations is cast as a Bayesian model calibration problem, which is solved using a variety of Bayesian techniques, such as the bias-enhanced Bayesian inference algorithm, which accounts for the intrinsic model deficiency, Polynomial Chaos Expansion to accelerate model evaluation and Markov Chain Monte Carlo sampling, and Karhunen-Lo %60 eve (KL) Expansion to reduce the dimensionality of the state space. We have established an atmospheric measurement site in Livermore, CA and are collect- ing continuous measurements of CO 2 , CH 4 and other species that are typically co-emitted with these GHGs. Measurements of co-emitted species can assist in attributing the GHGs to different emissions sectors. Automatic calibrations using traceable standards are performed routinely for the gas-phase measurements. We are also collecting standard meteorological data at the Livermore site as well as planetary boundary height measurements using a ceilometer. The location of the measurement site is well suited to sample air transported between the San Francisco Bay area and the California Central Valley.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Carbide Particle Redistribution in U-10Mo Alloy during Hot and Cold Rolling Processes: Synthetic Microstructure-Based Finite-Element Approach

The influence of initial carbide particle distribution features on particle fracture and the evolution of stringer distribution within an alloy of uranium with 10 weight percent molybdenum during cold rolling was investigated in this study. This work will provide guidance to casting optimization work and identify actionable process changes to eliminate the surface features directly attributed to carbides. In a previous study, various synthetic microstructure-based finite-element (FE) models were generated with consideration of different parameters (i.e., particle shape, particle volume fraction, and particle distribution), and used to examine the effects of hot rolling on the microstructure evolution. In this study, microstructures after ~80% hot rolling were used as the initial microstructures for the subsequent cold rolling simulation. Plane-strain compression loading was employed in the FE models to simulate cold rolling with reductions from ~80% to ~94%. Stringer analyses with the consideration of particle fracture were then performed with the simulated cold rolling results. The results of void evolution analyses show that microstructures with elliptical particles incorporate a much larger (almost double) void volume fraction than those with circular particles, indicating that the elliptical particles are fractured more easily than the circular ones. The results from the fracture mode analyses show that, under the adopted simple compression loading assumption, almost all the particles (i.e., >99%) were fractured by tensile separation. The results of stringer analyses show that, during cold rolling, the stringer volume fractions increase substantially due to the particle fracture, especially for the microstructures with the elliptical particles, and that the stringer volume fractions from stringers comprising many particles tend to increase faster. The maximum stringer lengths for the microstructures with circular particles tend to remain constant or decrease during cold rolling, whereas those for elliptical particles tend to slightly increase after a slight decrease. The number of particles at the surface increases very rapidly due to particle fracture during cold rolling.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Radiation behavior of high-entropy alloys for advanced reactors. Final report

In the first task, we have demonstrated the radiation damage and the recrystallization behaviors in multicomponent alloys through molecular-dynamics simulations. It is found that by alloying with atoms of different sizes, the atomic-level strain increases, and the propensity of the radiation-induced crystalline to amorphous transition increases as the defects cluster in the cascade body. Recrystallization of the radiation induced supercooled or glass regions show that by tuning the composition and the equilibrium temperature, the multicomponent alloys can be healed. The crystalline-amorphous-crystalline transitions predict the potential high radiation resistance in multicomponent alloys. In the second task, three types of high-entropy alloys (HEAs) were fabricated from AlCoCrFeNi and AlCuCrFeNi quinary alloys. Hardness and reduced contact modulus were measured using nanoindentation tests. Heavy ion irradiation were performed using 10 MeV gold and 5 MeV nickel to study radiation effects. Al 0.5 CrCuFeNi 2 shows phase separation upon the presence of copper. Both hardness and contact modulus exhibit the same trend as increasing the applied load, and it indicates that excessive free volume may alter the growth rate of the plastic zone. The as-cast Al 0.1 CoCrFeNi specimen undergone the hot isostatic pressing (HIP) process and steady cooling rate which mitigate the quenching effect. The swelling behavior was characterized by the atomic force microscopy (AFM), and the swelling rate is approximately 0.02% dpa. Selected area diffraction (SAD) patters show irradiation-induced amorphization throughout the ion projected range. Within the peak damage region, an amorpous ring is observed, and a mixture of amorphous/ crystalline structure at deeper depth is found. The Al 0.3 CoCrFeNi HEAs shows good radiation resistance up to 60 peak dpa. No voids or dislocations are observed. The crystal structures remain face-centered-cubic (FCC) before and after 5 MeV Ni irradiation. Higher dpa might be required to study defects formation mechanisms. In the third task, all the constituent binary and ternary systems of the Al-Co-Cr-Fe-Ni system were thermodynamically modeled within the whole composition range. Comparisons between the calculated phase diagrams and literature data are in good agreement. The multi-component thermodynamic database of the Al-Co-Cr-Fe-Ni system was then obtained via extrapolation. The current Al-Co-Cr-Fe-Ni thermodynamic database enables us to carry out the calculations of phase diagrams, which can be used as useful guidelines to identify the Al-Co-Cr-Fe-Ni HEAs with desirable microstructures. In the fourth task, we discuss how as-cast and homogenized phases can be identified, what phases are usually found in the as-cast and homogenized conditions, and what the thermodynamics and kinetics of phase transformations are in the AlCoCrFeNi HEA. The microstructure and phase composition were studied in as-cast and homogenized conditions. It showed the dendritrical structure in the as-cast condition consisting primarily of a nano-lamellar mixture of A2 [disorder body-centered-cubic (BCC)] and B2 (ordered BCC) phases, in addition to a very small amount of A1 [disorder face-centered-cubic (FCC)] phases. The homogenization heat treatment resulted in an increase in the volume fraction of the A1 phase and formation of a Sigma phase. Tensile properties in as-cast and homogenized conditions are reported at 700 °C. Thermodynamic modeling of non-equilibrium and equilibrium phase diagrams for the AlCoCrFeNi HEA gave good agreement with the experimental observations of the phase contents. The reasons for the improvement of ductility after the heat treatment are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Comment on: “Solving many-body Schrödinger equations with kinetic energy partition method”, Ann. Phys. 388 (2018) 54–68 by Y-H. Chen and S.D. Chao

We discuss a recent test of the performance of the kinetic energy partition method (KEP) through its application to two separable quantum-mechanical models. We argue that one of the benchmark models is exceedingly simple for testing any realistic approximate method and that almost any reasonable approach yields better results. In the second example our exact benchmark eigenvalues disagree considerably with those chosen by the authors for comparison, which casts doubts on the accuracy of their KEP approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗
DOE OSTI2022

Minimizing thickness variation in monolithic U-10Mo fuel foil and Zr interlayer during hot rolling: A microstructure-based finite element method analysis

Low-enriched uranium alloyed with 10 wt. % molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium fuel for the United States’ high performance research reactors. The monolithic U-10Mo fuel plate consists of a metallic U-10Mo fuel foil with a 25 µm Zr interlayer and a relatively thick cladding of aluminum alloy 6061. The Zr interlayer is typically applied during the hot co-rolling process, and this process dictates the uniformity of the Zr interlayer. Thickness variation observed in the U-10Mo and Zr interlayer has been attributed to several sources: the initial grain size of the U-10Mo castings, can materials, rolling temperature, inhomogeneous molybdenum content, and porosity in the cast U-10Mo. This thickness variation limits the ability to meet the dimensional specification; thus, a better understanding of the factors causing the nonuniform thickness is needed. In this work, we used a novel, microstructure-based finite element method to model the hot rolling process to address these concerns. Grain microstructures in U-10Mo were tessellated and explicitly considered in the finite element model. Each grain was assigned a random material property to mimic the grain strength variations induced by different grain orientations. Simulations were performed using six steel can thicknesses, four grain sizes, and with or without a Zr interlayer to investigate the influences of those variables on the thickness nonuniformity. The simulation results showed that a thinner steel can and finer U-10Mo grain size reduce thickness variations in both the U-10Mo fuel foil and Zr interlayer. The direct findings from the simulations and analysis can be used to optimize the hot rolling schedule, reduce fabrication defects, and meet the dimensional specifications. The proposed microstructure-based finite element model can be also coupled with experimental microstructure characterization data, images, and models to simulate multi-pass hot rolling.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A Comprehensive Analysis of Macrosegregation Formation During Twin-Roll Casting

A two-phase Eulerian–Eulerian volume-averaged model is used to predict the formation of macrosegregation during the twin-roll casting of inoculated Al-4 wt pct Cu alloys. For low solid fractions, the equiaxed crystals are modeled according to the submerged object approach. However, above a given transition limit, they are assumed to behave like a viscoplastic material. This means that the solid phase behaves as a coherent structure that can influence the motion of the liquid. Such approach allows one to take into account the flow dynamics arising from the occurrence of both solidification and hot rolling simultaneously, which usually occurs in twin-roll casting. Therefore, it is possible to explain the origin of the macrosegregation patterns obtained in the simulations based on the relative motion between the phases. Compression-induced expulsion of segregated melt is observed as a result of the deformation of the solidifying shells. Such occurrence leads to a negative macrosegregation region in the outer part of the as-cast strip. Then, because the solute-enriched melt is squeezed out toward adjacent regions, two positively segregated bands can be found near the center of the domain. Furthermore, it is shown how solidification-induced feeding weakens the absolute value of the negative and positive segregation bands.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The discrete-continuum connection in dislocation dynamics: I. Time coarse graining of cross slip

A recent continuum dislocation dynamics formalism (Xia and El-Azab 2015 Model. Simul. Mater. Sci. Eng. 23 055009) has been enriched by incorporating an improved cross slip model. 3D discrete dislocation dynamics simulations were used to collect cross slip rate data in the form of time series that were analysed to estimate the correlation time for cross slip, which was subsequently used as a time scale for local window averaging of the collected cross slip rate data. This time averaging filters out the cross slip rate fluctuations over time intervals less than the correlation time, thus resulting in relatively smoother time series for the cross slip rates. The coarse grained series were further cast in the form of smooth trends with superposed fluctuations and implemented in continuum dislocation dynamics simulations using a Monte Carlo scheme. Here, this approach resulted in a significant improvement of the predicted stress–strain response and a more realistic dislocation cell structure evolution. The similitude law for the average cell size evolution with inverse of stress, however, remains unaffected by the cross slip rates used in continuum dislocation dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

The Radial Distribution of Mono-metallicity Populations in the Galactic Disk as Evidence for Two-phase Disk Formation

Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [α/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin- and thick-disk dichotomy. The analysis of these observations led to the non-[α/Fe] enhanced populations splitting into MMPs with different surface densities according to their [Fe/H]. By contrast, [α/Fe] enhanced (i.e., old) populations show a homogeneous behavior. We analyze these results in the wider context of disk formation within non-isolated halos embedded in the Cosmic Web, resulting in a two-phase mass assembly. By performing hydrodynamical simulations in the context of the ΛCDM model, we have found that the two phases of halo mass assembly (an early fast phase, followed by a slow phase with low mass-assembly rates) are very relevant to determine the radial structure of MMP distributions, while radial mixing only plays a secondary role, depending on the coeval dynamical and/or destabilizing events. Indeed, while the frequent dynamical violent events occuring at high redshift remove metallicity gradients and imply efficient stellar mixing, the relatively quiescent dynamics after the transition keeps [Fe/H] gaseous gradients and prevents newly formed stars from suffering strong radial mixing. By linking the two-component disk concept with the two-phase halo mass-assembly scenario, our results set halo virialization (the event marking the transition from the fast to the slow phases) as the separating event that marks periods that are characterized by different physical conditions under which thick- and thin-disk stars were born.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

MCTP system model based on linear programming optimization of apertures obtained from sequencing patient image data maps

Purpose: The authors present a hybrid direct multileaf collimator (MLC) aperture optimization model exclusively based on sequencing of patient imaging data to be implemented on a Monte Carlo treatment planning system (MC-TPS) to allow the explicit radiation transport simulation of advanced radiotherapy treatments with optimal results in efficient times for clinical practice. Methods: The planning system (called CARMEN) is a full MC-TPS, controlled through aMATLAB interface, which is based on the sequencing of a novel map, called “biophysical” map, which is generated from enhanced image data of patients to achieve a set of segments actually deliverable. In order to reduce the required computation time, the conventional fluence map has been replaced by the biophysical map which is sequenced to provide direct apertures that will later be weighted by means of an optimization algorithm based on linear programming. A ray-casting algorithm throughout the patient CT assembles information about the found structures, the mass thickness crossed, as well as PET values. Data are recorded to generate a biophysical map for each gantry angle. These maps are the input files for a home-made sequencer developed to take into account the interactions of photons and electrons with the MLC. For each linac (Axesse of Elekta and Primus of Siemens) and energy beam studied (6, 9, 12, 15 MeV and 6 MV), phase space files were simulated with the EGSnrc/BEAMnrc code. The dose calculation in patient was carried out with the BEAMDOSE code. This code is a modified version of EGSnrc/DOSXYZnrc able to calculate the beamlet dose in order to combine them with different weights during the optimization process. Results: Three complex radiotherapy treatments were selected to check the reliability of CARMEN in situations where the MC calculation can offer an added value: A head-and-neck case (Case I) with three targets delineated on PET/CT images and a demanding dose-escalation; a partial breast irradiation case (Case II) solved with photon and electron modulated beams (IMRT + MERT); and a prostatic bed case (Case III) with a pronounced concave-shaped PTV by using volumetric modulated arc therapy. In the three cases, the required target prescription doses and constraints on organs at risk were fulfilled in a short enough time to allow routine clinical implementation. The quality assurance protocol followed to check CARMEN system showed a high agreement with the experimental measurements. Conclusions: A Monte Carlo treatment planning model exclusively based on maps performed from patient imaging data has been presented. The sequencing of these maps allows obtaining deliverable apertures which are weighted for modulation under a linear programming formulation. The model is able to solve complex radiotherapy treatments with high accuracy in an efficient computation time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Predicting Metallic Defects in Casting BrA9ZhZL Alloy Tool Electrodes

A method is proposed for modeling the lost-wax casting of tool electrodes by means of fusible BrA9ZhZL copper-based alloy. The functionality of the casting system is analyzed. The heating and cooling of a ceramic mold is taken into account in visualizing the temperature fields. Sections where shrinkage pores may form in the casting are identified.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the Correspondence Between Atmosphere-Only and Coupled Simulations for Radiative Feedbacks and Forcing From CO 2

Atmosphere-only experiments are widely used to investigate climate feedbacks simulated in more computationally expensive fully coupled global climate model simulations. We confirm that this remains a valid approach by comparing the radiative feedbacks and forcing between coupled and atmosphere-only simulations for the latest models taking part in the 6th phase of the Coupled Model Intercomparison Project (CMIP6). For global-mean cloud feedbacks, we find a better than previously known correspondence between these experiments, which applies even to the response of individual cloud properties (amount, altitude, and optical depth) and holds even when considering atmosphere-only simulations of only 1 yr duration. For regional cloud feedbacks, the correspondence between the two experiments is generally present at every geographic location except for the tropical Pacific but takes longer experiments to reveal. For the lapse rate and surface albedo feedbacks, the correspondence between the two experiments is weaker due to the non-uniform warming pattern and loss of sea ice in the coupled experiment. For the across-model relationship between 4xCO 2 effective radiative forcing and feedback, we find a different behavior across experiments in CMIP6 than in CMIP5, casting doubt on the physical significance of previous results that highlighted an anti-correlation between the two quantities. Overall, these results confirm the utility of atmosphere-only experiments particularly to study cloud feedbacks, which are the dominant source of inter-model spread in climate sensitivity.

58 GEOSCIENCES↗
DOE OSTI2016

Quantifying the Economic and Grid Reliability Impacts of Improved Wind Power Forecasting

Wind power forecasting is an important tool in power system operations to address variability and uncertainty. Accurately doing so is important to reducing the occurrence and length of curtailment, enhancing market efficiency, and improving the operational reliability of the bulk power system. This research quantifies the value of wind power forecasting improvements in the IEEE 118-bus test system as modified to emulate the generation mixes of Midcontinent, California, and New England independent system operator balancing authority areas. To measure the economic value, a commercially available production cost modeling tool was used to simulate the multi-timescale unit commitment (UC) and economic dispatch process for calculating the cost savings and curtailment reductions. To measure the reliability improvements, an in-house tool, FESTIV, was used to calculate the system's area control error and the North American Electric Reliability Corporation Control Performance Standard 2. The approach allowed scientific reproducibility of results and cross-validation of the tools. A total of 270 scenarios were evaluated to accommodate the variation of three factors: generation mix, wind penetration level, and wind fore-casting improvements. The modified IEEE 118-bus systems utilized 1 year of data at multiple timescales, including the day-ahead UC, 4-hour-ahead UC, and 5-min real-time dispatch. The value of improved wind power forecasting was found to be strongly tied to the conventional generation mix, existence of energy storage devices, and the penetration level of wind energy. The simulation results demonstrate that wind power forecasting brings clear benefits to power system operations.

17 WIND ENERGY↗