Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D temperature modeling”

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 415 records · Page 23

ElasTool v3.0: Efficient computational and visualization toolkit for elastic and mechanical properties of materials

Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. Here, the ElasTool v3.0 toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like elastic tensor, Young's modulus, bulk modulus, and Poisson's ratio. It introduces capabilities for computing minimum thermal conductivity, linear compressibility, rendering the Christoffel equation, and elastic energy density. Notably, it integrates advanced visualization tools, including compatibility with Plotly and Elate web platforms for interactive web-based property exploration. A key feature of ElasTool v3.0 is the implementation of second-order elastic constants (SOECs) for tubular 2D-based nanostructures and nanotubes. Leveraging high-efficiency strain-matrix sets (OHESS), the toolkit now facilitates efficient computation of elastic constants and mechanical properties at both zero and finite temperatures for 1D, 2D, and 3D dimensions. ElasTool is openly accessible on GitHub: https://github.com/gmp007/elastool.

1D, 2D, 3D, and tubular 2D nanostructure and nanot↗

Acceleration of the Fast Solar Wind by Solitary Waves in Coronal Holes

The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear, time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that solitary-like waves are generated in coronal holes nonlinearly by torsional Alfven waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature, and magnetic field geometry) that will become available from the recently launched SOHO spacecraft.

Ofman, Leon↗

Acceleration of the Fast Solar Wind by Solitary Waves in Coronal Holes

The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear. time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that nonlinear waves are generated in coronal holes by torsional Alfv\'{e}n waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the, fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature. and magnetic field geometry,) that will become available from the recently launched SOHO spacecraft.

Ofman, Leon↗

Experimentally Validated Computational Modeling of Creep and Creep-Cracking for Nuclear Concrete Structures

In a Nuclear Power Plant, one of the most important components is the concrete nuclear reactor cavity, which serves both a structural and protective function as the biological radiation shield. Given that creep has been identified as a major knowledge gap in the assessment of nuclear structures (NUREG/CR-7153), this work helps to further the understanding of creep behavior of massive concrete containment structures for decades to enable safe and long-term operation of these facilities. This project has developed a robust, experimentally validated model to predict creep in nuclear concrete structures for up 60 years using short-term creep data thereby enabling a longer service life of critical facilities and early detection of structural failure. The work presented in this report is a pairing of computational and experimental methods. For the first time, the time temperature superposition (TTS) principle was successfully used to generate a uniaxial creep compliance master curve to predict mortar creep response for up to 22,500 days (nearly 60 years) at a reference temperature of 20°C. These data were used as input into finite element analysis (FEA) codes that use highly realistic random, 3D concrete microstructures from reconstructed coarse limestone aggregates. Finite element analysis performed provides the ability to quickly upscale mortar viscoelastic behavior to long-term concrete creep/relaxation data. A master creep compliance curve, constructed from the TTS principle, spanning 27 years, was used to validate two and a half decades of simulated concrete creep. Concurrently, three different simulated wall specimens were designed to mimic the behavior of post-tensioned concrete nuclear containment facility vessel walls over time as a result of concrete creep. The specimens were designed with different thicknesses, transverse and longitudinal reinforcement ratios, and level of post-tensioning stress. Each specimen contained various instrumentation to measure internal concrete temperature, concrete strain, and post-tensioning strain hourly for over 3 years. The concrete creep model developed in this project, based on the FEA concrete simulations, was applied to simulate the structural-scale experiments of prestressed concrete walls conducted in this project using the Grizzly code. These models can represent the effects of reinforcing and prestressing. Although there are some discrepancies with the experimental data, the model can predict the overall trends of the creep response in these experiments. One of these experimental models was also applied to an extended time to demonstrate how the findings from this study can be used to predict the behavior of actual structures of interest that have been in service for extended periods of time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Wrought FeCrAl alloy (C26M) cladding behavior and burst under simulated loss-of-coolant accident conditions

Cladding burst experiments for FeCrAl cladding were performed in the Severe Accident Test Station facility at Oak Ridge National Laboratory. These experiments were simulated using the BISON fuel performance code to better understand the cladding plastic behavior and failure under simulated loss-of-coolant accident conditions. 3D cladding surface boundary conditions were generated using composite axial and azimuthal profiles from experiment thermocouple data. To improve the simulation analysis capabilities in BISON for cladding burst behavior, new thermal creep, plasticity, and failure stress models specific to C26M, a wrought FeCrAl alloy, were developed and implemented. Initial cladding burst results indicated a general underprediction in the failure temperature of the six cladding burst simulations versus the observed failure temperatures. Close investigation of the experiment timing versus the underlying tensile test data revealed that, compared with the tensile specimens, the cladding tubes did not experience the same long holding time at high temperatures. New tensile tests were performed at high temperatures using a temperature ramp similar to the simulated loss-of-coolant accident experiments. These new tensile curves showed an approximately 80% increase in the ultimate tensile strength of the C26M alloy, indicating that a holding time of 10 min at 700 °C and 800 °C allows annealing to change the material microstructure. Using the updated tensile properties, the burst temperatures and stresses from the simulations showed remarkable agreement with the experimental results. This study was then extended by varying the initial pressure to highlight the burst temperature difference between standard Zircaloy-4 and C26M cladding under equivalent conditions. The results show that C26M has a burst temperature that is approximately 70–130 K greater than that of Zircaloy-4. In conclusion, these modeling predictions can be further improved by collecting high-temperature tensile data for C26M beyond the temperature ranges used in this work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Control of Surface Morphology during the Growth of (110)-Oriented GaAs by Hydride Vapor Phase Epitaxy

We present a study of the effect of hydride vapor phase epitaxy (HVPE) growth conditions on the morphology of GaAs grown on vicinal and nominally exact (110) GaAs substrates. We evaluate epilayer surfaces using a combination of Nomarski microscopy and atomic force microscopy. The surface morphology strongly depends on the growth conditions employed, as well as the substrate orientation. On substrates offcut 3° toward (111)A, faceting generally develops under conditions in which the growth is Ga-limited, although Ga-limited growth is not the only requirement for faceting. Low growth temperatures lead to a faceted morphology, while higher growth temperatures favor smooth growth and facet-free morphology. We show that the surface morphology is kinetically controlled and that the tendency to facet correlates with the growth rate. 3D growth occurs on nominally exact substrates, resulting in hillock formation, but reduction of the group V precursor partial pressure results in smooth morphology. Overall, our results are consistent with models for step-bunching-induced surface faceting that invoke a negative Ehrlich–Schwoebel step-edge barrier.

(110) oriented GaAs↗

A Technique to Measure Coronal Electron Density, Temperature, and Velocity Above 2.5 R from Sun Center Using Polarized Brightness Spectrum

The current model for the polarized brightness (pB) spectrum has a decades-long history of progressively incorporating its dependence on electron density N(e), temperature T(e), and flow velocity in the radial direction V(e). The pBNe spectrum follows the exact shape of the photosphere spectrum, which is not smooth, which is expected from the thermal Doppler broadening of the photosphere spectrum due to the high coronal T(e). The pBN(e) spectrum is smooth, but the free coronal electrons remain static and unaffected by solar wind, and the pBN(e)T(e)V(e) spectrum is red-shifted by electrons seeing a red-shifted photosphere spectrum as they flow away from the Sun as solar wind, which takes a radial direction above 2.5 R from Sun center. In this article, we review the progress of the above three model pB spectra in describing the observations and highlight the differences, first by comparing the three model pB spectra against wavelength using a model for Ne and constant values for T(e) and V (e), and second by generating three model 2D pB maps by integrating over a selected wavelength region in the three model pB spectra along lines of sight passing through the 14 July 2000 (“Bastille Day”) coronal mass ejection (CME) model, which contains 3D information on N(e), T(e), and V e. In this regard, the COronal Diagnostic EXperiment (CODEX) on the International Space Station (ISS) in 2024 will measure N(e), T(e), and V(e) by matching the measured pB with modeled pBN(e)T(e)V(e) in selected wavelength regions using multiple filters.

Electron density↗

Computer-controlled multi-parameter mapping of 3D compressible flowfields using planar laser-induced iodine fluorescence

A computer-controlled technique, using planar laser-induced iodine fluorescence, for measuring complex compressible flowfields is presented. A new laser permits the use of a planar two-line temperature technique so that all parameters can be measured with the laser operated narrowband. Pressure and temperature measurements in a step flowfield show agreement within 10 percent of a CFD model except in regions close to walls. Deviation of near wall temperature measurements from the model was decreased from 21 percent to 12 percent compared to broadband planar temperature measurements. Computer-control of the experiment has been implemented, except for the frequency tuning of the laser. Image data storage and processing has been improved by integrating a workstation into the experimental setup reducing the data reduction time by a factor of 50.

Donohue, James M.↗

Engineering Graphene-Ceramic 3D Composite Foams by Freeze Drying

A 3D graphene foam produced by the chemical vapor deposition (CVD) technique is recognized as an effective nanofiller material, as it does not agglomerate in the matrix. Though CVD facilitates pristine graphene foam production, the method poses limitations in developing large-scale 3D graphene composite foams as advanced engineering materials. In this paper, a freeze-drying (FD) process is used to produce a composite 3D foam, which is a mixture of graphene nanoplatelet (GNP) and a low-temperature co-fired ceramic (LTCC). The freeze-dried GNP-LTCC reticulated 3D composite foam has an average pore size ranging from 70 to 100 μm. Pore size is varied by regulating the heat transfer rate during the freezing process using a thermally conductive aluminum (Al) mold and a thermally insulating acrylonitrile butadiene styrene (ABS) mold. Computational thermal modeling is used to visualize the heat transfer and its effect on foam pore size. Subsequently, a freeze-dried GNP-LTCC composite foam is embedded into the LTCC matrix to form a hierarchical assembly by a spark plasma technique without compromising the 3D structure of the FD foam. This study established that a simple, eco-friendly, and scalable processing methodology can produce advanced graphene-based 3D composite foams as the future tailorable nanofillers for designing multimatrix materials.

36 MATERIALS SCIENCE↗

Strong and fragile topological Dirac semimetals with higher-order Fermi arcs

Abstract Dirac and Weyl semimetals both exhibit arc-like surface states. However, whereas the surface Fermi arcs in Weyl semimetals are topological consequences of the Weyl points themselves, the surface Fermi arcs in Dirac semimetals are not directly related to the bulk Dirac points, raising the question of whether there exists a topological bulk-boundary correspondence for Dirac semimetals. In this work, we discover that strong and fragile topological Dirac semimetals exhibit one-dimensional (1D) higher-order hinge Fermi arcs (HOFAs) as universal, direct consequences of their bulk 3D Dirac points. To predict HOFAs coexisting with topological surface states in solid-state Dirac semimetals, we introduce and layer a spinful model of an s – d -hybridized quadrupole insulator (QI). We develop a rigorous nested Jackiw–Rebbi formulation of QIs and HOFA states. Employing ab initio calculations, we demonstrate HOFAs in both the room- ( α ) and intermediate-temperature ( α ″ ) phases of Cd 3 As 2 , KMgBi, and rutile-structure ( $$ \beta ^{\prime} $$ β ′ -) PtO 2 .

Science & Technology - Other Topics↗

Critical behavior in a chiral molecular model

Understanding the condensed-phase behavior of chiral molecules is important in biology as well as in a range of technological applications, such as the manufacture of pharmaceuticals. Here, we use molecular dynamics simulations to study a chiral four-site molecular model that exhibits a second-order symmetry-breaking phase transition from a supercritical racemic liquid into subcritical D-rich and L-rich liquids. We determine the infinite-size critical temperature using the fourth-order Binder cumulant, and we show that the finite-size scaling behavior of the order parameter is compatible with the 3D Ising universality class. We also study the spontaneous D-rich to L-rich transition at a slightly subcritical temperature of T = 0.985Tc, and our findings indicate that the free energy barrier for this transformation increases with system size as N2/3, where N is the number of molecules, consistent with a surface-dominated phenomenon. The critical behavior observed herein suggests a mechanism for chirality selection in which a liquid of chiral molecules spontaneously forms a phase enriched in one of the two enantiomers as the temperature is lowered below the critical point. Furthermore, the increasing free energy barrier with system size indicates that fluctuations between the L-rich and D-rich phases are suppressed as the size of the system increases, trapping it in one of the two enantiomerically enriched phases. Such a process could provide the basis for an alternative explanation for the origin of biological homochirality. We also conjecture the possibility of observing nucleation at subcritical temperatures under the action of a suitable chiral external field.

Chemistry↗

Thermodynamics and magnetism in the two-dimensional to three-dimensional crossover of the Hubbard model

The realization of antiferromagnetic (AF) correlations in ultracold fermionic atoms on an optical lattice is a significant achievement. Experiments have been carried out in one, two, and three dimensions, and have also studied anisotropic configurations with stronger tunneling in some lattice directions. Such anisotropy is relevant to the physics of cuprate superconductors and other strongly correlated materials. Moreover, this anisotropy might be harnessed to enhance AF order. Here we numerically investigate, using the determinant quantum Monte Carlo method, a simple realization of anisotropy in the three-dimensional (3D) Hubbard model in which the tunneling between planes, $t_⊥$, is unequal to the intraplane tunneling $\textit{t}$. This model interpolates between the three-dimensional isotropic ($t_⊥ = t$) and two-dimensional (2D;$t_⊥$ = 0) systems. In this work, we show that at fixed interaction strength to tunneling ratio ($\textit{U/t}$), anisotropy can enhance the magnetic structure factor relative to both 2D and 3D results. However, this enhancement occurs at interaction strengths below those for which the Néel temperature $T_{\text{Néel}}$ is largest, in such a way that the structure factor cannot be made to exceed its value in isotropic 3D systems at the optimal $\textit{U/t}$. We characterize the 2D-3D crossover in terms of the magnetic structure factor, real space spin correlations, number of doubly occupied sites, and thermodynamic observables. An interesting implication of our results stems from the entropy's dependence on anisotropy. As the system evolves from 3D to 2D, the entropy at a fixed temperature increases. Correspondingly, at fixed entropy, the temperature will decrease going from 3D to 2D. This suggests a cooling protocol in which the dimensionality is adiabatically changed from 3D to 2D.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Technique to Measure Coronal Electron Density, Temperature, and Velocity Above 2.5 R⨀ from Sun Center using Polarized Brightness Spectrum

The current model for the polarized brightness (pB) spectrum has a decades-long history of progressively incorporating its dependence on electron density Ne, temperature Te, and flow velocity in the radial direction V e . The pB N e spectrum follows the exact shape of the photosphere spectrum, which is not smooth, which is expected from the thermal Doppler broadening of the photosphere spectrum due to the high coronal T e ; the pB N e T e spectrum is smooth, but the free coronal electrons remain static and unaffected by solar wind, and the pB N e T e V e spectrum is red-shifted by electrons seeing a red-shifted photosphere spectrum as they flow away from the Sun as solar wind, which takes a radial direction above 2.5 R ⨀ from Sun center. In this article, we review the progress of the above three model pB spectra in describing the observations and highlight the differences, first by comparing the three model pB spectra against wavelength using a model for Ne and constant values for T e and V e , and second by generating three model 2D pB maps by integrating over a selected wavelengthregion in the three model pB spectra along lines of sight passing through the 14 July 2000 (“Bastille Day”) coronal mass ejection (CME) model, which contains 3D information on N e , T e , and V e . In this regard, the COronal Diagnostic EXperiment (CODEX) on the International Space Station (ISS) in 2024 will measure N e , T e , and V e by matching the measured pB with modeled pB N e T e V e in selected wavelength regions using multiple filters.

Electron density↗

HPC-Enabled Optimization of High Temperature Heat Exchangers (CRADA Final Report)

This project was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Materials Sciences, LLC, to develop a technology for design and optimization of heat exchangers using powerful desktop and laptop computers. The project was originally designated as a 12-month project, and consisted of 3# major tasks and the following 8# major deliverables: 1) CFD models of 3D heat exchangers based on existing and new geometry. 2) Validation against experimental data provided by MSC and published in the literature. 3) CFD models of 3D unit cells based on TPMS. 4) Surrogate models capable of delivering the gradients of the homogenized properties with respect to the parametrization. 5) 3D design methodology using TO algorithms. 6) Conventional reference and topology optimized designs. 7) 3D optimized designs stored in a 3D printer build format. 8) Verification of the improved performance. All of the deliverables for this project were successfully completed with two no-cost time extensions.

36 MATERIALS SCIENCE↗

HPC-Enabled Optimization of High Temperature Heat Exchangers (CRADA Final Report)

This project was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Materials Sciences, LLC, to develop a technology for design and optimization of heat exchangers using powerful desktop and laptop computers. The project was originally designated as a 12-month project, and consisted of 3# major tasks and the following 8# major deliverables: 1) CFD models of 3D heat exchangers based on existing and new geometry. 2) Validation against experimental data provided by MSC and published in the literature. 3) CFD models of 3D unit cells based on TPMS. 4) Surrogate models capable of delivering the gradients of the homogenized properties with respect to the parametrization. 5) 3D design methodology using TO algorithms. 6) Conventional reference and topology optimized designs. 7) 3D optimized designs stored in a 3D printer build format. 8) Verification of the improved performance. All of the deliverables for this project were successfully completed with two no-cost time extensions.

13 HYDRO ENERGY↗

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗