Search NASA⌕ Search

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.

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

At least 55 records · Page 3

Triple Plate Mold Final Report: Optimization of the Mold Design and Casting Parameters for a Thin U-10mo Fuel Casting

This work describes the experiments and modeling that have been performed to improve and try to optimize the simultaneous casting of three plates of U-10wt%Mo in a single coil vacuum induction melting (VIM) furnace. The plates of interest are 280 mm wide by 203 mm tall by 5 mm thick (11" x 8" x 0.2"). The initial mold design and processing parameters were supplied by Y-12. The mold and casting cavity were instrumented with a number of thermocouples, and the casting performed to determine the thermal history of the mold and casting. The resulting cast plates were radiographed and numerous defects identified. Metallography was performed to help identify the nature of the radiographically observed defects. This information was then used to validate a mold filling and solidification model of that casting. Based on the initial casting, good casting design practice, and process simulation of several design alternatives, a revised design was developed with the goal of minimizing casting defects such as porosity. The redesigned mold had a larger hot-top and had its long axis along the horizontal direction. These changes were to try to develop a strong thermal gradient conducive to good feeding and minimization of micro- and macroporosity in the cast plates. An instrumented casting was then performed with the revised mold design and a linear distributor. This design yielded cast plates with significantly less radiographically identified defects. Unfortunately, there was significant variation in plate weight and metal content in their hot-tops. Fluid flow simulations were then performed on this mold/distributor design. This helped identify the issue with this linear distributor design. Additional simulations were then performed on candidate distributor redesigns and a preferred distributor annular design was identified. This improved annular design was used to produce a third instrumented casting with favorable results. These refined designs and their radiographic characterization are compared to the initial design.

36 MATERIALS SCIENCE↗

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of Dendritic Evolution of Continuously Cast Steel Billet with Cellular Automaton

In order to predict the dendritic evolution during the continuous steel casting process, a simple mechanism to connect the heat transfer at the macroscopic scale and the dendritic growth at the microscopic scale was proposed in the present work. As the core of the across-scale simulation, a two-dimensional cell automaton (CA) model with a decentered square algorithm was developed and parallelized. Apart from nucleation undercooling and probability, a temperature gradient was introduced to deal with the columnar-to-equiaxed transition (CET) by considering its variation during continuous casting. Based on the thermal history, the dendritic evolution in a 4 mm × 40 mm region near the centerline of a SWRH82B steel billet was predicted. The influences of the secondary cooling intensity, superheat, and casting speed on the dendritic structure of the billet were investigated in detail. The results show that the predicted equiaxed dendritic solidification of Fe-5.3Si alloy and columnar dendritic solidification of Fe-0.45C alloy are consistent with in situ experimental results [Yasuda et al. Int J Cast Metals Res 22:15–21 (2009); Yasuda et al. ISIJ Int 51:402–408 (2011)]. Moreover, the predicted dendritic arm spacing and CET location agree well with the actual results in the billet. The primary dendrite arm spacing of columnar dendrites decreases with increasing secondary cooling intensity, or decreasing superheat and casting speed. Meanwhile, the CET is promoted as the secondary cooling intensity and superheat decrease. However, the CET is not influenced by the casting speed, owing to the adjusting of the flow rate of secondary spray water. Compared with the superheat and casting speed, the secondary cooling intensity can influence the cooling rate and temperature gradient in deeper locations, and accordingly exerts a more significant influence on the equiaxed dendritic structure.

36 MATERIALS SCIENCE↗

Transitional flow in thin tubes for space station freedom radiator

A two dimensional finite volume method is used to predict the film coefficients in the transitional flow region (laminar or turbulent) for the radiator panel tubes. The code used to perform this analysis is CAST (Computer Aided Simulation of Turbulent Flows). The information gathered from this code is then used to augment a Sinda85 model that predicts overall performance of the radiator. A final comparison is drawn between the results generated with a Sinda85 model using the Sinda85 provided transition region heat transfer correlations and the Sinda85 model using the CAST generated data.

Loney, Patrick↗

Application of a Near Infrared Imaging System for Thermographic Imaging of the Space Shuttle during Hypersonic Re-Entry

High resolution calibrated near infrared (NIR) imagery was obtained of the Space Shuttle s reentry during STS-119, STS-125, and STS-128 missions. The infrared imagery was collected using a US Navy NP-3D Orion aircraft using a long-range infrared optical package referred to as Cast Glance. The slant ranges between the Space Shuttle and Cast Glance were approximately 26-41 nautical miles at point of closest approach. The Hypersonic Thermodynamic Infrared Measurements (HYTHIRM) project was a NASA Langley led endeavor sponsored by the NASA Engineering Safety Center, the Space Shuttle Program Office and the NASA Aeronautics Research Mission Directorate to demonstrate a quantitative thermal imaging capability. HYTHIRM required several mission tools to acquire the imagery. These tools include pre-mission acquisition simulations of the Shuttle trajectory in relationship to the Cast Glance aircraft flight path, radiance modeling to predict the infrared response of the Shuttle, and post mission analysis tools to process the infrared imagery to quantitative temperature maps. The spatially resolved global thermal measurements made during the Shuttle s hypersonic reentry provides valuable flight data for reducing the uncertainty associated with present day ground-to-flight extrapolation techniques and current state-of-the-art empirical boundary-layer transition or turbulent heating prediction methods. Laminar and turbulent flight data is considered critical for the development of turbulence models supporting NASA s next-generation spacecraft. This paper will provide the motivation and details behind the use of an upgraded NIR imaging system used onboard a Navy Cast Glance aircraft and describe the characterizations and procedures performed to obtain quantitative temperature maps. A brief description and assessment will be provided of the previously used analog NIR camera along with image examples from Shuttle missions STS-121, STS-115, and solar tower test. These thermal observations confirmed the challenge of a long-range acquisition during re-entry. These challenges are due to unknown atmospheric conditions, image saturation, vibration etc. This provides the motivation for the use of a digital NIR sensor. The characterizations performed on the digital NIR sensor included radiometric, spatial, and spectral measurements using blackbody radiation sources and known targets. An assessment of the collected data for three Space Shuttle atmospheric re-entries, STS-119, STS-125, and STS-128, are provided along with a description of various events of interest captured using the digital NIR imaging system such as RCS firings and boundary layer transitions. Lastly the process used to convert the raw image counts to quantitative temperatures is presented along with comparisons to the Space Shuttle's onboard thermocouples.

Zalameda, Joseph N.↗

Distributed Sand-Grain Roughness Effects on Blunt Body Hypersonic Transition and Heating

An experimental investigation of distributed sand-grain surface roughness effects on boundary-layer transition and convective heating has been performed. Two representative entry vehicle geometries, a spherical-cap aeroshell and a sphere-cone aeroshell, were considered. Multiple cast ceramic wind tunnel models of each geometry were fabricated with various roughness heights to simulate an ablated thermal protection system. Wind tunnel testing was performed at Mach 6 over a range of Reynolds numbers sufficient to produce laminar, transitional, and turbulent flow. Aeroheating and boundary-layer transition onset data were obtained using global phosphor thermography. The experimental heating data are presented herein, as are comparisons to laminar and turbulent smooth-wall heat transfer distributions from computational flow field simulations.

Brian R. Hollis↗

Mechanism of Macrosegregation Formation in Continuous Casting Slab: A Numerical Simulation Study

Solidified shell bulging is supposed to be the main reason for slab center segregation, while the influence of thermal shrinkage rarely has been considered. In this article, a thermal shrinkage model coupled with the multiphase solidification model is developed to investigate the effect of the thermal shrinkage, solidification shrinkage, grain sedimentation, and thermal flow on solute transport in the continuous casting slab. In this model, the initial equiaxed grains contract freely with the temperature decrease, while the coherent equiaxed grains and columnar phase move directionally toward the slab surface. The results demonstrate that the center positive segregation accompanied by negative segregation in the periphery zone is mainly caused by thermal shrinkage. During the solidification process, liquid phase first transports toward the slab surface to compensate for thermal shrinkage, which is similar to the case considering solidification shrinkage, and then it moves opposite to the slab center near the solidification end. It is attributed to the sharp decrease of center temperature and the intensive contract of solid phase, which cause the enriched liquid to be squeezed out. With the effect of grain sedimentation and thermal flow, the negative segregation at the external arc side (zone A1) and the positive segregation near the columnar-to-equiaxed transition at the inner arc side (position B1) come into being. Besides, it is found that the grain sedimentation and thermal flow only influence solute transport before equiaxed grains impinge with each other, while the solidification and thermal shrinkage still affect solute redistribution in the later stage.

36 MATERIALS SCIENCE↗

Scale-Adaptive Simulation of Transient Two-Phase Flow in Continuous-Casting Mold

Scale-adaptive simulation (SAS) of the transient gas–liquid two-phase flow in a laboratory-scale continuous-casting mold is presented. The main objective is to investigate the applicability of the scale-adaptive unsteady Reynolds-averaged Navier–Stokes turbulent model (URANS SAS) for predicting the transient multiscale turbulent structures in a two-phase flow. Good quantitative agreements with the experimental data and the large eddy simulation (LES) results are obtained both for the time-averaged velocity field and for the transient turbulent characteristics. The introduction of the von Karman length-scale into the turbulence-scale equation allows the SAS model to dynamically adjust to the resolved turbulent structures. The LES-like pulsating behavior of the air gas and the large-scale liquid eddy magnitudes in the unsteady regions of flow field are captured by the SAS model. The classical − 5/3 law of power spectrum density (PSD) of the axial velocity is kept properly for the single-phase turbulent flow. For two-phase flow, the decay of PSD is too steep at the high-frequency region; the predicted PSD obtained with SAS is damped stronger than that estimated by LES. The SAS model offers an attractive alternative to the existing LES approach or to the other hybrid RANS/LES models for strongly unsteady flows.

36 MATERIALS SCIENCE↗

Modeling and Computation of Casting Process by Particle Method

Since particle-based computational methods, such as smoothed particle hydrodynamics (SPH) method, are suitable to treat free surfaces and moving boundaries, it is expected that the simulation with particle methods can represent more accurate molten metal flow. Hence, we apply the SPH method to some casting processes and validate the numerical results by comparing with experiments. While fluid behaviors of the transportation and the mold filling processes are expressed sufficiently, ladle pouring process of aluminum alloy is not reproduced by the existing scheme. So we developed a new numerical model of surface oxide film, and our numerical results show good agreement with the experimental data.

36 MATERIALS SCIENCE↗

Experimental Investigation of Block-TPS Fence/Gap Roughness Effects on Transition Onset and Turbulent Heating at Mach 6

An experimental investigation of boundary-layer transition and convective surface heating on wind tunnel models that simulated aeroshells comprised of multiple blocks of thermal protection system (TPS) material has been performed. Two representative entry vehicle geometries, a sphere-cone aeroshell and a sphere-cone aeroshell, were considered. Multiple cast ceramic models of each geometry were fabricated with varying fence heights or gap depths to represent gap fillers that either protrude above, or are recessed below, the TPS blocks due to differential ablation during reentry. Wind tunnel testing was performed at Mach 6 over a range of Reynolds numbers sufficient to produce laminar, transitional, and turbulent flow. Convective surface heating and boundary-layer transition onset data were obtained using global phosphor thermography. The experimental heating data are presented herein, as are comparisons to laminar and turbulent smooth-wall heat transfer distributions from computational flow field simulations.

Aerothermodynamics↗

Application of Al-Cu-W-Ta graded density impactors in dynamic ramp compression experiments

Graded density impactors (GDIs) are utilized to dynamically compress materials to extreme conditions. Two modifications to a previously developed Mg-Cu-W GDI are made in this work before using it in a dynamic compression experiment: Mg is replaced with Al and a Ta disk is glued to the back. The Mg phase is replaced by Al because FCC Al remains solid to higher pressure along its Hugoniot compared to Mg. The addition of the Ta disk creates a constant particle velocity regime and facilitates a definition of peak pressure states. Microstructure analysis, profilometry, and ultrasonic C-scans of the Al-Cu-W GDI all confirm excellent uniformity. We evaluated signal variation in the radial direction of a dynamically compressed Al-LiF bilayer target to evaluate the contribution of spatial nonuniformity to errors. Velocity traces from five photon Doppler velocimetry (PDV) probes located at different radial distances from the center of the target varied at most by 1.1% with a root mean square of 0.3% during the compression ramp, demonstrating low PDV measurement error over a relatively large experimental area. The experimental PDV data additionally agrees well with 1D simulations that use inputs from predictive characterization models developed for the material properties resulting from tape casting, laminating, and powder consolidation processes. Low measurement error during quasi-isentropic compression, leading to better precision, promises a robust platform to reach extreme compression and low-temperature recovery states and facilitates discovery via synthesis, quenching, and preservation of new high-pressure phases.

36 MATERIALS SCIENCE↗

Parallelizing Timed Petri Net simulations

The possibility of using parallel processing to accelerate the simulation of Timed Petri Nets (TPN's) was studied. It was recognized that complex system development tools often transform system descriptions into TPN's or TPN-like models, which are then simulated to obtain information about system behavior. Viewed this way, it was important that the parallelization of TPN's be as automatic as possible, to admit the possibility of the parallelization being embedded in the system design tool. Later years of the grant were devoted to examining the problem of joint performance and reliability analysis, to explore whether both types of analysis could be accomplished within a single framework. In this final report, the results of our studies are summarized. We believe that the problem of parallelizing TPN's automatically for MIMD architectures has been almost completely solved for a large and important class of problems. Our initial investigations into joint performance/reliability analysis are two-fold; it was shown that Monte Carlo simulation, with importance sampling, offers promise of joint analysis in the context of a single tool, and methods for the parallel simulation of general Continuous Time Markov Chains, a model framework within which joint performance/reliability models can be cast, were developed. However, very much more work is needed to determine the scope and generality of these approaches. The results obtained in our two studies, future directions for this type of work, and a list of publications are included.

Nicol, David M.↗

Experimental Investigation of Hexcomb-Pattern Roughness Effects on Transition Onset and Turbulent Heating Augmentation at Mach 6

An experimental investigation of hexcomb-patterned surface roughness effects on boundary-layer transition and turbulent heating has been performed. Two representative entry vehicle geometries, a spherical- cap aeroshell and a sphere-cone aeroshell, were considered. Multiple cast ceramic models of each geometry were fabricated with varying roughness pattern densities and depths that simulated an ablated hexcomb-structure thermal protection system. Wind tunnel testing was performed at Mach 6 over a range of Reynolds numbers sufficient to produce laminar, transitional, and turbulent flow. Aeroheating and boundary-layer transition onset data were obtained using global phosphor thermography. The experimental heating data are presented herein, as are comparisons to laminar and turbulent smooth-w ll heat transfer distributions from computational flow field simulations.

Hollis, Brian R.↗

Hexcomb Pattern Roughness Effects on Blunt Body Transition and Heating

An experimental investigation of hexcomb pattern surface roughness effects on boundary-layer transition and convective heating has been performed. Two representative entry vehicle geometries, a spherical-cap aeroshell and a sphere-cone aeroshell, were considered. Multiple cast ceramic models of each geometry were fabricated with various roughness pattern densities and depths that simulated an ablated hexcomb-structure thermal protection system. Wind tunnel testing was performed at Mach 6 over a range of Reynolds numbers sufficient to produce laminar, transitional, and turbulent flow. Aeroheating and boundary-layer transition onset data were obtained using global phosphor thermography. The experimental heating data are presented herein, as are comparisons to laminar and turbulent smooth-wall heat transfer distributions from computational flow field simulations. The experimental data were used to develop a unified boundary-layer transition correlation for both sandgrain distributed roughness and hexcomb pattern roughness.

Brian R Hollis↗

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aeroservoelastic modeling and applications using minimum-state approximations of the unsteady aerodynamics

Various control analysis, design, and simulation techniques for aeroelastic applications require the equations of motion to be cast in a linear time-invariant state-space form. Unsteady aerodynamics forces have to be approximated as rational functions of the Laplace variable in order to put them in this framework. For the minimum-state method, the number of denominator roots in the rational approximation. Results are shown of applying various approximation enhancements (including optimization, frequency dependent weighting of the tabular data, and constraint selection) with the minimum-state formulation to the active flexible wing wind-tunnel model. The results demonstrate that good models can be developed which have an order of magnitude fewer augmenting aerodynamic equations more than traditional approaches. This reduction facilitates the design of lower order control systems, analysis of control system performance, and near real-time simulation of aeroservoelastic phenomena.

Tiffany, Sherwood H.↗

Formation Mechanism of Surface Crack in Low Pressure Casting of A360 Alloy

A surface crack defect is normally found in low pressure castings of Al alloy with a sudden contraction structure. To further understand the formation mechanism of the defect, the mold filling process is simulated by a two-phase flow model. The experimental results indicate that the main reason for the defect deformation is the mismatching between the height of liquid surface in the mold and pressure in the crucible. In the case of filling, a sudden contraction structure with an area ratio smaller than 0.5 is obtained, and the velocity of the liquid front increases dramatically with the influence of inertia. Meanwhile, the pressurizing speed in the crucible remains unchanged, resulting in the pressure not being able to support the height of the liquid level. Then the liquid metal flows back to the crucible and forms a relatively thin layer solidification shell on the mold wall. With the increasing pressure in the crucible, the liquid level rises again, engulfing the shell and leading to a surface crack. As the filling velocity is characterized by the damping oscillations, surface cracks will form at different heights. The results shed light on designing a suitable pressurizing speed for the low pressure casting process.

36 MATERIALS SCIENCE↗

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗