SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE 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↗

NASA NTRS1993

Feedback control for unsteady flow and its application to the stochastic Burgers equation

The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.

Choi, Haecheon↗

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↗

NASA NTRS1984

Effect of mild atherosclerosis on flow resistance in a coronary artery casting of man

An in-vitro flow study was conducted in a mildly atherosclerotic main coronary artery casting of man using sugar-water solutions simulating blood viscosity. Steady flow results indicated substantial increases in pressure drop, and thus flow resistance at the same Reynolds number, above those for Poiseuille flow by 30 to 100 percent in the physiological Reynolds number range from about 100 to 400. Time-averaged pulsatile flow data showed additional 5 percent increases in flow resistance above the steady flow results. Both pulsatile and steady flow data from the casting were found to be nearly equal to those from a straight, axisymmetric model of the casting up to a Reynolds number of about 200, above which the flow resistance of the casting became gradually larger than the corresponding values from the axisymmetric model.

Back, L. H.↗

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↗

NASA NTRS2016

Method of Modeling and Simulation of Shaped External Occulters

The present invention relates to modeling an external occulter including: providing at least one processor executing program code to implement a simulation system, the program code including: providing an external occulter having a plurality of petals, the occulter being coupled to a telescope; and propagating light from the occulter to a telescope aperture of the telescope by scalar Fresnel propagation, by: obtaining an incident field strength at a predetermined wavelength at an occulter surface; obtaining a field propagation from the occulter to the telescope aperture using a Fresnel integral; modeling a celestial object at differing field angles by shifting a location of a shadow cast by the occulter on the telescope aperture; calculating an intensity of the occulter shadow on the telescope aperture; and applying a telescope aperture mask to a field of the occulter shadow, and propagating the light to a focal plane of the telescope via FFT techniques.

Lyon, Richard G.↗

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↗

NASA NTRS1997

A Study of the Role of Clouds in the Relationship Between Land Use/Land Cover and the Climate and Air Quality of the Atlanta Area

The goal of Project ATLANTA is to derive a better scientific understanding of how land cover changes associated with urbanization affect local and regional climate and air quality. Clouds play a significant role in this relationship. Using GOES images, we found that in a 63-day period (5 July-5 September 1996) there were zero days which were clear for the entire daylight period. Days which are cloud-free in the morning become partly cloudy with small cumulus clouds in the afternoon in response to solar heating. This result casts doubt on the applicability of California-style air quality models which run in perpetual clear skies. Days which are clear in the morning have higher ozone than those which are cloudy in the morning. Using the RAMS model, we found that urbanization increases the skin surface temperature by about 1.0-1.5 C on average under cloudy conditions, with an extreme of +3.5 C. Clouds cool the surface due to their shading effect by 1.5-2.0 C on average, with an extreme of 5.0 C. RAMS simulates well the building stage of the cumulus cloud field, but does poorly in the decaying phase. Next year's work: doing a detailed cloud climatology and developing improved RAMS cloud simulations.

Kidder, Stanley Q.↗

NASA NTRS1988

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗