Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Matrix”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 541 records · Page 30

Micrometer to Atomic Scale Characterisation of Primitive Astromaterials Using A Novel Method, Metis-Fa: A Coordinated Atom Probe Tomography, Transmission Electron Microscopy and NanoSIMS Approach

Introduction: Presolar grains preserve isotopic, chemical and microstructural records of physical and chemical processing, and formation mechanisms within a vast range of evolved stellar systems, the interstellar medium, solar nebula and their parent bodies. These evolutionary records are preserved at the micrometric to atomic scale, requiring coordinated studies to expand our understanding of evolutionary processes occurringthroughout ours and external stellar systems [1]. NanoSIMS enabled rapid in situ identification and isotopic characterisation of presolar grains and their stellar origins using 17O/16O and 18O/16O, and 13C/12C isotopic ratios [1]. Coordination with transmission electron microscopy (TEM) revealed crystallographic and localised contextual relationships and quantitively constrained their major and minor compositions [1]. However, trace elements cannot be quantified, the most sensitive geochemical tracers of environmental conditions, essential to unravelling the chemical record of their evolutionary pathway and parent stellar systems [2-3] . Furthermore, owing to the combination of technical limitations (only 5 – 7 isotopes can be measured per NanoSIMS run) and their small grain sizes of 100 nm < 3 μm (with rare exceptions in nanodiamonds (2 nm ≤) and SiC (< 40 μm)), the number of measurable isotopes per grain volume is limited [1,3] . Through more comprehensive isotopic studies of presolar grains, NanoSIMS studies have shown the importance of the latter, identifying Fe and Mg as important indicators of nuclear synthetic processing and their stellar origins, respectively [4- 5]. Coordination of NanoSIMS and Atom Probe Tomography (APT) revealed morphological signatures, and isotopic and chemical signatures at major to trace levels without requirements for preselection of elements [6]. However, crystallographic signatures in localized contextual relationships cannot be measured. Consequently, coordination of NanoSIMS, TEM and APT is essential to gain access to almost all contextual, structural and geochemical signatures within each presolar grain.Transmission electron microscopy requires a 100 nm thin lamella which is unstable in APT and would not produce any viable data. Atom probe tomography requires a needle-shaped specimen which when measured in TEM removes the local context, impacts the quality of the TEM diffraction images due to the shank angle of the needle, and can alter the chemistry of beam sensitive materials from the higher degree of surface exposure at the tip. To address these issues, we developed METIS-Fa (Multi-technical measurements of Electron Transparent materials using an Indium Sandwich - a FIB approach). A novel method which enables coordination of NanoSIMS, TEM and APT for generalized and targeted studies of individual grains, including beam sensitive materials, without compromising sample preparation requirements for TEM and APT. This method requires only indium and a Focus Ion Beam (FIB), minimizing the movement of fragile materials while still enabling preparation of TEM lamella into APT needles. Samples: Initial experimental development and testing of the method occurred at Astromaterials Research and Exploration Science (ARES), Johnson Space Centre (JSC), NASA and APT measurements and needle preparation occurred at JdLC, Curtin University. Synthetic silicate samples were used as analogs for presolar silicates when performing a trial run of the method. Samples were extracted from a polished thin section created at JSC, NASA, comprised of 38 wt.% Si, 17 wt.% FeO, 13 wt.% MgO, 12 wt.% Al, 11 wt.% Ca based on electron microprobe analysis (EMPA) [8] . Experimental details, pressure and temperature conditions were presented in [8] and references therein. Testing of the capability to target individual grains in mineral matrices using this method for acquisition in APT, measured matrix regions in meteoritic thin sections of primitive meteorites. These meteorites and their identified presolar grains for future targeted studies are detailed in [9]. Techniques: The TEM-FIB lamella were prepared using a FIB. An e-beam assisted pt deposition was used as a protective coating for the synthetic and meteoritic samples. When targeting individual grains, a secondary e-beam assisted pt deposition button is placed over the desired grain before the protective coating to denote its location. A JEOL 2500SE field-emission TEM was used for high-resolution imaging, energy-dispersive X-ray (EDX) and electron diffraction data.TheMETIS-Fa method was experimentally designed, tested and executed using a FIB at ARES, JSCNASA. Needles for APT were prepared using the Tescan Lyra3 GM Dual Beam Focus Ion Beam (FIB) Field Emission SEM (FE-SEM) at the JdLC, Curtin University. Atom probe tomography measurements were conducted using a CAMECA Local Electrode Atom Probe, LEAP 4000X HR. Two pure indium needles were analyzed initially to constraining acquisition parameters and stability under the beam. Manual acquisition was required to maintain evaporation of specimen’s at the apex and monitor interactions with measurement parameters. Experimental Design: Indium foil is pressed onto an Al stub with a pneumatic press and mounted into the FIB adjacent to the TEM-FIB lamella of interest. Using a FIB, two indium slices (5 μm x ~300 nm x 3 μm) are extracted from indium foil and aligned with the TEM-FIB lamella before touching the TEM-FIB lamella. Each slice is then attached through cold welding to the FIB-TEM lamella. This approach eliminates the need for chemical treatments and proved effective for aligning the Indium within the region of interest for APT, holding it in place for up to 4 days during testing.Once both indium slices are attached within their pre-determined region per grain targeting requirements, they are gradually melted onto the FIB-TEM lamella.When targeting a specific grain, measurements should be taken of the pt button and its distance from edge to edge of the lamella before and after sandwiching. A secondary button should be placed over the same region after the Indium slices have been attached to improve precision when preparing APT needles. Results: Figure 1 shows two indium slices melted onto a FIB-TEM lamella, adding additional bulk for preparation into APT needles as shown in Figure 2 [7] . The latter was essential so samples could be measured in TEM and APT without compromising sample preparation requirements and consequently data quality and acquisition stability. METIS-Fa proved effective forimproving geometry. Figure 3 shows a successful APTrun of the synthetic silicate. EMPA, TEM and APTshowed no chemical alterations. During targetingtesting, a solar silicate grain was successfully identifiedand measured in TEM, and prepared into an APTneedle. However, the indium was melted too long during sample preparation, causing expansion andformation of internal porosity leading to sample loss.Conclusion: METIS-Fa greatly expands the number of isotopic and chemical signatures measured per grain volume, and enables measurements of contextual, structural, crystallographic, isotopic and geochemical signatures within individual grains. Gaining access to such a vast range of evolutionary signatures required for expanding our understanding of external stellar and planetary systems and the evolution of our solar system. This method was designed for application to a vast range of phases including being sensitive materials and thus provides a way for coordination of NanoSIMS, TEM and APT not just for the study of presolar grains and by extension primitive astromaterials, but studies in a vast range of other fields including the geosciences and material sciences.Acknowledgments: Thankyou to ARES, JSC, NASA; JdLC Curtin University and Space Science Technology Centre for the use of laboratory facilities and funding [confirm].

Nicole D Nevill↗

Digital halftoning methods for selectively partitioning error into achromatic and chromatic channels

A method is described for reducing the visibility of artifacts arising in the display of quantized color images on CRT displays. The method is based on the differential spatial sensitivity of the human visual system to chromatic and achromatic modulations. Because the visual system has the highest spatial and temporal acuity for the luminance component of an image, a technique which will reduce luminance artifacts at the expense of introducing high-frequency chromatic errors is sought. A method based on controlling the correlations between the quantization errors in the individual phosphor images is explored. The luminance component is greatest when the phosphor errors are positively correlated, and is minimized when the phosphor errors are negatively correlated. The greatest effect of the correlation is obtained when the intensity quantization step sizes of the individual phosphors have equal luminances. For the ordered dither algorithm, a version of the method can be implemented by simply inverting the matrix of thresholds for one of the color components.

Mulligan, Jeffrey B.↗

Probabilistic sizing of laminates with uncertainties

A reliability based design methodology for laminate sizing and configuration for a special case of composite structures is described. The methodology combines probabilistic composite mechanics with probabilistic structural analysis. The uncertainties of constituent materials (fiber and matrix) to predict macroscopic behavior are simulated using probabilistic theory. Uncertainties in the degradation of composite material properties are included in this design methodology. A multi-factor interaction equation is used to evaluate load and environment dependent degradation of the composite material properties at the micromechanics level. The methodology is integrated into a computer code IPACS (Integrated Probabilistic Assessment of Composite Structures). Versatility of this design approach is demonstrated by performing a multi-level probabilistic analysis to size the laminates for design structural reliability of random type structures. The results show that laminate configurations can be selected to improve the structural reliability from three failures in 1000, to no failures in one million. Results also show that the laminates with the highest reliability are the least sensitive to the loading conditions.

Shah, A. R.↗

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↗

Application and development of the Schwinger multichannel scattering theory and the partial differential equation theory of electron-molecule scattering

One version of the multichannel theory for electron-target scattering based on the Schwinger variational principle, the SMC method, requires the introduction of a projection parameter. The role of the projection parameter a is investigated and it is shown that the principal-value operator in the SMC equation is Hermitian regardless of the value of a as long as it is real and nonzero. In a basis that is properly orthonormalizable, the matrix representation of this operator is also Hermitian. The use of such basis is consistent with the Schwinger variational principle because the Lippmann-Schwinger equation automatically builds in the correct boundary conditions. Otherwise, an auxiliary condition needs to be introduced, and Takatsuka and McKoy's original value of a is one of the three possible ways to achieve Hermiticity. In all cases but one, a can be uncoupled from the Hermiticity condition and becomes a free parameter. An equation for a based on the variational stability of the scattering amplitude is derived; its solution has an interesting property that the scattering amplitude from a converged SMC calculation is independent of the choice of a even though the SMC operator itself is a-dependent. This property provides a sensitive test of the convergence of the calculation. For a static-exchange calculation, the convergence requirement only depends on the completeness of the one-electron basis, but for a general multichannel case, the a-invariance in the scattering amplitude requires both the one-electron basis and the N plus 1-electron basis to be complete. The role of a in the SMC equation and the convergence property are illustrated using two examples: e-CO elastic scattering in the static-exchange approximation, and a two-state treatment of the e-H2 Chi(sup 1)Sigma(sub g)(+) yields b(sup 3)Sigma(sub u)(+) excitation.

Weatherford, Charles A.↗

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↗

Influence of Thermally-Grown Oxide (TGO) Layer on the Driving Forces Associated with Failure in Environmental Barrier Coating (EBC) Systems

Environmental barrier coatings (EBC) is an enabling technology for the successful application of ceramic matrix composites (CMCs) in air-breathing gas turbine engines. Spallation of environmental barrier coating (EBC) induced by thermally grown oxide (TGO) layer is a key EBC failure mode. The TGO layer, resulting from steam oxidation, grows either from a silicon bond coat layer (if present) or from the silicon carbide (SiC) based substrate itself. Critical thickness of the TGO layer for failure is in the range of 20-30 microns but it can vary due to exposure temperature, microstructure etc. Current work at NASA Glenn Research Center, under the Revolutionary Tools and Methods (RTM) project is aimed at addressing associated failure modes in EBC systems and developing robust analysis tools to aid in the design/analysis of these systems. The objective of the current work is to conduct a sensitivity study to examine the influence of uniformly and non-uniformly grown oxide layers with or without damage on the associated driving forces leading to spallation of the EBC when subjected to isothermal loading. Initial results indicate that the presence of damage (vertical cracks caused by in-plane stresses) enhances the stresses that are present due to non-uniformity. However, the presence of non-uniformity itself is still the main factor influencing the magnitude of peel and shear stresses in the TGO layer.

Mital, Subodh K.↗

An efficient design sensitivity analysis of eigenvectors

Subspace iteration has been a major advance in solving large eigen problems when only a subset of eigen-pairs is required. The essence of this method is a transformation from displacement coordinates of an n-th order eigensystem to generalized coordinates of a smaller q-th order. The eigenvalue problem is then solved in the reduced space. The method was first developed by Clint and Jennings for real symmetric systems and was then called 'simultaneous iteration'. The success of the method prompted further research along this line and there have been many improved algorithms developed. This approach has been widely used by structural engineers for extracting the most useful natural frequencies and mode shapes of large-scale dynamic systems. This paper exploits into a new direction which is in the form of iterative process for simultaneously calculating eigenvector derivatives of many eigenvectors with respect to multi-variables. The method fully uses all the available information from preceding eigenvalue solution and, thus, effectively economizes computational efforts. It iterates through two equations derived from the first variation of the two fundamental equations used in subspace iteration method. There is no expensive large matrix decomposition required and the process converges to acceptable solution in a finite number of iterations. Therefore, the procedure increases its efficiency superiority over the others as the system size or the number of interested eigenvectors become larger and larger.

Ting, T.↗

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↗

Vibration localization in mono- and bi-coupled bladed disks - A transfer matrix approach

A transfer matrix approach to the analysis of the dynamics of mistuned bladed disks is presented. The study focuses on mono-coupled systems, in which each blade is coupled to its two neighboring blades, and bi-coupled systems, where each blade is coupled to its four nearest neighbors. Transfer matrices yield the free dynamics, both the characteristic free wave and the normal modes - in closed form for the tuned assemblies. Mistuned assemblies are represented by random transfer matrices and an examination of the effect of mistuning on harmonic wave propagation yields the localization factor - the average rate of spatial wave amplitude decay per blade - in the mono-coupled assembly. Based on a comparison of the wave propagation characteristics of the mono- and bi-coupled assemblies, important conclusions are drawn about the effect of the additional coupling coordinate on the sensitivity to mistuning and the strength of mode localization predicted by a mono-coupled analysis.

Ottarsson, Gisli↗

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↗

Interpretation of satellite-measured bidirectional reflectance from Cirrus cloudy atmospheres

The interpretation of the observed bidirectional reflectance from cirrus cloudy atmospheres is presented. A theoretical model was developed for the computation of the transfer of solar radiation in an anisotropic medium with particular applications to oriented ice crystals in cirrus clouds. In this model, the adding principle for radiative transfer was used with modifications to account for the anisotropy of scattering particles and the associated scattering phase matrix. The single-scattering properties, including the phase function, single-scattering albedo, and extinction cross sections, were used for randomly and horizontally oriented hexagonal ice crystals in radiative transfer computations. The radiative transfer model developed for the cirrus clouds was modified to account for the scattering contributions from the atmosphere and the surface. In order to test the relevance and significance of the ice crystal model for the interpretation of observed bidirectional reflectance from satellites, visible radiances collected on the half hour by the GOES series were selected. The data are calibrated and corrected with the proper filter functions, then navigated to match selected landmark data. A number of clear and cloudy cases during the cirrus IFO of the Fire experiment were chosen for theoretical analyses. The cloud particle shape and size distributions that were taken during satellite overpasses are used in radiative transfer calculations. The sensitivities of the shape, orientation, and size distribution of ice crystals on the reflected intensities at the top of the atmosphere are investigated. Finally, the relative importance of these cloud microphysical properties in the interpretation of satellite bidirectional reflectance are assessed and presented.

Takano, Y.↗

The Density Matrix of H20 - N2 In the Coordinate Representation: A Monte Carlo Calculation of the Far-Wing Line Shape

The far-wing line shape theory within the binary collision and quasistatic framework has been developed using the coordinate representation. Within this formalism, the main computational task is the evaluation of multidimensional integrals whose variables are the orientational angles needed to specify the initial and final positions of the system during transition processes. Using standard methods, one is able to evaluate the 7-dimensional integrations required for linear molecular systems, or the 7-dimensional integrations for more complicated asymmetric-top (or symmetric-top) molecular systems whose interaction potential contains cyclic coordinates. In order to obviate this latter restriction on the form of the interaction potential, a Monte Carlo method is used to evaluate the 9-dimensional integrations required for systems consisting of one asymmetric-top (or symmetric-top) and one linear molecule, such as H20-N2. Combined with techniques developed previously to deal with sophisticated potential models, one is able to implement realistic potentials for these systems and derive accurate, converged results for the far-wing line shapes and the corresponding absorption coefficients. Conversely, comparison of the far-wing absorption with experimental data can serve as a sensitive diagnostic tool in order to obtain detailed information on the short-range anisotropic dependence of interaction potentials.

Ma, Q.↗

Characterization of a New Phase and Its Effect on the Work Characteristics of a Near-Stoichiometric Ni30Pt20Ti50 High-Temperature Shape Memory Alloy (HTSMA)

A new phase observed in a nominal Ni30Pt20Ti50 (at.%) high temperature shape memory alloy has been characterized using transmission electron microscopy and 3-D atom probe tomography. This phase forms homogeneously in the B2 austenite matrix by a nucleation and growth mechanism and results in a concomitant increase in the martensitic transformation temperature of the base alloy. Although the structure of this phase typically contains a high density of faults making characterization difficult, it appears to be trigonal (-3m point group) with a(sub o) approx. 1.28 nm and c(sub o) approx. 1.4 nm. Precipitation of this phase increases the microhardness of the alloy substantially over that of the solution treated and quenched single-phase material. The effect of precipitation strengthening on the work characteristics of the alloy has been explored through load-biased strain-temperature testing in the solution-treated condition and after aging at 500 C for times ranging from 1 to 256 hours. Work output was found to increase in the aged alloy as a result of an increase in transformation strain, but was not very sensitive to aging time. The amount of permanent deformation that occurred during thermal cycling under load was small but increased with increasing aging time and stress. Nevertheless, the dimensional stability of the alloy at short aging times (1-4 hours) was still very good making it a potentially useful material for high-temperature actuator applications.

Garg, A.↗

Closed-form analysis of fiber-matrix interface stresses under thermo-mechanical loadings

Closed form techniques for calculating fiber matrix (FM) interface stresses, using repeating square and diamond regular arrays, were presented for a unidirectional composite under thermo-mechanical loadings. An Airy's stress function micromechanics approach from the literature, developed for calculating overall composite moduli, was extended in the present study to compute FM interface stresses for a unidirectional graphite/epoxy (AS4/3501-6) composite under thermal, longitudinal, transverse, transverse shear, and longitudinal shear loadings. Comparison with finite element results indicate excellent agreement of the FM interface stresses for the square array. Under thermal and longitudinal loading, the square array has the same FM peak stresses as the diamond array. The square array predicted higher stress concentrations under transverse normal and longitudinal shear loadings than the diamond array. Under transverse shear loading, the square array had a higher stress concentration while the diamond array had a higher radial stress concentration. Stress concentration factors under transverse shear and longitudinal shear loadings were very sensitive to fiber volume fraction. The present analysis provides a simple way to calculate accurate FM interface stresses for both the square and diamond array configurations.

Naik, Rajiv A.↗

Processing and Properties of Fiber Reinforced Polymeric Matrix Composites: Processing Robustness of IM7/PETI Polyimide Composites - Part 2

The processability of a phenylethynyl terminated imide (PETI) resin matrix composite was investigated. Unidirectional prepregs were made by coating an N-methylpyrrolidone solution of the amide acid oligomer onto unsized IM7. Two batches of prepregs were used: one was made by NASA in-house, and the other was from an industrial source. The composite processing robustness was investigated with respect to the effect of B-staging conditions, the prepreg shelf life, and the optimal processing window. Rheological measurements indicated that PETI's processability was only slightly affected over a wide range of B-staging temperatures (from 250 C to 300 C). The open hole compression (OHC) strength values were statistically indistinguishable among specimens consolidated using various B-staging conditions. Prepreg rheology and OHC strengths were also found not to be affected by prolonged (i.e., up to 60 days) ambient storage. An optimal processing window was established using response surface methodology. It was found that IM7/PETI composite is more sensitive to the consolidation temperature than to the consolidation pressure. A good consolidation was achievable at 371 C/100 Psi, which yielded an OHC strength of 62 Ksi at room temperature. However, processability declined dramatically at temperatures below 350 C.

Hou, Tan-Hung↗

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗

Development of Three-Dimensional Multicellular Tissue-Like Constructs for Mutational Analysis Using Macroporous Microcarriers

A three-dimensional (3-D), tissue-like model was developed for the genotoxic assessment of space environment. In previous experiments, we found that culturing mammalian cells in a NASA-designed bioreactor, using Cytodex-3 beads as a scaffold, generated 3-D multicellular spheroids. In an effort to generate scaffold-free spheroids, we developed a new 3-D tissue-like model by coculturing fibroblast and epithelial cell in a NASA bioreactor using macroporous Cultispher-S(TradeMark) microcarriers. Big Blue(Registered Trademark) Rat 2(Lambda) fibroblasts, genetically engineered to contain multiple copies (>60 copies/cell) of the Lac I target gene, were cocultured with radio-sensitive human epithelial cells, H184F5. Over an 8-day period, samples were periodically examined by microscopy and histology to confirm cell attachment, growth, and viability. Immunohistochemistry and western analysis were used to evaluate the expression of specific cytoskeletal and adhesion proteins. Key cell culture parameters (glucose, pH, and lactate concentrations) were monitored daily. Controls were two-dimensional mono layers of fibroblast or epithelial cells cultured in T-flasks. Analysis of 3-D spheroids from the bioreactor suggests fibroblast cells attached to and completely covered the bead surface and inner channels by day 3 in the bioreactor. Treatment of the 3-day spheroids with dispase II dissolved the Cultisphers(TradeMark) and produced multicellular, bead-less constructs. Immunohistochemistry confirmed the presence of vi.mentin, cytokeratin and E-cadherin in treated spheroids. Examination of the dispase II treated spheroids with transmission electron microscopy (TEM) also showed the presence of desmosomes. These results suggest that the controlled enzymatic degradation of an artificial matrix in the low shear environment of the NASA-designed bioreactor can produce 3-D tissue-like spheroids. 2

Jordan, Jacqueline A.↗