Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermomechanical 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 73 records · Page 4

Characterization of Disconnection Content Associated with Grain Boundary Motion in Tungsten

The objective of this study is to develop a new computational methodology to predict the motion of grain boundaries in tungsten associated with recrystallization. This approach accounts recently discovered high-temperature grain boundary phases in tungsten and other metals and alloys. This paper reports the initial work in the project that developed a novel technique for molecular dynamics to characterize disconnections, dislocation-like defects in grain boundaries that mediate grain boundary motion. Grain boundary properties are needed for a thermomechanical model of recrystallization of tungsten for magnetic fusion applications.

36 MATERIALS SCIENCE↗

Long-Length SiC/SiC Tube Irradiation Experiment at LWR Temperatures Inserted into HFIR

This report describes the pre-irradiation characterization and successful assembly of a High Flux Isotope Reactor (HFIR) irradiation experiment designed to assess radiation-induced lateral bowing of silicon carbide fiber–reinforced, silicon carbide matrix composite (SiC/SiC) components under a radial fast neutron flux gradient at representative light-water reactor (LWR) temperatures of approximately 300 °C. The over-arching goal of this work is to validate thermomechanical models of SiC/SiC composite fuel assembly distortion used to predict the extent of lateral bowing and dose at which bowing reaches a maximum value. Ultimately, these models are used to ensure that fuel assembly distortion does not block coolant channels, interfere with control rod/blade movements, and/or impact local reactivity. The experiment contains six cladding tube specimens that are approximately 558.5 mm long with representative pressurized water reactor (PWR) outer and inner diameters of 9.5 mm and 7.5 mm, respectively. The specimens were thoroughly characterized prior to irradiation using traditional dimensional inspection and surface profilometry so that these measurements can later be compared with similar measurements that will be made post-irradiation to determine radiation-induced deformations. Furthermore, fine engraving markers were inscribed along the outer surfaces of the specimen and mapped using analytic photography and a three-dimensional stage. This technique allowed for accurate measurements of the marker spacings, which can be compared with similar measurements that will be made post-irradiation to provide local radiation-induced strain mapping. A light curtain optical micrometer was also employed to measure the pre-irradiation bowing of each specimen. The experiment was successfully assembled and inserted into position VXF-19 for HFIR cycle 516, which began May 12, 2026. The experiment will be removed after cycle 519, which is currently scheduled to end in December 2026.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Aperiodic behaviour of a non-linear oscillator.

The aperiodic behavior of the solution of the equation of motion derived previously (1966) when considering a model thermomechanical oscillator is examined. Periodic solutions of this equation are studied numerically and analytically. Conditions for the instability of the solutions are determined. This instability seems to be the cause of the observed aperiodicity.

Baker, N. H.↗

The role of a pressure-dependent rheology in the dynamics of mantle circulation

A thermomechanical model for upper mantle convection was constructed such that the thickness and the structure of the lithosphere are determined self consistently by the heat transported by convection. In this study of the interaction between the lithosphere and upper mantle, strongly temperature and pressure dependent rheologies for both Newtonian and non-Newtonian creep mechanisms are employed. For a strictly temperature dependent rheology an insignificant amount of heat, less than 12.5 mW/sq m, can be transported convectively for an interior viscosity, 0(10 sup 21 Pas), compatible with post glacial rebound. On the other hand, for similar values of the interior viscosity, steady heat fluxes between 20 and 40 mW/sq m are produced by introducing pressure dependence into the rheology. For the temperature and pressure dependent flow law the horizontally averaged interior temperature displays very little variation with the amount of heat evacuated, once all of the rheological parameters are fixed. This finding may have important ramifications for parameterized convection.

Yuen, D. A.↗

Sunrise/sunset thermal shock disturbance analysis and simulation for the TOPEX satellite

It is shown here that during normal on-orbit operations the TOPEX low-earth orbiting satellite is subjected to an impulsive disturbance torque caused by rapid heating of its solar array when entering and exiting the earth's shadow. Error budgets and simulation results are used to demonstrate that this sunrise/sunset torque disturbance is the dominant Normal Mission Mode (NMM) attitude error source. The detailed thermomechanical modeling, analysis, and simulation of this torque is described, and the predicted on-orbit performance of the NMM attitude control system in the face of the sunrise/sunset disturbance is presented. The disturbance results in temporary attitude perturbations that exceed NMM pointing requirements. However, they are below the maximum allowable pointing error which would cause the radar altimeter to break lock.

Dennehy, C. J.↗

Statistical discrete particle simulation studies

A particle simulation code was developed that is suited to Multiple Instruction Multiple Data multiprocessor computers. The resulting code is designed to be portable between a large class of multiprocessor computer architectures, including hypercubes (e.g., Intel iPSC/860), shared memory machines (e.g., Cray, SGI, DASH), and simple uniprocessors (e.g., SUN, VAX, IBM PC). This code provides a powerful general framework that can be adapted to specific applications, such as thermomechanical modeling, general 3-D geometry support, or visualization support via a codeveloped CPlot data reduction and visualization system. Chemistry models were implemented and validated with comparison to continuum solutions of thermochemically relaxing gas mixtures and with experimental results from the high speed flow about a circular cylinder. Furthermore, the exchange of energy between translational and internal modes was modeled. New models were developed, combining greater efficiency than earlier phenomenological models. They offer greater physical detail by addressing vibrational energy as a discretely distributed quantity.

Mcdonald, Jeffrey D.↗

Multidisciplinary tailoring of hot composite structures

A computational simulation procedure is described for multidisciplinary analysis and tailoring of multilayered multimaterial hot composite engine structural components subjected to simultaneous multiple discipline-specific thermal, structural, vibration, and acoustic loadings including the effect of aggressive environments. The simulation is based on a 3D finite element analysis technique in conjunction with structural mechanics codes, thermal/acoustic analysis methods, and tailoring procedures. The integrated multidisciplinary simulation procedure is general-purpose including the coupled effects of nonlinearities in structure geometry, material, loading, and environmental complexities. The composite material behavior is assessed at all composite scales, i.e., the laminate/ply/constituents (fiber/matrix), via a nonlinear material characterization hygro-thermomechanical model. Sample tailoring cases exhibiting nonlinear material/loading/environmental behavior of aircraft engine fan blades, are presented. The various multidisciplinary loadings lead to different tailored designs, even those opposite of each other, as in the case of minimum material cost versus minimum structure weight and in the case of minimum vibration frequency versus minimum acoustic noise.

Singhal, Surendra N.↗

Ice-age Ice-sheet Rheology: Constraints from the Last Glacial Maximum Form of the Laurentide Ice Sheet

State-ot-the-art thermomechanical models of the modern Greenland ice and the ancient Laurenticle ice sheet that covered Canada at the Last Glacial Maximum (LGM) are not able to explain simultaneously the observed forms of these cryospheric structures when the same, anisotropy-enhanced, version of the conventional Glen flow law is employed to describe their rheology. The LGM Laurenticle ice sheet. predicted to develop in response to orbital climate forcing, is such that the ratio of its thickness to its horizontal extent is extremely large compared to the aspect ratio inferred on the basis of surface-geomorphological and solid-earth-geophysical constraints. We show that if the Glen flow law representation of the rheology is replaced with a new rheology based upon very high quality laboratory measurements of the stress-strain-rate relation, then the aspect ratios of both the modern Greenland ice sheet and the Laurenticle ice sheet, that existed at the LGM, are simultaneously explained with little or no retuning of the flow law.

Peltier, W. Richard↗

Analysis of SMA Hybrid Composite Structures using Commercial Codes

A thermomechanical model for shape memory alloy (SMA) actuators and SMA hybrid composite (SMAHC) structures has been recently implemented in the commercial finite element codes MSC.Nastran and ABAQUS. The model may be easily implemented in any code that has the capability for analysis of laminated composite structures with temperature dependent material properties. The model is also relatively easy to use and requires input of only fundamental engineering properties. A brief description of the model is presented, followed by discussion of implementation and usage in the commercial codes. Results are presented from static and dynamic analysis of SMAHC beams of two types; a beam clamped at each end and a cantilevered beam. Nonlinear static (post-buckling) and random response analyses are demonstrated for the first specimen. Static deflection (shape) control is demonstrated for the cantilevered beam. Approaches for modeling SMAHC material systems with embedded SMA in ribbon and small round wire product forms are demonstrated and compared. The results from the commercial codes are compared to those from a research code as validation of the commercial implementations; excellent correlation is achieved in all cases.

Turner, Travis L.↗

Evolving Material Porosity on an Additive Manufacturing Simulation with the Generalized Method of Cells

The effect of material porosity on final part distortion and residual stresses in a selective laser sintering manufacturing simulation is presented here. A time-dependent thermomechanical model is used with the open-source FEA software CalculiX. Effective homogenized material properties for Inconel 625 are precomputed using NASA’s Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC). The evolving porosity of the material is estimated with each pass of the laser beam during simulation runtime. A comparison with a homogenous model and the evolving model shows that the evolving porous model predicts larger distortions with greater residual stresses.

Silva, Luis F.↗

Melt Inclusions in Zircon: A Window to Understanding the Structure and Evolution of the Magmatic System Beneath the Laguna del Maule Volcanic Field

Explosive silicic eruptions pose a significant threat to society, yet the development and destabilization of the underlying silicic magmatic systems are still controversial. Zircons provide simultaneous information on the trace element composition and age of silicic magmatic systems, while melt inclusions in quartz and plagioclase yield important constraints on their volatile content as well as storage depth. Melt inclusions in zircons (MIZs) combine these data from a single mineral grain, recording the age, storage depth, temperature, and composition of magmas, and thus provide unique constraints on the structure and evolution of silicic magmatic systems. We studied MIZs from the Laguna del Maule (LdM) volcanic field in the southern Andes that is among the most active Pleistocene-Holocene rhyolitic volcanic centers worldwide and a potentially hazardous system displaying inflation rates in excess of 25 cm/yr. The host zircon ages suggest that the LdM MIZ record extends to ~30 kyr before eruption, in contrast to the melt inclusions in LdM plagioclase and quartz crystals that formed only decades to centuries before eruption. The major element compositions of MIZs are minimally affected by post-entrapment crystallization, and agree well with the LdM rhyolitic whole rock data. The more evolved major element composition of rle MIZs than rdm MIZs, suggests a long-term deeper connection of the rdm crystal mush to a more primitive magma body than that of the rle. The evidence of slow H diffusion observed in MIZs suggest that their H 2 O contents are not significantly affected by diffusion of H through the host zircon. The storage depths of 1.1 to 2.8 kbars recorded by the H 2 O contents of rdm and rle MIZs are consistent with the optimal emplacement window (2.0 ± 0.5 kbar) of silicic magma reservoir growth, storage, and eruptibility based on thermomechanical modeling (Huber et al. 2019).

Kei Shimizu↗

Technical Performance of Refractory Liners for Molten Chloride Salt Thermal Energy Storage Systems

A chloride-based molten-salt system that uses a ternary blend of MgCl2/KCl/NaCl is investigated to provide higher temperature thermal energy storage capability. Despite higher thermal stability, molten chlorides present several unique challenges, including the design of internal refractory-ceramic liners to prevent the corrosion and thermal stress of alloy tank shells. This work discusses issues and potential solutions related to containment of molten chloride salt, specifically the optimization of the refractory material at the molten salt interface (hot face). The down-selected hot face candidate limits permeation of salt through the material and forms a highly stable secondary surface phase in equilibrium with the molten salt. A mortar is created using the corrosion resistant hot face brick. Brick and mortar composite are subjected to mechanical stress/strain analysis, in order to calculate composite material properties and better inform thermomechanical models. The U.S. Department of Energy Generation 3 (DOE Gen3) program seeks to develop higher efficiency CSP plants that can provide cost-competitive, flexible power in the U.S. electric grid. The proposed Gen3 Liquid Pathway CSP plant closely resembles the configuration of current nitrate salt power towers with two-tank storage (Gen2). The differences between Gen2 and Gen3 include the types of compatible materials used in salt storage tank construction. Stainless steel loses strength at Gen3 temperatures, and although nickel superalloys would be capable of withstanding sustained high temperatures, these materials are prohibitively expensive at scale. Uninsulated tank shells also pose a significant risk as common steels are highly susceptible to chemical attack from molten chloride salt. To address these concerns, refractory-ceramic based containment materials are proposed to line the inside of the hot and cold storage tanks. In doing so, stainless or carbon steel shells may be used in construction depending on the level of insulation provided. The composition of the internal liner requires careful consideration to maximize the efficacy of multiple parameters including corrosion resistance, strength at operating temperature, durability, and cost. This is particularly true for the material at the interface with the salt, known as the "hot face", which is responsible for protecting the insulating layers between the tank shell and the hot face brick layer. The molten salt in this system is superheated over 300 °C above its freeze temperature. Therefore, unlike other industrial processes which use refractory-lined vessels, it is not expected that a freeze plane will develop in the hot face. Therefore, the hot face must be designed to withstand chemical corrosion and inhibit permeation of molten salt into the insulating layers. A down selection was performed to identify a hot face candidate best equipped to maintain thermal, mechanical, and chemical integrity when exposed to molten salt over extended periods of time. Long-duration chemical capability experiments were conducted with the down selected hot face refractory fully immersed in molten chloride salt for up to 3000 hours. The average salt penetration does not exceed 100 microns. When extrapolated to 20 and 30 years of continuous exposure, the expected salt penetration depth is approximately 2.0 mm and 2.9 mm, respectively. A magnesium-rich secondary phases develops at the salt/refractory interface. X-ray diffraction identifies the material as forsterite (Mg2SiO4), which is reported to form synthetically in molten chloride salt solutions. These results suggest the selected hot face will adequately inhibit salt permeation. While there is optimism that the hot face brick will inhibit salt permeation, mortar joints are typically the weakest point of a refractory brick liner. From a thermochemical perspective, differing thermal expansion coefficients may result in the mortar and brick to grow independent of each other, creating gaps through which molten salt can penetrate. To address this issue, NREL has developed an in-house mortar composed of the down selected hot face brick that has been shown to be compatible with the salt. Compressive stress/strain analyses have been performed on the brick/ mortar composites to generate stress/strain curves. Modulus of elasticity and Poisson's ratio of the composite may be calculated from the stress/strain curves, in order to provide more representative data to finite element mechanical models for accurate approximation of stress on the tank shell and the amount of thermal expansion expected within the tank liner.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Thermomechanical analysis and modeling of a high-temperature light trapping planar cavity solar receiver

The development of durable particle-based high-temperature solar receivers is critical for advancing concentrating solar-thermal (CST) technologies to enable high-efficiency power generation and industrial process heat. Here, this study presents a computational framework to evaluate the thermomechanical performance of a proposed enclosed light-trapping planar cavity receiver designed for particle-based thermal energy systems. The receiver incorporates absorptive cavities and fluidized particle-bed channels to enhance heat capture and reduce thermal losses. Finite element analysis (FEA) is employed to assess stress, strain, and creep-fatigue behavior under concentrated solar flux using realistic thermal boundary conditions derived from coupled system models and experimental assembly parameters. The analysis investigates the influence of particle-to-wall heat transfer coefficients (HTC) ranging from 800 to 1800 W/m 2 .K on the thermomechanical response of six candidate high-temperature alloys: Alloy 740H, Alloy 282, Alloy 617, 316H, Alloy 230, and 800H. Results show that increasing HTC reduces thermal gradients, leading to lower stresses and strains and extended minimum predicted creep life. While all materials satisfy fatigue life requirements under the investigated conditions, significant differences in creep resistance are observed. Alloy 740H consistently exhibits the longest minimum predicted creep life and the most favorable durability margins, followed by Alloy 282, with the remaining materials showing reduced creep resistance under identical loading. The reported creep lives are conservative lower-bound estimates intended for comparative material evaluation. This framework highlights the critical roles of material selection and geometry optimization in improving mechanical durability and reliability of solar-thermal receivers, forming a foundation for future experimental validation and design optimization.

14 SOLAR ENERGY↗

Thermomechanical deformation testing and modeling in the presence of metallurgical instabilities

A number of viscoplastic constitutive models were developed to describe deformation behavior under complex combinations of thermal and mechanical loading. Questions remain, however, regarding the validity of procedures used to characterize these models for specific structural alloys. One area of concern is that the majority of experimental data available for this purpose are determined under isothermal conditions. This experimental study is aimed at determining whether viscoplastic constitutive theories characterized using an isothermal data base can adequately model material response under the complex thermomechanical loading conditions typical of power generation service. The approach adopted was to conduct a series of carefully controlled thermomechanical experiments on a nickel-based superalloy, Hastelloy Alloy X. Previous investigations had shown that this material experiences metallurgical instabilities leading to complex hardening behavior, termed dynamic strain aging. Investigating this phenomenon under full thermomechanical conditions leads to a number of challenging experimental difficulties which up to the present work were unresolved. To correct this situation, a number of advances were made in thermomechanical testing techniques. Advanced methods for dynamic temperature gradient control, phasing control and thermal strain compensation were developed and incorporated into real time test control software. These advances allowed the thermomechanical data to be analyzed with minimal experimental uncertainty. The thermomechanical results were evaluated on both a phenomenological and microstructural basis. Phenomenological results revealed that the thermomechanical hardening trends were not bounded by those displayed under isothermal conditions. For the case of Hastelloy Alloy X (and similar dynamic strain aging materials), this strongly suggests that some form of thermomechanical testing is necessary when characterizing a thermoviscoplastic deformation model. Transmission electron microscopy was used to study the microstructural physics, and analyze the unique phenomenological behavior.

Castelli, Michael G.↗

Improved techniques for thermomechanical testing in support of deformation modeling

The feasibility of generating precise thermomechanical deformation data to support constitutive model development was investigated. Here, the requirement is for experimental data that is free from anomalies caused by less than ideal equipment and procedures. A series of exploratory tests conducted on Hastelloy X showed that generally accepted techniques for strain controlled tests were lacking in at least three areas. Specifically, problems were encountered with specimen stability, thermal strain compensation, and temperature/mechanical strain phasing. The source of these difficulties was identified and improved thermomechanical testing techniques to correct them were developed. These goals were achieved by developing improved procedures for measuring and controlling thermal gradients and by designing a specimen specifically for thermomechanical testing. In addition, innovative control strategies were developed to correctly proportion and phase the thermal and mechanical components of strain. Subsequently, the improved techniques were used to generate deformation data for Hastelloy X over the temperature range, 200 to 1000 C.

Castelli, Michael G.↗

An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem

The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.

36 MATERIALS SCIENCE↗

Development of a Numerical Model for High-Temperature Shape Memory Alloys

A thermomechanical hysteresis model for a high-temperature shape memory alloy (HTSMA) actuator material is presented. The model is capable of predicting strain output of a tensile-loaded HTSMA when excited by arbitrary temperature-stress inputs for the purpose of actuator and controls design. Common quasi-static generalized Preisach hysteresis models available in the literature require large sets of experimental data for model identification at a particular operating point, and substantially more data for multiple operating points. The novel algorithm introduced here proposes an alternate approach to Preisach methods that is better suited for research-stage alloys, such as recently-developed HTSMAs, for which a complete database is not yet available. A detailed description of the minor loop hysteresis model is presented in this paper, as well as a methodology for determination of model parameters. The model is then qualitatively evaluated with respect to well-established Preisach properties and against a set of low-temperature cycled loading data using a modified form of the one-dimensional Brinson constitutive equation. The computationally efficient algorithm demonstrates adherence to Preisach properties and excellent agreement to the validation data set.

DeCastro, Jonathan A.↗