Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normalization”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Detailed mixing measurements from single-hole converging ECN injectors using Rayleigh scattering

Time resolved liquid and vapor fields of dodecane and oxymethylene ethers are measured from Spray A-3 and Spray D using high speed Rayleigh scattering and diffuse back illumination at the Engine Combustion Network (ECN) Spray A condition of 900 K and 22.8 kg/m 3 . Global quantities including mixture fraction, vapor and liquid penetration, as well as spreading angle are measured. The mixture fraction fields and vapor penetration profiles are well predicted by the 1-D Musculus-Kattke model. The mixture fraction field and vapor penetration from Spray A-3 are similar to those measured from Spray A in previous works. Spray D exhibits higher mixture fraction fields and vapor penetration due to its larger nozzle diameter. The quasi-steady mixture fraction fields from these injectors scale well when distance from the injector is normalized by the nozzle diameter. The turbulent dissipation structures were also analyzed based on the orientation, thickness, and magnitude of the mixing layers. Furthermore, the orientation and thickness are similar to other measurements in atmospheric gas jets, while the magnitude is lower. The thickness and magnitude are subject to uncertainties due to limitations in the imaging resolution of the system but still provide an order of magnitude as a useful reference for comparison against computational fluid dynamic simulations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

A machine-learning-aided data recovery approach for predicting multi-material thermal behaviors in advanced test reactor capsules

Instrumented experiments conducted at test reactors are essential to the deployment of new advanced reactor systems. Designing new experiments and generating data on specific reactor conditions require significant investments in terms of both time and cost. Finite element analysis software can be used to create high-fidelity models of experiment environments in order to support the actual experiments, but computation time remains a concern in terms of applying outcomes to real-time usage of data (e.g., a digital twin [DT]). Here, the present research proposes a machine-learning (ML) aided approach to making temperature and displacement predictions based on the thickness of the outer gas gap on the experimental capsule used for in-pile demonstration of a novel new thermal conductivity probe in the Advanced Test Reactor (ATR). This capsule consisted of U10Zr fuel, a rodlet, sodium, and inner and outer capsules. Gas gaps existed between the fuel and the rodlet, and between the inner and the outer capsule. The learning data pertained to an experimental capsule's radial distributions of temperature and displacement, as obtained based on Abaqus and the physical features. For the first step of ML sequence, the temperature was predicted using three positional parameters. Next, the displacement was predicted using seven additional parameters. Each physical feature was normalized in order to be both nondimensional and standardized. The temperature and displacement predictions showed good agreement with the simulation results in all cases involving interpolation and extrapolation. Furthermore, data similarity enhancement increased the similarity between the training and the target data, thereby increasing the predictive accuracy of the ML models. In certain extrapolation cases involving limited original ML model accuracy, data similarity enhancement and data recovery was able to somewhat improve this accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Thermal conductivity of 3D-printed block-copolymer-inspired structures

This study primarily focuses on examining the impact that geometric structure has on thermal conductivity of multi-phase constructs in different 3D-printed poly(lactic acid), PLA, samples. The investigated structures are inspired by morphologies formed by diblock copolymers: lamellae, hexagonally packed cylinders, and gyroid. This research also investigates how volume percentage and material combination influence the thermal conductivity of these structures. Further, the samples can be tailored to simulate various thermal management structures observed in practical applications, such as thermal interface materials in electronic devices. Thermal conductivity ratio is controlled using air, the least conductive material at 0.026 W/(m K), PLA at 0.136 W/(m K), and thermal paste at 5.11 W/(m K). Different models were tested against thermal conductivity measurements in order to capture the effect of material type (PLA-Air versus PLA-Thermal Paste), volume percentage, structure, and orientation. Simple, effective medium models were good predictions of thermal conductivity in lamellar structures, but it was necessary to develop models for conduction through cylindrical and gyroid structures. Finally, all results were normalized to find a universal model that is independent of structure and material. This approach provides a simple method to predict how to reduce or enhance transport properties and heat management capabilities of 3D printed objects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Alternating conduction and convection drying of paper – an experimental analysis with a continuous data acquisition approach

In conventional multi-cylinder drying of paper and board, both conductive drying from steam-heated dryer cylinders and convective drying by flowing air over the paper surface in the pockets are used. Conductive drying from steam-heated drying cylinders is a critical component in providing the necessary thermal energy to paper and board as they dry. Steam temperature and internal and external resistances at the contacting surface are critical process parameters influencing the conductive drying process. An experimental setup was developed to study the alternating conductive and convective drying of paper and board. Paper sheet moisture, temperature, and temperature distribution within the heated platen and the instantaneous heat flux as the sheet was being dried were measured. The instantaneous heat flux, contact heat transfer coefficient, and drying rates were determined as drying proceeds. Experimental results, as well as comparisons to literature and commercial data, are presented. The conductive heat transfer coefficients determined were compared to traditional correlations normally used in the modeling of paper drying. Similarly, the convective heat and mass transfer coefficients are also determined and compared to literature data. In addition to the evaluation of alternating conductive and convective drying characteristics of paper and board, the potential inclusion of auxiliary energy components will also be included. Experimental results from the conduction and convection drying system are presented. Furthermore, this data will be useful in process development, intensification of manufacturing processes, and modeling and simulation of paper drying processes.

42 ENGINEERING↗
DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗
DOE OSTI2020

Determination of the enthalpy of adsorption of hydrogen in activated carbon at room temperature

The development of high-performance materials for hydrogen storage by adsorption requires detailed understanding of the adsorbate-adsorbent interactions, e.g., the enthalpy of adsorption ΔH, which measures the interaction strength. The determination of ΔH for a weakly adsorbing gas such as hydrogen in a carbonaceous porous material is difficult experimentally, normally requiring measuring two cryogenic adsorption isotherms. Here we demonstrate a calculation of ΔH based on ca. room temperature adsorption isotherms at 273 K and 296 K using the Clausius-Clapeyron equation. This requires an estimation of the volume of the adsorbed film (~40%, ~12% of the total pore volume at 77 K, 296 K, respectively) obtained from fits of the excess adsorption isotherms to an Ono-Kondo model with the auxiliary use of a fixed point corresponding to the saturation film density (estimated as 100 ± 20 g/L) which appears to be remarkably sample and temperature independent, i.e., a property of the adsorbate. The calculated room temperature enthalpy of adsorption ΔH = 8.3 ± 0.4 kJ/mol is in excellent agreement with the low-coverage cryogenic determination of ΔH. The methodology hereby proposed facilitates reliable calculations of the enthalpy of adsorption at room temperatures for weakly-adsorbing gases.

08 HYDROGEN↗
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets—Part I: Mesoscale simulations

Recently a mesoscale model was developed that includes thermomechanically consistent constitutive equations with dissipation due to porous compaction, inelastic distortional deformation rate, and slipping on persistent dry joint surfaces. In this Part I of a two-part paper, we describe enhancements of the mesoscale model which include: normal compliance with hysteresis, non-persistent joints with finite area, joint sets with angular deviation, and the weakening effect of fluid pressure in the joints. Examples study the influences of these effects on average strength. These computationally expensive simulations are used to inspire functional forms in a continuum model, described in Part II, which is needed for large scale simulations of real applications.

42 ENGINEERING↗
DOE OSTI2020

Scaling of two-phase water-steam relative permeability and thermal fluxes in porous media

Two-phase water-steam flow conditions are frequently encountered in many engineering applications, including geothermal reservoirs. Although routine calculations are based on the multiphase Darcy’s law, the role of the topology of the flowing phases at the pore-scale is usually neglected in the estimation of relative permeabilities. Instead, the latter are frequently computed using empirical models like the Corey correlation. In this work, we first apply the model for relative permeabilities based on pore-scale flow regimes developed by Picchi and Battiato (2019), Relative permeability scaling from pore-scale flow regimes, Water Resour. Res. 55, 3215–3233, to scenarios typical of geothermal reservoirs and then extend it by deriving the scaling laws for the transmissibilities and the thermal properties as a function of temperature. First, we discuss the scaling behavior of normalized relative permeabilities in terms of viscosity ratio and capillary number of water-steam systems and, then, we provide a validation of the model against experimental data available in the literature. The model captures well the data trends collected in real 3D porous media. These results suggest that water-steam relative permeabilities follow the same scaling behavior of gas-liquid systems where the non-wetting phase is much less viscous than the wetting phase. Finally, we investigate the impact that relative permeabilities have on heat transfer rates at two-phase flow conditions and the scaling of mass and energy transmissibility and thermal properties of the mixture. An estimation of the exergy carried by the two-phase water-steam mixture is also included.

42 ENGINEERING↗
DOE OSTI2021

A critical dislocation velocity for serration mechanism transition in a nickel-chromium solid solution alloy

Here, the influence of strain rate across three orders of magnitude (1.70 × 10 –5 /s to 1.43 × 10 –2 /s) along with the effect of the plastic strain accumulation (up to 10%) on the serrated plastic flow were investigated in the nickel-chromium (Ni-Cr) solid solution alloy Nimonic 75 by performing constant-strain-rate tension testing at 600 °C. As the strain rate decreased, the critical strain for the onset of serrations transitioned from normal behavior to inverse behavior. The serrated flow was characterized as Type A+B serration at high strain rate (1.43 × 10 –2 /s). In the intermediate strain-rate regime (1.43 × 10 –3 /s and 1.45 × 10 –4 /s), Type B serrations were observed and followed by a transformation to Type C+B serrations. At the low strain rate (1.70 × 10 –5 /s), the plastic flow immediately displayed Type C serrations, which later evolved into Type C+B serrations. Regardless of the strain rate, plastic strain, or dislocation density, a critical dislocation velocity falling in the range of 1.2 × 10 –6 – 2.2 × 10 –6 m/s was identified to signify the onset of Type C serration, whereby the mobile dislocations break free from the solute cloud for short bursts of deformation. Finally, a novel model by solute rearrangement across dislocation cores was used to understand how the critical dislocation velocity is quantitatively determined by the rate at which solute atoms are able to hop across the glide plane as a partial dislocation core moves through the lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deformation behavior of CP-titanium under strain path changes: Experiment and crystal plasticity modeling

The deformation behavior of commercially pure rolled titanium subjected to strain path changes is studied using experiments and a crystal plasticity model. Four different loading combinations are performed at room temperature to study the activation of slip, twinning, de-twinning and double-twinning in hexagonal closed packed titanium. The strain paths considered are: rolling direction compression (RDC) followed by normal direction compression (NDC), RDC followed by transverse direction compression (TDC), NDC followed by RDC, and NDC followed by TDC. An EBSD-based analysis of the distribution of boundary misorientation angles before and after reload was developed to analyze the evolution of {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twins. This analysis supports the model results concerning the treatment of twin reorientation. A de-twinning and double-twinning model accounting for back stress effects, an important feature of strain path changes, is implemented within the framework of the visco-plastic self-consistent (VPSC) model along with a dislocation density (DD) based hardening scheme. In the model, plasticity is accommodated by prismatic $\langle$a$\rangle$, basal $\langle$a$\rangle$ and pyramidal $\langle$c+a$\rangle$ slip modes, and {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twinning modes. The VPSC model predicts the evolution of twinning, de-twinning and double-twinning processes for both tensile and compression twinning modes under strain path change. The model predicts macroscopic stress-strain response, texture evolution, and twin volume fraction that are in agreement with experimental observations. The evolution of texture is investigated in detail by separately analyzing the twinned domains, rather than the evolution of the global texture.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling the non-Schmid crystallographic slip in MAX phases

We present a crystal plasticity constitutive relation for the description of experimentally observed non-Schmid crystallographic slip in a class of ternary carbides and nitrides commonly referred to as MAX phases. In the constitutive relation, we assume that the evolution of the slip system strength in MAX phases has two components – a classical component that depends on the Taylor cumulative shear strain and a non-Schmid component that depends on the stress normal to the slip plane. The non-Schmid crystal plasticity constitutive relation is then used to carry out finite element simulations of micropillar compression of single crystals of two MAX phases, Ti 2 AlC and Ti 3 AlC 2 . The finite element simulations not only quantitatively predict the stress – strain response of a wide range of crystallographic orientations of the micropillars but also rationalize the non-uniform deformation and the deformed shape of the micropillars observed in the experiments for the two materials. As a result, parametric studies are also carried out to quantify the role of the non-Schmid effect and understand the effects of key experimental parameters on the stress – strain response of the micropillars of the two MAX phases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling inter- and intra-granular dislocation transport using crystal plasticity

Here, this work presents the development of a crystal plasticity material model that incorporates both dislocation transport within grains and dislocation transfer across grain boundaries. This model has been implemented in the open-source finite element code MOOSE. In addition, a novel geometry-based criterion is developed to determine the direction of dislocation transfer across grain boundaries. The transfer criterion incorporates the geometric features of the grain boundary, such as the grain boundary plane normal, and its misorientation, which is accounted for through the orientation of the incoming and outgoing slip systems. The model is tested with several cases, including a copper single crystal, bi-crystal, and polycrystal. The development of the transfer criterion, implementation of the model, and its application to these test cases are discussed in detail.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Active magnetic regenerative cooling with smaller magnets

Magnetocaloric heat pumping near room temperature, a.k.a. magnetic cooling, relies on the active regenerator cycle to achieve functional temperature spans and realize economic and societal benefits promised by this emerging solid-state cooling technology. The cycle itself depends upon synchronizing oscillating flow of a heat transfer fluid through a solid porous active material matrix, or the refrigerant, with periods when the refrigerant is in the highest available (field-on) and in nearly zero (field-off) magnetic fields to accomplish heat transfer. With this in mind, in this work we analyze varying flow and magnetic field wave forms and the timing between when the fluid is pumped and when the magnetic field is turned on and off. We demonstrate that the volume and the cost of permanent magnet generating the field changes can be cut nearly in half with little to no effect on the device temperature span and cooling power normalized by the refrigerant mass.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗
DOE OSTI2021

Experimental workflow to estimate model parameters for evaluating long term viscoelastic response of CO2 storage caprocks

Understanding the time-dependent behavior of reservoir and sealing formations is critical to assessing risks associated with geological carbon storage since time-dependent deformation strongly influences mechanical responses of some rock types. Many studies have evaluated the risk of CO2 leakage and induced seismicity by assuming poroelastic rheology in sealing formations. Few have considered viscoelastic or other time-dependent responses, where the existing literature adopts 1D models to represent long-term time-dependent responses. This is primarily because to date, the general form of a reasonable 3D time-dependent model for rocks remains unclear. In this paper, we address this unclear issue by proposing a new workflow to select constitutive modeling parameters to evaluate if a 3D viscoelastic model is reasonable using several-hour-long experimental data and a power-law response to extrapolate to the decades-long time frames of interest in geologic carbon storage. To provide experimental data, we conducted multi-level loading/unloading triaxial relaxation tests with four rock types. The experimental results showed that the maximum load relaxation observed is approximately 49%, with some rock types showing as little as 1.4%. Using a simple linear viscoelastic model, parameters were chosen such that a maximum deviation of 1.5 MPa in axial stress and 7 MPa in radial stress was attained with the extrapolated 30-year data. We found that a reasonable parameter range for the normalized elastic modulus is 0.1~2 for rocks with significant time-dependent responses and 0.01~0.06 for those with small time-dependent responses. No matter how significant time-dependent responses are for rocks considered, our results showed that the relaxation time has a general range of 1~10^10 s, whose time scale can be one or two orders higher than a time frame typically envisioned for CO2 injection projects.

Stress relaxation, 3D time-dependent model, viscoe↗
DOE OSTI2022

Critical review of the Mogi failure criterion based on true-triaxial laboratory data analysis and theoretical considerations

The influence of the intermediate principal stress, σ 2 , is disregarded when rock strength is analyzed using the Mohr-Coulomb or Hoek-Brown criteria. Several polyaxial strength criteria have been proposed to account for the influence of σ 2 . One of the methods that received considerable attention is the Mogi criterion. The criterion modifies the Nadai or Drucker-Prager approaches by exchanging the normal octahedral stress with what is introduced as the mean effective stress, σ m,2 . This change is based on experimental evidence that the influence of the intermediate principal stress on strength is less pronounced than the influence of the minimum principal stress. However, the Mogi criterion was derived somewhat empirically and consequently it is essential to analyze its potential theoretical deficiencies. In the current note, we analyzed the Mogi failure criterion, describing its characteristics in both τ oct –σ m,2 and σ 1 -σ 2 spaces. We investigated the relation between Mohr-Coulomb and Mogi criteria, which marks the basis for the commonly used Mogi-Coulomb approach. Here, we compared the performance of the Mogi criterion against available polyaxial strength data using three fitting approaches: Mogi-Coulomb based on triaxial compression results only, linear Mogi based on all polyaxial strength results, and power law Mogi also based on all data. The analysis reveals several serious deficiencies of the Mogi approach which all suggest that the Mogi criterion cannot properly represent the σ 2 dependence of rock strength in many cases.

58 GEOSCIENCES↗