Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Reduced Order 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.

608 records · Page 2

Thrust Chamber Thermal Barrier Coating Techniques

Methods for applying thermal barrier coatings to the hot-gas side wall of rocket thrust chambers in order to significantly reduce the heat transfer in high heat flux regions has been the focus of technology efforts for many years. A successful technique developed by NASA-Lewis that starts with the coating on a mandrel and then builds the thrust chamber around it by electroforming appropriate materials is described. This results in a smooth coating with exceptional adherence, as was demonstrated in hot fire rig tests. The low cycle fatigue life of chambers with coatings applied in this manner was increased dramatically compared to uncoated chambers.

Richard J Quentmeyer

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda

Ground Slope Effects on Lander Radiator Performance

Two sensitivity studies were performed to examine the drop in radiator heat rejection from slopes. Study #1 used a 5.5m tall reduced order asset to represent a small-lander/rover. Study #2 used a 20m asset represent a tall lander. The assets had body mounted horizontal, 45º inward-tilted, and vertical radiator orientations. The impact of slopes on different radiator orientations and on radiators pointing different directions (i.e., slope verses sun) was investigated. The studies were conducted at a polar location of -82.5ºS to capture the incident solar flux at the northmost, and therefore the hottest, range of potential Artemis landing sites (Mons Mouton Plateau).

Lunar Slope

Spray Combustion Stability

The central purpose of this project is the improvement of liquid-fueled rocket motor design technology in order to assist the establishment of economical commercial access to space through the development of engines with enhanced performance and reliability. Specific research effort is focused on spray physics and associated combustion instability phenomena. Results concerning high pressure droplet gasification model, droplet turbulent dispersion model, and spray atomization model will contribute to the development of new computational tools for design of stable liquid propellant rocket engines.

Pak-Yan Liang

Technology Assessment and Modelling of Three Terminal Tandem Solar Cells for In-Space Utilization

NASA seeks a sustained human presence on the lunar surface. High-efficiency and resilient energy systems are critical to supporting habitats, scientific outposts, and lunar surface operations. Photovoltaics have long been the backbone of space power systems, transitioning from silicon solar cells to multi-junction III-V solar cells, which now dominate state-of-the-art (SOA) technology. Multi-junction cells achieve high efficiency by layering semiconductors, each absorbing a specific portion of the solar spectrum. However, efficiency gains are increasingly constrained by device physics and manufacturing complexity. Additionally, higher-order junctions pose challenges for accurate modelling and performance characterization due to difficulties in replicating the AM0 spectrum.

solar cells

Laser Beam Welding Benchmark Experiments Performed in Reduced Gravity and Vacuum

Laser beam welding (LBW) is affected by the extreme temperatures, reduced pressure, and reduced gravity present in space environments. Gravity and pressure especially influence its melt pool and solidification dynamics. A compact, modular vacuum chamber adaptable to flight platforms from parabolic to orbital currently hosts an experiment to investigate the combined influence of reduced gravity and pressure on LBW. A swappable cartridge contains a rotating platen on which customizable workpieces can be welded under vacuum, greatly increasing experimental throughput. Instrumentation includes weld and thermal cameras observing the process, thermocouples placed on workpieces, accelerometers, and vacuum sensors. Experimental data gathered during the welding process will be combined with post-flight nondestructive evaluation, metallography, and mechanical testing to provide validation datasets for computational modeling. Phase I of this effort involves a parabolic flight campaign in low gravity while an anticipated Phase II would proceed to in-space demonstration to access extended duration microgravity.

in-space welding

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster

Spacecraft Radiator Protection from Ionizing Radiation, Dust, and Excessive Heat Loss

Under a Phase II SBIR project funded by NASA Johnson Space Center (Contract No. 80NSSC25C0088), Analytical Scientific Products LLC (ASP) has been developing an actively controlled louver to protect spacecraft radiators from degradation due to exposure to various types of environmental conditions. Of particular interest are ionizing radiation during spacecraft transit through the Van Allen belts, dust during spacecraft landing and surface operations on the moon, and excessive heat loss during the long lunar night especially near the poles where the local ambient temperatures can drop below -200°C. Exposure to ionizing radiation and dust can degrade the optical properties of the radiator coating that in turn reduces its ability to reject excess heat from the spacecraft into the environment. Exposure to the extremely low temperature conditions during the lunar night near the poles can freeze the radiator fluids that can compromise the integrity of the radiator. Passive louvers constructed using thick aluminum vanes are used currently to protect spacecraft radiators from some of the above effects, but they have several drawbacks: (i) their high aerial density makes it impractical to scale them to protect the much larger spacecraft radiators needed for future manned missions to the moon and beyond, (ii) the bimetallic actuators used to open and close the vanes in passive louvers rely on external temperature alone and so the louver cannot offer protection from dust and ionizing radiation when the ambient temperatures are high, and (iii) the bimetallic actuators need time scales of the order of hours to open and close. We have addressed the above limitations of passive louvers through the development of a low aerial density and rapidly actuating actively controlled louver. It is constructed from a custom alloy that offers much higher levels of protection against the ionizing radiation, dust and excessive heat loss at a fraction of weight compared to passive louvers while providing opening and closing time scales of the order of a second. Our modular design makes it easy to scale the system up or down depending on the application. We have recently constructed a 31 inch × 31 inch module of this louver and tested its functionality and effectiveness. These tests have shown that the louver can be opened and closed in less than 2 s. It reduces the transmission of ionizing radiation by 78% to 100% (depending on the radiation source), dust transmission by more than 93%, and heat loss by more than 97%. We are currently preparing to evaluate this device at the Johnson Space Center’s cryogenic vacuum chamber where it can be subjected to simulated lunar surface conditions. This paper provides a detailed discussion of the test designs as well as the data from tests that demonstrate the effectiveness in reducing the transmission of ionizing radiation, dust, and heat under laboratory conditions.

Radiator

Spacecraft Radiator Protection from Ionizing Radiation, Dust, and Excessive Heat Loss

Under a Phase II SBIR project funded by NASA Johnson Space Center (Contract No. 80NSSC25C0088), Analytical Scientific Products LLC (ASP) has been developing an actively controlled louver to protect spacecraft radiators from degradation due to exposure to various types of environmental conditions. Of particular interest are ionizing radiation during spacecraft transit through the Van Allen belts, dust during spacecraft landing and surface operations on the moon, and excessive heat loss during the long lunar night especially near the poles where the local ambient temperatures can drop below -200°C. Exposure to ionizing radiation and dust can degrade the optical properties of the radiator coating that in turn reduces its ability to reject excess heat from the spacecraft into the environment. Exposure to the extremely low temperature conditions during the lunar night near the poles can freeze the radiator fluids that can compromise the integrity of the radiator. Passive louvers constructed using thick aluminum vanes are used currently to protect spacecraft radiators from some of the above effects, but they have several drawbacks: (i) their high aerial density makes it impractical to scale them to protect the much larger spacecraft radiators needed for future manned missions to the moon and beyond, (ii) the bimetallic actuators used to open and close the vanes in passive louvers rely on external temperature alone and so the louver cannot offer protection from dust and ionizing radiation when the ambient temperatures are high, and (iii) the bimetallic actuators need time scales of the order of hours to open and close. We have addressed the above limitations of passive louvers through the development of a low aerial density and rapidly actuating actively controlled louver. It is constructed from a custom alloy that offers much higher levels of protection against the ionizing radiation, dust and excessive heat loss at a fraction of weight compared to passive louvers while providing opening and closing time scales of the order of a second. Our modular design makes it easy to scale the system up or down depending on the application. We have recently constructed a 31 inch × 31 inch module of this louver and tested its functionality and effectiveness. These tests have shown that the louver can be opened and closed in less than 2 s. It reduces the transmission of ionizing radiation by 78% to 100% (depending on the radiation source), dust transmission by more than 93%, and heat loss by more than 97%. We are currently preparing to evaluate this device at the Johnson Space Center’s cryogenic vacuum chamber where it can be subjected to simulated lunar surface conditions. This paper provides a detailed discussion of the test designs as well as the data from tests that demonstrate the effectiveness in reducing the transmission of ionizing radiation, dust, and heat under laboratory conditions.

Radiator

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson

Analysis of Developing Laminar Flows in Circular Pipes Using a Higher-Order Finite-Difference Technique

A higher-order finite-difference technique is developed to calculate the developing-flow field of steady incompressible laminar flows in the entrance regions of circular pipes. Navier-Stokes equations governing the motion of such a flow field are solved by using this new finite-difference scheme. This new technique can increase the accuracy of the finite-difference approximation, while also providing the option of using unevenly spaced clustered nodes for computation such that relatively fine grids can be adopted for regions with large velocity gradients. The velocity profile at the entrance of the pipe is assumed to be uniform for the computation. The velocity distribution and the surface pressure drop of the developing flow then are calculated and compared to existing experimental measurements reported in the literature. Computational results obtained are found to be in good agreement with existing experimental correlations and therefore, the reliability of the new technique has been successfully tested.

Herbert J Gladden

Model Spread in Tropical Low Cloud Feedback Tied to Overturning Circulation Response to Warming

Among models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), here we show that the magnitude of the tropical low cloud feedback, which contributes considerably to uncertainty in estimates of climate sensitivity, is intimately linked to tropical deep convection and its effects on the tropical atmospheric overturning circulation. First, a reduction in tropical ascent area and an increased frequency of heavy precipitation result in high cloud reduction and upper-tropospheric drying, which increases longwave cooling and reduces subsidence weakening, favoring low cloud reduction (Radiation-Subsidence Pathway). Second, increased longwave cooling decreases tropospheric stability, which also reduces subsidence weakening and low cloudiness (Stability-Subsidence Pathway). In summary, greater high cloud reduction and upper-tropospheric drying (negative longwave feedback) lead to a more positive cloud feedback among CMIP6 models by contributing to a greater reduction in low cloudiness (positive shortwave feedback). Varying strengths of the two pathways contribute considerably to the intermodel spread in climate sensitivity.

Atmospheric dynamics

Effect of Turbulence Models on Criticality Conditions in Swirling Flows

The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon , RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Robert E Spall

An Overview of CMIP5 and CMIP6 Simulated Cloud Ice, Radiation Fields, Surface Wind Stress, Sea Surface Temperatures and Precipitation over Tropical and Subtropical Oceans

The potential links between ice water path (IWP), radiation, circulation, sea surface temperature (SST) and precipitation over the Pacific and Atlantic Oceans resulting from the falling ice radiative effects (FIREs) are examined from present day model outputs of CMIP5 and CMIP6. The latter is divided into two subsets with (SON6) and without FIREs (NOS6) as more models with FIREs are included in CMIP6 than in CMIP5. Improvement in floating cloud ice (~20 g m-2) is noticeable over convective regions in CMIP6 relative to CMIP5. The inclusion of FIREs in SON6 subset may contribute to reduce biases of overestimated outgoing longwave radiation and downward surface shortwave and overestimated reflected shortwave at the top of the atmosphere (TOA) by magnitudes of 4?8 W m-2 over convective regions against CERES, compared to NOS6 subset. The reduced biases in radiative fluxes in convective regions stabilize the atmosphere and lead to circulation, SST, cloud and precipitation changes over the trade-wind regions, as seen from improved radiative fluxes (4?15 W m-2), surface wind stress biases, SST (0.2?0.8 K) and precipitation (1 mm day-1) biases. The significant improvement from NOS6 to SON6 leads to improved multi-model means for CMIP6 relative to CMIP5 for radiation fields over the trade wind regions but the degradation over convective zones is attributed to NOS6 subset. The results suggest that other sources of uncertainty and deficiencies in climate models may play significant roles for reducing discrepancies although FIREs, via radiation-circulation coupling, may be one of the factors that help to reduce regional biases.

Jui-Lin F Li

Behavior of Ullage Bubbles During Blowdown in Low-g Experiment (BUBBLE): Overview of a Cryogenic Tank Depressurization Test

Behavior of Ullage Bubbles during Blowdown in Low-g Experiment (BUBBLE) is a cryogenic tank depressurization drop tower experiment currently being designed at NASA Glenn Research Center. There is a need for liquid level rise data while venting a cryogenic propellant tank below the liquid saturation pressure in a reduced gravity environment. When a cryogenic tank is vented to a vapor pressure below the liquid saturation pressure, bubbles become entrained in the liquid causing the bulk liquid-vapor interface to rise. In reduced gravity, bubble rise velocity is diminished leading to a larger liquid level rise compared to a 1g environment. The purpose of the experiment is to gain further understanding of tank depressurization fluid physics and obtain high-fidelity data for model validation in 1g and reduced gravity environments. On-orbit cryogenic propellant tanks operating at high fill levels must efficiently manage venting operations to avoid the risk of liquid entrainment in the vent line, which could lead to asymmetric control thruster loads, freezing and clogging, and loss of liquid propellant. Validated models could be used to design settling and venting profiles to reduce risk and increase efficiency for cryogenic storage and transfer operations.

tank venting

A Simplified Model of VIPER Thermal Management System. Part I: Loop Heat Pipe

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) uses a combination of loop heat pipes (LHPs) and heaters as the primary hardware to modulate the temperature of the electronic boxes and payloads instruments. LHPs have been designed and instrumented to be capable of passive and active shutdown. The goal of the passive shutdown, using a thermal control valve (TCV) to limit the flow of the working fluid, is to preserve thermal energy and avoid unnecessary heat leaks during lunar night. The goal of the active shutdown on the other hand, is to forcefully increase the vehicle internal temperature to bakeout the payload instruments. Several numerical modeling tools are available to analyze LHP performance under the complex set of operational requirements and environments levied on the VIPER Thermal Management System (TMS). Yet, while those numerical modeling tools are capable of detailed performance analysis, they tend to take large computational resources and long computational time. In this paper, a model of a single LHP developed based on control volume approach is presented. By using the control volume approach, the complex geometry of the LHP, and the components attached to it, are reduced to thermal lump capacitances, thus reducing the computation effort. This paper also describes the effort to correlate the proposed LHP model using thermal vacuum (TVAC) test data.

Loop Heat Pipe

A Simplified Model of VIPER Thermal Management System. Part I: Loop Heat Pipe

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) uses a combination of loop heat pipes (LHPs) and heaters as the primary hardware to modulate the temperature of the electronic boxes and payloads instruments. LHPs have been designed and instrumented to be capable of passive and active shutdown. The goal of the passive shutdown, using a thermal control valve (TCV) to limit the flow of the working fluid, is to preserve thermal energy and avoid unnecessary heat leaks during lunar night. The goal of the active shutdown on the other hand, is to forcefully increase the vehicle internal temperature to bakeout the payload instruments. Several numerical modeling tools are available to analyze LHP performance under the complex set of operational requirements and environments levied on the VIPER Thermal Management System (TMS). Yet, while those numerical modeling tools are capable of detailed performance analysis, they tend to take large computational resources and long computational time. In this paper, a model of a single LHP developed based on control volume approach is presented. By using the control volume approach, the complex geometry of the LHP, and the components attached to it, are reduced to thermal lump capacitances, thus reducing the computation effort. This paper also describes the effort to correlate the proposed LHP model using thermal vacuum (TVAC) test data.

Loop Heat Pipe

Diffusion Quantum Monte Carlo Calculation of the Austenite and Martensite Phases of NiTi

NiTi is a promising material for smart and active technologies due to its exhibition of the shape memory effect, superelasticity, and biocompatibility. The shape memory effect is tied to the reversible transition between the austenite and martensite phases. A major research direction is to alloy NiTi with Zr, Hf, Pd, Pt, etc., in order to tune the martensitic transition temperature (MTT). Modeling the MTT from first principles is challenging because the lattice dynamics is complicated by anharmonicity and various low-energy structures. Using density functional theory, the energy difference between the austenite and martensite phases of NiTi varies by up to 100 meV/atom depending on the choice of density functional, which is of the same order of the energy difference itself. Consequently, free energy calculations with different functionals can result in estimates of the MTT that vary by several hundred K. Using diffusion quantum Monte Carlo, we calculated the energy difference between the B2 and B19' structures of NiTi to be 70.9 +- 2.5 meV/atom.

Kevin K Ly