Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interior methods”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Direct Simulation Monte Carlo Studies of the Gas Sampling for the VATMOS-SR Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept which would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere, and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus, with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a full picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of atmospheric volatile signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process in order to determine how measured isotopic ratios of noble gases in the sample compare to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept, and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed at Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield a comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters as well as the modeling parameters will be performed.

direct simulation Monte Carlo↗
NASA NTRS2000

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. The most direct approach to this problem is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing energy and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following questions about these proteins: (1) How do small proteins self-organize into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides form membrane-spanning structures (e.g. channels)? (3) By what mechanisms do such structures perform their functions? The simulations are performed using the molecular dynamics method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated in real space while long-range forces are evaluated in reciprocal space, using a particle-mesh algorithm which is of order O(NInN). With a time step of 2 femtoseconds, problems occurring on multi-nanosecond time scales (10(exp 6)-10(exp 8) time steps) are accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude, which requires algorithmic improvements and codes scalable to a large number of processors. Work in this direction is in progress. Two series of simulations are discussed. In one series, it is shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and fold into helical structures (see Figure). Once in the membrane, the peptides exhibit orientational flexibility with changing conditions, which may have provided a mechanism of transmitting signals between the protocell and its environment. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations identify the mechanisms by which protons move through the gate.

Pohorille, Andrew↗
DOE OSTI2025

Exploring strain rate effects upon 3D materials using high speed in situ X-ray tomoscopy

Cellular materials are ubiquitous in our modern society. They may be stochastic gas-blown foams (e.g., polyurethane), foamed starches (e.g., cereals), or, in this case, 3D printed microlattices. Failure in these materials is often driven by surface or sub-surface defects, which may be nucleated at a surface roughness, an interior void, or inclusion interfaces that may not be typically observable. Obfuscating our understanding further, bulk materials are known to exhibit strain-rate-dependent mechanical response, making a subsurface understanding of damage even more critical. For the first time, an in situ uniaxial mechanical loading stage that simultaneously rotates specimens up to 18 Hz was fielded at a synchrotron for 3D tomographic imaging. This capability opens a plethora of materials science opportunities to explore strain rate effects in materials and examining deformation, fracture, and delamination’s (in composites) for a complete 3D picture (movie) of material response. We demonstrate the deformation of 3D printed polymer lattice structures, of three different material types, at 0.25, 1.1, and 2.2 s −1 strain rates. We successfully imaged the 3D deformation of these materials and can directly compare the same printed structure to the three material types at three strain rates, all in 3D. Material point method simulations were applied to one of the materials to better understand the role of voids on the 3D printed structure’s performance.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Acquisition and Retaining Granular Samples via a Rotating Coring Bit

This device takes advantage of the centrifugal forces that are generated when a coring bit is rotated, and a granular sample is entered into the bit while it is spinning, making it adhere to the internal wall of the bit, where it compacts itself into the wall of the bit. The bit can be specially designed to increase the effectiveness of regolith capturing while turning and penetrating the subsurface. The bit teeth can be oriented such that they direct the regolith toward the bit axis during the rotation of the bit. The bit can be designed with an internal flute that directs the regolith upward inside the bit. The use of both the teeth and flute can be implemented in the same bit. The bit can also be designed with an internal spiral into which the various particles wedge. In another implementation, the bit can be designed to collect regolith primarily from a specific depth. For that implementation, the bit can be designed such that when turning one way, the teeth guide the regolith outward of the bit and when turning in the opposite direction, the teeth will guide the regolith inward into the bit internal section. This mechanism can be implemented with or without an internal flute. The device is based on the use of a spinning coring bit (hollow interior) as a means of retaining granular sample, and the acquisition is done by inserting the bit into the subsurface of a regolith, soil, or powder. To demonstrate the concept, a commercial drill and a coring bit were used. The bit was turned and inserted into the soil that was contained in a bucket. While spinning the bit (at speeds of 600 to 700 RPM), the drill was lifted and the soil was retained inside the bit. To prove this point, the drill was turned horizontally, and the acquired soil was still inside the bit. The basic theory behind the process of retaining unconsolidated mass that can be acquired by the centrifugal forces of the bit is determined by noting that in order to stay inside the interior of the bit, the frictional force must be greater than the weight of the sample. The bit can be designed with an internal sleeve to serve as a container for granular samples. This tube-shaped component can be extracted upon completion of the sampling, and the bottom can be capped by placing the bit onto a corklike component. Then, upon removal of the internal tube, the top section can be sealed. The novel features of this device are: center dot A mechanism of acquiring and retaining granular samples using a coring bit without a closed door. center dot An acquisition bit that has internal structure such as a waffle pattern for compartmentalizing or helical internal flute to propel the sample inside the bit and help in acquiring and retaining granular samples. center dot A bit with an internal spiral into which the various particles wedge. center dot A design that provides a method of testing frictional properties of the granular samples and potentially segregating particles based on size and density. A controlled acceleration or deceleration may be used to drop the least-frictional particles or to eventually shear the unconsolidated material near the bit center.

Bar-Cohen, Yoseph↗
NASA NTRS1988

Impact-induced devolatilization and hydrogen isotopic fractionation of serpentine: Implications for planetary accretion

Impact-induced devolatilization of porous serpentine was investigated using two independent experimental methods, the gas recovery and the solid recovery method, each yielding nearly identical results. For shock pressures near incipient devolatilization, the hydrogen isotopic composition of the evolved H2O is very close to that of the starting material. For shock pressures at which up to 12 percent impact-induced devolatilization occurs, the bulk evolved gas is significantly lower in deuterium than the starting material. There is also significant reduction of H2O to H2 in gases recovered at these higher shock pressures, probably caused by reaction of evolved H2O with the metal gas recovery fixture. Gaseous H2O-H2 isotopic fractionation suggests high temperature isotopic equilibrium between the gaseous species, indicating initiation of devolatilization at sites of greater than average energy deposition. Bulk gas-residual solid isotopic fractionations indicate nonequilibrium, kinetic control of gas-solid isotopic ratios. Impact-induced hydrogen isotopic fractionation of hydrous silicates during accretion can strongly affect the long-term planetary isotopic ratios of planetary bodies, leaving the interiors enriched in deuterium. Depending on the model used for extrapolation of the isotopic fractionation to devolatilization fractions greater than those investigated experimentally can result from this process.

Tyburczy, James A.↗
NASA NTRS2016

Shrink Tape Technique for Heat Forming Aluminum Substrates for Thin Foil X-Ray Mirrors and the Neutron Star Interior Composition Explorer X-Ray Concentrators

Consistent improvements in the design and fabrication of thin-foil, epoxy-replicated x-ray mirrors for astronomical telescopes have yielded increasingly higher quality and more precise astrophysical data. The Neutron Star Interior Composition Explorer (NICER) x-ray timing mission optics continues this tradition and introduces design elements that promise even more accurate measurements and precise astrophysical parameters. The singly reflecting concentrators have a curved axial profile to improve photon concentration and a sturdy full shell structure for enhanced module stability. These design elements introduced the challenge of reliably forming mirror substrates at an acceptable production rate. By developing a technique using heat shrink tape to compress and conform thin aluminum mirror substrates to shaping mandrels, production rate improved with successful fabrication. The technique's efficiency was analyzed by measuring hundreds of substrate profiles postforming, performance testing completely assembled concentrators composed of every size substrate, and comparing the results to simulated fabrication scenarios. On average, the profiles were copied within 4.6 plus or minus 3.7 percent. These measurements and the overall success of NICER's optics, via ground calibration, have shown that the heat-shrink tape method is reliable, repeatable, and could be used in future missions to increase production rate and improve performance.

Balsamo, Erin↗
DOE OSTI2024

On a Spectral Method for β -particle Bound Excitation Collisions in Kilonovae

The interaction of β-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova (KN) transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti following the seminal formulation of Bethe, in a spectral solver of the Vlasov–Maxwell–Boltzmann equation. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon–matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric KN ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for β-particles in the KN is a bulk velocity parameter u = 0, a thermal velocity parameter α = 0.5c, and a 9 × 9 × 9 mode velocity basis set (Hermite orders of 0–8 in each dimension). For interior-ejecta sample zones, we estimate that the ratio of thermalization from large-angle (≳2fdg5) bound excitation scattering to total thermalization is ~0.002–0.003.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

The determination of the mass and mean density of Enceladus from its observed shape

Application of limb-fitting methods to the 11 best Voyager 2 images of Enceladus has shown that the shape of this satellite is closely represented by a triaxial ellipsoid. The observed ratio of the differences of the principal axes F = (b - c)/(a - c) is 0.23 (sup +0.04 sub - 0.01), consisting with the value F = 0.23 expected for a synchronously rotating satellite in hydrostatic equilibrium. We also deduce from the Voyager observations, after allowing for limb topography, that the mean radius of the satellite is 249.4 +/- 0.2 km. For satellites of known mass, measurement of the size and shape leads to a determination of the satellite's mean density and moment of inertial. We have used this method to determine the moments of inertia of Mimas (1988) and Tethys (1991). Enceladus appears to be hydrostatically relaxed, making it an ideal candidate for this type of analysis. However, none of the Pioneer or Voyager spacecraft had a close encounter with this satellite and thus its mass is effectively unknown. Enceladus is trapped in a 2:1 orbit-orbit resonance with Dione, but the amplitudes of libration are too small to allow a useful mass determination. Using the observed shape alone, without any other assumptions other than that the satellite is in hydrostatic equilibrium at its present orbital radius, we place an upper bound on the mean density of 1.12 +/- 0.05 g/cu cm. Thus, the mean density of Enceladus is probably little more than that of water-ice and we conclude that this satellite is markedly deficient in rock. If the mass of a satellite is unknown, but the satellite is differentiated and has a deep mantle of known composition, then we show that measurement of the shape alone can lead to a determination of the satellite's mass, mean density, and moment of inertia. Application of this method to Enceladus, assuming that the satellite has a deep mantle of water-ice of density 0.93 g/cu cm, gives the result that the mean density of the satellite is 1.00 +/- 0.03 g/cu cm. This result fills the one remaining gap in our knowledge of the structure of the Saturnian satellite system. We now know the mean densities of all the primary Saturnian satellites in the sequence from the coorbital satellites, Janus and Epimetheus, through to the outer satellite Iapetus (the densities of the small, secondary satellites in Trojan-type orbits are still unknown). The Saturnian system possesses two striking features. (1) Because of significant porosity, the mean material densities of the satellites Janus, Epimetheus, and Mimas could be substantially greater than the apparent mean densities of these satellies. (2) The densities of the satellites are not correlated with their distances from the planet; in particular, the satellites Enceladus and Tethys have lower mean densities than their interior and exterior neighbors, Mimas and Dione. This may be the result of gross postformation redistribution of rock and ice, possibly due to satellite disruptions are suggested by Smith et al. (1982).

Dermott, Stanley F.↗
DOE OSTI2021

Deployment of salt sample extraction system at an engineering-scale electrorefiner

The goal of the salt sampling program at Argonne is to develop and deploy automated molten salt sampling approaches for interfacing relevant unit operations with salt analysis to improve the timeliness of sampling-based accountancy measurements. Two technologies under development in support of this goal are a vacuum sampling loop module and a high-throughput pneumatic sample generator module. Compared to traditional point sampling approaches (i.e., dip probes), the vacuum sampling loop facilitates the collection of a larger cross-section of the bulk salt in order to collect more representative samples. The vacuum sampling approach also eliminates the risk of dross contamination of samples and avoids the use of moving parts in the salt. The pneumatic sample generator module is used to facilitate high-throughput sample analysis to improve the measurement precision of any given analytical technique by averaging out random sampling and measurement errors. In FY21, two methods for integrating these two modules were tested including direct fluidic coupling and coupling using a solid salt transfer mechanism. Solid salt transfer was ultimately selected over fluidic coupling, primarily to enable the transport of samples over longer distances to support automated at-line integration with high-precision techniques (such as microcalorimetry) that cannot withstand the conditions near an electrorefining process. To facilitate rapid solid salt coupling, new mechanisms were developed for rapidly charging and discharging salt sample tubes at the vacuum sampling loop and pneumatic sample generator modules, respectively. While the charging mechanism will be deployed in FY22, the discharge mechanism was tested in FY21 and is described here. The solid salt tube transfer method was deployed at one of Argonne’s engineering-scale electrorefiners to implement at-line high-throughput pneumatic micro-sample generation capabilities. The approach was used to generate precise uranium- and lanthanide-bearing electrorefiner micro-samples with the specific dimensions requested by researchers at Los Alamos National Laboratory for use in testing their novel microcalorimeter x-ray techniques. The solid salt transfer mechanism proved not only to be an effective means of integrating the precision sample generator with vacuum sampling, but also improved the performance of the sampler generator. To discharge salt from the sample tubes at the sampler generator, tube segments were inserted directly into the sample generator’s Helmholtz chamber and pressure pulse actuations were used to generate precision molten salt samples directly from the tube segments. The direct insertion of sample tubes into the sample generator enabled rapid loading of the salt and prevented salt from contacting most of the interior surfaces of the sample generator, which eliminated cross-contamination between runs. The vacuum sampling-loop tube charging mechanism will support high-throughput tube sampling operations by employing a dynamic vacuum filling process to fill short charge tubes that are configured to be rapidly connected and disconnected from the loop. The dynamic vacuum sampling operation will be automated, and sample tube handling can be executed with simple overhead actuation. Because the modular sampling approach described here eliminates the need for new high-radiation sample handling capabilities, salt-wetted seals, salt-wetted moving parts, and heated transfer lines outside the electrorefiner, it will address most of the remaining technical challenges for the at-line deployment of high-precision analytical techniques which will enable significant reductions in the time delay for sampling-based high-precision accountancy measurements.

42 ENGINEERING↗
NASA NTRS1989

Tissue simulating gel for medical research

A tissue simulating gel and a method for preparing the tissue simulating gel are disclosed. The tissue simulating gel is prepared by a process using water, gelatin, ethylene gylcol, and a cross-linking agent. In order to closely approximate the characteristics of the type of tissue being simulated, other material has been added to change the electrical, sound conducting, and wave scattering properties of the tissue simulating gel. The result of the entire process is a formulation that will not melt at the elevated temperatures involved in hyperthermia medical research. Furthermore, the tissue simulating gel will not support mold or bacterial growth, is of a sufficient mechanical strength to maintain a desired shape without a supporting shell, and is non-hardening and non-drying. Substances were injected into the tissue simulating gel prior to the setting-up thereof just as they could be injected into actual tissue, and the tissue simulating gel is translucent so as to permit visual inspection of its interior. A polyurethane spray often used for coating circuit boards can be applied to the surface of the tissue simulating gel to give a texture similar to human skin, making the tissue simulating gel easier to handle and contributing to its longevity.

Companion, John A.↗
NASA NTRS2006

Development of a Liner Design Methodology and Relevant Results of Acoustic Suppression in the Farfield for Mixer-Ejector Nozzles

We have developed a process to predict noise field interior to the ejector and in the farfield for any liner design for a mixer-ejector of arbitrary scale factor. However, a number of assumptions, not verified for the current application, utilized in this process, introduce uncertainties in the final result, especially, on a quantitative basis. The normal impedance model for bulk with perforated facesheet is based on homogeneous foam materials of low resistivity. The impact of flow conditions for HSCT application as well as the impact of perforated facesheet on predicted impedance is not properly accounted. Based on the measured normal impedance for deeper bulk samples (i.e., 2.0 in.) the predicted reactance is much higher compared to the data at frequencies above 2 kHz for T-foam and 200 ppi SiC. The resistance is under predicted at lower frequencies (below 4 kHz) for these samples. Thus, the use of such predicted data in acoustic suppression is likely to introduce inaccuracies. It should be noted that the impedance prediction methods developed recently under liner technology program are not utilized in the studies described in this report due to the program closeout. Acoustic suppression prediction is based on the uniform flow and temperature conditions in a two-sided treated constant area rectangular duct. In addition, assumptions of equal energy per mode noise field and interaction of all frequencies with the treated surface for the entire ejector length may not be accurate. While, the use of acoustic transfer factor minimizes the inaccuracies associated with the prediction for a known test case, the assumption of the same factor for other liner designs and with different linear scale factor ejectors seems to be very optimistic. As illustrated in appendix D that the predicted noise suppression for LSM-1 is lower compared to the measured data is an indication of the above argument. However, the process seems to be more reliable when used for the same scale models for different liner designs as demonstrated for Gen. 1 mixer-ejectors.

Salikuddin, M.↗
NASA NTRS1991

Tissue simulating gel for medical research

A tissue simulating gel and a method for preparing the tissue simulating gel are disclosed. The tissue simulating gel is prepared by a process using water, gelatin, ethylene glycol, and a cross-linking agent. In order to closely approximate the characteristics of the type of tissue being simulated, other material has been added to change the electrical, sound conducting, and wave scattering properties of the tissue simulating gel. The result of the entire process is a formulation that will not melt at the elevated temperatures involved in hyperthermia medical research. Furthermore, the tissue simulating gel will not support mold or bacterial growth, is of a sufficient mechanical strength to maintain a desired shape without a supporting shell, and is non-hardening and non-drying. Substances have been injected into the tissue simulating gel prior to the setting-up thereof just as they could be injected into actual tissue, and the tissue simulating gel is translucent so as to permit visual inspection of its interior. A polyurethane spray often used for coating circuit boards can be applied to the surface of the tissue simulating gel to give a texture similar to human skin, making the tissue simulating gel easier to handle and contributing to its longevity.

Companion, John A.↗
NASA NTRS2023

Using Satellite Data to Characterize Land Surface Processes in Morocco

This study endeavors to produce a comprehensive land cover map for Morocco, addressing the absence of such a detailed map in the country. Our research encompasses ecological and climatic aspects specific to Morocco, while the methods used can be adapted to various regions and countries, considering their unique climatic conditions and land cover types. A combination of MODIS and Landsat datasets was employed to create a 5 km resolution Land Use and Land Cover (LULC) map for the entire nation. The process involved the aggregation and advanced processing of these datasets using surface processes algorithms. The resulting LULC map is the first of its kind for Morocco, shedding light on land cover distribution nationwide. It shows that approximately 13.5% of the country is covered by forests, predominantly in the Atlas and Rif mountains, Rabat–Sale, and the southern regions. Grasslands occupy over 16% of the study area, mainly in the north-east and west. Urban areas, including major cities like Casablanca, Rabat, and Marrakech, span nearly 3400 km². Moreover, large areas of shrublands and bare lands are evident across the country, while agricultural lands account for almost 20% of the national territory, mainly in the interior plains and north-western Atlantic coast. This study forms a crucial basis for ecological and climatic research in Morocco and serves as a valuable reference for various disciplines such as agriculture, natural resource management, and climate modeling. The mapping of biophysical parameters for each land cover class is a key feature of our research, and these parameters will be instrumental in a subsequent study examining the impact of urban development on surface climate in Morocco. Overall, our study underscores the importance of understanding biophysical parameters in addressing environmental and societal challenges.

Mohammed Thaiki↗
NASA NTRS2015

A Fully Nonmetallic Gas Turbine Engine Enabled by Additive Manufacturing of Ceramic Composites: Additive Manufacturing and Characterization of Ceramic Composites - Part III

This publication is the third part of a three part report of the project entitled "A Fully Nonmetallic Gas Turbine Engine Enabled by Additive Manufacturing" funded by NASA Aeronautics Research Institute (NARI). The objective of this project was to conduct additive manufacturing to produce ceramic matrix composite materials and aircraft engine components by the binder jet process. Different SiC powders with median sizes ranging from 9.3 to 53.0 microns were investigated solely and in powder blends in order to maximize powder packing. Various infiltration approaches were investigated to include polycarbosilane (SMP-10), phenolic, and liquid silicon. Single infiltrations of SMP-10 and phenolic only slightly filled in the interior. When the SMP-10 was loaded with sub-micron sized SiC powders, the infiltrant gave a much better result of filling in the interior. Silicon carbide fibers were added to the powder bed to make ceramic matrix composite materials. Microscopy showed that the fibers were well distributed with no preferred orientation on the horizontal plane and fibers in the vertical plane were at angles as much as 45deg. Secondary infiltration steps were necessary to further densify the material. Two to three extra infiltration steps of SMP-10 increased the density by 0.20 to 0.55 g/cc. However, the highest densities achieved were 2.10 to 2.15 g/cc. Mechanical tests consisting of 4 point bend tests were conducted. Samples from the two CMC panels had higher strengths and strains to failure than the samples from the two nonfiber reinforced panels. The highest strengths were from Set N with 65 vol% fiber loading which had an average strength of 66 MPa. Analysis of the fracture surfaces did not reveal pullout of the reinforcing fibers. Blunt fiber failure suggested that there was not composite behavior. The binder jet additive manufacturing method was used to also demonstrate the fabrication of turbine engine vane components of two different designs and sizes. The binder jet method has benefits over the conventional manufacturing of CMCs in that prototype and production parts can be fabricated quickly and economically with no tooling and extensive hand layup.

Additive Manufacturing↗
NASA NTRS1980

Numerical generation of two-dimensional grids by the use of Poisson equations with grid control at boundaries

A method for generating boundary-fitted, curvilinear, two dimensional grids by the use of the Poisson equations is presented. Grids of C-type and O-type were made about airfoils and other shapes, with circular, rectangular, cascade-type, and other outer boundary shapes. Both viscous and inviscid spacings were used. In all cases, two important types of grid control can be exercised at both inner and outer boundaries. First is arbitrary control of the distances between the boundaries and the adjacent lines of the same coordinate family, i.e., stand-off distances. Second is arbitrary control of the angles with which lines of the opposite coordinate family intersect the boundaries. Thus, both grid cell size (or aspect ratio) and grid cell skewness are controlled at boundaries. Reasonable cell size and shape are ensured even in cases wherein extreme boundary shapes would tend to cause skewness or poorly controlled grid spacing. An inherent feature of the Poisson equations is that lines in the interior of the grid smoothly connect the boundary points (the grid mapping functions are second order differentiable).

Sorenson, R. L.↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Robust Feedback Control of Flow Induced Structural Radiation of Sound

A significant component of the interior noise of aircraft and automobiles is a result of turbulent boundary layer excitation of the vehicular structure. In this work, active robust feedback control of the noise due to this non-predictable excitation is investigated. Both an analytical model and experimental investigations are used to determine the characteristics of the flow induced structural sound radiation problem. The problem is shown to be broadband in nature with large system uncertainties associated with the various operating conditions. Furthermore the delay associated with sound propagation is shown to restrict the use of microphone feedback. The state of the art control methodologies, IL synthesis and adaptive feedback control, are evaluated and shown to have limited success for solving this problem. A robust frequency domain controller design methodology is developed for the problem of sound radiated from turbulent flow driven plates. The control design methodology uses frequency domain sequential loop shaping techniques. System uncertainty, sound pressure level reduction performance, and actuator constraints are included in the design process. Using this design method, phase lag was added using non-minimum phase zeros such that the beneficial plant dynamics could be used. This general control approach has application to lightly damped vibration and sound radiation problems where there are high bandwidth control objectives requiring a low controller DC gain and controller order.

Heatwole, Craig M.↗
NASA NTRS2004

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗