Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix power”

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,279 records · Page 71

Hypersonic Second-Mode Instability Response to Shaped Roughness

An experimental campaign was conducted on a 7-degree half-angle cone in the NASA Langley Research Center 20-Inch Mach 6 Wind Tunnel to examine the influence of arrays of regularly spaced roughness elements on instability growth and transition. The primary element shape was a pair of elliptical planform ramps that were inclined at equal and opposite angles with respect to the local streamwise direction. The element shapes were designed to induce transient growth disturbances that would lead to sustained azimuthal modulation of the boundary layer flow while limiting the nearfield disturbances to avoid an immediate, i.e., effective tripping of the boundary layer. The bulk of the run matrix consisted of testing different element height sat free stream unit Reynolds numbers ranging from 9.8 to 13.1 million per meter. Other element shapes previously designed for tripping hypersonic boundary layers were also implemented. The model was instrumented with surface mounted Kulite ® and PCB ® pressure transducers and thermocouples. Spectra from the PCBs ® indicated clear suppression of the second-mode instability; however, neither the PCB ® spectra nor the heat transfer data presented strong evidence for delayed turbulent flow. Complementary stability computations likewise demonstrated second-mode reduction, particularly just downstream of the roughness, but also revealed a rise in first mode (streak-instability) amplitudes from the baseline that was likely responsible for the earlier transition observed for taller roughness cases.

Hypersonics↗

Application of modern control theory to the design of optimum aircraft controllers

The synthesis procedure presented is based on the solution of the output regulator problem of linear optimal control theory for time-invariant systems. By this technique, solution of the matrix Riccati equation leads to a constant linear feedback control law for an output regulator which will maintain a plant in a particular equilibrium condition in the presence of impulse disturbances. Two simple algorithms are presented that can be used in an automatic synthesis procedure for the design of maneuverable output regulators requiring only selected state variables for feedback. The first algorithm is for the construction of optimal feedforward control laws that can be superimposed upon a Kalman output regulator and that will drive the output of a plant to a desired constant value on command. The second algorithm is for the construction of optimal Luenberger observers that can be used to obtain feedback control laws for the output regulator requiring measurement of only part of the state vector. This algorithm constructs observers which have minimum response time under the constraint that the magnitude of the gains in the observer filter be less than some arbitrary limit.

Power, L. J.↗

Evaluation of Joint Cyber/Safety Risk in Nuclear Power Systems

This report presents an analysis of the Emergency Core Cooling System (ECCS) for a generic Boiling Water Reactor (BWR)-4 NPP. The Electric Power Research Institute (EPRI) developed Hazards and Consequences Analysis for Digital Systems (HAZCADS) process is applied to the ECCS and its subsystems to identify unsafe control actions (UCAs) which act as possible cyber events of concern. The analysis is performed for two design basis events: Small-break Loss of Coolant Accident (SLOCA) and general transients (TRANS), such as unintended reactor trip. In previous work, HAZCADS UCAs were combined with other cyber-attack analysis to develop a risk-informed approach; however, this was for a single system. This report explores advanced systems engineering modeling approaches to model the interactions between digital assets across multiple systems which may be targeted by cyber adversaries. The complex and interdependent design of digital systems has the potential to introduce emergent cyber properties that are generally not covered by hazard analyses nor formal nuclear Probabilistic Risk Assessment (PRA). The R&D and supporting analysis presented here explores approaches to predict and manage how interdependent system properties effect risk. To show the potential impact of a successful cyber-attack to formal PRA event tree probabilities, HAZCADS analysis was also used. HAZCADS was also used to model the automatic depressurization system (ADS) automatic actuation. This analysis extended to an integrated system analysis for common-cause failure (CCF). In this aspect, the HAZCADS analysis continued by analyzing plant design details for system connectivity in support of critical plant functions. A dependency matrix was developed to depict the integrated functionality of the interconnected systems. Areas of potential CCF are indicated. Future work could include adversary attack development to show how CCF could be caused, resulting in PRA events. Across the multiple systems that comprise the ECCS, the analysis shows that the change in such probabilities was very different between systems. This indicates that some systems have a larger potential risk impact from successful cyber-attack or digital failure, which indicates a need for these systems to have a higher priority for design and defensive measures. Furthermore, we were able to establish that a risk analysis using any arbitrary threat model establishes an ordering of components with regard to cyber-risk. This ordering can be used to influence the overall system design with an eye to lowering risk, or as a way to understand real-time risk to operational systems based on a current threat landscape. Expert knowledge of both the analysis process and the system being analyzed is required to perform a HAZCADS analysis. The need for a tiered risk analysis is demonstrated by the results of this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗

Measure this, not that: Optimizing the cost and model-based information content of measurements

Model-based design of experiments (MBDoE) is a powerful framework for selecting and calibrating science-based mathematical models from data. Here, this work extends popular MBDoE workflows by proposing a convex mixed integer (non)linear programming (MINLP) to optimize the selection of measurements. The solver MindtPy is modified to support calculating the D-optimality objective and its gradient via an external package, scipy, using the grey-box module in Pyomo. The new approach is demonstrated in two case studies: estimating highly correlated kinetics from a batch reactor and estimating transport parameters in a large-scale rotary packed bed for CO 2 capture. Both case studies show how examining the Pareto optimal trade-offs between information content measured by A- and D-optimality versus measurement budget offers practical guidance for selecting measurements for scientific experiments.

97 MATHEMATICS AND COMPUTING↗

Thermodynamics of quantum measurement and Maxwell's demon's arrow of time

Here, we discuss the thermodynamic aspects of a single qubit based device, powered by weak quantum measurements, and feedback controlled by a quantum Maxwell's demon. We discuss both discrete and time-continuous operation of the measurement based device at finite temperature of the reservoir. In the discrete example where a demon acquires information via discrete weak measurements, we find that the thermodynamic variables including the heat exchanged, extractable work, and the entropy produced are completely determined by an information theoretic measure of the demon's perceived arrow of time. We also discuss a realistic time-continuous operation of the device where the feedback is applied after a sequence of weak measurements. In the time-continuous limit, we derive the exact finite-time statistics of work, heat and entropy changes along individual quantum trajectories of the quantum measurement process, and relate them to the demon's arrow of time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hypersonic Second-Mode Instability Response to Shaped Roughness

An experimental campaign was conducted on a 7-degree half-angle cone in the NASA Langley Research Center 20-Inch Mach 6 Wind Tunnel to examine the influence of arrays of regularly spaced roughness elements on instability growth and transition. The primary element shape was a pair of elliptical planform ramps that were inclined at equal and opposite angles with respect to the local streamwise direction. The element shapes were designed to induce transient growth disturbances that would lead to sustained azimuthal modulation of the boundary layer flow while limiting the nearfield disturbances to avoid an immediate, i.e., effective tripping of the boundary layer. The bulk of the run matrix consisted of testing different element height sat free stream unit Reynolds numbers ranging from 9.8 to 13.1 million per meter. Other element shapes previously designed for tripping hypersonic boundary layers were also implemented. The model was instrumented with surface mounted Kulite ® and PCB ® pressure transducers and thermocouples. Spectra from the PCBs ® indicated clear suppression of the second-mode instability; however, neither the PCB ® spectra nor the heat transfer data presented strong evidence for delayed turbulent flow. Complementary stability computations likewise demonstrated second-mode reduction, particularly just downstream of the roughness, but also revealed a rise in first mode (streak-instability) amplitudes from the baseline that was likely responsible for the earlier transition observed for taller roughness cases.

Hypersonics↗

Improved Underwater Excitation-Emission Matrix Fluorometer

A compact, high-resolution, two-dimensional excitation-emission matrix fluorometer (EEMF) has been designed and built specifically for use in identifying and measuring the concentrations of organic compounds, including polluting hydrocarbons, in natural underwater settings. Heretofore, most EEMFs have been designed and built for installation in laboratories, where they are used to analyze the contents of samples collected in the field and brought to the laboratories. Because the present EEMF can be operated in the field, it is better suited to measurement of spatially and temporally varying concentrations of substances of interest. In excitation-emission matrix (EEM) fluorometry, fluorescence is excited by irradiating a sample at one or more wavelengths, and the fluorescent emission from the sample is measured at multiple wavelengths. When excitation is provided at only one wavelength, the technique is termed one-dimensional (1D) EEM fluorometry because the resulting matrix of fluorescence emission data (the EEM) contains only one row or column. When excitation is provided at multiple wavelengths, the technique is termed two-dimensional (2D) EEM fluorometry because the resulting EEM contains multiple rows and columns. EEM fluorometry - especially the 2D variety - is well established as a means of simultaneously detecting numerous dissolved and particulate compounds in water. Each compound or pool of compounds has a unique spectral fluorescence signature, and each EEM is rich in information content, in that it can contain multiple fluorescence signatures. By use of deconvolution and/or other mixture-analyses techniques, it is often possible to isolate the spectral signature of compounds of interest, even when their fluorescence spectra overlap. What distinguishes the present 2D EEMF over prior laboratory-type 2D EEMFs are several improvements in packaging (including a sealed housing) and other aspects of design that render it suitable for use in natural underwater settings. In addition, the design of the present 2D EEMF incorporates improvements over the one prior commercial underwater 2D EEMF, developed in 1994 by the same company that developed the present one. Notable advanced features of the present EEMF include the following: 1) High sensitivity and spectral resolution are achieved by use of an off-the-shelf grating spectrometer equipped with a sensor in the form of a commercial astronomical- grade 256 532-pixel charge-coupled-device (CCD) array. 2) All of the power supply, timing, control, and readout circuits for the illumination source and the CCD, ancillary environmental monitoring sensors, and circuitry for controlling a shutter or filter motor are custom-designed and mounted compactly on three circuit boards below a fourth circuit board that holds the CCD (see figure). 3) The compactness of the grating spectrometer, CCD, and circuit assembly makes it possible to fit the entire instrument into a compact package that is intended to be maneuverable underwater by one person. 4) In mass production, the cost of the complete instrument would be relatively low - estimated at approximately $30,000 at 2005 prices.

Moore, Casey↗

Microscopic Origins of the Nonlinear Behavior of Particle-Filled Rubber Probed with Dynamic Strain XPCS

The underlying microscopic response of filler networks in reinforced rubber to dynamic strain is not well understood due to the experimental difficulty of directly measuring filler network behavior in samples undergoing dynamic strain. This difficulty can be overcome with in situ X-ray photon correlation spectroscopy (XPCS) measurements. The contrast between the silica filler and the rubber matrix for X-ray scattering allows us to isolate the filler network behavior from the overall response of the rubber. This in situ XPCS technique probes the microscopic breakdown and reforming of the filler network structure, which are responsible for the nonlinear dependence of modulus on strain, known in the rubber science community as the Payne effect. These microscopic changes in the filler network structure have consequences for the macroscopic material performance, especially for the fuel efficiency of tire tread compounds. In this work, we elucidate the behavior with in situ dynamic strain XPCS experiments on industrially relevant, vulcanized rubbers filled (13 vol %) with novel air-milled silica of ultrahigh-surface area (UHSA) (250 m 2 /g). The addition of a silane coupling agent to rubber containing this silica causes an unexpected and counterintuitive increase in the Payne effect and decrease in energy dissipation. For this rubber, we observe a nearly two-fold enhancement of the storage modulus and virtually equivalent loss tangent compared to a rubber containing a coupling agent and conventional silica. Interpretation of our in situ XPCS results simultaneously with interpretation of traditional dynamic mechanical analysis (DMA) strain sweep experiments reveals that the debonding or yielding of bridged bound rubber layers is key to understanding the behavior of rubber formulations containing the silane coupling agent and high-surface area silica. These results demonstrate that the combination of XPCS and DMA is a powerful method for unraveling the microscale filler response to strain which dictates the dynamic mechanical properties of reinforced soft matter composites. With this combination of techniques, we have elucidated the great promise of UHSA silica when used in concert with a silane coupling agent in filled rubber. Such composites simultaneously exhibit large moduli and low hysteresis under dynamic strain.

36 MATERIALS SCIENCE↗

Hypersonic Second Mode Instability Response to Shaped Roughness

An experimental campaign was conducted on a 7-degree half-angle cone in the NASA Langley Research Center 20-Inch Mach 6 Wind Tunnel to examine the influence of arrays of regularly spaced roughness elements on instability growth and transition. The primary element shape was a pair of elliptical planform ramps that were inclined at equal and opposite angles with respect to the local streamwise direction. The element shapes were designed to induce transient growth disturbances that would lead to sustained azimuthal modulation of the boundary layer flow while limiting the nearfield disturbances to avoid an immediate, i.e., effective tripping of the boundary layer. The bulk of the run matrix consisted of testing different element height sat free stream unit Reynolds numbers ranging from 9.8 to 13.1 million per meter. Other element shapes previously designed for tripping hypersonic boundary layers were also implemented. The model was instrumented with surface mounted Kulite ® and PCB ® pressure transducers and thermocouples. Spectra from the PCBs ® indicated clear suppression of the second-mode instability; however, neither the PCB ® spectra nor the heat transfer data presented strong evidence for delayed turbulent flow. Complementary stability computations likewise demonstrated second-mode reduction, particularly just downstream of the roughness, but also revealed a rise in first mode (streak-instability) amplitudes from the baseline that was likely responsible for the earlier transition observed for taller roughness cases.

Hypersonics↗

Static and Fatigue Characterization of Large Composite T-Bolt Connections in Marine Hygrothermal Environments

Fiber-reinforced polymer composites have been highlighted as ideal candidates for structural applications in marine renewable energy devices, such as tidal turbines and wave energy converters. It is well understood that harsh marine environments can cause strength degradation of composite laminates, which has been extensively researched at the coupon scale; however, no research has investigated how this translates into larger-scale composite structures. This paper presents a subcomponent-scale study which investigates the effects of hygrothermal aging and subsequent static and fatigue characterization of thick composite T-bolt connections as part of a large, multilaboratory materials research effort. Of the glass-reinforced epoxy and vinylester-epoxy matrix composites tested, both showed measurable static strength degradation (4–36%) after being hygrothermally aged, even though the composite specimens were only partially saturated with water. Under tension–tension fatigue loading, the epoxy specimens performed very well in their dry states but exhibited significant degradation after hygrothermal aging. In comparison, the vinylester-epoxy specimens had much shorter fatigue lives in their dry states but exhibited no degradation after hygrothermal aging. Overall, this research demonstrates that hygrothermal aging can have significant effects on the ultimate strengths and fatigue lives of even partially saturated thick composite T-bolt connections, indicating that degradation of the outer plies on thick composite laminates can have pronounced effects on the whole structure. It discusses the challenges of building an understanding of the effects of harsh marine environments in large-scale composite structures.

16 TIDAL AND WAVE POWER↗

Development of Novel Ferritic-Martensitic Steels with Superior Creep Properties for Power Plant Applications

A creep resistant martensitic steel, CPJ-7, was developed with an operating temperature approaching 650°C. Subsequently, another creep resistant martensitic steel, JMP, was designed with the potential to operate at, or slightly above, 650°C. This report describes the development of these alloys from early iterations of CPJ steels to the final CPJ-7 formulation as well as the more heat resistant JMP steel formulation. The design originated from computational modeling for phase stability and precipitate strengthening using fifteen constituent elements. Approximately forty heats of CPJ and ten of JMP, each weighing ~7 kg, were vacuum induction melted. A computationally optimized heat treatment schedule was developed to homogenize the ingots prior to hot forging and rolling prior to final normalization and tempering. Overall, wrought and cast versions of CPJ-7 present superior creep properties when compared to wrought and cast versions of COST alloys for steam turbine and wrought and cast versions of P91/92 for boiler applications. For instance, the Larson Miller Parameter curve for CPJ-7 at 650°C almost coincides with that of COST E at 620°C. The prolonged creep life was attributed to slowing down the process of the destabilization of the MX and M 23 C 6 precipitates at 650°C. On the other hand, the cast version of CPJ-7 also revealed superior mechanical performance, well above commercially available cast 9% Cr martensitic steel or derivatives, especially for creep life at 650°C. The casting process employed slow cooling to simulate the conditions of a thick wall full-size steam turbine casing but utilized a separate homogenization step prior to final normalization and tempering. To advance the development of CPJ-7 for commercial applications, a process was used to scale up the production of the alloy using vacuum induction melting (VIM) and electroslag remelting (ESR), which underlined the importance of melt processing control of intentionally designed minor and trace elements in these advanced alloys. Following the work on CPJ-7, the JMP steels were designed with higher Co for increased solid solution strengthening, Si for oxidation resistance and increased W (with low Mo content) for matrix strength and stability as well as solid solution strengthening. The JMP steels showed increases in creep life compared to CPJ-7 between 118 to 150% at 650°C for testing at various stresses between 138 MPa and 207 MPa. On a Larson-Miller plot, the performance of the JMP steels surpasses that of state-of-the-art MARBN and other MARBN-type steels. The influence of various elements within the composition of the alloys on the microstructure and mechanical properties are discussed. This report presents approximately 420,000 h of in-house creep testing, the equivalent of almost 50 years of cumulative creep tests.

36 MATERIALS SCIENCE↗

Microcracking in composite laminates under thermal and mechanical loading

Composites used in space structures are exposed to both extremes in temperature and applied mechanical loads. Cracks in the matrix form, changing the laminate thermoelastic properties. The goal of the present investigation is to develop a predictive methodology to quantify microcracking in general composite laminates under both thermal and mechanical loading. This objective is successfully met through a combination of analytical modeling and experimental investigation. In the analysis, the stress and displacement distributions in the vicinity of a crack are determined using a shear lag model. These are incorporated into an energy based cracking criterion to determine the favorability of crack formation. A progressive damage algorithm allows the inclusion of material softening effects and temperature-dependent material properties. The analysis is implemented by a computer code which gives predicted crack density and degraded laminate properties as functions of any thermomechanical load history. Extensive experimentation provides verification of the analysis. AS4/3501-6 graphite/epoxy laminates are manufactured with three different layups to investigate ply thickness and orientation effects. Thermal specimens are cooled to progressively lower temperatures down to -184 C. After conditioning the specimens to each temperature, cracks are counted on their edges using optical microscopy and in their interiors by sanding to incremental depths. Tensile coupons are loaded monotonically to progressively higher loads until failure. Cracks are counted on the coupon edges after each loading. A data fit to all available results provides input parameters for the analysis and shows them to be material properties, independent of geometry and loading. Correlation between experiment and analysis is generally very good under both thermal and mechanical loading, showing the methodology to be a powerful, unified tool. Delayed crack initiation observed in a few cases is attributed to a lack of preexisting flaws assumed by the analysis. Some interactions between adjacent ply groups are attributed to local stress concentrations. These two effects are not captured by the analysis due to its global nature. The analysis is conservative in these cases and agrees well with data after the observed onset of cracking.

Maddocks, Jason R.↗

Investigation of Ga 2 O 3 as a new transparent conductive oxide for photovoltaics applications

This small innovation project intended to leverage recent activity and advances in ultrawide-bandgap oxide power electronics toward the investigation of Ga 2 O 3 as a transparent conductive oxide (TCO) for use in photovoltaics (PV). Ga 2 O 3 ’s theoretical advantage over incumbent TCOs is its very large bandgap of 4.8 eV, ensuring optical transparency of photons with λ ≥ 260 nm, effectively the full terrestrial solar spectrum. At the time of proposal writing, literature on Ga 2 O 3 as a TCO/optical material was relatively sparse — what little existed was focused on for UV sensors and/or emitters, with none related to PV. As such, this project was intended to help determine the potential of Ga 2 O 3 as a PV-oriented TCO by investigating its deposition using tools common to the PV and electronics industries — atomic layer deposition (ALD) and RF sputtering (RS), both of which are already used in PV manufacturing — and its resultant optical and electronic properties. Final project goals were to test its application to Si and III-V solar cells. Controllable film thickness and excellent uniformity was established for both methods deposition methods, with ALD providing higher precision for thinner films and RS more effective for growing thicker films. The as-deposited films were found to be amorphous in nature. Spectroscopic ellipsometry (SE) confirmed bandgaps of at least 4.8 eV and non-parasitic absorption across both AM0 and AM1.5G spectra; a refractive index approaching the expected value of 1.8 was also observed, with some degree of tunability based on process parameters. Using this initial optical data, a transfer matrix model, which interfaces with in-house EQE and LIV models, was developed to simulate Ga 2 O 3 optical effects on various solar cells, included potential in antireflection coatings (ARCs). Despite promising optical properties, the resultant resistivity / conductivity metrics did not meet expectations, regardless of deposition process. Although thick RS-deposited films, using both undoped and 1 at% Ge (n-type) doped sintered targets — demonstrated high net carrier concentrations (via C-V), and low specific resistance Ohmic contacts were demonstrated, transmission line measurements (TLM) showed very high resistivity. Further analysis indicated high vertical film conductivity, but very low lateral conductivity. Deeper characterization revealed a relatively high density of nano/microcrystalline inclusions that were ostensibly the source of the vertical conductivity, with the amorphous matrix serving as an effective insulator (likely due to high concentrations of electronic trap states). Despite extensive work to increase film polycrystallinity — demonstrated via high temperature annealing — sufficient lateral conductivity for TCO use was not achieved. The final phase of the project shifted focus toward investigation of ALD-deposited Ga 2 O 3 ’s potential as a passivant and/or passivating contact for both Si and III-V solar cells, with expectation of performance similar to Al 2 O 3 . However, initial rounds of testing using our baseline ALD process yielded minimal (but not quite zero) passivation of both GaInP and Si surfaces. Dielectric passivation is well-known to be highly process sensitive; given the unoptimized nature of the process used, this work was deemed inconclusive. The final outlook with respect to feasibility of Ga 2 O 3 as a PV-oriented TCO is, ultimately, still uncertain. The work performed in this project confirmed the optical properties and deposition methods, but the achieved electrical properties do not meet technological needs; literature reports in recent years are somewhat inconsistent and potentially untrustworthy, but generally appear to be in line with our results. At the very least it is clear that, should Ga 2 O 3 still be under consideration for PV TCO and/or selective contact applications, a significant amount of optimization and study is still needed.

14 SOLAR ENERGY↗

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Performance SiC/SiC Ceramic Composite Systems Developed for 1315 C (2400 F) Engine Components

As structural materials for hot-section components in advanced aerospace and land-based gas turbine engines, silicon carbide (SiC) ceramic matrix composites reinforced by high performance SiC fibers offer a variety of performance advantages over current bill-of-materials, such as nickel-based superalloys. These advantages are based on the SiC/SiC composites displaying higher temperature capability for a given structural load, lower density (approximately 30- to 50-percent metal density), and lower thermal expansion. These properties should, in turn, result in many important engine benefits, such as reduced component cooling air requirements, simpler component design, reduced support structure weight, improved fuel efficiency, reduced emissions, higher blade frequencies, reduced blade clearances, and higher thrust. Under the NASA Ultra-Efficient Engine Technology (UEET) Project, much progress has been made at the NASA Glenn Research Center in identifying and optimizing two highperformance SiC/SiC composite systems. The table compares typical properties of oxide/oxide panels and SiC/SiC panels formed by the random stacking of balanced 0 degrees/90 degrees fabric pieces reinforced by the indicated fiber types. The Glenn SiC/SiC systems A and B (shaded area of the table) were reinforced by the Sylramic-iBN SiC fiber, which was produced at Glenn by thermal treatment of the commercial Sylramic SiC fiber (Dow Corning, Midland, MI; ref. 2). The treatment process (1) removes boron from the Sylramic fiber, thereby improving fiber creep, rupture, and oxidation resistance and (2) allows the boron to react with nitrogen to form a thin in situ grown BN coating on the fiber surface, thereby providing an oxidation-resistant buffer layer between contacting fibers in the fabric and the final composite. The fabric stacks for all SiC/SiC panels were provided to GE Power Systems Composites for chemical vapor infiltration of Glenn designed BN fiber coatings and conventional SiC matrices. Composite panels with system B were heat treated at Glenn, and the pores that remained open were filled by silicon melt infiltration (MI). Panels with system A and the other SiC/SiC systems were not heat treated, and remaining open pores in these systems were filled with SiC slurry and silicon MI.

DiCarlo, James A.↗

Fluid-like shock compression of dilute polymer nanocomposites

The shock Hugoniot of heterogenous mixtures of discrete particles has been experimentally investigated for porous agglomerates and fluids, while the study of full-density solids has been primarily limited to compressed powders and high-dimensional composites. By dispersing ceria nanoparticles in a polyethylene matrix, we are able to examine the hydrodynamic behavior of a nonporous, heterogenous solid in thermal equilibrium during weak shock compression. Phase-driven discontinuities in the Hugoniot particle velocity–shock velocity (u−D) relationship of pure polyethylene are replicated in the nanocomposites but are shifted to lower velocity and to higher pressure with higher particle concentration. The results are explained using an isothermal, two-velocity fluid model under the hydrodynamic approximation. The model, which assumes a theoretical equation-of-state for ceria and either a low-order or high-order fit to the measured polyethylene Hugoniot, reasonably predicts the Hugoniot for two different polyethylene/ceria nanocomposites. Using the model, the mixture Hugoniot is shown to be insensitive to the Hugoniot of the stiffer constituent when the moduli are sufficiently disparate, while dependence on particle density and volume fraction is preserved through fluid-like motion.

Moore, Nathan W. [Sandia National Laboratories (SN↗

Initial development of an RIA envelope for dispersed nuclear fuel

A reactivity insertion accident (RIA) is a design basis accident in which reactivity is rapidly injected as a result of a control rod ejection or blade drop scenario. The resultant power increase can result in fuel rod failure and subsequent release of radioactive material into the reactor’s primary system. For light water reactors, failure under RIAs is typically a result of fuel melt, pellet-cladding interaction, or rod over-pressurization. Fuel melt occurs when the fuel system is unable to transport heat out of the fuel system. Pellet-cladding interaction occurs when an aggressive fuel pellet expansion pushes on an embrittled cladding beyond its strain limit. Rod over-pressurization occurs when a fuel rod exceeds the critical heat flux and departs from nucleate boiling, causing the rod internal pressure to rapidly increase, exceeding the systems pressure, and balloon until it bursts. To prevent failure, regulators have developed safety requirements that limit the injected enthalpy based on the state of the fuel system. However, dispersed fuel could enable the removal or relaxation of regulator-imposed safety criteria. Dispersed fuel embeds fuel particles in a highly conductive metal. Dispersed particles can have a particle radius of 1 µm to greater than 100 µm. Traditional fuel systems are plagued by poor thermal conductivity, especially at higher burnups, thus reducing the ability of the fuel system to respond to an RIA. However, dispersed fuel could improve the ability of the fuel system to transport injected energy away from the fissile material and into the coolant more efficiently. Furthermore, dispersing the fuel particles mitigates hard contact that typically occurs in the traditional zirconium-uranium dioxide fuel system, thus removing pellet-cladding mechanical interaction as a potential failure mechanism. This paper evaluates fuel particles dispersed in a metal matrix using the BISON fuel performance code to develop an initial failure threshold based on melt temperature of the fuel particle and/or cladding material as a function of fission density for beginning-of-life and end-of-life conditions. BISON results show that (1) reducing the size of the particles allows fission density to increase as a function of time, (2) particle-to-particle proximity must be considered to evaluate the limiting conditions, and (3) improving the fuel particle thermal conductivity improves thermal performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗