Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “material flow analysis”

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 1,099 records · Page 61

Finite Element Modeling of a Cylindrical Contact Using Hertzian Assumptions

The turbine blades in the high-pressure fuel turbopump/alternate turbopump (HPFTP/AT) are subjected to hot gases rapidly flowing around them. This flow excites vibrations in the blades. Naturally, one has to worry about resonance, so a damping device was added to dissipate some energy from the system. The foundation is now laid for a very complex problem. The damper is in contact with the blade, so now there are contact stresses (both normal and tangential) to contend with. Since these stresses can be very high, it is not all that difficult to yield the material. Friction is another non-linearity and the blade is made out of a Nickel-based single-crystal superalloy that is orthotropic. A few approaches exist to solve such a problem and computer models, using contact elements, have been built with friction, plasticity, etc. These models are quite cumbersome and require many hours to solve just one load case and material orientation. A simpler approach is required. Ideally, the model should be simplified so the analysis can be conducted faster. When working with contact problems determining the contact patch and the stresses in the material are the main concerns. Closed-form solutions for non-conforming bodies, developed by Hertz, made out of isotropic materials are readily available. More involved solutions for 3-D cases using different materials are also available. The question is this: can Hertzian1 solutions be applied, or superimposed, to more complicated problems-like those involving anisotropic materials? That is the point of the investigation here. If these results agree with the more complicated computer models, then the analytical solutions can be used in lieu of the numerical solutions that take a very long time to process. As time goes on, the analytical solution will eventually have to include things like friction and plasticity. The models in this report use no contact elements and are essentially an applied load problem using Hertzian assumptions to determine the contact patch dimensions.

Knudsen, Erik↗

Liquid Bismuth Propellant Flow Sensor

Quantifying the propellant mass flow rate in liquid bismuth-fed electric propulsion systems has two challenging facets. First, the flow sensors must be capable of providing a resolvable measurement at propellant mass flow rates on the order of 10 mg/see with and uncertainty of less that 5%. The second challenge has to do with the fact that the materials from which the flow sensors are fabricated must be capable of resisting any of the corrosive effects associated with the high-temperature propellant. The measurement itself is necessary in order to properly assess the performance (thrust efficiency, Isp) of thruster systems in the laboratory environment. The hotspot sensor[I] has been designed to provide the bismuth propellant mass flow rate measurement. In the hotspot sensor, a pulse of thermal energy (derived from a current pulse and associated joule heating) is applied near the inlet of the sensor. The flow is "tagged" with a thermal feature that is convected downstream by the flowing liquid metal. Downstream, a temperature measurement is performed to detect a "ripple" in the local temperature associated with the passing "hotspot" in the propellant. By measuring the time between the upstream generation and downstream detection of the thermal feature, the flow speed can be calculated using a "time of flight" analysis. In addition, the system can be calibrated by measuring the accumulated mass exiting the system as a-function of time and correlating this with the time it takes the hotspot to convect through the sensor. The primary advantage of this technique is that it doesn't depend on an absolute measurement of temperature but, instead, relies on the observation of thermal features. This makes the technique insensitive to other externally generated thermal fluctuations. In this paper, we describe experiments performed using the hotspot flow sensor aimed at quantifying the resolution of the sensor technology. Propellant is expelled onto an electronic scale to provide an independent measure of the propellant mass flow rate as a function of time. In addition, two separate detection schemes are employed. The first uses a thermocouple to directly measure temperature in the fluid. The second involves the ,use of a fiber optic coupled to a photodiode allowing for detection of an increase in light emission from the fluid as the hotspot passes. the detection location.

Polzin, Kurt A.↗

Cooling Duct Analysis for Transpiration/Film Cooled Liquid Propellant Rocket Engines

The development of a low cost space transportation system requires that the propulsion system be reusable, have long life, with good performance and use low cost propellants. Improved performance can be achieved by operating the engine at higher pressure and temperature levels than previous designs. Increasing the chamber pressure and temperature, however, will increase wall heating rates. This necessitates the need for active cooling methods such as film cooling or transpiration cooling. But active cooling can reduce the net thrust of the engine and add considerably to the design complexity. Recently, a metal drawing process has been patented where it is possible to fabricate plates with very small holes with high uniformity with a closely specified porosity. Such a metal plate could be used for an inexpensive transpiration/film cooled liner to meet the demands of advanced reusable rocket engines, if coolant mass flow rates could be controlled to satisfy wall cooling requirements and performance. The present study investigates the possibility of controlling the coolant mass flow rate through the porous material by simple non-active fluid dynamic means. The coolant will be supplied to the porous material by series of constant geometry slots machined on the exterior of the engine.

Micklow, Gerald J.↗

Investigation of Lunar-Inspired Geopolymer Concrete Formulations Mixed and Cured in Microgravity on the International Space Station (ISS)

The research outlined in this presentation investigates the use of various lunar regolith simulants in geopolymer lunar concrete mixes mixed and cured on the International Space Station (ISS). The motivation for this work is to study the effects of gravity on the microstructure of alkali-activated materials cured with heat, and to develop materials for the construction of long-term infrastructure on the lunar surface with in-situ resource utilization (ISRU). ISRU for construction materials reduces the cost and mass of payloads related to lunar construction. The advantage of geopolymer concrete as opposed to traditional portland cement concrete is that water acts as a medium for the polymerization reaction and leaves the system throughout the process, reducing its demand. Twelve samples of lunar regolith simulant and a solution composed of sodium hydroxide and sodium silicate were sent to the ISS. The three simulants were OPRH2N, OPRL2N, and JSC-1AF, using only particles less than 53 µm in diameter to increase reactivity of the simulant. Simulant to solution ratios were determined by workability while mixing. The simulant and solution were sealed Burst Pouches® along with 2 other sealed bags to prevent material from leaking. Crew member F-14 conducted testing on the ISS by introducing the solution to the simulant in the Burst Pouch®, mixing the sample with a spatula, and then clamping the specimen in the fresh state to prevent flow inside the Burst Pouch®. These specimens were then put in a thermos heated to 80C via sealed drinking water bags to cure for 24 hours with a temperature logger. The cured specimens remained in microgravity for at least 28 days and were returned from the ISS in February 2025. The specimens were then brought to the NASA Marshall Space Flight Center (MSFC) to analyze. Material characterization consisted of conducting Micro-CT tests of entire samples in their sealed apparatus to a resolution of 25µm. 2D image slices were saved in each orthogonal direction of each specimen at a 0.03 mm step size from the 3D model to conduct analytical porosity calculations. Representative samples from each specimen were sampled to perform helium gas pycnometery and were then mounted in resin for SEM imaging, EDS, and nanoindentation. Porosity was analyzed analytically using micromechanics modelling with the assistance of the NASA Multiscale Analysis Tool (NASMAT), as well as the NASA Advanced Supercomputing (NAS) servers (V. Saseendran & N. Yamamoto, 2024). Density was measured using helium gas pycnometery and was then compared to the theoretical density for experimental porosity calculation. Due to the samples’ non-uniform shape being cured in a pouch, traditional compression and tensile strength testing could not be performed. Nanoindentation was conducted at Clarkson University to determine the microhardness and reduced modulus of elasticity. Results from flight samples can be compared to ground samples currently in DLR’s possession to determine the effect on microstructure from being mixed and cured in microgravity. This study gives further insight and understanding of geopolymer lunar concrete and its viability as a lunar construction material with ISRU.

Adam Johnson↗

A Tabletop Tool for Modeling Life Support Systems

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations.

Ramachandran, N.↗

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗

Space Laboratory on a Table Top: A Next Generative ECLSS design and diagnostic tool

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale-time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations. Preliminary experimental and modeling work in this area will be presented. This involves testing of a single inlet-exit model with detailed 3-D flow visualization and quantitative diagnostics and computational modeling of the system.

Ramachandran, N.↗

Single Degree-of-Freedom Modeling of SLS Liquid Hydrogen Pre-Valve Flow Guide to Enable Rapid Transient Analysis

A unique single degree-of-freedom approximation technique has been developed to enable rapid application of a temporally-defined multi-spectral semi-narrow-band loading for generation of realistic stress/cycle values compared to a resonant analysis. The technique uses the harmonic analysis at resonance of a high-fidelity finite element model to produce a transfer function, which is then used to calibrate the response of the SDOF model. A standard numerical ordinary differential equation solver is then used to obtain the temporal response, and its histogram is used in a fatigue/fracture model. This technique is related to other SDOF methods used widely in industry, such as Miles' Equation and the Shock Response Spectra, but it is unique in that it produces a realistic time history of the response. The most obvious error in the process, which is the effect of closely-spaced modes, was also assessed using the parallel application of several SDOF models, and the error is shown to be small. The application of this unique and tractable reduced-order methodology has enabled the SLS program to avoid substantial cost and schedule penalties if a redesign or change of material were required. It has also enabled quick analysis of a number of other structures undergoing the same or similar excitation fields, and quick assessment when the excitation and structural configuration has been altered due to design changes in the system.

Brown, Andrew M.↗

First lithium experiments in HIDRA and evidence of helium retention during quasi-steady-state stellarator plasma operations

Recent experiments in Hybrid Illinois Device for Research and Applications (HIDRA) have had operational discharges between t discharge = 60 and 1000 s using electron cyclotron resonant heating (ECRH) of the plasma. This means that quasi-steady-state plasma discharges reach conditions to study long-pulse plasma material interactions (PMIs). The newly commissioned HIDRA-Material Analysis Test-stand PMI diagnostic is used to place a drop of lithium onto a heated tungsten surface, transfer the sample in-vacuo and expose it in a helium plasma. Helium is of interest as there is an open question to whether lithium will be able to remove helium ash in real fusion devices. The introduction of the W-Li sample in HIDRA resulted in evaporation of lithium into the helium plasma during a 600 s pulse and caused a reduction of over 90% in neutral pressure during the discharge. It was also observed that the plasma density and temperature increased by over 2.5 times. Using spectroscopy and a helium collisional radiative model, the peak temperature and density of the helium plasma can be monitored during the discharge. During lithium evaporation, as significant lithium ionization occurs, there is a 85% drop in the HIDRA vessel neutral pressure, despite a constant flow rate of He gas. This reduction in neutral pressure is supported by spectroscopy data with corresponding reductions in He I line intensities (587 nm, 667 nm, 706 nm, and 728 nm), as well as those of other impurities. At one point in the discharge a lithium plasma is created, as indicated by an increase in Li + emission and a complete reduction in He + emission, but the electron density jumps from ne = 3 × 10 18 m –3 to over ne = 8 × 10 18 m –3 while the core temperature stays relatively constant between Te = 16 eV and 20 eV. Once lithium has completely evaporated from the sample and the majority of the ionized lithium has diffused from the plasma to the vessel walls, pressure and spectroscopy data paired with He collisional radiative model calculations shows a re-establishment of a helium plasma in a low recycling regime. In this regime, the density drops down to ne = 2 × 10 18 m –3 and the electron temperature increases from T e = 20 eV to over T e = 50 eV indicating an increase in helium heating efficiency. This is also indicated by the He+ emission re-establishing and having a higher intensity. Here, we show the results from the first lithium campaign in HIDRA. In the presence of lithium, and in particular when lithium ions are present, the helium disappears from the plasma via an as of yet unknown complex relationship that needs to be further studied. The most likely explanation is that the lithium ions are distributed around the vessel and able to trap helium to the surface turning HIDRA into a large gettering surface. These results have potential implications on future plasma facing component design using liquid lithium for impurity and recycling control using limiters and divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aerothermodynamic assessment of corrugated panel thermal protection systems

The results of an analysis of wind-tunnel data for the effect of wave shape on local boundary layer flows and heat transfer are used to assess the feasibility of using corrugated panels on a Single-Stage-To-Orbit (SSTO) vehicle. In the flight assessment, available engineering design procedures are used to compute the inviscid and boundary layer flow for a given flight trajectory. The flight vehicle flow field calculations are then used in conjunction with wind-tunnel corrugated panel heating data correlations to examine the heating penalties associated with using corrugated panels at a number of locations on the vehicle. Severe heating penalties, as a result of the corrugations, were found over large areas of the windward surface. Maximum surface temperature penalties ranged from 100K to 500K and very large structural advantages for the corrugations would be required to override the heating penalties. Computed peak corrugation temperatures were within allowable limits of the equivalent smooth wall material for most leeward and side fuselage locations and corrugated panels will be most likely limited to these locations.

Brandon, H. J.↗

Technology Development and Infusion by NASA's Entry Systems Modeling Project

This paper describes recent development of modeling and simulation technologies for entry systems and their infusion into NASA's exploration missions. Technology development is organized and prioritized using a system-level perspective, resulting in four broad technical areas of investment: (1) Thermal protection material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. The paper will highlight key contributions from each of these areas, their impacts from a spacecraft and mission design perspective, and discuss planned future investment. Aspects of each technical area are only briefly summarized here. Thermal protection material modeling is geared toward high-fidelity, predictive models capable of optimizing design performance, post-flight reconstruction, and quantifying thermal protection system reliability. New computational tools and experimental techniques have been applied to Orion, MSL/Mars 2020, Mars InSight, and Mars Sample Return missions. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. In both cases, substantial improvements in model uncertainty have directly impacted the development of mission margin policies, flight instrumentation design and analysis (Orion and Mars 2020), and have even revealed the importance of neglected phenomena like mid-wave infrared radiation of CO2. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics affecting Orion, Commercial Crew, and Mars programs; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate augmentation of aerothermal environments at Mars during dust storms; and studies of roughness-induced heating augmentation relevant to tiled (Orion, Mars 2020) and woven (Mars Sample Return) thermal protection systems. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent, with precision landing requirements driven by Mars human exploration goals.

Barnhardt, Michael D.↗

Time-Resolved X-ray Operando Observations of Lithiation Gradients across the Cathode Matrix and Individual Oxide Particles during Fast Cycling of a Li-Ion Cell

Lithiated transition metal oxides serve as active materials in the positive electrode (cathode) of lithium-ion cells. During electrochemical cycling, lithium ions intercalate and deintercalate into these oxide particles. This behavior causes two types of lithiation gradients to emerge: (i) a bulk gradient across the depth of the cathode matrix (averaged over individual oxide particles) and (ii) a microscopic gradient across the particles themselves, which also depends on their location in the electrode. Here we show how both gradients can be studied using operando X-ray diffraction during 4C charge and 4C discharge. The oxide (de)lithiation is estimated from the unit cell parameters by indexing the X-ray diffraction spectra. By fitting the lithiation profiles with orthogonal polynomials, the bulk gradients across the electrode thickness are quantified. These gradients develop as the current flows through the cell and dissipate during open-circuit and potentiostatic-hold periods. Further details of lithiation dynamics can be obtained through shape analysis of the Bragg peaks. In particular, from electrochemical model simulations, we show that the width and skewness of the (003) peak track (de)lithiation fronts moving across the individual oxide particles.

(003) peak↗

Meeting the Grand Challenge of Protecting Astronauts Health: Electrostatic Active Space Radiation Shielding for Deep Space Missions

This report describes the research completed during 2011 for the NASA Innovative Advanced Concepts (NIAC) project. The research is motivated by the desire to safely send humans in deep space missions and to keep radiation exposures within permitted limits. To this end current material shielding, developed for low earth orbit missions, is not a viable option due to payload and cost penalties. The active radiation shielding is the path forward for such missions. To achieve active space radiation shielding innovative large lightweight gossamer space structures are used. The goal is to deflect enough positive ions without attracting negatively charged plasma and to investigate if a charged Gossamer structure can perform charge deflections without significant structural instabilities occurring. In this study different innovative configurations are explored to design an optimum active shielding. In addition, to establish technological feasibility experiments are performed with up to 10kV of membrane charging, and an electron flux source with up to 5keV of energy and 5mA of current. While these charge flux energy levels are much less than those encountered in space, the fundamental coupled interaction of charged Gossamer structures with the ambient charge flux can be experimentally investigated. Of interest are, will the EIMS remain inflated during the charge deflections, and are there visible charge flux interactions. Aluminum coated Mylar membrane prototype structures are created to test their inflation capability using electrostatic charging. To simulate the charge flux, a 5keV electron emitter is utilized. The remaining charge flux at the end of the test chamber is measured with a Faraday cup mounted on a movable boom. A range of experiments with this electron emitter and detector were performed within a 30x60cm vacuum chamber with vacuum environment capability of 10-7 Torr. Experiments are performed with the charge flux aimed at the electrostatically inflated membrane structure (EIMS) in both charged and uncharged configurations. The amount of charge shielding behind and around the EIMS was studied for different combinations of membrane structure voltages and electron energies. Both passive and active shielding were observed, with active shielding capable of deflecting nearly all incoming electrons. The pattern of charge distribution around the structure was studied as well as the stability of the structures in the charge flow. The charge deflection experiments illustrate that the EIMS remain inflated during charge deflection, but will experience small amplitude oscillations. Investigations were performed to determine a potential cause of the vibrations. It is postulated these vibrations are due to the charge flux causing local membrane charge distribution changes. As the membrane structure inflation pressure is changed, the shape responds, and causes the observed sustained vibration. Having identified this phenomenon is important when considering electrostatically inflated membrane structures (EIMS) in a space environment. Additionally, this project included a study of membrane material impacts, specifically the impact of membrane thickness. Extremely thin materials presented new challenges with vacuum preparation techniques and rapid charging. The thinner and lighter membrane materials were successfully inflated using electrostatic forces in a vacuum chamber. However, care must be taken when varying the potentials of such lighter structures as the currents can cause local heating and melting of the very thin membranes. Lastly, a preliminary analysis is performed to study rough order of magnitude power requirements for using EIMS for radiation shielding. The EIMS power requirement becomes increasingly more challenging as the spacecraft voltage is increased. As a result, the emphasis is on the deflection of charges away from the spacecraft rather than totally stopping them. This significantly alleviates the initial power requirements. With modest technological development(s) active shielding is emerging to be a viable option.

Radiation↗

NRC Multiphysics Analysis Capability Deployment (FY2021 - Part 2)

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” Task 4b was completed for this report. INL developed a sample problem showing how to use Serpent 2 to calculate macroscopic cross sections for use in Griffin for a typical sodium fast reactor Unprotected Loss of Flow (ULOF) transient. The complicating factor is the relatively large axial gradient of the coolant (and hence fuel) temperature and its evolution during the transient. A 3D Griffin model of the Advanced Burner Test Reactor (ABTR) is coupled to the System Analysis Module (SAM) to perform the ULOF transient. The model includes various explicit feedback mechanisms for fast reactors, including Doppler, radial expansion from the displacement of the support plate, and axial expansion from the displacement of the fuel pins. The results for the ULOF transient are consistent with published values. We observed small effects in the transient results that arise from the superhomogenization equivalence correction of uniform and nonuniform temperature data sets, but they diminish as more dominant reactivity mechanisms are added to the model. Potential improvements to the neutronics model include: adding the differential control rod and sodium feedback mechanisms, implementing cross section re-homogenization based on the partial volume of the various materials, and enhancing the kinetics parameters. Furthermore, the fluids model only includes four subchannels to represent the entire active core due to performance limitations in SAM when attempting to model one subchannel per assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗

Comparing energy and water use of aqueous and gas-based metalworking fluids

Gas-based metalworking fluids (MWFs) have been proposed as alternative coolants and lubricants in machining operations to mitigate concerns surrounding water use and pollution, industrial hygiene, occupational health, and performance limitations associated with water-based (aqueous) MWFs that are ubiquitously used in the metals manufacturing industry. This study compares the primary energy and water use associated with the consumptive use, delivery, and disposal of aqueous MWFs with three gas-based MWFs in the literature-minimum quantity lubricant-in-compressed air (MQL), liquid/gaseous N-2, and liquid/supercritical CO2. The comparison accounts for reported differences in machining performance in peer-reviewed experimental studies across several machining processes and materials. The analysis shows that despite the reported improvement in tool life with N-2 and CO2-based MWFs, the electricity- and water-intensive separation and purification processes for N-2 and CO2 lead to their higher primary energy and water use per volume of material machined relative to water-based MWFs. Although MQL is found to have lower primary energy use, significant consumptive water use associated with the vegetable oil commonly used with this MWF leads to higher overall water use than aqueous MWF, which is operated in a recirculative system. Gas-based MWFs thus shift the water use upstream of the manufacturing plant. Primary energy and water use of gas-based MWFs could be reduced by focusing on achieving higher material removal rates and throughput compared to aqueous MWF instead of solely targeting improvements in tool life. Additionally, the consumptive use of CO2 and N-2 MWFs could be minimized by optimizing their flow rates and delivery to precisely meet the cooling and lubrication needs of specific machining processes instead of flooding the tool and workpiece with these gases. This article met the requirements for a gold-gold JIE data openness badge described at http://jie.click/badges.

cutting fluids↗

A knowledge gap analysis for transient CHF prediction within RELAP5-3D

Here we compare several experimental pool and flow boiling critical heat flux experiments with corresponding computational models developed with RELAP5-3D. Encompassed within this paper are several sensitivity analyses for the boiling heat transfer correlations, critical heat flux (CHF) prediction within RELAP5-3D and the thermal properties of the cladding for experiments in the Transient Reactor Test Loop (TRTL) and Transient Reactor Test (TREAT) facilities. The sensitivity analyses were consistent in determining that there is a large contribution from the critical heat flux multiplier which in this application, accounts for the effects of power transient critical heat flux. Comparison to experimental results of both the TREAT CHF-SERTTA experimental campaign and a University of New Mexico (UNM) pool boiling experimental campaign demonstrate both the under-prediction of transient CHF and nucleate boiling heat transfer. To address the gaps in knowledge regarding the transient heating effect for both boiling heat transfer and CHF prediction, we present an experimental test matrix for the TRTL facility that will isolate the power transient effect under PWR relevant operational conditions. The test matrix will incorporate 6 power pulses with a constant energy deposition and a timescale that spans over several orders of magnitude, including a steady state test with an exponential period of 10 s.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗