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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shrinkage-induced deformations and creep of structural concrete: 1-year measurements and numerical prediction

Highlights: • Extensive experimental study on drying shrinkage, creep and microcracking of concrete • All specimens prepared from a single batch of ordinary-strength structural concrete • 1st year of measurements of (not only) non-uniformly drying beams with span up to 3-m • The database is downloadable from free-to-use research data repository • Modified MPS model for concrete creep used in blind prediction of all experiments The material models for creep and shrinkage operating on the material point level in FEM are usually intended for challenging complex applications and structures, where the average cross-sectional approach does not suffice. The identification of the growing number of material parameters induced by increasing model capabilities relies on very specific and narrow-oriented yet interconnected experiments which are scarce. The presented comprehensive experiments aim to provide a clearer image of the complicated interaction among the basic phenomena: drying, shrinkage, creep, and microcracking. The cornerstone of this ongoing research is a unique set of 30 partially-sealed unreinforced concrete beams with span 1.75–3.0 m subjected to drying. To minimize material variation, all specimens in this study were cast from a single batch of ordinary strength structural concrete with slag-blended binder. The resulting experimental database will be suitable both for validation and development of the constitutive models.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Constitutive modeling for isotropic materials

The first year of progress on a NASA-Lewis contract with the General Electric Co is documented. The purpose of this contract (NAS3-23927) is to develop and evaluate unified constitutive equations for applications to hot-path components of aircraft gas turbine engines such as high pressure turbine blades and vanes. To accomplish this goal, uniaxial, notched, and multiaxial specimens made of conventionally cast Rene 80 are being tested under conditions that simulate engine operating conditions. To reduce the raw data, automated data reduction techniques are being developed that produce computer files containing the information needed to analyze proposed constitutive theories. Described are the analytical methods being developed to determine the parameters for these nonlinear unified theories by using the reduced data files. In another activity, a dedicated finite-element computer code is being developed to use unified theories in the structural analysis of hot-section components. This code was extensively verified for one such theory by successfully predicting the strain histories measured experimentally at the notch root of complex specimens taken from complex laboratory specimens.

Ramaswamy, V. G.↗
DOE OSTI2020

Imaging transient solidification behavior

Solidification processing offers the first opportunity to control microstructure, properties, and performance in metallic alloy components. Until recently, microstructural evaluations were limited to post-solidification characterization by destructive methods. We review the development of time-resolved, in situimaging techniques capable of capturing solid-liquid interfacial evolution in metallic alloys with high spatial and temporal resolution underdiverse solidification conditions relevant for applications ranging from conventional directional solidification, crystal growth, and casting,to welding and additive manufacturing. These experiments enable direct visualization of transient behaviors that would otherwise remain unknown, uniquely providinginsights into the physics that impact microstructure and defect development and strategies for microstructuralcontrol and defect mitigation. Understanding microstructural evolution and the characteristics that form under various solidification conditions is essential for the development of multiscale, experimentallyinformed predictive modeling. This is highlightedby solidification simulations that utilize in situmeasurements of solidification dynamics from state-of-the-art experimental techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Characterization of compositional variation and solidification condition on rapid solidification behavior of beta stabilized titanium

Rapid solidification (RS) is a field of study that has increased in importance due to the rise of additive manufacturing (AM) techniques. The breakdown of equilibrium at the solid-liquid interface results in an inability to use traditional solidification models that have been applied to great success in casting. As such, there is a current need for a high-throughput testing method to screen alloys at rapid solidification conditions without the significant cost for powder atomization. In this work, twin hammer splat quenching (SQ) is utilized and theoretically developed as a technique uniquely placed to test alloy systems for AM applications. Using analytical thermal spray models and numerical heat transfer simulations, the cooling rate and solidification conditions of splat quenching are compared across several parameters to evaluate the key parameters for generating cooling rates relevant to AM process. Once a framework for the cooling rates and relevancy of SQ are developed, SQ was used to process alloy compositions in the Ti5553 alloy space (AMS 7026) to study the formation of detrimental secondary phases occurring upon rapid quenching from the melt. A martensitic phase is found to occur at the lower end of the standard, but the composition listed in the standard were found to be safe. Finally, the rapid solidification behavior of Ti-5553 is evaluated by comparing splat quenching, laser welding, and plasma atomization that offer a wide window of solidification conditions. The segregation and solidification structure are analyzed and compared to models to provide a deeper understanding of the solidification behavior of beta stabilized titanium alloys.

Williamson, Christopher J.↗
DOE OSTI2019

The Influence of Segregation of Mn on the Recrystallization Behavior of C-Mn Steels

The influence of Mn segregation that occurs during casting on recrystallization kinetics has been explored for a C-Mn automotive steel. A homogenization heat treatment was used to remove Mn segregation, while maintaining a similar initial austenite grain size to the segregated condition, to provide base-line comparison data. Deformation trials were carried out in a Gleeble thermomechanical simulator to 0.3 strains at 900 °C, and significant differences in both recrystallization times and final grain sizes were seen between the two conditions. The segregated sample incurred a longer recrystallization time and a wider recrystallized grain size distribution than the homogenized sample. A finite element model was used to predict the local deformation strains based on variations in material properties caused by segregation; the model used data from XRF elemental mapping of Mn to define the local mechanical properties based on differences in solid solution strengthening. Local strain variations ranging from 0.19 to 0.49 were predicted in the segregated sample for an overall applied 0.3 strain (compared to 0.26 to 0.35 in the homogenized sample); using these strains to predict the recrystallized grain sizes gave a much better prediction for the grain size distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design of Silicide-Strengthened Nb–Si–Cr–(Mo) alloys for additive manufacturing

Three high–intermetallic volume Nb–Si–Cr–(Mo) alloys were designed using CALPHAD modeling with the goal of identifying high–specific strength, oxidation-resistant alloys that can be additively manufactured using powder bed fusion. The silicides Nb 5 Si 3 and Nb 9 Si 2 Cr3 were targeted as the primary strengthening phases, and the addition of Cr promoted the NbCr 2 phase. These alloys were cast and surface-processed with electron beam welding at different speeds to simulate additive manufacturing, and the phases and microstructures of both cast and welded regions were characterized. The weld processing was found to produce fine-grained microstructures in each alloy with fine-scale intermetallics uniformly distributed among a body-centered cubic Nb matrix. Microstructural refinement and hardness were found to increase with weld velocity; one alloy reached its highest hardness of approximately 16 GPa before the brittleness at higher velocities became detrimental. One alloy was found to be qualitatively the least brittle while also attaining a hardness of 13 GPa and was therefore identified as a good candidate for additive manufacturing.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Probabilistic analysis for fatigue strength degradation of materials

This report presents the results of the first year of a research program conducted for NASA-LeRC by the University of Texas at San Antonio. The research included development of methodology that provides a probabilistic treatment of lifetime prediction of structural components of aerospace propulsion systems subjected to fatigue. Material strength degradation models, based on primitive variables, include both a fatigue strength reduction model and a fatigue crack growth model. Linear elastic fracture mechanics is utilized in the latter model. Probabilistic analysis is based on simulation, and both maximum entropy and maximum penalized likelihood methods are used for the generation of probability density functions. The resulting constitutive relationships are included in several computer programs, RANDOM2, RANDOM3, and RANDOM4. These programs determine the random lifetime of an engine component, in mechanical load cycles, to reach a critical fatigue strength or crack size. The material considered was a cast nickel base superalloy, one typical of those used in the Space Shuttle Main Engine.

Royce, Lola↗
NASA NTRS1994

Shoemaker-Levy 9 and the tidal disruption of comets

The break-up of Periodic Comet Shoemaker-Levy 9 into multiple pieces following its grazing encounter with Jupiter in July 1992 can be used to study tidally-induced fracture in comets. This spectacular event allows us not only to set limits on the size, strength and density of Shoemaker-Levy 9 itself, but provides invaluable guidance to numerical modeling of such encounters. In an extensive treatment of tidal breakup which assumed self-gravitating, homogeneous, perfectly elastic bodies, Dobrovolskis derived simple analytical expressions for the tidally-induced surface and central stresses. Both can be cast in such a way that Poisson's ratio is the only material dependent constant entering these expressions. Whether both surface and central failure must be initiated as a criterion for breakup, or either one of them is sufficient, remains a subject of disagreement. To resolve this debate, we model the details of cometary breakup using a three-dimensional Smooth Particle Hydrodynamics (SPH) code modified to simulate fracture in small solid objects. At lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution.

Benz, W.↗
DOE OSTI2021

Improved Gamma-Ray Point Source Quantification in Three Dimensions by Modeling Attenuation in the Scene

We discuss using a series of detector measurements taken at different locations to localize a source of radiation as a well-studied problem. The source of radiation is sometimes constrained to a single point-like source, in which case the location of the point source can be found using techniques such as maximum likelihood. Recent advancements have shown the ability to locate point sources in 2D and even 3D, but few have studied the effect of intervening material on the problem. In this work we examine gamma-ray data taken from a freely moving system and develop voxelized 3-D models of the scene using data from the onboard LiDAR. Ray casting is used to compute the distance each gamma ray travels through the scene material, which is then used to calculate attenuation assuming a single attenuation coefficient for solids within the geometry. Parameter estimation using maximum likelihood is performed to simultaneously find the attenuation coefficient, source activity, and source position that best match the data. Using a simulation, we validate the ability of this method to reconstruct the true location and activity of a source, along with the true attenuation coefficient of the structure it is inside, and then we apply the method to measured data with sources and find good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Current-driven plasma acceleration versus current-driven energy dissipation. III - Anomalous transport

In the present paper the linear stability description and weak turbulence theory are used to develop a second order description of wave-particle transport and anomalous dissipation. The goal is to arrive at anomalous transport coefficients that can be readily included in fluid flow codes. In particular, expressions are derived for the heating rates of ions and electrons by the unstable waves and for the electron-wave momentum exchange rate that controls the anomalous resistivity effect. Comparative calculations were undertaken assuming four different saturation models: ion trapping, electron trapping, ion resonance broadening, and thermodynamic bound. A foremost finding is the importance of the role of electron Hall parameter in scaling the level of anomalous dissipation for the parameter range of the MPD thruster plasma. Polynomial expressions of the relevant transport coefficients cast solely in terms of macroscopic parameters are also obtained for inclusion in plasma fluid codes for the self-consistent numerical simulation of real thruster flows including microturbulent effects.

Choueiri, Edgar Y.↗
DOE OSTI2019

Computational Modeling of Fracture Toughness of Al-Si and Al-Zn-Mg-Cu Alloys with Detected Porosity

This article presents the results of an experimental and numerical analysis of the fracture resistance of cast aluminum alloys subjected to variable loads over time based on the fracture mechanics method. The subjects of the analysis were the AlSi9 and AlZn6MgCu alloys (usually used for plastic processing) in which the presence of porosities in the form of separate or dispersed voids was revealed experimentally. The material for the experimental procedures was obtained directly from a casting or from casting samples. Depending on the type of the defect detected, different numerical models were proposed. The results of the numerical calculations were verified experimentally. Static tensile tests, fatigue low-cycle tests, fracture mechanics tests, fatigue crack propagation tests, a microscopic examination of the sample fracture area, and an x-ray tomographic analysis were all performed. A correlation of the results between the numerical analysis and the experimental work has been found. On the basis of the proposed calculation models, the fatigue life of castings of the analyzed alloys can be predicted in the case of concentrated and dispersed defects.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗
NASA NTRS1984

Heat transfer distributions around nominal ice accretion shapes formed on a cylinder in the NASA Lewis icing research tunnel

Local heat transfer coefficients were obtained on irregular cylindrical shapes which typify the accretion of ice on circular cylinders in cross flow. The ice shapes were grown on a 5.1 cm (2.0 in.) diameter cylinder in the NASA Lewis Icing Research Tunnel. The shapes were 2, 5, and 15 min accumulations of glaze ice and 15 min accumulation of rime ice. Heat transfer coefficients were also measured around the cylinder with no ice accretion. These icing shapes were averaged axially to obtain a nominal shape of constant cross section for the heat transfer tests. Heat transfer coefficients around the perimeter of each shape were measured with electrically heated copper strips embedded in the surface of the model which was cast from polyurethane foam. Each strip contained a thermocouple to measure the local surface temperature. The models were run in a 15.2 x 68.6 cm (6 x 27 in.) wind tunnel at several velocities. Background turbulence in the wind tunnel was less than 0.5 percent. The models were also run with a turbulence producing grid which gave about 3.5 percent turbulence at the model location with the model removed. The effect of roughness was also simulated with sand grains glued to the surface. Results are presented as Nusselt number versus angle from the stagnation line for the smooth and rough models for both high and low levels of free stream turblence. Roughness of the surface in the region prior to flow separation plays a major role in determining the heat transfer distribution.

Vanfossen, G. J.↗
NASA NTRS1985

A three-dimensional incompressible flow simulation method and its application to the Space Shuttle main engine. II Turbulent flow

An implicit finite difference code cast in general curvilinear coordinates is further developed for three-dimensional incompressible turbulent flows. The code is based on the method of pseudocompressibility and utilizes the Beam and Warming implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is further extended to include composite-grids to overcome the excessive computer memory required for solving turbulent flows in complex three-dimensional geometries. A simple turbulence model is proposed for internal flows. The code is being used for the Space Shuttle Main Engine (SSME) internal flow analyses.

Chang, J. L. C.↗
DOE OSTI2017

Model Uncertainties for Valencia RPA Effect for MINERvA

This technical note describes the application of the Valencia RPA multi-nucleon effect and its uncertainty to QE reactions from the GENIE neutrino event generator. The analysis of MINERvA neutrino data in Rodrigues et al. PRL 116 071802 (2016) paper makes clear the need for an RPA suppression, especially at very low momentum and energy transfer. That published analysis does not constrain the magnitude of the effect; it only tests models with and without the effect against the data. Other MINERvA analyses need an expression of the model uncertainty in the RPA effect. A well-described uncertainty can be used for systematics for unfolding, for model errors in the analysis of non-QE samples, and as input for fitting exercises for model testing or constraining backgrounds. This prescription takes uncertainties on the parameters in the Valencia RPA model and adds a (not-as-tight) constraint from muon capture data. For MINERvA we apply it as a 2D ($q_0$,$q_3$) weight to GENIE events, in lieu of generating a full beyond-Fermi-gas quasielastic events. Because it is a weight, it can be applied to the generated and fully Geant4 simulated events used in analysis without a special GENIE sample. For some limited uses, it could be cast as a 1D $Q^2$ weight without much trouble. This procedure is a suitable starting point for NOvA and DUNE where the energy dependence is modest, but probably not adequate for T2K or MicroBooNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1989

Some Navier-Stokes calculations for the CAST-10 airfoil

Results of the simulation of the viscous flow past the CAST-10 airfoil were shown for different flow conditions. Since the experiments provide only surface pressures and force coefficients, the comparison to the numerical results relies on these. Good agreement of these results is found for the lower Mach number cases except for the shock position. As numerical experiments indicate, this seems to be due to the turbulent shock boundary layer interaction which is not correctly modelled by the algebraic turbulence model employed. For the lower Mach number case the influence of the transition location was investigated, too. Changing the transition location at the lower surface has much more influence on the pressure distribution than changing it on the upper side. For the higher Mach number case the double shock structure found in the experiment for the lower Reynolds numbers was not reproduced by the numerical solutions. The reason for this is unknown though it may be due to the turbulence modelling. For the higher Reynolds number a better resolution of the boundary layer is needed in the computation in order to recover the experimental pressure plateau; but then the shock position is still found downstream of the experimental one.

Schwamborn, D.↗