Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface structure prediction”

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 343 records · Page 19

Variational Coupled Loads Analysis using the Hybrid Parametric Variation Method

Time-domain coupled loads analysis (CLA)is used to determine the response of a launch vehicle and payload system to transient forces, such as liftoff, engine ignitions and shutdowns, jettison events, and atmospheric flight loads, such as buffet. CLA, using Hurty/Craig-Bampton (HCB)component models, is the accepted method for the establishment of design-level loads for launch systems. However, uncertainty in the component models flows into uncertainty in predicted system results. Uncertainty in the structural responses during launch is a significant concern because small variations in launch vehicle and payload mode shapes and their interactions can result in significant variations in system loads. Uncertainty quantification (UQ)is used to determine statistical bounds on prediction accuracy based on model uncertainty. In this paper uncertainty is treated at the HCB component-model level. In an effort to account for model uncertainties and statistically bound their effect on CLA predictions, this work combines CLA with UQ in a process termed variational coupled loads analysis (VCLA). The modeling of uncertainty using a parametric approach, in which input parameters are represented by random variables, is common, but its major drawback is the resulting uncertainty is limited to the form of the nominal model. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. Model-form uncertainty can be represented using a nonparametric approach based on random matrix theory (RMT). In this work, UQ is performed using the hybrid parametric variation (HPV)method, which combines parametric with nonparametric uncertainty at the HCB component model level. The HPV method requires the selection of dispersion values for the HCB fixed-interface (FI)eigenvalues, and the HCB mass and stiffness matrices. The dispersions are based upon component test-analysis modal correlation results. During VCLA, random component models are assembled into an ensemble of random systems using a Monte Carlo (MC)approach. CLA is applied to each of the ensemble members to produce an ensemble of system-level responses for statistical analysis. The proposed methodology is demonstrated through its application to a buffet loads analysis of NASA’s Space Launch System (SLS)during the transonic regime fifty seconds after liftoff. Core stage (CS)section shears and moments are recovered, and statistics are computed.

Uncertainty Quantification↗

CIGS Technology Advancement via Fundamental Modeling of Defect/Impurity Interactions (Final Technical Report)

The primary goals of the proposed work were to provide modeling tools (and the associated insight which comes along with model development) for design and optimization of CuIn x Ga 1-x Se 2 (CIGS) and CdSeTe (CST) solar cell manufacturing processes and to establish the foundation for comprehensive end-to-end predictive modeling tools to enable optimization of thin film photovoltaic technology for performance, cost, yield, and reliability. The initial focus of efforts within this project was to develop coupled process/optical/device models for CIGS PV technology and to work with Siva Power to apply that TCAD (technology computer-aided design) system to improve the efficiency and reduce manufacturing costs for CIGS solar cells. Our approach to that end was to generate an extensive database of DFT calculations and to use those calculations via statistical thermodynamics methods and Monte Carlo simulation to develop and characterize models for the behavior of native defects as well as intentional and unintentional impurities, including the redistribution of the primary components of CIGS films. Increased effort went toward coupling those models for defect behavior and composition evolution to the performance of multicrystalline CIGS solar cells via prediction of doping level and recombination lifetime as function of manufacturing process. In the second budget period, the project pivoted to developing a similar system for the CdSeTe system, focused especially on understanding the role of Se/Te alloy concentration. Execution of the project resulted in the successful development of TCAD systems for both CIGS and CdSeTe thin film PV within the Synopsys Sentaurus framework by utilizing the Alagator interface. In the first budget period of the project, we developed quantitative models for the major components of CIGS PV and implemented them within a framework that couples process, optical, and device simulation. From the insights we have gained, we identified novel opportunities for enhancing CIGS solar cell performance and have laid the groundwork to further optimize the layer structure, composition profile, and thermal cycles for substantially improved efficiency and lower manufacturing costs. For the CIGS system, process changes to achieve greater than 1% absolute enhancement in efficiency were identified, but testing of those approaches was stymied by lack of a domestic CIGS manufacturing partner after the closure of Siva Power as well as Miasole. For CdSeTe, a fully capable TCAD system only became ready to apply near the end of the project period, so substantial opportunities remain to apply those models to enhance the leading thin film PV technology.

14 SOLAR ENERGY↗

Sensitivity Analysis of Multidisciplinary Rotorcraft Simulations

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. An unstructured sensitivity-enabled Navier-Stokes solver, FUN3D, and a nonlinear flexible multibody dynamics solver, DYMORE, are coupled to predict the aerodynamic loads and structural responses of helicopter rotor blades. A discretely-consistent adjoint-based sensitivity analysis available in FUN3D provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute DYMORE structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Numerical results verify accuracy of the FUN3D/DYMORE system by conducting simulations for a benchmark rotorcraft test model and comparing solutions with established analyses and experimental data. Complex-variable implementation of sensitivity analysis of DYMORE and the coupled FUN3D/DYMORE system is verified by comparing with real-valued analysis and sensitivities. Correctness of adjoint formulations for FUN3D/DYMORE interfaces is verified by comparing adjoint-based and complex-variable sensitivities. Finally, sensitivities of the lift and drag functions obtained by complex-variable FUN3D/DYMORE simulations are compared with sensitivities computed by the multidisciplinary sensitivity analysis, which couples adjoint-based flow and grid sensitivities of FUN3D and FUN3D/DYMORE interfaces with complex-variable sensitivities of DYMORE structural responses.

Wang, Li↗

Microstructural Development of Directionally Solidified Hg(1-x)Zn(x)Se Alloys

Hg(1-x)Zn(x)Se alloys have been studied as an alternative to HgCdTe for the detection of electromagnetic radiation, because the shorter ZnSe and HgSe bonds were predicted to improve lattice stability. Previous studies showed that optical properties were stable with time and that the microhardness was higher than for HgTe based alloys. However, for this material to be commercially viable, compositional variations must be reduced, and the distribution of defects induced by the growth process must be understood. The basic structural properties of bulk materials were characterized, and dislocation etch pit densities were measured to determine whether the addition of Zn to the HgSe lattice reduced the potential for dislocation formation. Hg(0.9)Zn(0.1)Se alloys were directionally solidified using a modified Bridgman-Stockbarger method and in an applied magnetic field. Radial compositional variations were greatly reduced when solidification occurred in an applied magnetic field. The combination of a convex liquid-solid interface shape and a method for reduced wetting produced boules that were a single crystal after reaching full diameter. Observed surface features indicated ampoule wetting could be eliminated using a graphite getter. Microstructural characteristics were greatly improved over HgCdTe alloys, In six boules, a total of only one twin was observed, and no inclusions were found. A method for polishing and producing dislocation etch pits was developed for these alloys, revealing dislocation etch pit densities one to two orders of magnitude less than HgTe based alloys. A kink in the thermal profile during processing of one boule generated more dislocations than did lattice mismatch due to compositional variations. The results of this investigation indicate this alloy has improved microstructural properties and resistance to dislocation formation compared with similar II-VI alloys, and should be further investigated for the detection of electromagnetic radiation.

Cobb, S. D.↗

Experimental and Numerical Analysis of Skin-Stiffener Separation Using a Seven-Point Bend Configuration

Skin-stiffener separation in stiffened composite panels consists of a complex interaction between multiple scales of progressive damage and failure mechanisms. This work used a superposed cohesive element method of modeling resistance curve effects that represents the structural interface of a unidirectional tape skin and a fabric stiffener. Finite element models using cohesive elements with input properties obtained from tape-to-fabric interface characterization tests were compared to experimental results of a stable skin-stiffener separation characterization test using a seven-point bend fixture. This fixture deformed the stiffened panel specimens into a buckled configuration which induced mixed-mode interlaminar stress states at the skin-stiffener interface. The advantage of this fixture was the potential for more stable damage initiation and delamination growth compared to a stringer-stiffened panel under axial compression. Use of the superposed cohesive elements showed promise, but the characterization of interface material properties as input to the cohesive elements remains a crucial component to be investigated to ensure accurate prediction of skin-stiffener separation.

Kosztowny, Cyrus J. R.↗

High‐Resolution National‐Scale Water Modeling Is Enhanced by Multiscale Differentiable Physics‐Informed Machine Learning

Abstract The National Water Model (NWM) is a key tool for flood forecasting, planning, and water management. Key challenges facing the NWM include calibration and parameter regionalization when confronted with big data. We present two novel versions of high‐resolution (∼37 km 2 ) differentiable models (a type of hybrid model): one with implicit, unit‐hydrograph‐style routing and another with explicit Muskingum‐Cunge routing in the river network. The former predicts streamflow at basin outlets whereas the latter presents a discretized product that seamlessly covers rivers in the conterminous United States (CONUS). Both versions use neural networks to provide a multiscale parameterization and process‐based equations to provide a structural backbone, which were trained simultaneously (“end‐to‐end”) on 2,807 basins across the CONUS and evaluated on 4,997 basins. Both versions show great potential to elevate future NWM performance for extensively calibrated as well as ungauged sites: the median daily Nash‐Sutcliffe efficiency of all 4,997 basins is improved to around 0.68 from 0.48 of NWM3.0. As they resolve spatial heterogeneity, both versions greatly improved simulations in the western CONUS and also in the Prairie Pothole Region, a long‐standing modeling challenge. The Muskingum‐Cunge version further improved performance for basins >10,000 km 2 . Overall, our results show how neural‐network‐based parameterizations can improve NWM performance for providing operational flood predictions while maintaining interpretability and multivariate outputs. The modeling system supports the Basic Model Interface (BMI), which allows seamless integration with the next‐generation NWM. We also provide a CONUS‐scale hydrologic data set for further evaluation and use.

Song, Yalan [Civil and Environmental Engineering T↗

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution, electric propulsion system, and heat rejection system. Specific mass, or α (kg/kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass predictions for each individual component. To inform technology maturation planning activities, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively-cooled heat pipe reactor concepts. Each concept requires specific geometries, fluids, and power conversion interface conditions (temperature, pressure, flow rate) to meet desired performance and mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support concept down-selection along with related technology development planning. The reactor and shield α are a function of several PCS and heat rejection system design choices, and reactor scaling with these parameters must be considered to enable an informed decision on an optimal reactor geometry and working fluid combination.

Nuclear Electric Propulsion↗

An Investigation of the Grain Size and Adhesion Strength of High-Speed Impact Ice

Ice accumulation on aircraft from supercooled water droplets (impact icing) can increase drag and decrease lift, leading to a decrease in fuel efficiency and a potential for unsafe situations. Anti-icing materials and coatings are of high interest to prevent aircraft icing. Because no anti-icing material suitable for in-flight situations has been realized, icing models aim to predict the accumulation and shape of the ice to understand its effects. In both cases, understanding the physical mechanism behind ice adhesion and shedding is necessary. A recent literature review has shown that the apparent ice adhesion strength to any specific material varies over several orders of magnitude across studies and test methods. One of the reasons for this high variance is likely to be the varying microstructure of the ice used for different studies. Although the microstructure of impact ice has been studied in the past, those conditions were not relevant to that of in-flight aircraft icing. The microstructure of the ice formed will determine both its bulk and surface properties, yet it has been overlooked in the ice adhesion community. Here, we use optical microscopy to study the non-uniform ice crystal grain size in impact ice accreted in NASA Glenn Research Center’s Icing Research Tunnel at various icing conditions. Similar to other environmental conditions studied in the literature, we show an increase in grain size at larger distances from the growth interface; however, the data shows substantially larger grains than what was expected from the prior literature. Ice accrued at lower temperatures has a smaller grain size at all distances from the growth interface. No change in grain size was noticed when varying only the wind speed. A correlational analysis shows that smaller grain sizes at the growth interface corresponded to an increase in the apparent ice adhesion strength when the adhesive failure was interfacial, but not when the failure was cohesive. These results may explain the overall trend in the literature of increasing apparent adhesion strength with decreasing temperature. Incorporating microstructure trends for various environmental conditions in icing models may increase the certainty in predicting ice adhesion strength and shedding, and will provide a better basis for comparison between adhesive tests on ice produced in different facilities and at different storage times.

• Ice Fabric↗

Universal Workflow Language and Software Enable Geometric Learning and FAIR Scientific Protocol Reporting

Written language and conventional data structures for representing scientific procedures suffer from low process detail, often fail to accurately represent protocols, and lack universality. New strategies for the handling of experimental data are needed to provide viable process information for both humans and machines. In this work, we present the universal workflow language (UWL) and interface (UWLi). UWL is a findable, accessible, interoperable, and reusable (FAIR)-compatible, graph-based data architecture that can capture arbitrary scientific procedures through workflow representation, and UWLi is an accompanying software package for building, manipulating, and interpreting UWL entries. The UWL format was found to be highly effective in identifying deficiencies in the reported process details of high-impact, peer-reviewed scientific journals, and in simulated scenarios, the graph format was shown to be more effective than conventional methods in predictively modeling the outcome of diverse scientific protocols. Implementation of UWL could enable more accurate scientific communication and more impactful process datasets.

14 SOLAR ENERGY↗

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗

Unraveling nonperturbative QCD with transverse momentum hadronic structures

The main purpose of this thesis is to investigate the theoretical foundations of the factorization theorems that involve transverse quark and gluon momentum distributions, and to better understand the interface between foundational questions and phenomenological applications. Many of these applications involve using transverse parton momentum dependent (TMD) correlation functions as tools for probing the complex nonperturbative structures of the hadrons. We first introduce the concepts of collinear and TMD factorization in a toy-model theory and, by leveraging the comparative simplicity of such a model, we test the range of validity and limitations of the factorization approach. At the same time, we introduce the conventional approaches that are adopted for the phenomenological extractions of the TMD distributions, and catalog the advantages and disadvantages. Following the general theoretical framework of factorization, we then develop a different approach to TMD phenomenology that is able to systematically integrate the perturbative and non perturbative information in parametrizations of TMD distributions, while being completely consistent with their operator definitions and the predictions of QCD, including their nontrivial evolution. The name Hadron Structure Oriented approach, or HSO, is coined in order to emphasize the central role of the nonperturbative nature of QCD and its relation to the partonic structure of hadrons. Finally, we provide a first phenomenological implementation of the HSO approach by extracting TMD distributions relying on low-to-moderate energy Drell-Yan data. The quality and robustness of the results are then compared against higher energy processes like high-Q2 Z0 boson production. Our analysis demonstrates the feasibility of this novel approach and points the way toward the improvements it can achieve in future applications.

Rainaldi, Tommaso [Old Dominion Univ., Norfolk, VA↗

Molecular Dynamics of a Water-Lipid Bilayer Interface

We present results of molecular dynamics simulations of a glycerol 1-monooleate bilayer in water. The total length of analyzed trajectories is 5ns. The calculated width of the bilayer agrees well with the experimentally measured value. The interior of the membrane is in a highly disordered fluid state. Atomic density profile, orientational and conformational distribution functions, and order parameters indicate that disorder increases toward the center of the bilayer. Analysis of out-of-plane thermal fluctuations of the bilayer surfaces occurring at the time scale of the present calculations reveals that the distribution of modes agrees with predictions of the capillary wave model. Fluctuations of both bilayer surfaces are uncorrelated, yielding Gaussian distribution of instantaneous widths of the membrane. Fluctuations of the width produce transient thinning defects in the bilayer which occasionally span almost half of the membrane. The leading mechanism of these fluctuations is the orientational and conformational motion of head groups rather than vertical motion of the whole molecules. Water considerably penetrates the head group region of the bilayer but not its hydrocarbon core. The total net excess dipole moment of the interfacial water points toward the aqueous phase, but the water polarization profile is non-monotonic. Both water and head groups significantly contribute to the surface potential across the interface. The calculated sign of the surface potential is in agreement with that from experimental measurements, but the value is markedly overestimated. The structural and electrical properties of the water-bilayer system are discussed in relation to membrane functions, in particular transport of ions and nonelectrolytes across membranes.

Wilson, Michael A.↗

Experimental and theoretical investigation of the crystalline surface, film, and interface properties of antiperovskite Mn 3 GaN grown by molecular beam epitaxy on MgO(001)

Here, we present a study of the epitaxial growth, characterization, and theoretical modeling of thin film antiperovskite Mn 3 GaN, an antiferromagnetic material with kagome structure which is grown on MgO (001) substrates using N-plasma-assisted molecular beam epitaxy. Reflection high energy electron diffraction is used to assess the in-plane evolution of the film structure during growth, and the surface is investigated in-situ using scanning tunneling microscopy and Auger electron spectroscopy. These results are combined with precision measurements done ex-situ determining the film lattice constants using a combination of x-ray diffraction with reciprocal space mapping and scanning transmission electron microscopy. Overall, a uniform, homogeneous film with an atomically smooth vacuum surface and atomically sharp substrate interface is found having very small in-plane tensile strain and mild out-of-plane compressive strain. First-principles theoretical calculations are applied in order to ascertain the lowest energy models for both the Mn 3 GaN surface and the Mn 3 GaN/MgO film/substrate interface. Models including MnGa versus MnN surface layers and MnGa versus MnN interfacial layers are considered as functions of both the Mn and Ga chemical potentials. The predictions are discussed in comparison to the experimental results. The overall findings suggest that Mn 3 GaN on MgO(001) is a viable epitaxial film which can be further explored in connection with antiferromagnetic spintronics.

Density functional theory↗

CELFE: Coupled Eulerian-Lagrangian Finite Element program for high velocity impact. Part 1: Theory and formulation

A 3-D finite element program capable of simulating the dynamic behavior in the vicinity of the impact point, together with predicting the dynamic response in the remaining part of the structural component subjected to high velocity impact is discussed. The finite algorithm is formulated in a general moving coordinate system. In the vicinity of the impact point contained by a moving failure front, the relative velocity of the coordinate system will approach the material particle velocity. The dynamic behavior inside the region is described by Eulerian formulation based on a hydroelasto-viscoplastic model. The failure front which can be regarded as the boundary of the impact zone is described by a transition layer. The layer changes the representation from the Eulerian mode to the Lagrangian mode outside the failure front by varying the relative velocity of the coordinate system to zero. The dynamic response in the remaining part of the structure described by the Lagrangian formulation is treated using advanced structural analysis. An interfacing algorithm for coupling CELFE with NASTRAN is constructed to provide computational capabilities for large structures.

Lee, C. H.↗

Large format imaging arrays for the Atacama Cosmology Telescope

We describe progress in the fabrication, characterization, and production of detector arrays for the Atacama Cosmology Telescope (ACT). The completed ACT instrument is specified to image simultaneously at 145, 225, and 265 GHz using three 32x32 filled arrays of superconducting transition edge sensors (TES) read out with time-division-multiplexed SQUID amplifiers. We present details of the pixel design and testing including the optimization of the electrical parameters for multiplexed readout. Using geometric noise suppression and careful tuning of operation temperature and device bias resistance, the excess noise in the TES devices is balanced with detector speed for interfacing with the ACT optics. The design also accounts for practical tolerances such as transition temperature gradients and scatter that occur in the production of multiple wafers to populate fully the kilopixel cameras. We have developed an implanted absorber layer compatible with our silicon-on-insulator process that allows for tunable optical resistance with requisite on-wafer uniformity and wafer-to-wafer reproducibility. Arrays of 32 elements have been tested in the laboratory environment including electrical, optical, and multiplexed performance. Given this pixel design, optical tests and modeling are used to predict the performance of the filled array under anticipated viewing conditions. Integration of the filled array of pixels with a tuned backshort and dielectric plate in front of the array maximize absorption and the focal plane and suppress reflections. A mechanical design for the build of the full structure is completed and we report on progress toward the construction of a prototype array for first light on the ACT.

Chervenak, J. A.↗

Linearized Aeroelastic Solver Applied to the Flutter Prediction of Real Configurations

A fast-running unsteady aerodynamics code, LINFLUX, was previously developed for predicting turbomachinery flutter. This linearized code, based on a frequency domain method, models the effects of steady blade loading through a nonlinear steady flow field. The LINFLUX code, which is 6 to 7 times faster than the corresponding nonlinear time domain code, is suitable for use in the initial design phase. Earlier, this code was verified through application to a research fan, and it was shown that the predictions of work per cycle and flutter compared well with those from a nonlinear time-marching aeroelastic code, TURBO-AE. Now, the LINFLUX code has been applied to real configurations: fans developed under the Energy Efficient Engine (E-cubed) Program and the Quiet Aircraft Technology (QAT) project. The LINFLUX code starts with a steady nonlinear aerodynamic flow field and solves the unsteady linearized Euler equations to calculate the unsteady aerodynamic forces on the turbomachinery blades. First, a steady aerodynamic solution is computed for given operating conditions using the nonlinear unsteady aerodynamic code TURBO-AE. A blade vibration analysis is done to determine the frequencies and mode shapes of the vibrating blades, and an interface code is used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor is used to interpolate the mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh. Then, LINFLUX is used to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. Finally, a post-processor uses the unsteady pressures to calculate the generalized aerodynamic forces, eigenvalues, an esponse amplitudes. The eigenvalues determine the flutter frequency and damping. Results of flutter calculations from the LINFLUX code are presented for (1) the E-cubed fan developed under the E-cubed program and (2) the Quiet High Speed Fan (QHSF) developed under the Quiet Aircraft Technology project. The results are compared with those obtained from the TURBO-AE code. A graph of the work done per vibration cycle for the first vibration mode of the E-cubed fan is shown. It can be seen that the LINFLUX results show a very good comparison with TURBO-AE results over the entire range of interblade phase angle. The work done per vibration cycle for the first vibration mode of the QHSF fan is shown. Once again, the LINFLUX results compare very well with the results from the TURBOAE code.

Reddy, Tondapu S.↗

Atomistic Modeling of the Interphase Structure between Alumina and Al12Fe2Cr Quasicrystal Approximant

Molecular dynamics simulations are performed to study the stability of an fcc metallic interlayer phase between an iron aluminide quasicrystal approximant (QCA) phase (84.2 Al, 5.7 Fe, 5.4 Ni, and 4.7 at.% Cr) and aluminum oxides (α-Al 2 O 3 and γ-Al 2 O 3 ). The presence of the interlayer phase was experimentally observed in some regions of the inner diameter surface of a TPBAR cladding tube. The simulations employed existing many-body (MEAM and EAM) interatomic potentials to describe atomic interactions in the interlayer and QCA phases, and Buckingham+Coulomb pair interactions for the alumina and between the alumina and the QCA (or fcc). Interface systems of QCA/alumina, fec/alumina, and QCA/fcc are modeled and the fracture energy of the interfaces is calculated. Based on the fracture energies, for the γ-Al 2 O 3 , both MEAM and EAM suggest that the fcc interlayer phase would not form. For the α-Al 2 O 3 , the MEAM predicts that the fcc phase would form, however, the EAM predicts that it forms only if the α-Al 2 O 3 is Al-terminated at the interface. For the O-terminated α-Al 2 O 3 , the EAM predicts the fcc phase would not form, consistent with the results for γ-Al 2 O 3 .

36 MATERIALS SCIENCE↗