Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response simulation”

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 1,711 records · Page 95

Revisiting Multi-Material Composite Structures with Homogenized Composite Properties

Composite structures inherently develop residual stresses during their curing process. Driven predominately by mismatched thermal strains between differing materials or ply orientations, but also affected by curing process phenomena like polymer shrinkage, these residual stresses can lead to failure within composite structures. There are several methods varying in complexity that can be used to model the development of residual stresses, all of which are capable of capturing sufficient detail to understand the residual stress state at the ply level. However, explicitly modeling all plies of a layup in a composite structure can be prohibitively expensive based on the number of plies, structure size, and required element size. The computational cost can be reduced through the homogenization of the composite layup without losing much fidelity of the overall response of the structure. The homogenization process reduces the many plies of a laminate to a single lamina that reduces complexity and increases the mesh size where a single element can span multiple plies. This report focuses on verification and validation efforts for a homogenization process using a suite of finite element simulations rather than an analytic solution derived from classical laminate theory. Initial verification using representative element volumes indicated there was minimal error in the homogenization process; however, this compounded to a small, but acceptable error in strip and split ring experimental composite structures. The error does under predict the residual stress state in the strip and split ring and should be accounted for when simulating composite structures with homogenized properties.

36 MATERIALS SCIENCE↗

GC/MS Method Development for Separating Lunar Volatile Ice Simulant Headspace Gases

Various investigators propose the lunar surface contains widely distributed volatiles, especially water- like species, i.e. OH and H2O. Surface volatiles are theorized to exist as a hydrated regolith layer, concentrated in extremely cold polar permanently shadowed regions (PSR), and/or solar wind implantation reservoirs in lunar glasses. The proposed sources of lunar surface volatiles range from cometary impacts, solar wind, or a supply present during moon formation. Future Artemis missions aim to collect and return the samples containing volatiles collected near lunar polar craters or PSRs. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving returned sample integrity as much as possible. Pristine volatile-bearing samples are invaluable to the scientific community seeking to unravel the history of the solar system. Realistically, a sample will experience alteration during collection, transportation back to earth, and storage. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on high-fidelity volatile-containing regolith simulants, the foundation for the future of cold curation. This abstract outlines the separation, identification, and quantification of headspace gases over volatile ice feed stock material using gas chromatography/mass spectrometry (GC/MS). Preliminary objectives concentrated on sample handling, reproducibility, and understanding the elution characteristics for each analyte. Initial GC/MS method development experiments utilized diluted static headspace sample preparation. Diluted samples were used because sampling headspace gases directly from a vial containing liquid analyte resulted in overloading of the column and detector. Overloading is evident based on chromatogram peak shapes and instrument contamination, or carry over, between experiments. A mixture of three alcohols were used for a majority of the sample handling and reproducibility studies. Reproducibility was tested via multiple users, calibration curves, and check standards. Stock solutions of condensed lunar volatile analytes included methanol, ammonia in methanol, hydrogen sulfide in water, and an equal volume mixture of methanol, ethanol, and isopropanol. Current samples use room air as the headspace sample matrix, however future experiments will incorporate an inert purge gas, such as argon or nitrogen. Three mL of each analyte solution were capped in separate 20 mL crimp top GC vials. Dilutions were carried out by removing an aliquot of headspace gases with a calibrated 1 mL gastight syringe and immediately transferring to a 20 mL capped crimp top vial. The GC/MS is a Thermo Fisher Trace 1310/ISQ 7000 with a TriPlus RSH autosampler and split/splitless injector module. The experiments outlined in this abstract use the following hardware: a 2.5 mL gastight headspace syringe tool, 1 mm ID x 78.5 mm length ultra-inert straight injection liner, and a TG-BondQ 30 m × 0.32 mm × 10 μm column. Various parameters, such as hardware selection and the temperature, pressure, and split ratio set points, continue to evolve as the overall experiment is refined. Diluted headspace chromatograms were collected for the individual stock solutions. Retention times, peak shapes, and mass spectra were evaluated and added to the data processing method for each molecule of interest. Figure 1 shows the total ion chromatograms for the three major lunar volatile simulant stock solutions: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. Tailing peak shapes for ammonia (2.98 min rt) and water (4.06 min rt) indicate the molecules are not properly eluting from the selected column with the current separation method. Additionally, hydrogen sulfide and ammonia have overlapping peak windows, which could impact quantification. Ongoing experiments aim to address the peak shape and overlapping via the separation method and hardware selection. Sample preparation reproducibility experiments used stock solution containing equal volumes of a non- interactive mixture of methanol, ethanol, and isopropanol. Mass spectrum ion traces were used to identify and quantify all three alcohols. Peaks were automatically detected, identified, and integrated through the mass spectra detection and processing parameters. Calibration response curves and check standards were used to evaluate the validity of the sample preparation procedure. Figure 2 shows the methanol chromatogram peak area versus total headspace dilution volume transferred from the alcohol mixture vial. The calibration response curves and check standards validate sample preparation procedure. Continuing data analysis efforts are working towards correlating the peak area and instrument response factor to the headspace analyte concentration and condensed phase composition. Static headspace gas chromatography theory relies on Dalton’s law, Raoult’s law, Henry’s Law, and the Kolb and Ettre equation to associate peak area to the analyte composition in a non-ideal solution. Equation 1 is a simplified expression derived from the aforementioned theories. Future experiments involve liquid injections of the individual stock solutions, liquid and headspace analysis of various stock solution combinations, and the addition of regolith simulants to the mixtures. Temperature is another variable expected to affect reaction rates and will be explored.

Cecilia L. Amick↗

A Machine Learning Framework for Modeling Ensemble Properties of Atomically Disordered Materials

Atomic disorder can strongly influence material properties such as charge transport, optical response, and catalytic activity. However, efficiently modeling these disorder effects remains challenging for first-principles methods due to the cost of sampling large configurational spaces and computing complex physical quantities. Recent advances of machine learning techniques, particularly graph neural networks (GNNs), has enabled the efficient and accurate predictions of complex material properties, offering promising tools for studying disordered systems. In this work, we present a general machine-learning-assisted computational framework that integrates equivariant GNNs with Monte Carlo simulations to compute the thermodynamic and ensemble-averaged functional properties of disordered materials. Using the surface-termination-disordered MXene monolayer Ti 3 C 2 T 2–x as a representative system, we find that electrical conductivity exhibits an emergent peak near the order–disorder phase transition temperature due to the interplay between electron scattering and doping. In contrast, optical conductivity remains largely insensitive to local atomic disorder and reflects the global surface chemical composition. These results highlight the role of atomic disorder in affecting material properties and demonstrate the potential of our approach for statistically modeling disorder effects in a wide range of materials such as high-entropy alloys and spin liquids.

MXene↗

Negative capacitance regime in antiferroelectric PbZrO 3

The potential of antiferroelectrics to exhibit a negative capacitance regime has been largely overlooked as all the attention focused on their polar counterparts, ferroelectrics. We use nonequilibrium first-principles-based molecular dynamics to probe a negative capacitance regime in prototypical antiferroelectric PbZrO 3 . Simulations predict that this antiferroelectric can exhibit a negative capacitance/susceptibility regime in response of polarization to an internal electric field, which is a superposition of applied and residual depolarizing fields. Consequently, the regime emerges when the polarization surface charge in the polar phase of antiferroelectric is not fully screened, as is often the case in thin films and nanostructures. The negative capacitance regime occurs below the Curie temperature and disappears in the paraelectric phase. We find that the time the material spends in the negative capacitance regime is proportional to the time needed to complete antipolar–polar (or its reverse) transition and shortens as the frequency of the applied field increases. Furthermore, a negative susceptibility value exhibits strong dependence on the quality of surface charge screening with the largest in magnitude values occurring in the vicinity of the transition into a negative capacitance regime.

36 MATERIALS SCIENCE↗

Damage Resistance of Titanium Aluminide Evaluated

As part of the aviation safety goal to reduce the aircraft accident rate, NASA has undertaken studies to develop durable engine component materials. One of these materials, g-TiAl, has superior high-temperature material properties. Its low density provides improved specific strength and creep resistance in comparison to currently used titanium alloys. However, this intermetallic is inherently brittle, and long life durability is a potential problem. Of particular concern is the material s sensitivity to defects, which may form during the manufacturing process or in service. To determine the sensitivity of TiAl to defects, a team consisting of GE Aircraft Engines, Precision Cast Parts, and NASA was formed. The work at the NASA Glenn Research Center at Lewis Field has concentrated on the fatigue response to specimens containing defects. The overall objective of this work is to determine the influence of defects on the high cycle fatigue life of TiAl-simulated low-pressure turbine blades. Two types of defects have been introduced into the specimens: cracking from impact damage and casting porosity. For both types of defects, the cast-to-size fatigue specimens were fatigue tested at 650 C and 100 Hz until failure.

Lerch, Bradley A.↗

Development of a Tool to Recreate the Mars Science Laboratory Aerothermal Environment

The Mars Science Laboratory will enter the Martian atmosphere in 2012 with multiple char depth sensors and in-depth thermocouples in its heatshield. The aerothermal environment experienced by MSL may be computationally recreated using the data from the sensors and a material response program, such as the Fully Implicit Ablation and Thermal (FIAT) response program, through the matching of the char depth and thermocouple predictions of the material response program to the sensor data. A tool, CHanging Inputs from the Environment of FIAT (CHIEF), was developed to iteratively change different environmental conditions such that FIAT predictions match within certain criteria applied to an external data set. The computational environment is changed by iterating on the enthalpy, pressure, or heat transfer coefficient at certain times in the trajectory. CHIEF was initially compared against arc-jet test data from the development of the MSL heatshield and then against simulated sensor data derived from design trajectories for MSL. CHIEF was able to match char depth and in-depth thermocouple temperatures within the bounds placed upon it for these cases. Further refinement of CHIEF to compare multiple time points and assign convergence criteria may improve accuracy.

Beerman, A. F.↗

Thermal Analysis of Small Re-Entry Probe

The Small Probe Reentry Investigation for TPS Engineering (SPRITE) concept was developed at NASA Ames Research Center to facilitate arc-jet testing of a fully instrumented prototype probe at flight scale. Besides demonstrating the feasibility of testing a flight-scale model and the capability of an on-board data acquisition system, another objective for this project was to investigate the capability of simulation tools to predict thermal environments of the probe/test article and its interior. This paper focuses on finite-element thermal analyses of the SPRITE probe during the arcjet tests. Several iterations were performed during the early design phase to provide critical design parameters and guidelines for testing. The thermal effects of ablation and pyrolysis were incorporated into the final higher-fidelity modeling approach by coupling the finite-element analyses with a two-dimensional thermal protection materials response code. Model predictions show good agreement with thermocouple data obtained during the arcjet test.

Agrawal, Parul↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗

NREL/NASA Internal Short-Circuit Instigator in Lithium Ion Cells

NREL has developed a device to test one of the most challenging failure mechanisms of lithium-ion (Li-ion) batteries -- a battery internal short circuit. Many members of the technical community believe that this type of failure is caused by a latent flaw that results in a short circuit between electrodes during use. As electric car manufacturers turn to Li-ion batteries for energy storage, solving the short circuit problem becomes more important. To date, no reliable and practical method exists to create on-demand internal shorts in Li-ion cells that produce a response that is relevant to the ones produced by field failures. NREL and NASA have worked to establish an improved ISC cell-level test method that simulates an emergent internal short circuit, is capable of triggering the four types of cell internal shorts, and produces consistent and reproducible results. Internal short circuit device design is small, low-profile and implantable into Li-ion cells, preferably during assembly. The key component is an electrolyte-compatible phase change material (PCM). The ISC is triggered by heating the cell above PCM melting temperature (presently 40 degrees C -- 60 degrees C). In laboratory testing, the activated device can handle currents in excess of 300 A to simulate hard shorts (< 2 mohms). Phase change from non-conducting to conducting has been 100 percent successful during trigger tests.

25 ENERGY STORAGE↗

Predicting Local Material Thickness from Steady-State Ultrasonic Wavefield Measurements Using a Convolutional Neural Network

Acoustic steady-state excitation spatial spectroscopy (ASSESS) is a full-field, ultrasonic non-destructive evaluation (NDE) technique used to locate and characterize defects in plate-like structures. ASSESS generates a steady-state, single-tone ultrasonic excitation in a structure and a scanning laser Doppler Vibrometer (LDV) measures the resulting full-field surface velocity response. Traditional processing techniques for ASSESS data rely on wavenumber domain analysis. This paper presents the alternative use of a convolutional neural network (CNN), trained using simulated ASSESS data, to predict the local plate thickness at every pixel in the wavefield measurement directly. The defect detection accuracy of CNN-based thickness predictions are shown to improve for defects of greater size, and for defects with higher thickness reductions. The CNN demonstrates the ability to predict thickness accurately in regions where Lamb wave dispersion relations are complex or unknown, such as near the boundaries of a test specimen, so long as the CNN is trained on data that accounts for these regions. The CNN also shows generalizability to ASSESS experimental data, despite an entirely simulated training dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dynamic System Scaling Application to Accelerated Nuclear Fuel Testing

The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phase-field simulation of domain size effect on dielectric and piezoelectric responses in K 0.5 Na 0.5 NbO 3 epitaxial thin films with superdomain structures

Size effects of mesoscale ferroelectric domains on the macroscopic dielectric and piezoelectric responses in domain-engineered bulk piezocrystals have been extensively studied for more than two decades. However, less is known about the domain size effects in ferroelectric epitaxial thin films, especially for films consisting of low-symmetry ferroelectric phases and exhibits hierarchical superdomain structures. Herein, using phase-field simulations, we systemically evaluate the effective out-of-plane dielectric and piezoelectric coefficients, κ$^{*}_{33}$ and d$^{*}_{33}$, as a function of the domain periods for two types of superdomain structures in ferroelectric K 0.5 Na 0.5 NbO 3 epitaxial thin films. In one type of the superdomain structures, we find that more than 70% increase of κ$^{*}_{33}$ and nearly 20% increase of d$^{*}_{33}$ can be achieved by tuning the domain period by a few tens of nanometers. Dissimilar behaviors are found when the domain period varies along different lateral directions, suggesting anisotropic domain size effects in thin films. By analyzing the local dielectric and piezoelectric responses from each domain variants and domain walls, we reveal that the domain size effect is governed by the variation of out-of-plane polarization inside the domains. Moreover, we also demonstrate enhanced domain size effects by modulating the misfit strains and temperature to approach the polymorphic phase boundaries, suggesting tunability of the size effect by doping and strain engineering. Furthermore, our results reveal multimodal domain size dependence of dielectric and piezoelectric responses in low-symmetry ferroelectric epitaxial thin films, implying that domain size engineering can be used to tune macroscopic properties of thin-film ferroelectrics, similar to their bulk counterparts.

36 MATERIALS SCIENCE↗

Analysis of tow architecture variability in biaxially braided composite tubes

Spatial deviations in composite braid unit cell dimensions limit mechanical performance and are challenging to model. Full-field measurements of the irregular tow placement in biaxially braided composite tubes were obtained by digital image correlation. Analysis of the reconstructed tow paths showed long-range defects along the hoop direction with short-range defects dominating the biased tow direction. A route toward numerical simulation of the composite response is established by developing a 3D model of the composite using the measured tow trajectories. Finally, a comparison to previously-published triaxial braid results indicates axial tows work to constrain unit cells dimensions along the hoop direction.

3-Dimensional reinforcement↗

The origin of the X-ray emission from the high-velocity cloud MS30.7–81.4–118

A soft X-ray enhancement has recently been reported toward the high-velocity cloud MS30.7–81.4–118 (MS30.7), a constituent of the Magellanic Stream. In order to investigate the origin of this enhancement, we have analyzed two overlapping XMM-Newton observations of this cloud. We find that the X-ray enhancement is ∼6' or ∼100 pc across, and is concentrated to the north and west of the densest part of the cloud. We modeled the X-ray enhancement with a variety of spectral models. A single-temperature equilibrium plasma model yields a temperature of (3.69{sub −0.44}{sup +0.47})×10{sup 6} K and a 0.4-2.0 keV luminosity of 7.9 × 10{sup 33} erg s{sup –1}. However, this model underpredicts the on-enhancement emission around 1 keV, which may indicate the additional presence of hotter plasma (T ≳ 10{sup 7} K), or that recombination emission is important. We examined several different physical models for the origin of the X-ray enhancement. We find that turbulent mixing of cold cloud material with hot ambient material, compression or shock heating of a hot ambient medium, and charge exchange reactions between cloud atoms and ions in a hot ambient medium all lead to emission that is too faint. In addition, shock heating in a cool or warm medium leads to emission that is too soft (for reasonable cloud speeds). We find that magnetic reconnection could plausibly power the observed X-ray emission, but resistive magnetohydrodynamical simulations are needed to test this hypothesis. If magnetic reconnection is responsible for the X-ray enhancement, the observed spectral properties could potentially constrain the magnetic field in the vicinity of the Magellanic Stream.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of Soret effect on hydrogen and helium retention in PFC tungsten under ELM-like conditions

Abstract In our previous work, we have demonstrated using nonequilibrium molecular-dynamics simulations that the fluxes of helium and self-interstitial atoms in the presence of a thermal gradient in tungsten are directed opposite to the heat flux, indicating that species transport is governed by a Soret effect, namely, thermal-gradient-driven diffusion, characterized by a negative heat of transport that drives species transport uphill, i.e. from the cooler to the hot regions of the tungsten sample. In this work, the findings of our thermal and species transport analysis have been implemented in our cluster-dynamics code, Xolotl, which has been used to compute temperature and species profiles over spatiotemporal scales representative of plasma-facing component (PFC) tungsten under typical reactor operating conditions, including extreme heat loads at the plasma-facing surface characteristic of plasma instabilities that induce edge localized modes (ELMs). We demonstrate that the steady-state species profiles, when properly accounting for the Soret effect, vary significantly from those where temperature-gradient-driven transport is not accounted for and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten. Although our cluster-dynamics simulations do not yet include self-clustering of helium or hydrogen blister formation, our simulation results show that the Soret effect substantially reduces helium and hydrogenic species retention inside PFC tungsten.

36 MATERIALS SCIENCE↗

Stripe antiferromagnetism and disorder in the Mott insulator NaFe 1₋ x Cu x As ( x ≲0.5)

Neutron scattering measurements have demonstrated that the heavily Cu doped NaFe 1₋ x Cu x As compound behaves like a Mott insulator exhibiting both real-space Fe-Cu stripes and antiferromagnetism below a Néel temperature for x ≲ 0.5. In this work, we have investigated the evolution of structural and magnetic ordering using 23 Na and 75 As NMR for single crystals ( x = 0.39 and 0.48), confirming antiferromagnetism in the form of magnetic stripes. We show that end-chain defects in these stripes are the principal source of magnetic disorder and are responsible for cluster spin-glass transitions in both compounds, in the latter case coexistent with antiferromagnetism. Aided by our numerical simulation of the 75 As spectra, we show that a staggered magnetization at the Fe sites is induced by nonmagnetic Cu dopants.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simulator for an Accelerator-Driven Subcritical Fissile Solution System

LANL has developed a process to generate a progressive family of system models for a fissile solution system. This family includes a dynamic system simulation comprised of coupled nonlinear differential equations describing the time evolution of the system. Neutron kinetics, radiolytic gas generation and transport, and core thermal hydraulics are included in the DSS. Extensions to explicit operation of cooling loops and radiolytic gas handling are embedded in these systems as is a stability model. The DSS may then be converted to an implementation in Visual Studio to provide a design team the ability to rapidly estimate system performance impacts from a variety of design decisions. This provides a method to assist in optimization of the system design. Once design has been generated in some detail the C++ version of the system model may then be implemented in a LabVIEW user interface to evaluate operator controls and instrumentation and operator recognition and response to off-normal events. Taken as a set of system models the DSS, Visual Studio, and LabVIEW progression provides a comprehensive set of design support tools.

42 ENGINEERING↗