Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Ground-Ground Data Communication-Assisted Planning and Coordination: Shorter Verbal Communications

A human-in-the-loop simulation was conducted to investigate the operational feasibility, technical requirements, and potential improvement in airspace efficiency of adding a Multi-Sector Planner position. A subset of the data from that simulation is analyzed here to determine the impact, if any, of ground-ground data communication (Data Comm) on verbal communication and coordination for multi-sector air traffic management. The results suggest that the use of Data Comm significantly decreases the duration of individual verbal communications. The results also suggest that the use of Data Comm, as instantiated in the current simulation, does not obviate the need for accompanying voice calls.

Kessell, Angela Mary↗
NASA NTRS1995

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗
NASA NTRS1995

Feasibility Study of SSTO Base Heating Simulation in Pulsed-Type Facilities

A laboratory simulation of the base heating environment of the proposed reusable Single-Stage-To-Orbit vehicle during its ascent flight was proposed. The rocket engine produces CO2 and H2, which are the main combustible components of the exhaust effluent. The burning of these species, known as afterburning, enhances the base region gas temperature as well as the base heating. To determine the heat flux on the SSTO vehicle, current simulation focuses on the thermochemistry of the afterburning, thermophysical properties of the base region gas, and ensuing radiation from the gas. By extrapolating from the Saturn flight data, the Damkohler number for the afterburning of SSTO vehicle is estimated to be of the order of 10. The limitations on the material strengths limit the laboratory simulation of the flight Damkohler number as well as other flow parameters. A plan is presented in impulse facilities using miniature rocket engines which generate the simulated rocket plume by electric ally-heating a H2/CO2 mixture.

Park, Chung Sik↗
NASA NTRS1992

The use of typed lambda calculus for comprehension and construction of simulation models in the domain of ecology

We are concerned with two important issues in simulation modelling: model comprehension and model construction. Model comprehension is limited because many important choices taken during the modelling process are not documented. This makes it difficult for models to be modified or used by others. A key factor hindering model construction is the vast modelling search space which must be navigated. This is exacerbated by the fact that many modellers are unfamiliar with the terms and concepts catered to by current tools. The root of both problems is the lack of facilities for representing or reasoning about domain concepts in current simulation technology. The basis for our achievements in both of these areas is the development of a language with two distinct levels; one for representing domain information, and the other for representing the simulation model. Of equal importance, is the fact that we make formal connections between these two levels. The domain we are concerned with is ecological modelling. This language, called Elklogic, is based on the typed lambda calculus. Important features include a rich type structure, the use of various higher order functions, and semantics. This enables complex expressions to be constructed from relatively few primitives. The meaning of each expression can be determined in terms of the domain, the simulation model, or the relationship between the two. We describe a novel representation for sets and substructure, and a variety of other general concepts that are especially useful in the ecological domain. We use the type structure in a novel way: for controlling the modelling search space, rather than a proof search space. We facilitate model comprehension by representing modelling decisions that are embodied in the simulation model. We represent the simulation model separately from, but in terms of a domain mode. The explicit links between the two models constitute the modelling decisions. The semantics of Elklogic enables English text to be generated to explain the simulation model in domain terms.

Uschold, Michael↗
NASA NTRS2009

Return to Flying Duties Following Centrifuge or Vibration Exposures

Introduction: In an effort to determine the human performance limits for vibration in spacecraft being developed by NASA, astronauts were evaluated during a simulated launch profile in a centrifuge/vibration environment and separate vibration-only simulation. Current USAF and Army standards for return to flight following centrifuge exposures require 12-24 hours to pass before a crewmember may return to flying duties. There are no standards on vibration exposures and return to flying duties. Based on direct observation and provocative neurological testing of the astronauts, a new standard for return to flying duties following centrifuge and/or vibration exposures was established. Methods: 13 astronaut participants were exposed to simulated launch profiles in a + 3.5 Gx bias centrifuge/vibration environment and separately on a vibration table at the NASA-Ames Research Center. Each subject had complete neurological evaluations pre- and post-exposure for the centrifuge/vibration runs with the NASA neurological function rating scale (NFRS). Subjects who participated in the vibration-only exposures had video oculography performed with provocative maneuvers in addition to the NFRS. NFRS evaluations occurred immediately following each exposure and at 1 hour post-run. Astronauts who remained symptomatic at 1 hour had repeat NFRS performed at 1 hour intervals until the crewmember was asymptomatic. Results: Astronauts in the centrifuge/vibration study averaged a 3-5 point increase in NFRS scores immediately following exposure but returned to baseline 3 hours post-run. Subjects exposed to the vibration-only simulation had a 1-3 point increase following exposure and returned to baseline within 1-2 hours. Pre- and post- vibration exposure video oculography did not reveal any persistent ocular findings with provocative testing 1 hour post-exposure. Discussion: Based on direct observations and objective measurement of neurological function in astronauts following simulated launch profiles, asymptomatic individuals are allowed to return to flying duties within 3 hours following centrifuge/vibration and 2 hours after vibration-only exposures.

Scheuring, Richard A.↗
DOE OSTI2025

Accurate and efficient predictions of keyhole dynamics in laser materials processing using machine learning-aided simulations

The keyhole phenomenon has been widely observed in laser materials processing, including laser welding, remelting, cladding, drilling, and additive manufacturing. Keyhole-induced defects, primarily pores, dramatically affect the performance of final products, impeding the broad use of these laser-based technologies. The formation of these pores is typically associated with the dynamic behavior of the keyhole. So far, the accurate characterization and prediction of keyhole features, particularly keyhole depth, as a function of time, has been a challenging task. In situ characterization of keyhole dynamic behavior using the synchrotron X-ray technique is informative but complicated and expensive. Current simulations are generally hindered by their poor accuracy and generalization abilities in predicting keyhole depths due to the lack of accurate laser absorptance data. In this study, we develop a machine learning-aided simulation method that accurately predicts keyhole dynamics, especially in keyhole depth fluctuations, over a wide range of processing parameters. In two case studies involving titanium and aluminum alloys, we achieve keyhole depth prediction with a mean absolute percentage error of 10 %, surpassing those simulated using the ray-tracing method with an error margin of 30 %, while also reducing computational time. This exceptional fidelity and efficiency empower our model to serve as a cost-effective alternative to synchrotron experiments. Our machine learning-aided simulation method is affordable and readily deployable for a large variety of materials, opening new doors to eliminate or reduce defects for a wide range of laser materials processing techniques.

Computational fluid dynamics↗
DOE OSTI2021

Entanglement and many-body effects in collective neutrino oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings, such as supernovae, where the neutrino density is large. In this regime, neutrino-neutrino interactions are important and simulations in the mean-field approximation show evidence for collective oscillations occurring at timescales much shorter than those associated with vacuum oscillations. In this work, we study the out-of-equilibrium dynamics of a corresponding spin model using matrix product states and show how collective bipolar oscillations can be triggered by many-body correlations if appropriate initial conditions are present. We find entanglement entropies scaling at most logarithmically in the system size suggesting that classical tensor network methods could be efficient in describing collective neutrino dynamics more generally. These observation provide a clear path forward, not only to increase the accuracy of current simulations, but also to elucidate the mechanism behind collective flavor oscillations without resorting to the mean-field approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Modeling and Simulation to Support Current and Advanced Reactors : Thermal Hydraulics

Presentation on thermohydraulic modeling and simulation capabilities in support of current and advanced reactors. The presentation includes information on the capabilities of RELAP5-3D, MOOSE and the experimental testbed for this purpose. This presentation is intended for a Taiwan delegation that will be visiting INL in August 2024.

42 ENGINEERING↗
NASA NTRS2013

Reduction of Non-uniform Beam Filling Effects by Vertical Decorrelation: Theory and Simulations

Algorithms for estimating precipitation rates from spaceborne radar observations of apparent radar reflectivity depend on attenuation correction procedures. The algorithm suite for the Ku-band precipitation radar aboard the Tropical Rainfall Measuring Mission satellite is one such example. The well-known problem of nonuniform beam filling is a source of error in the estimates, especially in regions where intense deep convection occurs. The error is caused by unresolved horizontal variability in precipitation characteristics such as specific attenuation, rain rate, and effective reflectivity factor. This paper proposes the use of vertical decorrelation for correcting the nonuniform beam filling error developed under the assumption of a perfect vertical correlation. Empirical tests conducted using ground-based radar observations in the current simulation study show that decorrelation effects are evident in tilted convective cells. However, the problem of obtaining reasonable estimates of a governing parameter from the satellite data remains unresolved.

radar reflectivity↗
DOE OSTI2023

Exploration of non-resonant divertor features on the Compact Toroidal Hybrid

Abstract Non-resonant divertors (NRDs) separate the confined plasma from the surrounding plasma facing components (PFCs). The resulting striking field line intersection pattern on these PFCs is insensitive to plasma equilibrium effects. However, a complex scrape-off layer (SOL), created by chaotic magnetic topology in the plasma edge, connects the core plasma to the PFCs through varying magnetic flux tubes. The Compact Toroidal Hybrid (CTH) serves as a test-bed to study this by scanning across its inductive current. Simulations observe a significant change of the chaotic edge structure and an effective distance between the confined plasma and the instrumented wall targets. The intersection pattern is observed to be a narrow helical band, which we claim is a resilient strike line pattern. However, signatures of finger-like structures, defined as heteroclinic tangles in chaotic domains, within the plasma edge connect the island chains to this resilient pattern. The dominant connection length field lines intersecting the targets are observed via heat flux modeling with EMC3-EIRENE. At low inductive current levels, the excursion of the field lines resembles a limited plasma wall scenario. At high currents, a private flux region is created in the area where the helical strike line pattern splits into two bands. These bands are divertor legs with distinct SOL parallel particle flow channels. The results demonstrate the NRD strike line pattern resiliency within CTH, but also show the underlying chaotic edge structure determining if the configuration is diverted or limited. This work supports future design efforts for a mechanical structure for the NRD.

Physics↗
DOE OSTI2022

Improved Vertical Carrier Transport for Green III-Nitride LEDs Using ( In , Ga ) N Alloy Quantum Barriers

We report on experimental and simulation-based results using ( In , Ga ) N alloy quantum barriers in c -plane green light-emitting diode (LED) structures as a means to improve vertical carrier transport and reduce forward voltage ( V F ) . Three-dimensional device simulations that include random alloy fluctuations are used to understand carrier behavior in a disordered potential. The simulated current density–voltage ( J - V ) characteristics and modified electron-hole overlap | F mod | 2 indicate that increasing the indium fraction in the ( In , Ga ) N quantum barriers leads to a reduced polarization discontinuity at the interface between the quantum barrier and quantum well, thereby reducing V F and improving | F mod | 2 . Maps of electron and hole current through the device show a relatively homogenous distribution in the X Y plane for structures using Ga N quantum barriers; in contrast, preferential pathways for vertical transport are identified in structures with ( In , Ga ) N barriers as regions of high and low current. A positive correlation between hole (electron) current in the p -side ( n -side) barrier and indium fraction reveals that preferential pathways exist in regions of high indium content. Furthermore, a negative correlation between the strain ε z z and indium fraction shows that high indium content regions have reduced strain-induced piezoelectric polarization in the Z direction due to the mechanical constraint of the surrounding lower indium content regions. Experimentally, multiple quantum well green LEDs with ( In , Ga ) N quantum barriers exhibit lower V F and blue-shifted wavelengths relative to LEDs with Ga N quantum barriers, consistent with simulation data. These results can be used to inform heterostructure design of low V F , long-wavelength LEDs and provide important insight into the nature of carrier transport in III-nitride alloy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Noise Prediction of Multi-stream Internally Mixed Jets with External and Internal Plugs using Large-Eddy Simulation

This paper presents results from ongoing research on jet noise prediction using wall-modeled large eddy simulations (WMLES) conducted within the Launch Ascent and Vehicle Aerodynamics (LAVA) computational framework. The study primarily focuses on the aeroacoustic implications of multi-stream nozzle configurations with internal mixing and both internal and external plugs, at a Reynolds number of 1×10 6 based on the nozzle exit diameter. While internal mixing nozzles have long been considered for jet noise reduction, the complexities of their impact on overall noise levels remain insufficiently understood. This research applies established best practices for WMLES to these complex nozzle configurations, aiming to assess their efficacy and identify limitations in accurately predicting noise behaviors. Through detailed comparisons with experimental data obtained from NASA’s Glenn Research Center, initial findings underscore the challenges inherent in current simulation practices when confronted with intricate geometric and operational conditions. Responding to these challenges, the study explores innovative computational approaches to rectify discrepancies noted between experimental outcomes and CFD predictions. Once example is the introduction of a resonating sound source in the simulation environment to mimic potential unobserved acoustic phenomena.A second example is the strategic modifications to the geometry of external plugs to account for real-world deformations caused by heating and gravity. These novel strategies aim to enhance the accuracy and reliability of noise predictions from supersonic jets, advancing our understanding and capability to effectively reduce jet noise in commercial supersonic aircraft

CST↗
NASA NTRSDate not supplied

Human Factors Impact of Different Ramp Controller Scheduling Advisories for ATD-2 Surface Metering in a Human-in-the-Loop Simulation

A human-in-the-loop (HITL) simulation took place in 2019 at NASA Ames Research Center’s Future Flight Central, an airport simulator with a full-scale 360-degree view that showed real-time airport operations at Dallas/Fort Worth International Airport (DFW). The HITL was part of the NASA/FAA Airspace Technology Demonstration 2 (ATD-2) project which involves an on-going field demonstration at Charlotte Douglas International Airport (CLT). In the field demonstration, metering on the airport surface takes place with one of the goals being to reduce runway queue length by holding aircraft at the gate during busy times, thereby reducing fuel burn and CO2 emissions. This goal has been successfully accomplished at CLT by showing Ramp Controllers when to release aircraft at the gate so as to reduce aircraft waiting in the queue. Among the questions that the current simulation at DFW was designed to address were, "Which types of advisories to the Ramp Controllers work best, those that show a ±2 minute window for a release time from the gate, or a ±5 minute window for an arrival time at the spot (where Air Traffic Control takes over) or both?" and "How do these different types of advisories affect Ramp Controller workload and situation awareness?" Another question was, “Does a ramp that is less restrictive than the one at CLT affect compliance with these advisories?”

Bonny K Parke↗
NASA NTRS1993

The giant impact produced a precipitated Moon

The author's current simulations of Giant Impacts on the protoearth show the development of large hot rock vapor atmospheres. The Balbus-Hawley mechanism will pump mass and angular momentum outwards in the equatorial plane; upon cooling and expansion the rock vapor will condense refractory material beyond the Roche distance, where it is available for lunar formation. During the last seven years, the author together with several colleagues has carried out a series of numerical investigations of the Giant Impact theory for the origin of the Moon. These involved three-dimensional simulations of the impact and its aftermath using Smooth Particle Hydrodynamics (SPH), in which the matter in the system is divided into discrete particles whose motions and internal energies are determined as a result of the imposed initial conditions. Densities and pressures are determined from the combined overlaps of the particles, which have a bell-shaped density distribution characterized by a smoothing length. In the original series of runs all particle masses and smoothing lengths had the same values; the matter in the colliding bodies consisted of initial iron cores and rock (dunite) mantles. Each of 41 runs used 3,008 particles, took several weeks of continuous computation, and gave fairly good representations of the ultimate state of the post-collision body or bodies but at best crude and qualitative information about individual particles in orbit. During the last two years an improved SPH program was used in which the masses and smoothing lengths of the particles are variable, and the intent of the current series of computations is to investigate the behavior of the matter exterior to the main parts of the body or bodies subsequent to the collisions. These runs are taking times comparable to a year of continuous computation in each case; they use 10,000 particles with 5,000 particles in the target and 5,000 in the impactor, and the particles thus have variable masses and smoothing lengths (the latter are dynamically adjusted so that a particle typically overlaps a few tens of its neighbors). Since the matter in the impactor provides the majority of the mass left in orbit after the collision, and since the masses of the particles that originated in the impactor are smaller than those in the target, the mass resolution in the exterior parts of the problem is greatly improved and the exterior particles properly simulate atmospheres in hydrostatic equilibrium.

Cameron, A. G. W.↗
DOE OSTI2024

Multiphysics Demonstration of Temperature-Driven Assembly Bowing in SFRs using MOOSE-Based Codes

Core bowing is an important passive safety mechanism in liquid metal cooled fast reactors. When the core restraint system is properly designed, temperature and flux gradients influence assemblies in the core to bow into less reactive configurations during accident scenarios, resulting in negative reactivity feedback. Prediction of core bowing involves complex interplay of radiation transport, impacts of fluid flow and heat transfer on duct temperature, and mechanical responses to the induced temperature and flux gradients. Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program [1], an integrated multiphysics approach is being developed to model the core bowing phenomena in liquid metal-cooled fast reactors with the Multiphysics Object Oriented Simulation Environment (MOOSE) [2]. In this methodology, the MOOSE-based reactor physics code Griffin [3] will solve the neutron transport equation and determine the power distribution. With the detailed power distribution from Griffin, the subchannel analysis codes MOOSE-Subchannel [4] and Pronghorn [5] are utilized to calculate the assembly temperature distribution. MOOSE’s Solid Mechanics [6] and Contact [7] Modules are leveraged to calculate the thermal expansion and duct bowing displacement with the duct wall temperature from thermal hydraulics calculation. In this work, an initial one-way coupling demonstration of the integrated multiphysics approach has been performed on a seven-assembly problem based on the sodium-cooled fast reactor ABR-1000 design [8]. The neutronics calculation with Griffin is not yet involved in the current simulation. MOOSE-Subchannel and Pronghorn evaluate fluid and solid temperature based on a fixed power distribution. In addition, one-way coupling is utilized in this coupled calculation, via Pronghorn passing the duct temperature data to the MOOSE Solid Mechanics calculation. An assessment of the Solid Mechanics module was performed in parallel to verify duct bowing behavior with duct-to-duct contact phenomenon [9]. The displacement from MOOSE Solid Mechanics is not yet transferred back and utilized in the Pronghorn and MOOSE-Subchannel calculation. This model will be available on the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) repository [10]. Future stages of this work will involve solving problems of increasing complexity as well as adding more physics (e.g. reactor physics) to the integrated workflow to reach the end goal of modeling the core bowing phenomenon with an integrated multiphysics workflow.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2006

Comanche Helmet-Mounted Display Symbology Simulation

The Aeroflightdynamics Directorate (AMRDEC) conducted a simulation to examine the performance of the Comanche Contact Analog world-referenced symbology displayed on the Comanche HIDSS when compared with a compressed symbology design similar to that specified by the former MIL-STD 1295. Six experimental test pilots flew one modified ADS-33 maneuver (hover turn, bob-up), an unusual attitude recovery, and two terrain flight tactical tasks in the NASA Vertical Motion Simulator (VMS). Analysis of the pilot objective performance data and subjective data showed the following results. Objective test results showed that 1295 symbology yielded more rapid maneuvering in the hover turn bob-up than Contact Analog symbology. The average margin of difference in the time to complete the maneuver was approximately two seconds, which was statistically significant. There were no significant differences measured between symbology sets with respect to altitude or position performance measures for all other maneuvers. The NOE target ID task data showed improved accuracy in determining heading to target when using Contact Analog over MIL-STD-1295. Subjective test results, including handling qualities ratings (HQRs) and NASA-TLX workload ratings, showed small but consistent advantages of 1295 symbology over Contact Analog for most parameters. For the bob-up maneuver, 1295 symbology handling qualities were rated Desired for lateral position error and time to complete whereas Contact Analog was rated adequate. The average HQRs for all other maneuvers were rated the same for both symbology sets. Pilot comments and the results of an online questionnaire more strongly favored 1295 over Contact Analog. Repeated comments from all six pilots led to a focus on design issues with six Contact Analog symbols. Those symbols were the heading tape, horizon line, radar altitude six-second predictor, the position of the torque symbol, the absence of a hover position cue, and the widespread positioning of the symbology to the outer edges of the display field-of-view. Pilots rated the present design of three of the six symbols as having safety-of-flight implications. Those symbols were the heading tape, horizon line, and six-second predictor. To summarize, test results showed no objective data that would warrant restricting experimental test pilots from flying constrained tasks. Small but consistent advantages were recorded for the MILSTD-1295 symbology design when executing ADS-33 constrained tasks. No general performance differences were recorded for the operational maneuvers (NOE, contour flight modes) except for the azimuth-to-target task, which favored the Contact Analog heading tape design. There were consistent and strong pilot comments supporting the MIL-STD-1295 design over Contact Analog in both this current simulation and in Comanche Sim I. Three Contact Analog symbols warrant modification and further evaluation to mitigate safety-of-flight implications noted by participating test pilots. A plan for symbology redesign and testing was developed. However, the Comanche program was cancelled by Department of Army in February 2004.

Turpin, Terry S.↗
DOE OSTI2026

Observable optimization for precision theory: machine learning energy correlators

The practice of collider physics typically involves the marginalization of multi-dimensional collider data to uni-dimensional observables relevant for some physics task. In many cases, such as classification or anomaly detection, the observable can be arbitrarily complicated, such as the output of a neural network. However, for precision measurements, the observable must correspond to something computable systematically beyond the level of current simulation tools. In this work, we demonstrate that precision-theory-compatible observable space exploration can be systematized by using neural simulation-based inference techniques from machine learning. We illustrate this approach by exploring the space of marginalizations of the energy 3-point correlator to optimize sensitivity to the top quark mass. We first learn the energy-weighted probability density from simulation, then search in the space of marginalizations for an optimal triangle shape. Although simulations and machine learning are used in the process of observable optimization, the output is an observable definition which can be then computed to high precision and compared directly to data without any memory of the computations which produced it. We find that the optimal marginalization is isosceles triangles on the sphere with a side ratio approximately $1 : 1 : \sqrt{2}$ (i.e. right triangles) within the set of marginalizations we consider.

Jets and Jet Substructure↗
DOE OSTI2021

Ab Initio Modeling of fcc Fe-Co-Cr-Ni High Entropy Alloys with Full Composition Range

The Fe-Co-Cr-Ni quaternary system has been studied extensively in the past decade, not only because of the superior properties achieved like high tensile ductility and fracture toughness, but also it is the foundation for the development of FeCrCoNi-based HEAs. However, most of the investigations are mainly focusing on the equiatomic and semi-equiatomic compositions. The physical properties of non-equiatomic compositions are barely explored. In the current work, the ab initio approach is adopted to predict the lattice parameter, structure stabilities, elastic properties, and enthalpy of formations of fcc Fe-Cr-Co-Ni HEA single crystals with the special quasirandom structure (SQS) method. To expand the design of Fe-Co-Cr-Ni HEAs into a full composition range, the current simulations cover all the binaries, ternaries, and the whole quaternary system. The 2-D composition-property contour diagrams are created to show the physical properties in the ternary and quaternary systems. Furthermore, our investigation shows that the fcc alloys' lattice parameter increases with the addition of Fe and Cr and decreases with Ni and Co. In addition to the enthalpy of formation at o K, the shape distortion rate is proposed to be another essential parameter to evaluate crystal stability. The results show the addition of Cr would destabilize the fcc lattice and cause symmetry breaking. Finally, the result of Pugh’s ratio shows that most of the fcc alloys show ductile behavior (P > 1.75), especially for the alloys around Fe 3 Ni, FeCo, and Fe 2 NiCo (P > 2.5). The brittle behavior (P < 1.75) is located around Fe 3 Cr and Fe 2 CoCr. Meanwhile, Cr 2 NiCo and Cr 2 FeNiCo are considered promising compositions to be tried and verified experimentally because of the higher bulk modulus and moderate shear modulus.

36 MATERIALS SCIENCE↗