Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix completion”

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.

557 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Techniques Optimized for Reducing Instabilities in Advanced Nickel-Base Superalloys for Turbine Blades

The High-Speed Research (HSR) Airfoil Alloy program developed fourth-generation single-crystal superalloys with up to an 85 F increase in creep rupture capability over current production airfoil alloys. Recent results have been generated at the NASA Glenn Research Center on these fourth-generation alloys, but in coated form, for subsonic turbine blade applications under NASA's Ultra-Efficient Engine Technology (UEET) Program. One goal for UEET is to optimize the airfoil alloy/thermal barrier coating system for 3100 F turbine inlet temperatures. The state-of-the art turbine blade airfoil system consists of a superalloy single crystal that provides the basic mechanical performance of the airfoil. A thermal barrier coating is used to reduce the temperature of the base superalloy, and a bondcoat is deposited between the base material and the thermal barrier coating. The bondcoat improves the oxidation and corrosion resistance of the base superalloy and improves the spallation resistance of the thermal barrier coating. A commercial platinum aluminide bondcoat was applied to the HSR-developed alloys, and a diffusion zone developed as a result of interaction between the bondcoat and the superalloy. Optimized strength is obtained for superalloys when the refractory element content is high and the limits of microstructural stability are approached or exceeded slightly. For fourthgeneration alloys, instability leads to the formation of topologically close packed (TCP) phases, which form internally in the superalloy, and a secondary reaction zone (SRZ), which forms under the diffusion zone. There was a concern that excessive quantities of either TCP or SRZ might decrease the mechanical properties of the superalloy, with SRZ thought to be particularly detrimental and its formation unpredictable. Thus, an SRZreduction effort was initiated in the NASA UEET Program so that methods developed during the HSR project could be optimized further to reduce or eliminate the SRZ. An SRZ is a three-phase constituent composed of TCP and stringers of gamma phase in a matrix of gamma prime. An incoherent grain boundary separates the SRZ from the gammagamma prime microstructure of the superalloy. The SRZ is believed to form as a result of local chemistry changes in the superalloy due to the application of the diffusion aluminide bondcoat. Locally high surface stresses also appear to promote the formation of the SRZ. Thus, techniques that change the local alloy chemistry or reduce surface stresses have been examined for their effectiveness in reducing SRZ. These SRZ-reduction steps are performed on the test specimen or the turbine blade before the bondcoat is applied. Stressrelief heat treatments developed at NASA Glenn have been demonstrated to reduce significantly the amount of SRZ that develops during subsequent high-temperature exposures. Stress-relief heat treatments reduce surface stresses by recrystallizing a thin surface layer of the superalloy. However, in alloys with very high propensities to form SRZ, stress relief heat treatments alone do not eliminate SRZ entirely. Thus, techniques that modify the local chemistry under the bondcoat have been emphasized and optimized successfully at Glenn. One such technique is carburization, which changes the local chemistry by forming submicron carbides near the surface of the superalloy. Detailed characterizations have demonstrated that the depth and uniform distribution of these carbides are enhanced when a stress relief treatment and an appropriate surface preparation are employed in advance of the carburization treatment. Even in alloys that have the propensity to develop a continuous SRZ layer beneath the diffusion zone, the SRZ has been completely eliminated or reduced to low, manageable levels when this combination of techniques is utilized. Now that the techniques to mitigate SRZ have been established at Glenn, TCP phase formation is being emphasized in ongoing work under the UEET Program. The limitsf stability of the fourth-generation alloys with respect to TCP phase formation are currently being defined along with high-temperature creep rupture properties. In addition, a regression model is being developed at Glenn for the prediction of the presence of TCP phase in the microstructure and SRZ under the diffusion zone. The model is based on a design-of-experiments methodology with emphasis on the potential synergistic effects of alloying elements.

MacKay, Rebecca A.↗
NASA NTRS2003

Characterization of a Ultra-high Temperature Ceramic Composite

Ultra-high temperature ceramics (UHTC) are of interest for hypersonic vehicle leading edge applications. Monolithic UHTCs are of concern because of their low fracture toughness and brittle behavior. UHTC composites (UHTCC) are being investigated as a possible approach to overcome these deficiencies. In this study a small sample of a UHTCC was evaluated by limited mechanical property tests, furnace oxidation exposures, and oxidation exposures in a flowing environment. The composite was prepared from a carbon fiber perform using ceramic particulates and a preceramic polymer. The as-received composite plate was non-uniform from front to back surface. Plate dimensions were 150 x 150 x 6 mm. The back surface had a fibrous, uniform appearance; XRD analysis revealed the presence of Sic and C. The front surface was smooth and non-uniform in appearance with evidence of a coarse grain structure produced by a liquid phase; XRD analysis revealed the presence of HfB2. Microcracks were present throughout the thickness as one might expect from a carbon fiber reinforced composite with attendant large thermal expansion mismatch between the matrix phases and the fibers. The HfB2 phase on the front surface was comparable in thickness to a fiber ply or about 0.6 mm, and surface microcracks were evident. Limited four point flexural tests were carried out at span to depth ratios of approximately 14 and 16 with markedly different results. Tests were run with the front or the back surface in tension. At the shorter span to depth failures occurred under a loading pin for both orientations. At a span to depth of 16 failures occurred in the center of the span with fracture clearly initiating from a tensile failure. Ultimate flexural strength, strain at ultimate stress, stress and strain at deviation from linear elastic behavior are reported. Strains at ultimate stress ranged from about 0.6 to 0.7 % for the back surface in tension, and 0.4 to 0.6 for the front surface in tension. At constant span to depth the strain at ultimate stress was about 0.2% greater for the back surface in tension and the ultimate strength was also higher. Strengths were in line with predictions from theory. Furnace oxidation studies were carried out at 1627 and 1927OC in a static furnace environment using ten minute cycles and one, five, and ten cycles. Limited oxidation studies were also carried out in a flowing oxyacetylene torch environment. Specimens were photographed, and weight and dimensional changes were determined. XRD and SEM characterizations were performed. Weight losses were attributed primarily to carbon fiber oxidation. The composite survived the torch test with little visible distress. Further details will be determined once metallographic studies are completed.

Levine, Stanley R.↗
NASA NTRS2011

Should Pruning be a Pre-Processor of any Linear System?

There are many real-world problems whose mathematical models turn out to be linear systems Ax = b, where A is an m x n matrix. Each equation of the linear system is an information. An information, in a physical problem, such as 4 mangoes, 6 bananas, and 5 oranges cost $10, is mathematically modeled as an equation 4x(sub 1) + 6x(sub 2) + 5x(sub 3) = 10 , where x(sub 1), x(sub 2), x(sub 3) are each cost of one mango, that of one banana, and that of one orange, respectively. All the information put together in a specified context, constitutes the physical problem and need not be all distinct. Some of these could be redundant, which cannot be readily identified by inspection. The resulting mathematical model will thus have equations corresponding to this redundant information and hence are linearly dependent and thus superfluous. Consequently, these equations once identified should be better pruned in the process of solving the system. The benefits are (i) less computation and hence less error and consequently a better quality of solution and (ii) reduced storage requirements. In literature, the pruning concept is not in vogue so far although it is most desirable. It is assumed that at least one information, i.e. one equation is known to be correct and which will be our first equation. In a numerical linear system, the system could be slightly inconsistent or inconsistent of varying degree. If the system is too inconsistent, then we should fall back on to the physical problem (PP), check the correctness of the PP derived from the material universe, modify it, if necessary, and then check the corresponding mathematical model (MM) and correct it. In nature/material universe, inconsistency is completely nonexistent. If the MM becomes inconsistent, it could be due to error introduced by the concerned measuring device and/or due to assumptions made on the PP to obtain an MM which is relatively easily solvable or simply due to human error. No measuring device can usually measure a quantity with an accuracy greater that 0.005% or, equivalently with a relative error less than 0.005%. Hence measurement error is unavoidable in a numerical linear system when the quantities are continuous (or even discrete with extremely large number). Assumptions, though not desirable, are usually made when we find the problem sufficiently difficult to be solved within the available means/tools/resources and hence distort the PP and the corresponding MM. The . error thus introduced in the system could (not always necessarily though) make the system somewhat inconsistent. If the inconsistency (contradiction) is too much then one should definitely not proceed to solve the system in terms of getting a least-squares solution or the minimum-norm least-squares solution. All these solutions will be invariably of no real-world use. If, on the other hand, inconsistency is reasonably low, i.e. the system is near-consistent or, equivalently, has near-linearly-dependent rows, then the foregoing solutions are useful. Pruning in such a near-consistent system should be performed based on the desired accuracy and on the definition of near-linear dependence. In this article, we discuss pruning over various kinds of linear systems and strongly suggest its use as a pre-processor or as a part of an algorithm. Ideally pruning should (i) be a part of the solution process (algorithm) of the system, (ii) reduce both computational error and complexity of the process, and (iii) take into account the numerical zero defined in the context. These are precisely what we achieve through our proposed O(mn2) algorithm presented in Matlab, that uses a subprogram of solving a single linear equation and that has embedded in it the pruning.

Sen, Syamal K.↗
NASA NTRS2024

Hybrid Thermally Efficient Core (HyTEC) HyTEC Phase 1 – Combustor Final Report

The objective of the HyTEC – Combustor Technology project is to develop technology for a compact, low emissions, rich-burn combustor that maintains a high-level of durability. To accomplish this, the combustor will incorporate CMC liners and a CMC Dome to enable improved air utilization for mixing and NOx reduction. In particular, three technology areas were matured: metallic swirler attachment to the CMC dome, Nickel diffusion in CMC, and EBC advancements with improved durability. The maturation of these technology areas is critical to successfully designing a combustor for an engine intended to demonstrate increased thermal efficiency with integrated high-power density-core engine technologies. As part of the project, multiple swirler attachment architectures were matured through the design process, manufactured, and tested in simulated engine conditions. These architectures utilized two different approaches to attachment: a pressure loaded approach and a clamped approach. Each approach required testing at simulated engine conditions to mature the technology. These tests were developed to ensure that the probable failure mechanisms for each architecture were assessed relative to its long-term durability. After successfully testing both approaches, the pressure loaded approach was chosen based on its simplicity and its similarity to the attachment approach used on lean burn combustors. Finally, this successful testing resulted in the swirler attachment technology achieving TRL 4 status. Another aspect of concern with the metallic swirler attachment is the potential for nickel diffusion into the CMC to cause durability shortfalls. A series of tests were completed to understand both the amount of nickel that is expected to diffuse into the CMC and the effect this diffusion will have on mechanical properties of the CMC. After successfully completing this testing, the reduction in local mechanical properties was found to not have an impact on the overall durability of the CMC Dome. To further improve the overall durability, three different EBC coating architectures were also studied. One of these coatings being the current state of the art for CMC liners with the other two being proposed improvements. Coupons were produced using each coating and subjected to a series of laboratory testing. One coating was found to have an overall superior durability in laboratory testing and survived rig testing without indication of degradation. This successful testing resulted in the EBC technology achieving TRL 5 status. Overall, a design approach was successfully down-selected that achieves both the TPM requirements and the need for TRL 4+ status. This design incorporates a pressure loaded swirler attachment with an improved EBC coating (Architecture 1). This combination allows the continued use of a CMC dome and results in an improved durability for the proposed compact core combustor.

CMC↗
NASA NTRS1995

The glycan keratan sulfate in inner ear crystals

The otoconial matrix (OM) of chicks (Gallus domesticus) inner ear was analyzed. Histochemically the OM was reacted with phosphotungstic acid (PTA) and immunohistochemically with the monoclonal antibody antikeratan sulfate (antiKS). The OM was digested with the enzyme endo-beta-galactosidase (E beta Galase) or separated by 1D and 2D gel electrophoresis. PTA which reacts with glycoproteins precipitated the OM, suggesting that the OM contains glycoproteins. A central core in each crystal had no PTA staining, suggesting that the core lacked glycoproteins. Anti KS antibody stained the OM with increased density in older embryos as determined by color thresholding. E beta Galase, which cleaves the lactosamine repeating units in KS, decreased the immunostain by 30% in the OM and by 20% in the cartilage. The OM from the utricle, saccule and macula lagena contained similar molecular weight bands. Five dense bands in the OM were less dense in tissue and blood controls, suggesting that such bands are enriched in the OM. Isoelectric focusing of the OM showed a negatively charged high molecular weight smear not present in blood and faint in tissue controls. The high affinity of the OM for the cationic PTA stain, the strong immunohistochemical reaction of the OM with anti KS antibody and high molecular weight negative smear in 2D gels taken together suggest that: a) the OM contains large amounts of glycoproteins and glycans, one of which is keratan sulfate, because its immuno stain with antiKS antibody was decreased by the enzyme E beta Galase, b) the utricle, saccule and macula lagena may have similar composition, and c) the concentration of KS may increase gradually until complete mineralization of the OM is reached.

Non-NASA Center↗
NASA NTRS1997

Pilot Non-Conformance to Alerting System Commands

This research project examined the effects of consonance between cockpit displays and alerting system as a technique to encourage pilots to conform to alerting system commands. An experiment used the task of collision avoidance during closely spaced parallel approaches as a case study, building upon previous experiments which identified instances of non-conformance and conflicts between the alerting criteria preferred by pilots, compared to that used by alerting systems. Using a workstation based, part-task simulator, each of 45 subjects completed 45 experiment runs. In each run, the subjects were told they were flying an approach. Their primary task was to keep their wings level despite turbulence through the use of a sidestick. The sidestick commands did not affect the path of the aircraft, however, so that consistent approach paths were be followed. Their secondary task was to indicate when an aircraft on a parallel approach is blundering towards them, as evidenced by the traffic display. Subjects were asked to press different buttons indicating whether they feel an avoidance maneuver is required by the traffic situation or not. At the completion of each run, subjects were asked to rate their confidence in their decision and, if appropriate, to rate the timeliness of automatic alerts when had been given. Three different automatic alert conditions were tested. The "No Automatic Alerts Given" condition is self-explanatory. In the "Automatic Alerts Based on NTZ Criteria" condition, an automatic alert was given when the NTZ criteria was triggered; this criteria is consistent with subject reactions in other studies, in which subjects were found to react, on average, when the other aircraft was 1350 min to the side of the own aircraft. In the "Automatic Alerts Based on MIT Criteria" condition, an automatic alert was given when the MIT criteria was triggered; this criteria was developed by Carpenter and Kuchar for parallel approaches to have better performance, at the cost of increased complexity and higher sensitivity to pilot non-conformance. Three displays were tested. The "Baseline Display" is based on a moving map display, with a top-down view, track-up orientation, iconic presentation of the other aircraft's positions and a text presentation of the other aircraft's altitude. The "NTZ Alert Criteria Shown" display added to the Baseline Display an explicit presentation on the display of the boundary delineating the NTZ criteria, allowing for quick comparison of the other aircraft's position relative to this alert criteria. Likewise, the "MIT Alert Criteria Shown" display added to the Baseline Display an explicit presentation on the display of the boundary delineating the MIT criteria, allowing for quick comparison of the other aircraft's position relative to this alert criteria The test matrix for this experiment was three dimensional, varying displays, alerts and traffic conflict scenarios.

Pritchett, Amy↗
NASA NTRS2017

Constraints on Mechanisms of Chondrule Formation from Chondrule Precursors and Chronology of Transient Heating Events in the Protoplanetary Disk

The mineralogy, petrography, and oxygen-isotope compositions of porphyritic chondrules—dominant chondrule type in most chondrite groups—suggest formation by incomplete melting of isotopically diverse precursors during localized transient heating events in dust-rich regions of the protoplanetary disk characterized by 16O-poor compositions (D17Odust+gas~ –7‰ to +4‰) relative to the inferred Sun’s value (D17O ~ –28 ± 2‰). The chondrule precursors included Ca,Al-rich inclusions (CAIs), amoeboid olivine aggregates (AOAs), chondrules of earlier generations, fine-grained matrix-like material, and possibly fragments of pre-existing planetesimals. Like porphyritic chondrules, igneous CAIs formed by melting of isotopically diverse precursors during transient heating events, but in an isotopically distinct, solar-like reservoir of the protoplanetary disk (D17O dust + gas ~ –24‰), probably near the proto Sun. Based on a narrow range of the initial 26Al/27Alratios inferred from the internal Al-Mg isochrons in igneous CAIs, their melting started at the very beginning of Solar System formation (t0), defined by the CV CAIs with U-corrected Pb-Pb age of 4567.3 ± 0.16 Ma and the canonical 26Al/27Al ratio of (5.25 ± 0.02) ¥ 10–5, and lasted at least 0.3 Ma. The U-corrected Pb-Pb absolute and 26Al-26Mg relative ages of porphyritic chondrules from type 3 ordinary, CO, CV, and CR carbonaceous chondrites (assuming uniform distribution of 26Al in the disk at the canonical level) suggest chondrule formation started at t0 and lasted for about 4 Ma. These observations may preclude formation of the majority of porphyritic chondrules by splashing of differentiated planetesimals and by collisions between planetesimals; instead, they are consistent with melting of dust balls by bow shocks or magnetized turbulence in the disk. Some porphyritic chondrules in equilibrated (petrologic type 4–6) ordinary chondrites contain relict fragments of coarse-grained chromite, ilmenite, phosphates, and albitic plagioclase. The similar mineral assemblage is commonly observed in type 4–6 ordinary chondrites, but is absent in type 3 chondrites, suggesting these chondrules formed by incomplete melting of thermally metamorphosed ordinary chondrite material, possibly by impacts. The CB metal-rich carbonaceous chondrites contain exclusively magnesian non-porphyritic chondrules crystallized from complete melts. These chondrules formed in a gas-melt plume generated by a hypervelocity (≥20 km/s) collision between planetesimals ~4.8 Ma after t0 in a transition or a debris disk. One of the colliding bodies was probably differentiated. The CH metal-rich carbonaceous chondrites contain chondrules formed by different mechanisms. The magnesian non-porphyritic chondrules formed in the CB impact plume ~4.8 Ma after t0. The chemically diverse (magnesian, ferroan, and Al-rich) porphyritic chondrules formed by incomplete melting of isotopically diverse precursors in the protoplanetary disk, most likely prior the CB impact plume event. We conclude that there are multiple mechanisms of chondrule formation that operated over the entire life-time of the disk.

Chondrules↗
NASA NTRS1997

Large Scale Turbulent Structures in Supersonic Jets

Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations(DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those, of a spatially evolving jet, a temporal jet problem was solved, using periodicity ill the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible application to active noise suppression. In addition, the data generated can be used to compute various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved Direct Numerical Simulation (DNS) of a turbulent jet flow is produced.

Rao, Ram Mohan↗
NASA NTRS1997

Large Scale Turbulent Structures in Supersonic Jets

Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations (DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those of a spatially evolving jet, a temporal jet problem was solved, using periodicity in the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible appli(,a- tion to active noise suppression. In addition, the data generated can be used to compute, various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved DNS of a turbulent jet flow is produced.

Rao, Ram Mohan↗
NASA NTRS1993

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.↗
NASA NTRS2011

Method for Impeding Degradation of Porous Silicon Structures

This invention relates to surface modification of porosified silicon (pSi) structures with poly(alkylene) glycols for the purpose of controlled degradation of the silicon matrix and tailored release of encapsulated substances for biomedical applications. The pSi structures are currently used in diverse biomedical applications including bio-molecular screening, optical bio-sensoring, and drug delivery by means of injectable/orally administered carriers and implantable devices. The size of the pores and the surface chemistry of the pSi structure can be controlled during the microfabrication process and thereafter. A fine regulation of the degradation kinetics of mesoporous silicon structures is of fundamental importance. Polyethylene glycols (PEGs) represent the major category of surface modifying agents used in classical drug delivery systems and in pharmaceutical dosage forms. PEGylation enables avoidance of RES uptake, thus prolonging circulation time of intravenously injectable nanovectors. PEG molecules demonstrate little toxicity and immunogenicity, and are cleared from the body through the urine (molecular weight, MW less than 30 kDa) or in the feces (MW greater than 30kDa). The invention focuses on the possibility of finely tuning the degradation kinetics of the pSi nanovectors and other structures through surface conjugation of PEGs with various backbone lengths/MWs. To prove the concept, pSi nanovectors were covalently conjugated to seven PEGs with MW from 245 to 5,000 Da and their degradation kinetics in physiologically relevant media (phosphate buffer saline, PBS pH7.4, and fetal bovine serum) was assessed by the elemental analysis of the Si using inductive coupled plasma atomic emission spectroscopy (ICP-AES). The conjugation of the PEG with lowest MW to the nanovectors surface did not induce any change in the degradation kinetics in serum, but inhibited degradation and consequently the release of orthosilicic acid into buffer. When PEGs with the longer chains were evaluated, Si mass loss from the nanovectors was slowed down, and the PEGylated structures were almost fully degraded within 18 24 hours in serum and within 48 hours in PBS. The most dramatic effect was observed for high MW PEGs 3,400 and 5,000 Da, which prominently inhibited the degradation of the systems, with complete degradation achieved only after four days. For these PEGs, during the early stages of the degradation, there was a lag period of little or no Si mass loss from the nanovector. The obtained profiles were in agreement with the erosion of the nanovector surface as observed by scanning electron microscopy.

Vilentchouk, Biana Godin↗
NASA NTRS2007

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↗
NASA NTRS1997

High Temperature Degradation Mechanisms in Polymer Matrix Composites

Polymer matrix composites are increasingly used in demanding structural applications in which they may be exposed to harsh environments. The durability of such materials is a major concern, potentially limiting both the integrity of the structures and their useful lifetimes. The goal of the current investigation is to develop a mechanism-based model of the chemical degradation which occurs, such that given the external chemical environment and temperatures throughout the laminate, laminate geometry, and ply and/or constituent material properties, we can calculate the concentration of diffusing substances and extent of chemical degradation as functions of time and position throughout the laminate. This objective is met through the development and use of analytical models, coupled to an analysis-driven experimental program which offers both quantitative and qualitative information on the degradation mechanism. Preliminary analyses using coupled diffusion/reaction model are used to gain insight into the physics of the degradation mechanisms and to identify crucial material parameters. An experimental program is defined based on the results of the preliminary analysis which allows the determination of the necessary material coefficients. Thermogravimetric analyses are carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. Powdered samples are used to eliminate diffusion effects. Tests in both inert and oxidative environments allow the separation of thermal and oxidative contributions to specimen mass loss. The concentration dependency of the oxidative reactions is determined from the tests in pure oxygen. Short term isothermal tests at different temperatures are carried out on neat resin and unidirectional macroscopic specimens to identify diffusion effects. Mass loss, specimen shrinkage, the formation of degraded surface layers and surface cracking are recorded as functions of exposure time. Geometry effects in the neat resin, and anisotropic diffusion effects in the composites, are identified through the use of specimens with different aspect ratios. The data is used with the model to determine reaction coefficients and effective diffusion coefficients. The empirical and analytical correlations confirm the preliminary model results which suggest that mass loss at lower temperatures is dominated by oxidative reactions and that these reaction are limited by diffusion of oxygen from the surface. The mechanism-based model is able to successfully capture the basic physics of the degradation phenomena under a wide range of test conditions. The analysis-based test design is successful in separating out oxidative, thermal, and diffusion effects to allow the determination of material coefficients. This success confirms the basic picture of the process; however, a more complete understanding of some aspects of the physics are required before truly predictive capability can be achieved.

Cunningham, Ronan A.↗
NASA NTRS1996

DCTune Perceptual Optimization of Compressed Dental X-Rays

In current dental practice, x-rays of completed dental work are often sent to the insurer for verification. It is faster and cheaper to transmit instead digital scans of the x-rays. Further economies result if the images are sent in compressed form. DCTune is a technology for optimizing DCT (digital communication technology) quantization matrices to yield maximum perceptual quality for a given bit-rate, or minimum bit-rate for a given perceptual quality. Perceptual optimization of DCT color quantization matrices. In addition, the technology provides a means of setting the perceptual quality of compressed imagery in a systematic way. The purpose of this research was, with respect to dental x-rays, 1) to verify the advantage of DCTune over standard JPEG (Joint Photographic Experts Group), 2) to verify the quality control feature of DCTune, and 3) to discover regularities in the optimized matrices of a set of images. We optimized matrices for a total of 20 images at two resolutions (150 and 300 dpi) and four bit-rates (0.25, 0.5, 0.75, 1.0 bits/pixel), and examined structural regularities in the resulting matrices. We also conducted psychophysical studies (1) to discover the DCTune quality level at which the images became 'visually lossless,' and (2) to rate the relative quality of DCTune and standard JPEG images at various bitrates. Results include: (1) At both resolutions, DCTune quality is a linear function of bit-rate. (2) DCTune quantization matrices for all images at all bitrates and resolutions are modeled well by an inverse Gaussian, with parameters of amplitude and width. (3) As bit-rate is varied, optimal values of both amplitude and width covary in an approximately linear fashion. (4) Both amplitude and width vary in systematic and orderly fashion with either bit-rate or DCTune quality; simple mathematical functions serve to describe these relationships. (5) In going from 150 to 300 dpi, amplitude parameters are substantially lower and widths larger at corresponding bit-rates or qualities. (6) Visually lossless compression occurs at a DCTune quality value of about 1. (7) At 0.25 bits/pixel, comparative ratings give DCTune a substantial advantage over standard JPEG. As visually lossless bit-rates are approached, this advantage of necessity diminishes. We have concluded that DCTune optimized quantization matrices provide better visual quality than standard JPEG. Meaningful quality levels may be specified by means of the DCTune metric. Optimized matrices are very similar across the class of dental x-rays, suggesting the possibility of a 'class-optimal' matrix. DCTune technology appears to provide some value in the context of compressed dental x-rays.

Watson, Andrew B.↗
NASA NTRS1996

High Temperature Degradation Mechanisms in Polymer Matrix Composites

Polymer matrix composites are increasingly used in demanding structural applications in which they may be exposed to harsh environments. The durability of such materials is a major concern, potentially limiting both the integrity of the structures and their useful lifetimes. The goal of the current investigation is to develop a mechanism-based model of the chemical degradation which occurs, such that given the external chemical environment and temperatures throughout the laminate, laminate geometry, and ply and/or constituent material properties, we can calculate the concentration of diffusing substances and extent of chemical degradation as functions of time and position throughout the laminate. This objective is met through the development and use of analytical models, coupled to an analysis-driven experimental program which offers both quantitative and qualitative information on the degradation mechanism. Preliminary analyses using a coupled diffusion/reaction model are used to gain insight into the physics of the degradation mechanisms and to identify crucial material parameters. An experimental program is defined based on the results of the preliminary analysis which allows the determination of the necessary material coefficients. Thermogravimetric analyses are carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. Powdered samples are used to eliminate diffusion effects. Tests in both inert and oxidative environments allow the separation of thermal and oxidative contributions to specimen mass loss. The concentration dependency of the oxidative reactions is determined from the tests in pure oxygen. Short term isothermal tests at different temperatures are carried out on neat resin and unidirectional macroscopic specimens to identify diffusion effects. Mass loss, specimen shrinkage, the formation of degraded surface layers and surface cracking are recorded as functions of exposure time. Geometry effects in the neat resin, and anisotropic diffusion effects in the composites, are identified through the use of specimens with different aspect ratios. The data is used with the model to determine reaction coefficients and effective diffusion coefficients. The empirical and analytical correlations confirm the preliminary model results which suggest that mass loss at lower temperatures is dominated by oxidative reactions and that these reaction are limited by diffusion of oxygen from the surface. The mechanism-based model is able to successfully capture the basic physics of the degradation phenomena under a wide range of test conditions. The analysis-based test design is successful in separating out oxidative, thermal, and diffusion effects to allow the determination of material coefficients. This success confirms the basic picture of the process; however, a more complete understanding of some aspects of the physics are required before truly predictive capability can be achieved.

Cunningham, Ronan A.↗
NASA NTRSDate not supplied

An Optimization Approach to Support Science Decision Making for Lunar Surface Exploration

Introduction: Scientific exploration is one of the three pillars of NASA’s Moon2Mars architecture, with crew surface extra vehicular activities (EVA) serving a critical enabling function. Development of surface EVA operational planning and execution, specifically integrating science and flight control teams (FCT), is currently being explored through analog scenarios. This integration, exercised, for example, through the Joint EVA and Hu-man Surface Mobility Test Team (JETT), allows for science input on EVA activities in near real-time through a Science Evaluation Room (SER), or Arte-mis science backroom, which integrates with the broader FCT through the Science Officer. The SER works within the FCT to support dynamic EVA planning in response to changes in operational constraints as well as science opportunities and re-prioritization, increasing the mission science return and accelerating the accomplishment of the Moon2Mars science objectives. The SER works within the FCT to provide recommendations to traverse execution in near real-time. One challenge is the requirement to deliver SER inputs to the FCT on operationally relevant timelines. Failure to do so may result in suboptimal execution of science exploration EVAs or even loss of key science objectives. To close this gap, we present a network optimization tool to allow the SER to provide rapid input to the FCT in response to changes in operational constraints or science opportunities. Inputs are predicated on approved science objectives, and clear rationale must be provided to the FCT for any requested change. Accordingly, this tool incorporates the Science Traceability Matrix (STM), SER prioritization scheme, and station characterization and action planning with operational constraints such as duration, traverse speed, and distance to maximize science objectives based on SER priorities, consistent with FCT operational requirements. Method: As a proof of concept, we used an existing linear programing software package used to simulate optimal routes through cellular metabolism. We built a Demonstrative Model with three STM objectives and four stations on a region of the Moon. The objectives were given an arbitrary prioritization and mapped to the stations through four possible crew actions. (Figs. 1 and 2). This station to STM mapping is consistent with the method used by the JETT5 Science Team to develop analog surface EVA science planning. We used a grid system with the landing site at the origin and the four stations placed across the positive x,y quadrant. Actions were assigned to each station and the accomplishment of those actions resulted in a numerical “reward” based on the ability of that action to achieve science objectives. The aggregate reward from each individual STM objective contributes to a global score (Science Yield), weighted by its priority. Operational constraints included a requirement to start and end at the landing site, 5 minutes each for initial station characterization and “clean up,” and variable total EVA time, traverse rate (fixed to 0.5 meters per second in our example), and time to perform each action (10, 5, 7, and 15 min for actions 1, 2, 3, and 4, respectively). Additional constraints and variables will be added in the future (e.g., sample mass, number of stations, traverse route constraints, illumination). Optimization. We converted the connections (arcs) between these stations (nodes) into a mixed integer linear programming optimization problem (arcs = constraints, nodes = variables) with the objective to maximize Science Yield. For any action, the Science Yield is equal to the relevance of that action to an STM objective [3, 2, and 1 point(s) for High, Med., and Low relevance, respectively], multiplied by the STM Objective Priority [3, 2, and 1 point(s) for High, Med., and Low priority, respectively]. This resulted in a model that computes the optimal station and action combination to maximize the Science Yield. These weightings can be adjusted by the SER as desired. Results: We explored three test cases for the Demonstrative Model. First, we set the maximum EVA duration to 120 minutes and computed the optimal route (Fig. 3A). The model suggested per-forming Actions 1 and 2 at Station P01, followed by Actions 1 and 2 at Station P02, and finally Actions 1 and 3 at Station P04 before returning to the Landing Site. Second, we adjusted the STM Objective Priori-ty order and computed the new optimal route (Fig. 3B). Under this situation, the model suggested per-forming all Actions at Station P02 followed by all Actions at Station P03. The previous test cases were relevant to SER planning activities. Next, we explored providing mid-EVA replanning input to the FCT. Scenario: While executing the Route in Fig. 3A the crew finishes at Station P01 and FCT decides that the EVA needs to finish in 45 minutes back at the Landing Site. FCT asks SER to recommend changes to the plan to accommodate this operation-al change. Using the model and incorporating these new constraints (start at Station P01, max. time of 45 min), the model suggested performing Actions 2 and 4 at Station P03 (Fig. 4), requiring 41 minutes to complete and return to the Landing Site. Interestingly, Station 3 was not part of the original route. Using the model, we determined the EVA would need 66 minutes, instead of 45, in order for the original Station P04 to yield a larger Science Yield than Station P03. The parametrization and simulation was per-formed in less than a minute, demonstrating the operational relevance of the approach. Future Efforts: The results from the Demonstrative Model suggest this tool can accelerate SER decision making on operationally relevant timelines. Use in analog activities, such as JETT5 or follow-ons, which have over a dozen stations for a crew to explore and over a dozen actions per station, will provide needed validation of the utility of this tool for planning EVAs, replanning mid-EVA, or planning follow-on EVAs based on previous results. Further integration with FCT execution monitoring tools may provide additional efficiency gains, al-lowing rapid and iterative exploration of operation-al and science decision space by the FCT and SER.

Science Operations↗
NASA NTRSDate not supplied

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗