Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

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 847 records · Page 47

Magnetic properties of a staggered S = 1 chain with an alternating single-ion anisotropy direction

Materials composed of spin-1 antiferromagnetic (AFM) chains are known to adopt complex ground states that are sensitive to the single-ion-anisotropy (SIA) energy ( D ), and intrachain ( J 0 ) and interchain ( J 1 , 2 ′ ) exchange energy scales. While theoretical and experimental studies have extended this model to include various other energy scales, the effect of the lack of a common SIA axis is not well explored. Here we investigate the magnetic properties of Ni ( pyrimidine ) ( H 2 O ) 2 ( NO 3 ) 2 , a chain compound where the tilting of Ni octahedra leads to a twofold alternation of the easy-axis directions along the chain. Muon-spin relaxation measurements indicate a transition to long-range order at T N = 2.3 K and the magnetic structure is initially determined to be antiferromagnetic and collinear using elastic neutron diffraction experiments. Inelastic neutron scattering measurements were used to find J 0 = 5.107 ( 7 ) K , D = 2.79 ( 1 ) K , J 1 ′ = 0.00 ( 5 ) K , J 2 ′ = 0.18 ( 3 ) K , and a rhombic anisotropy energy E = 0.19 ( 9 ) K . Mean-field modeling reveals that the ground state structure hosts spin canting of ϕ ≈ 6 . 5 ∘ , which is not detectable above the noise floor of the elastic neutron diffraction data. Monte Carlo simulation of the powder-averaged magnetization, M ( H ) , is then used to confirm these Hamiltonian parameters, while single-crystal M ( H ) simulations provide insight into features observed in the data. Published by the American Physical Society 2025

Vaidya, S. (ORCID:000000033975497X)↗

Impact of direct exchange coupling via the insulator on the magnetic tunnel junction based molecular spintronics devices with competing molecule induced inter-electrode coupling

The magnetic tunnel junction (MTJ) based molecular spintronics device (MTJMSD) approach is suitable for mass production. This approach provides solutions to fabrication difficulties related to reliably connecting molecular device elements to the ferromagnets (FMs). To producing MTJMSD, the molecular channels are bridged across the insulator of an MTJ testbed with exposed side edges. In an MTJMSD, two FMs are simultaneously connected by an insulator film and the molecular channels along the exposed sides. In our prior experimental studies, we observed that molecules could produce strong coupling between ferromagnets in the presence of the competing coupling via an insulator. In this paper, our Monte Carlo Simulation (MCS) was used to study the impact of coupling variation via insulator (a.k.a. Ji) on the magnetic properties of an MTJMSD. We studied the effect of Ji while varying the molecule induced antiferromagnetic exchange coupling. The ferromagnetic or antiferromagnetic nature and magnitude of Ji determined the resultant effect. Antiferromagnetic Ji enhanced the pre-existing antiferromagnetic molecular coupling effect. Ferromagnetic Ji competed with the opposite nature of antiferromagnetic molecular coupling. Our MCS may help to understand the impact of insulator thickness and defects on the molecular spintronics device performance and design process.

36 MATERIALS SCIENCE↗

Formation and coarsening of epitaxially-supported metal nanoclusters

This mini-review describes developments over the last ~30 years in characterizing the nucleation & growth of epitaxially-supported metal nanoclusters (NCs) or islands during vapor deposition, as well as their post-deposition coarsening. A beyond-mean-field treatment for homogeneous nucleation & growth corrects the deficiencies of traditional treatments in describing, e.g., the island size distribution, but also necessitates consideration of the spatial distribution of islands and their capture zones. We discuss advances in modeling capabilities, including those based upon on an ab-initio level treatment of periphery diffusion kinetics, for description of the non-equilibrium growth shapes of these NCs, focusing on 2D NCs. For post-deposition coarsening of arrays of NCs, there is generally a competition between Ostwald Ripening (OR) and Smoluchowski Ripening (SR). SR is also known as Particle Migration & Coalescence. For 2D NCs in homoepitaxial systems, conventional OR is observed on pristine fcc(111) surfaces, dramatically enhanced OR in the presence of even trace amounts of chalcogens for Cu(111) and Ag(111), and anomalous OR on anisotropic fcc(110) surfaces. The unexpected discovery of SR for fcc(100) homoepitaxial systems prompted extensive analysis of the underlying diffusivities of 2D NCs as a function of size, as well as of NC coalescence dynamics. A comprehensive understanding of these processes is now available. Self-assembly of 3D NCs during deposition, issues related to heterogeneous nucleation, directed assembly, NC growth structure selection, and coarsening are addressed. Finally, for SR of 3D epitaxial NCs, recent insights into the size-dependence of diffusivity are described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dust and gas absorption in the high mass X-ray binary IGR J16318-4848

With an absorption column density on the order of 10 24 cm -2 , IGR J16318-4848 is one of the most extreme cases of a highly obscured high mass X-ray binary. In addition to the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. We aim at a self-consistent description of the continuum absorption, the absorption edges, and the fluorescence lines to constrain the properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescence emission from dust grains. We used XMM-Newton and NuSTAR spectra to first empirically constrain the incident continuum and fluorescence lines. Next we used XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In the third step we modeled the absorption and fluorescence in a dusty olivine absorber and employed both a simple analytical model for the fluorescence line emission and a Monte Carlo simulation of radiative transfer that generates line fluxes, which are very close to the observational data. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is able to produce the observed continuum absorption and edges. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318-4848. In particular, for low optical depth of individual grains, which is the case for typical volume densities and grain size distribution models, the dust will contribute significantly to the fluorescence emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Pairing Correlations across the Superfluid Phase Transition in the Unitary Fermi Gas

In the two-component Fermi gas with a contact interaction, a pseudogap regime can exist at temperatures between the superfluid critical temperature $T_c$ and a temperature $T^\ast > T_c$. This regime is characterized by pairing correlations without superfluidity. However, in the unitary limit of infinite scattering length, the existence of this regime is still debated. To help address this, we have applied finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods to study the thermodynamics of the superfluid phase transition and signatures of the pseudogap in the spin-balanced homogeneous unitary Fermi gas. We present results at finite filling factor $\textit{ν}$ ≃ 0.06 for the condensate fraction, an energy-staggering pairing gap, the spin susceptibility, and the heat capacity, and compare them to experimental data when available. In contrast to previous AFMC simulations, our model space consists of the complete first Brillouin zone of the lattice, and our calculations are performed in the canonical ensemble of fixed particle number. The canonical ensemble AFMC framework enables the calculation of a model-independent gap, providing direct information on pairing correlations without the need for numerical analytic continuation. We use finite-size scaling to estimate $T_c$ at the corresponding filling factor. In this work, we find that the energy-staggering pairing gap vanishes above $T_c$, showing no pseudogap effects, and that the spin susceptibility shows a substantially reduced signature of a spin gap compared to previously reported AFMC simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Waveform-dependent air fluorescence from neutral and ionic nitrogen molecules

Laser-induced air fluorescence in the ultraviolet regime is primarily attributed to transitions between the C and B states in excited neutral N 2 molecules and between the B and X states in N$^+_2$ ions. However, the mechanism underlying the former remains contentious, as direct population to the C state by light fields is forbidden by electron spin constraints. In this work, we investigate the mechanism of air fluorescence from excited neutral N 2 molecules by carrier-envelope phase–stabilized sub–4 femtosecond pulses. Our results show that fluorescence from N$^+_2$ ions reaches a maximum with cosine-like pulses, while fluorescence from excited neutral N 2 molecules peaks with sine-like pulses. In addition, by scanning the chirp of the driving pulse, we find that ionic fluorescence is maximized with chirp-free pulses, whereas neutral fluorescence favors negatively chirped pulses. These observations, supported by classical trajectory Monte Carlo simulations, support the mechanism of intersystem crossing from excited spin-singlet states, which are populated via recollision-induced strong-field excitation.

Science & Technology - Other Topics↗

Induced phase transformation in ionizable colloidal nanoparticles

Abstract Acid–base equilibria directly influence the functionality and behavior of particles in a system. Due to the ionizing effects of acid–base functional groups, particles will undergo charge exchange. The degree of ionization and their intermolecular and electrostatic interactions are controlled by varying the pH and salt concentration of the solution in a system. Although the pH can be tuned in experiments, it is hard to model this effect using simulations or theoretical approaches. This is due to the difficulty in treating charge regulation and capturing the cooperative effects in a colloidal suspension with Coulombic interaction. In this work, we analyze a suspension of ionizable colloidal particles via molecular dynamics (MD) simulations, along with Monte Carlo simulations for charge regulation (MC-CR) and derive a phase diagram of the system as a function of pH. It is observed that as pH increases, particles functionalized with acid groups change their arrangement from face-centered cubic (FCC) packing to a disordered state. We attribute these transitions to an increase in the degree of charge polydispersity arising from an increase in pH. Our work shows that charge regulation leads to amorphous solids in colloids when the mean nanoparticle charge is sufficiently high. Graphical abstract

36 MATERIALS SCIENCE↗

Space Launch System Implementation of Adaptive Augmenting Control

Given the complex structural dynamics, challenging ascent performance requirements, and rigorous flight certification constraints owing to its manned capability, the NASA Space Launch System (SLS) launch vehicle requires a proven thrust vector control algorithm design with highly optimized parameters to robustly demonstrate stable and high performance flight. On its development path to preliminary design review (PDR), the stability of the SLS flight control system has been challenged by significant vehicle flexibility, aerodynamics, and sloshing propellant dynamics. While the design has been able to meet all robust stability criteria, it has done so with little excess margin. Through significant development work, an adaptive augmenting control (AAC) algorithm previously presented by Orr and VanZwieten, has been shown to extend the envelope of failures and flight anomalies for which the SLS control system can accommodate while maintaining a direct link to flight control stability criteria (e.g. gain & phase margin). In this paper, the work performed to mature the AAC algorithm as a baseline component of the SLS flight control system is presented. The progress to date has brought the algorithm design to the PDR level of maturity. The algorithm has been extended to augment the SLS digital 3-axis autopilot, including existing load-relief elements, and necessary steps for integration with the production flight software prototype have been implemented. Several updates to the adaptive algorithm to increase its performance, decrease its sensitivity to expected external commands, and safeguard against limitations in the digital implementation are discussed with illustrating results. Monte Carlo simulations and selected stressing case results are shown to demonstrate the algorithm's ability to increase the robustness of the integrated SLS flight control system.

VanZwieten, Tannen S.↗

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modernizing the active neutron collar poison rod correction using simulations to enhance PWR fresh fuel verification

The declared inventory of 235 U in fresh, low enriched uranium, nuclear fuel assemblies is routinely verified by nuclear safeguards inspectorates using the UNCL (Uranium Neutron Collar – Light Water Reactor Fuel) nondestructive assay instrument. The trend in modern fuels is towards higher initial enrichment, which in turn requires a larger number of burnable poison pins with higher Gadolinia concentration to hold down the initial reactivity. UNCL assay error for modern fuel assemblies typically exceeds 10%. The traditional algorithm used to correct the response for burnable poison content needs revision to achieve accuracy comparable to the $\sigma_\text{R}$ of 4.5% expected for poison-free assemblies. Here we review the problem and use a large set of 287 Monte Carlo simulations based on an experimentally benchmarked model corresponding to the 16×16 PWR array used in Brazilian Angra type II and III fuel. Simulated relative responses are used to evaluate the functional form of the poison correction and determine parameters optimal for 16×16 PWR assemblies. This update is conducted using between 4–24 Gadolinia pins ranging 2–11 wt% Gd 2 O 3 with assemblies having mean enrichments ranging 2.5%–5%. Experimental validation of the updated coefficients using seven measurements shows bias is reduced by an order of magnitude and $\sigma_\text{R}$R reduced from 10% to less than 2%. Two updated sets of coefficients are provided, the most accurate of which is directly useable in the existing analysis code (INCC) in use by inspectorate at fuel fabrication facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

Bidirectional reflectance of oceanic waters: A comparison of modeled and measured upward radiance fields

The bidirectional reflectance of oceanic waters is conveniently described in a normalized way by forming the ratio of the upwelling irradiance E(sub u) to any upwelling radiance L(sub u)(theta prime, phi). This ratio, Q (theta prime, theta(sub 0), (phi(sub 0) - phi), where theta prime, phi are the nadir and azimuth angles for the upward radiance and theta(sub 0), phi(sub 0) are the zenith and azimuth angles of the Sun, has been determined from measurements at sea and computed via Monte Carlo simulations using the inherent optical properties measured in the field and appropriate boundary conditions (clear sky, no wind, varying Sun angle). Experimental ad computed Q values are in excellent agreement. This successful comparison confirms the importance of the bidirectional character of ocean reflectance, already pointed out from a purely numerical approach without field validation, and corroborates the extended range of the Q variations. The later point is of importance when interpreting the marine signals detected by an ocean color satellite-borne sensor. The validation is extended by considering the historical data for the radiance distributions in Lake Pend Oreille determined at various depths. The closure issue in ocean optics is examined by solving the direct problem of radiative transfer and through a model-data comparison in terms of radiance field.

Morel, Andre↗

Energy Calibration of the Scintillating Optical Fiber Calorimeter Chamber (SOFCAL)

The Scintillating Optical Fiber Calorimeter (SOFCAL) detector is designed to make direct measures of the primary cosmic ray spectrum from -200 GeV/amu - 20 TeV/amu. The primary particles are resolved into groups according to their charge (p, He, CNO, Medium Z, Heavy Z) using both active and passive components integrated into the detector. The principal part of SOFCAL is a thin ionization calorimeter that measures the electromagnetic cascades that result from these energetic particles interacting in the detector. The calorimeter is divided into two sections: a thin passive emulsion/x-ray film calorimeter, and a fiber calorimeter that uses crossing layers of small scintillating optical fibers to sample the energy deposition of the cascades. The energy determination is made by fitting the fiber data to transition curves generated by Monte Carlo simulations. The fiber data must first be calibrated using the electron counts from the emulsion plates in the calorimeter for a small number of events. The technique and results of this calibration will be presented together with samples of the data from a balloon flight.

Christl, M. C.↗

Monte Carlo simulation of the passage of $\gamma$-rays and $\alpha$-particles in CsI

In this work, theoretical and computational methods for simulating the creation of ionization tracks by fast ions in solids were applied to the passage of a-particles in CsI, an inorganic scintillator commonly used for radiation detection. The methods were implemented in a Monte Carlo program to simulate the interaction of $\alpha$-particles, with incident energies of up to 1 MeV, with CsI. The simulations followed the fate of individual electron-hole pairs and thus allowed for a detailed description of the microscopic structure of ionization tracks created by incident radiation. Simulations were also performed with $\gamma$-rays of the same energy to compare and contrast the ionization tracks obtained with both types of particle. Intrinsic properties such as the mean energy per electron-hole pair, Fano factor, maximum theoretical light yield, and spatial distributions of electron-hole pairs were computed for both $\alpha$-particles and $\gamma$-rays. $\alpha$-particles created cylindrical tracks that were initially aligned with the incident direction and with initial radii of a few nanometers, whereas $\gamma$-rays showed significant scattering, resulting in probability distributions with lower intensities and much greater radial extents.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

CTRANS: A Monte Carlo program for radiative transfer in plane parallel atmospheres with imbedded finite clouds: Development, testing and user's guide

The program called CTRANS is described which was designed to perform radiative transfer computations in an atmosphere with horizontal inhomogeneities (clouds). Since the atmosphere-ground system was to be richly detailed, the Monte Carlo method was employed. This means that results are obtained through direct modeling of the physical process of radiative transport. The effects of atmopheric or ground albedo pattern detail are essentially built up from their impact upon the transport of individual photons. The CTRANS program actually tracks the photons backwards through the atmosphere, initiating them at a receiver and following them backwards along their path to the Sun. The pattern of incident photons generated through backwards tracking automatically reflects the importance to the receiver of each region of the sky. Further, through backwards tracking, the impact of the finite field of view of the receiver and variations in its response over the field of view can be directly simulated.

Source record↗

Growth morphology with anisotropic surface kinetics

The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.

Xiao, Rong-Fu↗

Presentation, Analysis, and Simulation of Active Alignment Strategies for the James Webb Space Telescope

This paper presents three characteristics in the simulated active alignment strategy of the James Webb Space Telescope. The first includes the analysis and comparison of a baseline active alignment strategy with a damped least squares strategy. This baseline utilizes prior knowledge by means of direct human operator interaction to engage sets of telescope compensators to target specific aberration signatures. The baseline is compared to a damped least-squares strategy that utilizes simultaneous engagement of all telescope compensators without explicit human operator inter<actiori to achieve a least-squares telescope compensation. Second, we discuss how the active alignment of the JWST is encapsulated in a linear optical model developed at the Space Telescope Science Institute. This linear optical model provides a framework for an efficient and robust description of the optical control properties of the JWST and clearly articulates the necessity for having a multi-instrument multi-field wavefront sensing strategy to overcome control system non independence and the effects of non-common path errors in the main wavefront sensing camera. Finally. we present analytical results that explicitly map the telescope wavefront responses to the telescope control modes, and we present Monte-Carlo optical performance simulation results that demonstrate the efficacy of the damped least-squares active alignment and the prior-knowledge active alignment schemes.

Upton, Robert S.↗

Alloy Design Workbench-Surface Modeling Package Developed

NASA Glenn Research Center's Computational Materials Group has integrated a graphical user interface with in-house-developed surface modeling capabilities, with the goal of using computationally efficient atomistic simulations to aid the development of advanced aerospace materials, through the modeling of alloy surfaces, surface alloys, and segregation. The software is also ideal for modeling nanomaterials, since surface and interfacial effects can dominate material behavior and properties at this level. Through the combination of an accurate atomistic surface modeling methodology and an efficient computational engine, it is now possible to directly model these types of surface phenomenon and metallic nanostructures without a supercomputer. Fulfilling a High Operating Temperature Propulsion Components (HOTPC) project level-I milestone, a graphical user interface was created for a suite of quantum approximate atomistic materials modeling Fortran programs developed at Glenn. The resulting "Alloy Design Workbench-Surface Modeling Package" (ADW-SMP) is the combination of proven quantum approximate Bozzolo-Ferrante-Smith (BFS) algorithms (refs. 1 and 2) with a productivity-enhancing graphical front end. Written in the portable, platform independent Java programming language, the graphical user interface calls on extensively tested Fortran programs running in the background for the detailed computational tasks. Designed to run on desktop computers, the package has been deployed on PC, Mac, and SGI computer systems. The graphical user interface integrates two modes of computational materials exploration. One mode uses Monte Carlo simulations to determine lowest energy equilibrium configurations. The second approach is an interactive "what if" comparison of atomic configuration energies, designed to provide real-time insight into the underlying drivers of alloying processes.

Abel, Phillip B.↗