SEARCH · Engineering Papers

Results for “geometry optimization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Investigation of the Rocket Induced Flow Field in a Rectangular Duct

Rocket-Based Combined Cycle (RBCC) concepts attempt to improve the performance of launch vehicles at all points in the launch trajectory and make highly reusable launch vehicles a reality. The Aerojet Strutjet RBCC concept consists of a variable geometry duct with internal, vertical struts that functions in ducted rocket, ramjet, scramjet, and pure rocket modes. These struts have rocket and turbine exhaust nozzles imbedded within them. The rocket flows create an ejector effect with the ingested air at subsonic flight velocities. In ramjet and scramjet modes, the fuel rich nozzle flows react with the ingested air producing an afterburner effect. Under a NASA Marshall Space Flight Center contract, the UAH Propulsion Research Center (PRC) has designed and built a Strutjet simulation facility. A scale model of a single strut has been built and is undergoing cold-flow testing to investigate the mixing of the rocket and turbine exhausts with the ingested air. A complementary experimental program is also underway to examine the induced flow-field generated by rocket nozzles confined in a rectangular duct. Characterizing the induced flow behavior is critical to understanding and optimizing the performance of future Strutjet-based RBCC propulsion systems. The proposed paper will present results from the rocket induced flow investigation.

Landrum, D. Brian↗

NASA NTRS1993

LOITA: Lunar Optical/Infrared Telescope Array

LOITA (Lunar Optical/Infrared Telescope Array) is a lunar-based interferometer composed of 18 alt-azimuth telescopes arranged in a circular geometry. This geometry results in excellent uv coverage and allows baselines up to 5 km long. The angular resolution will be 25 micro-arcsec at 500 nm and the main spectral range of the array will be 200 to 1100 nm. For infrared planet detection, the spectral range may be extended to nearly 10 mu m. The telescope mirrors have a Cassegrain configuration using a 1.75 m diameter primary mirror and a 0.24 m diameter secondary mirror. A three-stage (coarse, intermediate, and fine) optical delay system, controlled by laser metrology, is used to equalize path lengths from different telescopes to within a few wavelengths. All instruments and the fine delay system are located within the instrument room. Upon exiting the fine delay system, all beams enter the beam combiner and are then directed to the various scientific instruments and detectors. The array instrumentation will consist of CCD detectors optimized for both the visible and infrared as well as specially designed cameras and spectrographs. For direct planet detection, a beam combiner employing achromatic nulling interferometry will be used to reduce star light (by several orders of magnitude) while passing the planet light. A single telescope will be capable of autonomous operation. This telescope will be equipped with four instruments: wide field and planetary camera, faint object camera, high resolution spectrograph, and faint object spectrograph. These instruments will be housed beneath the telescope. The array pointing and control system is designed to meet the fine pointing requirement of one micro-arcsec stability and to allow precise tracking of celestial objects for up to 12 days. During the lunar night, the optics and the detectors will be passively cooled to 70-80 K temperature. To maintain a continuous communication with the earth a relay satellite placed at the L4 libration point will be used in conjunction with the Advanced Tracking and Data Relay Satellite System (ATDRSS). Electrical power of about 10 kW will be supplied by a nuclear reactor based on the SP-100 technology. LOITA will be constructed in three phases of six telescopes each. The total mass of the first operational phase is estimated at 58,820 kg. The cost of the fully operational first phase of the observatory is estimated at $8.9 billion. LOITA's primary objectives will be to detect and characterize planets around nearby stars (up to ten parsec away), study physics of collapsed stellar objects, solar/stellar surface features and the processes in nuclear regions of galaxies and quasars. An interferometric array such as LOITA will be capable of achieving resolutions three orders of magnitude greater than Hubble's design goal. LOITA will also be able to maintain higher signal to noise ratios than are currently attainable due to long observation times available on the moon.

Source record↗

NASA NTRS1995

Micro-mechanics of micro-composites

The Structural Dynamics branch at NASA LaRC is working on developing an active passive mount system for vibration control. Toward this end a system utilizing piezoelectric actuators is currently being utilized. There are limitations to the current system related to space applications under which it is desired to eliminate deformations in the actuators associated with thermal effects. In addition, a material that is readily formable into complex shapes and whose mechanical properties can be optimized with regards to vibration control would be highly desirable. Microcomposite material are currently under study to service these needs. Microcomposite materials are essentially materials in which particles on the scale of microns are bound together with a polyimide (LaRC Si) that has been developed at LaRC. In particular a micro-composite consisting of LaRC Si binder and piezoelectric ceramic particles shows promise in satisfying the needs of the active passive mount project. The LaRC/ Si microcomposite has a unique combination of piezoelectric properties combined with a near zero coefficient of thermal expansion and easy machinability. The goal of this ASEE project is to develop techniques to analytically determine important material properties necessary to characterize the dynamic properties of actuators and mounts made from the LaRC Si / ceramic microcomposite. In particular, a generalized method of cells micromechanics originally developed at NASA Lewis is employed to analyze the microstructural geometry of the microcomposites and predict the overall mechanical properties of the material. A testing program has been established to evaluate and refine the GMC approach to these materials. In addition, a theory of mixtures analysis is being developed that utilizes the GMC micromechanics information to analyze complex behavior of the microcomposite material which has a near zero CTE.

Donovan, Richard P.↗

NASA NTRS1994

Fabric geometry distortion during composites processing

Waviness and tow misalignment are often cited as possible causes of data scatter and lower compression stiffness and strength in textile composites. Strength differences of as much as 40 percent have been seen in composites that appear to have the same basic material and structural properties -- i.e., yarn orientation, yarn size, interlacing geometry. Fabric geometry distortion has been suggested as a possible reason for this discrepancy, but little quantitative data or substantial evidence exists. The focus of this research is to contribute to the present understanding of the causes and effects of geometric distortion in textile composites. The initial part of the study was an attempt to gather qualitative information on a variety of textile structures. Existing and new samples confirmed that structures with a significant direction presence would be more susceptible to distortion due to the compaction process. Thus, uniweaves (fiber vol frac: 54-72 percent) biaxial braids (vf: 34-58 percent) demonstrated very little fabric geometry distortion. In stitched panels, only slight buckling of z-direction stitches was observed, primarily near the surface. In contrast, for structures with high compaction ratios -- e.g., large cylindrical yarns (2.5:1) orpowder towpreg (4:1) -- there were visible distortions where previously smooth and periodic undulations were transformed to abrupt changes in direction. A controlled study of the effect of forming pressure on distortion was conducted on type 162 glass plain weave fabrics. Panels (6 x 6 in) were produced via a resin infusion type setup, but with an EPON 815 epoxy resin. Pressures ranging from hand layup to 200 psi were used (vf: 34-54 percent). Photomicrographs indicated that at pressures up to 50 psi, large changes in thickness were due primarily to resin squeeze out. At higher pressures, when intimate contact was made between the layers, there was some tow flattening and in-plane shifting to optimize nesting. However, even at 200 psi the period and amplitude of the tow undulation remained constant, suggesting that for this relatively fine fabric, distortions from compaction were not a problem. Because of the interest in using larger tows (to reduce cost) and more complex structures, tests were also run on 2D triaxial glass braid (113 yd/lb at 0, 225 yd/lb at +/- 45). Forming pressures of 20, 50, 200, and 500 psi were used, and short block compression tests were run. The 500 psi specimen had a 10 percent decrease in modulus and an almost 50 percent decrease in strength (vs. 20 psi). Because the total fiber wgt/panel was kept constant, the thickness varied from 0.32 to 0.22 in (49-70 percent vf). Yet, the strength value is clearly below what would be expected, even with the decrease in thickness. Photomicrographs of these samples will be taken to determine if more fabric distortion exists in the 500 psi specimens. Finally, because the ultimate goal is to be able to predict and control distortion in a variety of textile structures, a model compaction test was developed to directly measure the deformation of the tows during compaction. Layers of dry glass fabric were placed in a mold with a clear plexiglass window. The yarn amplitude and period was then calculated using image analysis of the videotaped deformation. Preliminary tests demonstrated the feasibility of this technique for simple fabrics with large tows.

Chen, Julie↗

NASA NTRSDate not supplied

Surface meteorology and Solar Energy

The Release 5.1 Surface meteorology and Solar Energy (SSE) data contains parameters formulated for assessing and designing renewable energy systems. Parameters fall under 11 categories including: Solar cooking, solar thermal applications, solar geometry, tilted solar panels, energy storage systems, surplus product storage systems, cloud information, temperature, wind, other meteorological factors, and supporting information. This latest release contains new parameters based on recommendations by the renewable energy industry and it is more accurate than previous releases. On-line plotting capabilities allow quick evaluation of potential renewable energy projects for any region of the world. The SSE data set is formulated from NASA satellite- and reanalysis-derived insolation and meteorological data for the 10-year period July 1983 through June 1993. Results are provided for 1 degree latitude by 1 degree longitude grid cells over the globe. Average daily and monthly measurements for 1195 World Radiation Data Centre ground sites are also available. [Mission Objectives] The SSE project contains insolation and meteorology data intended to aid in the development of renewable energy systems. Collaboration between SSE and technology industries such as the Hybrid Optimization Model for Electric Renewables ( HOMER ) may aid in designing electric power systems that employ some combination of wind turbines, photovoltaic panels, or diesel generators to produce electricity. [Temporal_Coverage: Start_Date=1983-07-01; Stop_Date=1993-06-30] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180].

Stackhouse, Paul W.↗

NASA NTRS2010

Optimal Codes for the Burst Erasure Channel

Deep space communications over noisy channels lead to certain packets that are not decodable. These packets leave gaps, or bursts of erasures, in the data stream. Burst erasure correcting codes overcome this problem. These are forward erasure correcting codes that allow one to recover the missing gaps of data. Much of the recent work on this topic concentrated on Low-Density Parity-Check (LDPC) codes. These are more complicated to encode and decode than Single Parity Check (SPC) codes or Reed-Solomon (RS) codes, and so far have not been able to achieve the theoretical limit for burst erasure protection. A block interleaved maximum distance separable (MDS) code (e.g., an SPC or RS code) offers near-optimal burst erasure protection, in the sense that no other scheme of equal total transmission length and code rate could improve the guaranteed correctible burst erasure length by more than one symbol. The optimality does not depend on the length of the code, i.e., a short MDS code block interleaved to a given length would perform as well as a longer MDS code interleaved to the same overall length. As a result, this approach offers lower decoding complexity with better burst erasure protection compared to other recent designs for the burst erasure channel (e.g., LDPC codes). A limitation of the design is its lack of robustness to channels that have impairments other than burst erasures (e.g., additive white Gaussian noise), making its application best suited for correcting data erasures in layers above the physical layer. The efficiency of a burst erasure code is the length of its burst erasure correction capability divided by the theoretical upper limit on this length. The inefficiency is one minus the efficiency. The illustration compares the inefficiency of interleaved RS codes to Quasi-Cyclic (QC) LDPC codes, Euclidean Geometry (EG) LDPC codes, extended Irregular Repeat Accumulate (eIRA) codes, array codes, and random LDPC codes previously proposed for burst erasure protection. As can be seen, the simple interleaved RS codes have substantially lower inefficiency over a wide range of transmission lengths.

Hamkins, Jon↗

NASA NTRS2015

Evaluation of Variable-Depth Liner Configurations for Increased Broadband Noise Reduction

This paper explores the effects of variable-depth geometry on the amount of noise reduction that can be achieved with acoustic liners. Results for two variable-depth liners tested in the NASA Langley Grazing Flow Impedance Tube demonstrate significant broadband noise reduction. An impedance prediction model is combined with two propagation codes to predict corresponding sound pressure level profiles over the length of the Grazing Flow Impedance Tube. The comparison of measured and predicted sound pressure level profiles is sufficiently favorable to support use of these tools for investigation of a number of proposed variable-depth liner configurations. Predicted sound pressure level profiles for these proposed configurations reveal a number of interesting features. Liner orientation clearly affects the sound pressure level profile over the length of the liner, but the effect on the total attenuation is less pronounced. The axial extent of attenuation at an individual frequency continues well beyond the location where the liner depth is optimally tuned to the quarter-wavelength of that frequency. The sound pressure level profile is significantly affected by the way in which variable-depth segments are distributed over the length of the liner. Given the broadband noise reduction capability for these liner configurations, further development of impedance prediction models and propagation codes specifically tuned for this application is warranted.

Jones, M. G.↗

NASA NTRS2014

Guided Wave Delamination Detection and Quantification With Wavefield Data Analysis

Unexpected damage can occur in aerospace composites due to impact events or material stress during off-nominal loading events. In particular, laminated composites are susceptible to delamination damage due to weak transverse tensile and inter-laminar shear strengths. Developments of reliable and quantitative techniques to detect delamination damage in laminated composites are imperative for safe and functional optimally-designed next-generation composite structures. In this paper, we investigate guided wave interactions with delamination damage and develop quantification algorithms by using wavefield data analysis. The trapped guided waves in the delamination region are observed from the wavefield data and further quantitatively interpreted by using different wavenumber analysis methods. The frequency-wavenumber representation of the wavefield shows that new wavenumbers are present and correlate to trapped waves in the damage region. These new wavenumbers are used to detect and quantify the delamination damage through the wavenumber analysis, which can show how the wavenumber changes as a function of wave propagation distance. The location and spatial duration of the new wavenumbers can be identified, providing a useful means not only for detecting the presence of delamination damage but also allowing for estimation of the delamination size. Our method has been applied to detect and quantify real delamination damage with complex geometry (grown using a quasi-static indentation technique). The detection and quantification results show the location, size, and shape of the delamination damage.

Tian, Zhenhua↗

NASA NTRS1997

Validation of an Accurate Three-Dimensional Helical Slow-Wave Circuit Model

The helical slow-wave circuit embodies a helical coil of rectangular tape supported in a metal barrel by dielectric support rods. Although the helix slow-wave circuit remains the mainstay of the traveling-wave tube (TWT) industry because of its exceptionally wide bandwidth, a full helical circuit, without significant dimensional approximations, has not been successfully modeled until now. Numerous attempts have been made to analyze the helical slow-wave circuit so that the performance could be accurately predicted without actually building it, but because of its complex geometry, many geometrical approximations became necessary rendering the previous models inaccurate. In the course of this research it has been demonstrated that using the simulation code, MAFIA, the helical structure can be modeled with actual tape width and thickness, dielectric support rod geometry and materials. To demonstrate the accuracy of the MAFIA model, the cold-test parameters including dispersion, on-axis interaction impedance and attenuation have been calculated for several helical TWT slow-wave circuits with a variety of support rod geometries including rectangular and T-shaped rods, as well as various support rod materials including isotropic, anisotropic and partially metal coated dielectrics. Compared with experimentally measured results, the agreement is excellent. With the accuracy of the MAFIA helical model validated, the code was used to investigate several conventional geometric approximations in an attempt to obtain the most computationally efficient model. Several simplifications were made to a standard model including replacing the helical tape with filaments, and replacing rectangular support rods with shapes conforming to the cylindrical coordinate system with effective permittivity. The approximate models are compared with the standard model in terms of cold-test characteristics and computational time. The model was also used to determine the sensitivity of various circuit parameters including typical manufacturing dimensional tolerances and support rod permittivity. By varying the circuit parameters of an accurate model using MAFIA, these sensitivities can be computed for manufacturing concerns, and design optimization previous to fabrication, thus eliminating the need for costly experimental iterations. Several variations were made to a standard helical circuit using MAFIA to investigate the effect that variations on helical tape and support rod width, metallized loading height and support rod permittivity, have on TWT cold-test characteristics.

Kory, Carol L.↗

NASA NTRS2013

JWST Integrated Science Instrument Module Alignment Optimization Tool

During cryogenic vacuum testing of the James Webb Space Telescope (JWST) Integrated Science Instrument Module (ISIM), the global alignment of the ISIM with respect to the designed interface of the JWST optical telescope element (OTE) will be measured through a series of optical characterization tests. These tests will determine the locations and orientations of the JWST science instrument projected focal surfaces and entrance pupils with respect to their corresponding OTE optical interfaces. If any optical performance non-compliances are identified, the ISIM will be adjusted to improve its performance. In order to understand how to manipulate the ISIM's degrees of freedom properly and to prepare for the ISIM flight model testing, a series of optical-mechanical analyses have been completed to develop and identify the best approaches for bringing a non-compliant ISIM element into compliance. In order for JWST to meet its observatory-level optical requirements and ambitious science goals, the ISIM element has to meet approximately 150 separate optical requirements. Successfully achieving many of those optical requirements depends on the proper alignment of the ISIM element with respect to the OTE. To verify that the ISIM element will meet its optical requirements, a series of cryogenic vacuum tests will be conducted with an OTE Simulator (OSIM). An optical Ray Trace and Geometry Model tool was developed to help solve the multi-dimensional alignment problem. The tool allows the user to determine how best to adjust the alignment of the JWST ISIM with respect to the ideal telescope interfaces so that the approximately 150 ISIM optical performance requirements can be satisfied. This capability has not existed previously.

Bos, Brent↗

NASA NTRS1993

Parallel algorithm for dominant points correspondences in robot binocular stereo vision

This paper presents an algorithm to find the correspondences of points representing dominant feature in robot stereo vision. The algorithm consists of two main steps: dominant point extraction and dominant point matching. In the feature extraction phase, the algorithm utilizes the widely used Moravec Interest Operator and two other operators: the Prewitt Operator and a new operator called Gradient Angle Variance Operator. The Interest Operator in the Moravec algorithm was used to exclude featureless areas and simple edges which are oriented in the vertical, horizontal, and two diagonals. It was incorrectly detecting points on edges which are not on the four main directions (vertical, horizontal, and two diagonals). The new algorithm uses the Prewitt operator to exclude featureless areas, so that the Interest Operator is applied only on the edges to exclude simple edges and to leave interesting points. This modification speeds-up the extraction process by approximately 5 times. The Gradient Angle Variance (GAV), an operator which calculates the variance of the gradient angle in a window around the point under concern, is then applied on the interesting points to exclude the redundant ones and leave the actual dominant ones. The matching phase is performed after the extraction of the dominant points in both stereo images. The matching starts with dominant points in the left image and does a local search, looking for corresponding dominant points in the right image. The search is geometrically constrained the epipolar line of the parallel-axes stereo geometry and the maximum disparity of the application environment. If one dominant point in the right image lies in the search areas, then it is the corresponding point of the reference dominant point in the left image. A parameter provided by the GAV is thresholded and used as a rough similarity measure to select the corresponding dominant point if there is more than one point the search area. The correlation is used as a final decision tool when there is still more than one point in the search area. If there is no dominant point in the search area of if the points in the search area are below a correlation threshold, then the dominant point in the reference image is occluded and can not be corresponded. The algorithm has been modeled, implemented and shown to be fast, robust and parallel. The parallelism is created from three main features: locality of the operators; a memory optimization scheme; and the ability to fully parallelize the extraction phase which is the most computational intensive task in the algorithm. The last feature is achieved by performing the extraction phase on the two images simultaneously.

Al-Tammami, A.↗

NASA NTRS2026

Micromechanical Design of Carbon Nanotube Ribbon Reinforced Polymer Composite Materials

Lightweight materials are an important component of the design of aerospace structures. Carbon nanotube materials have been considered for this purpose due to the strength and stiffness of individual nanotubes, and the commercial availability of bulk formats such as fibers. These fibers can have a ribbon cross section which results in a different design space for their composites relative to traditional reinforcements which have a round cross section. This work applies brick-and-mortar micromechanical models and classical lamination theory with an inverse approach to investigate the design space of these composites. Using this approach, the influences of fiber geometry and axial and transverse mechanical properties are mapped. Finally, a sensitivity study is performed and the relative impacts of ±10% variations in the constituent material and geometric properties are ranked. Lamina axial moduli were found to range from a maximum of 3x to 1x minimum relative to a target quasi-isotropic laminate modulus depending on the anisotropy and shear modulus of the lamina. The fiber targets depended strongly on the fiber volume fraction in the lamina and the fiber axial modulus target was found to range from 4.8x to 3.2x the quasi-isotropic laminate target. The sensitivity analysis found that the largest driver of performance was the volume fraction, followed by the fiber axial modulus. While bio-based brick-and-mortar composites, such as nacre, can benefit from reinforcement aspect ratios above 10, for carbon nanotube ribbon (or carbon fiber)/polymer composites the sensitivity study indicated that the optimal cross-sectional aspect ratio was relatively smaller, potentially less than three.

Micromechanics↗

NASA NTRS2011

Laser Truss Sensor for Segmented Telescope Phasing

A paper describes the laser truss sensor (LTS) for detecting piston motion between two adjacent telescope segment edges. LTS is formed by two point-to-point laser metrology gauges in a crossed geometry. A high-resolution (<30 nm) LTS can be implemented with existing laser metrology gauges. The distance change between the reference plane and the target plane is measured as a function of the phase change between the reference and target beams. To ease the bandwidth requirements for phase detection electronics (or phase meter), homodyne or heterodyne detection techniques have been used. The phase of the target beam also changes with the refractive index of air, which changes with the air pressure, temperature, and humidity. This error can be minimized by enclosing the metrology beams in baffles. For longer-term (weeks) tracking at the micron level accuracy, the same gauge can be operated in the absolute metrology mode with an accuracy of microns; to implement absolute metrology, two laser frequencies will be used on the same gauge. Absolute metrology using heterodyne laser gauges is a demonstrated technology. Complexity of laser source fiber distribution can be optimized using the range-gated metrology (RGM) approach.

Liu, Duncan T.↗

NASA NTRS2013

Integration of Mirror Design with Suspension System using NASA's New Mirror Modeling Software

Advances in mirror fabrication is making very large space based telescopes possible. In the many applications, only monolithic mirrors meet the performance requirements. The existing and near-term planned heavy launch vehicles place a premium on lowest possible mass. Again, available and planned payload shroud size limits near term designs to 4 meter class mirror. Practical 8 meter and beyond designs could encourage planners to include larger shrouds if it can be proven that such mirrors can be manufactured. These two factors lower mass and larger mirrors, presents the classic optimization problem. There is a practical upper limit to how large a mirror can be supported by a purely kinematic mount system and be launched. This paper shows how the design of the suspension system and mirror blank needs to be designed simultaneously. We will also explore the concepts of auxiliary support systems, which act only during launch and disengage on orbit. We will define required characteristics of these systems and show how they can substantially reduce the mirror mass. The AMTD project is developing and maturing the processes for future replacements for HUBBLE, creating the design tools, validating the methods and techniques necessary to manufacture, test and launch extremely large optical missions. This paper will use the AMTD 4 meter "design point" as an illustration of the typical use of the modeler in generating the multiple models of mirror and suspension systems used during the conceptual design phase of most projects. The influence of Hexapod geometry, mirror depth, cell size and construction techniques (Exelsis Deep Core Low Temperature Fusion (c) versus Corning Frit Bonded (c) versus Schott Pocket Milled Zerodur (c) in this particular study) are being evaluated. Due to space and time consideration we will only be able to present snippets of the study in this paper. The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond.

Arnold, William↗

NASA NTRSDate not supplied

System Modeling of a Lunar Molten Regolith Electrolysis Plant

Introduction: In-Situ Resource Utilization (ISRU) is the process of extracting local resources to produce commodities for propulsion, life support systems, and off-planet construction rather than transporting consumables from Earth. Molten Regolith Electrolysis (MRE) is a novel ISRU method of extracting oxygen gas and metal alloy from lunar regolith. The MRE process involves placing lunar regolith between two electrodes, through which current is passed, to melt the regolith and reduce the metal oxide constituents by direct electrolysis (e.g. FeO, SiO2, MgO, Al2O3) into oxygen gas and metal alloys. The oxygen is liquefied and used as propellant for landers, while the metals (e.g. Ferro-alloys) are further processed and used in structural building materials and parts manufacturing. A system model was developed that accounted for the major processes of an MRE plant (from excavation of raw materials to storage of products) to assess the feasibility of a lunar MRE plant. The System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team utilized its previously documented system sizing model, the Mission Analysis and Integration Tool (MAIT) [1] as framework of the system model. MAIT uses MATLAB/Simulink to integrate subsystem models into a complete system model of the MRE plant. Total mass, volume, and power requirements were computed for numerous iterations of a MRE plant. System Model: Figure 1: MRE Plant Block Diagram The regolith excavation model determines the mass and power needed to excavate sufficient regolith. The preheating auger initiates the regolith heating process before regolith enters the MRE re-actor to reduce the energy required to turn the solid into a molten liquid. The MRE reactor is modeled in COMSOL Multiphysics and based on the research by Dominguez, Sibille, and Schreiner [2, 3, 4]. This preliminary reactor model provides an accurate calculation of thermal equilibrium during electrochemical operation of the reactor system to assess the optimal mass and power required to process the inlet flow of regolith. The model also computes the outlet flowrates of oxygen and molten products. For this analysis, the primary components of the metal alloy considered were iron and silicon. The oxygen is then purified using an Yttrium Stabilized Zirconia (YSZ) electrode, followed by liquefaction using a 90K cryocooler to be stored as liquid oxygen in insulated cylindrical tanks. In future iterations of the system model, the molten metal tapped from the MRE reactor will undergo additional processing or refinement. However, downstream handling of metals is currently a technology gap that is missing a high TRL subsystem model. Therefore, for this analysis, the accumulated metal alloy stream terminates after leaving the MRE reactor. Study Goals: This analysis investigates multiple input variables to the system to determine the sensitivity of a (near) complete plant at full-scale. This preliminary investigation ran parametric sweeps on the MRE reactor geometry, electrical current supply, layers of multi-layer insulation (MLI) on the reactor, size of the electrodes in the oxygen purification model, and regolith composition (based on landing site location). Three production targets of oxygen (1,000, 10,000, and 50,000 kg/yr) were investigated for this analysis. The parametric sweeps conducted in this analysis provide valuable insight into the expected impact of the various model inputs on plant size. This information can be used to identify the most critical components of the plant and guide future decisions on allocating funding for research and development, providing subsystem developers with appropriate interfaces with downstream and upstream processes, and assessing the overall feasibility of MRE when compared to other ISRU plants. References: [1] Carlson, A. et al. (2024) ICES, ICES-2024-53. [2] Dominguez, D.A., and Sibille, L. (2011) AIAA, AIAA-2011-700. [3] Schreiner, S.S. (2015) MIT, Dissertation. [4] Schreiner, S.S. et al. (2016) ASR, 57(7), pp.1585-1603.

ISRU↗

NASA NTRS2008

Joint Test Plan for Gas Dynamic Spray Technology Demonstration

Air Force Space Command (AFSPC) and NASA have similar missions, facilities, and structures located in similar harsh environments. Both are responsible for a number of facilities/structures with metallic structural and non-structural components in highly and moderately corrosive environments. Regardless of the corrosivity of the environment, all metals require periodic maintenance activity to guard against the insidious effects of corrosion and thus ensure that structures meet or exceed design or performance life. The standard practice for protecting metallic substrates in atmospheric environments is the use of an applied coating system. Current coating systems used across AFSPC and NASA contain volatile organic compounds (VOCs) and hazardous air pollutants (HAPs). These coatings are subject to environmental regulations at the Federal and State levels that limit their usage. In addition, these coatings often cannot withstand the high temperatures and exhaust that may be experienced by AFSPC and NASA structures. In response to these concerns, AFSPC and NASA have approved the use of thermal spray coatings (TSCs). Thermal spray coatings are extremely durable and environmentally friendly coating alternatives, but utilize large cumbersome equipment for application that make the coatings difficult and time consuming to repair. Other concerns include difficulties coating complex geometries and the cost of equipment, training, and materials. Gas Dynamic Spray (GDS) technology (also known as Cold Spray) will be evaluated as a smaller, more maneuverable repair method as well as for areas where thermal spray techniques are not as effective. The technology can result in reduced maintenance and thus reduced hazardous materials/wastes associated with current processes. Thermal spray and GDS coatings also have no VOCs and are environmentally preferable coatings. To achieve a condition suitable for the application of a coating system, including GDS coatings, the substrate must undergo some type of surface preparation and/or depainting operation to ensure adhesion of the new coating system. The GDS unit selected for demonstration has a powder feeding system that can be used for surface preparation or coating application. The surface preparation feature will also be examined. The primary objective of this effort is to demonstrate GDS technology as a repair method for TSCs. The project will also determine the optimal GDS coating thickness for acceptable performance. Successful completion of this project will result in approval of GDS technology as a repair method for TSCs at AFSPC and NASA installations and will improve corrosion protection at critical systems, facilitate easier maintenance activity, extend maintenance cycles, eliminate flight hardware contamination, and reduce the amount of hazardous waste generated.

Lewis, Pattie↗

NASA NTRS2012

Extension of the Time-Spectral Approach to Overset Solvers for Arbitrary Motion

Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.

Leffell, Joshua Isaac↗

NASA NTRSDate not supplied

Simulations of a Turbulent Flow Subjected to Favorable and Adverse Pressure Gradients

This paper reports the results from a direct numerical simulation of an initially turbulent boundary layer passing over a wall-mounted “speed bump” geometry. The speed bump, represented in the form of a Gaussian distribution profile, generates a favorable pressure gradient region over the upstream half of the geometry, followed by an adverse pressure gradient over the downstream half. The boundary layer approaching the bump undergoes strong acceleration in the favorable pressure gradient region before experiencing incipient or very weak separation within the adverse pressure gradient region. These types of flows have proven to be particularly challenging to predict using lower-fidelity simulation tools based on various turbulence modeling approaches and warrant the use of the highest-fidelity simulation techniques. Simulation results are utilized to examine the key phenomena present in the flowfield, such as relaminarization/stabilization in the strong acceleration region succeeded by retransition to turbulence near the onset of adverse pressure gradient, incipient/weak separation, and development of internal layers where the sense of streamwise pressure gradient changes at the foot, apex and tail of the bump. The present direct numerical simulation is performed using a flow solver developed exclusively for graphics processing units, which is found to provide a significant speedup compared to an earlier solver optimized for central processing unit architectures.

Ali Uzun↗