SEARCH · Engineering Papers

Results for “Classical mechanics”

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

SourceResearch connectionFrom yearThrough year

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

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Experimental Identification and Simulation of Time and/or Rate Dependent Reversible and Irreversible Deformation Regions for both a Titanium and Nickel Alloy

In this paper time and/or rate dependent deformation regions are experimentally mapped out as a function of temperature. It is clearly demonstrated that the concept of a threshold stress (a stress that delineate reversible and irreversible behavior) is valid and necessary at elevated temperatures and corresponds to the classical yield stress at lower temperatures. Also the infinitely slow modulus, (Es) i.e. the elastic modulus of the material if it was loaded at an infinitely slow strain rate, and the "dynamic modulus", modulus, Ed, which represents the modulus of the material if it is loaded at an infinitely fast rate are used to delineate rate dependent from rate independent regions. As demonstrated at elevated temperatures there is a significant difference between the two modulus values, thus indicating both significant time-dependence and rate dependence. In the case of the nickel-based super alloy, ME3, this behavior is also shown to be grain size specific. Consequently, at higher temperatures viscoelastic behavior exist below k (i.e., the threshold stress) and at stresses above k the behavior is viscoplastic. Finally a multi-mechanism, stress partitioned viscoelastic model, capable of being consistently coupled to a viscoplastic model is characterized over the full temperature range investigated for Ti-6-4 and ME3.

Arnold, Steven M.↗

NASA NTRS2022

Electroactive ZnO: Mechanisms, Conductivity, and Advances in Zn Alkaline Battery Cycling

Zinc oxide is of great interest for advanced energy devices because of its low cost, wide direct bandgap, non-toxicity, and facile electrochemistry. In zinc alkaline batteries, ZnO plays a critical role in electrode passivation, a process that hinders commercialization and remains poorly understood. Here, novel observations of an electroactive type of ZnO formed in Zn-metal alkaline electrodes are disclosed. The electrical conductivity of battery-formed ZnO is measured and found to vary by factors of up to 104, which provides a first-principles-based understanding of Zn passivation in industrial alkaline batteries. Simultaneous with this conductivity change, protons are inserted into the crystal structure and electrons are inserted into the conduction band in quantities up to ≈1020 cm−3 and ≈1 mAh gZnO−1. Electron insertion causes blue electrochromic coloration with efficiencies and rates competitive with leading electrochromic materials. The electroactivity of ZnO is evidently enabled by rapid crystal growth, which forms defects that complex with inserted cations, charge-balanced by the increase of conduction band electrons. This property distinguishes electroactive ZnO from inactive classical ZnO. Knowledge of this phenomenon is applied to improve cycling performance of industrial-design electrodes at 50% zinc utilization and the authors propose other uses for ZnO such as electrochromic devices.

Brendan E. Hawkins↗

NASA NTRS1990

Shape optimization and CAD

Structural optimization has attracted the attention since the days of Galileo. Olhoff and Taylor have produced an excellent overview of the classical research within this field. However, the interest in structural optimization has increased greatly during the last decade due to the advent of reliable general numerical analysis methods and the computer power necessary to use them efficiently. This has created the possibility of developing general numerical systems for shape optimization. Several authors, eg., Esping; Braibant & Fleury; Bennet & Botkin; Botkin, Yang, and Bennet; and Stanton have published practical and successful applications of general optimization systems. Ding and Homlein have produced extensive overviews of available systems. Furthermore, a number of commercial optimization systems based on well-established finite element codes have been introduced. Systems like ANSYS, IDEAS, OASIS, and NISAOPT are widely known examples. In parallel to this development, the technology of computer aided design (CAD) has gained a large influence on the design process of mechanical engineering. The CAD technology has already lived through a rapid development driven by the drastically growing capabilities of digital computers. However, the systems of today are still considered as being only the first generation of a long row of computer integrated manufacturing (CIM) systems. These systems to come will offer an integrated environment for design, analysis, and fabrication of products of almost any character. Thus, the CAD system could be regarded as simply a database for geometrical information equipped with a number of tools with the purpose of helping the user in the design process. Among these tools are facilities for structural analysis and optimization as well as present standard CAD features like drawing, modeling, and visualization tools. The state of the art of structural optimization is that a large amount of mathematical and mechanical techniques are available for the solution of single problems. By implementing collections of the available techniques into general software systems, operational environments for structural optimization have been created. The forthcoming years must bring solutions to the problem of integrating such systems into more general design environments. The result of this work should be CAD systems for rational design in which structural optimization is one important design tool among many others.

Rasmussen, John↗

NASA NTRS1980

The effect of radiation pressure on the restricted three-body problem

The classical restricted three-body problem is extended to account for the effects of radiation pressure and the Poynting-Robertson effect on a particle with a surface-to-volume relation such that radiation pressure forces are nonnegligible moving in the gravitational field of two larger, orbiting bodies. The positions of the L4 and L5 Lagrangian points are found as functions of the ratio of radiational to gravitational forces. The inclusion of the Poynting-Robertson effect, however, is found to render the L4 and L5 points unstable on a time scale which is long compared to the rotation period of the two massive bodies. The results imply that space colonies with solar facing adjustable sails may be able to remain stationary at any heliocentric distance between a planet and the sun, and may provide a mechanism for producing asymmetries in the interplanetary dust cloud.

Schuerman, D. W.↗

NASA NTRS1999

A Single-Lap Joint Adhesive Bonding Optimization Method Using Gradient and Genetic Algorithms

A natural process for any engineer, scientist, educator, etc. is to seek the most efficient method for accomplishing a given task. In the case of structural design, an area that has a significant impact on the structural efficiency is joint design. Unless the structure is machined from a solid block of material, the individual components which compose the overall structure must be joined together. The method for joining a structure varies depending on the applied loads, material, assembly and disassembly requirements, service life, environment, etc. Using both metallic and fiber reinforced plastic materials limits the user to two methods or a combination of these methods for joining the components into one structure. The first is mechanical fastening and the second is adhesive bonding. Mechanical fastening is by far the most popular joining technique; however, in terms of structural efficiency, adhesive bonding provides a superior joint since the load is distributed uniformly across the joint. The purpose of this paper is to develop a method for optimizing single-lap joint adhesive bonded structures using both gradient and genetic algorithms and comparing the solution process for each method. The goal of the single-lap joint optimization is to find the most efficient structure that meets the imposed requirements while still remaining as lightweight, economical, and reliable as possible. For the single-lap joint, an optimum joint is determined by minimizing the weight of the overall joint based on constraints from adhesive strengths as well as empirically derived rules. The analytical solution of the sin-le-lap joint is determined using the classical Goland-Reissner technique for case 2 type adhesive joints. Joint weight minimization is achieved using a commercially available routine, Design Optimization Tool (DOT), for the gradient solution while an author developed method is used for the genetic algorithm solution. Results illustrate the critical design variables as a function of adhesive properties and convergences of different joints based on the two optimization methods.

Smeltzer, Stanley S., III↗

NASA NTRS1985

Desorption from interstellar grains

Different desorption mechanisms from interstellar grains are considered to resolve the conflict between the observed presence of gaseous species in molecular clouds and their expected depletion onto grains. The physics of desorption is discussed with particular reference to the process of grain heating and the specific heat of the dust material. Impulsive heating by X-rays and cosmic rays is addressed. Spot heating of the grains by cosmic rays and how this can lead to desorption of mantles from very large grains is considered. It is concluded that CO depletion on grains will be small in regions with A(V) less than five from the cloud surface and n(H) less than 10,000, in agreement with observations and in contrast to expectations from pure thermal equilibrium. Even in very dense and obscured regions and in the absence of internal ultraviolet sources, the classical evaporation of CO or N2 and O2-rich mantles by cosmic rays is important.

Leger, A.↗

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

NASA NTRS1994

Coupling of microprocesses and macroprocesses due to velocity shear: An application to the low-altitude ionosphere

Recent observations indicate that low-altitude (below 1500 km) ion energization and thermal ion upwelling are colocated in the convective flow reversal region. In this region the convective velocity V(sub perpendicular) is generally small but spatial gradients in V(sub perpendicular) can be large. As a result, Joule heating is small. The observed high level of ion heating (few electron volts or more) cannot be explained by classical Joule heating alone but requires additional heating sources such as plasma waves. At these lower altitudes, sources of free energy are not obvious and hence the nature of ion energization remains ill understood. The high degree of correlation of ion heating with shear in the convective velocity (Tsunoda et al., 1989) is suggestive of an important role of velocity shear in this phenomenon. We provide more recent evidence for this correlation and show that even a small amount of velocity shear in the transverse flow is sufficient to excite a large-scale Kelvin-Helmholtz mode, which can nonlinearly steepen and give rise to highly stressed regions of strongly sheared flows. Futhermore, these stressed regions of strongly sheared flows may seed plasma waves in the range of ion cyclotron to lower hybrid frequencies, which are potential sources for ion heating. This novle two-step mechanism for ion energization is applied to typical observations of low-altitude thermal ion upwelling events.

Ganguli, G.↗

NASA NTRS1981

Momentum and energy balance in late-type stellar winds

Observations at ultraviolet and X-ray wavelengths indicate that the classical picture of a static stellar atmosphere containing a radiative equilibrium temperature distribution is inapplicable to the majority of late type stars. Mass loss and the presence of atmospheric regions characterized by gas temperatures in excess of the stellar effective temperature appear to be almost ubiquitous throughout the HR diagram. Evidence pertaining to the thermal and dynamical structure of the outer envelopes of cool stars is summarized. These results are compared with the predictions of several theoretical models which were proposed to account for mass loss from latetype stars. Models in which the outflow is thermally radiatively, or wave driven are considered for identification of the physical processes responsible for the observed wind properties. The observed variation of both the wind, thermal and dynamical structure as one proceeds from the supergiant branch toward the main sequence in the cool portion of the HR diagram give consideration to potential mechanisms for heating and cooling the flow from low gravity stars.

Macgregor, K. B.↗

NASA NTRS2004

The Search for Young Planetary Systems And the Evolution of Young Stars

The Space Interferometer Mission (SIM) will provide a census of planetary systems by con- ducting a broad survey of 2,000 stars that will be sensitive to the presence of planets with masses as small as approx. 15 Earth masses (1 Uranus mass) and a deep survey of approx. 250 of the nearest, stars with a mass limit of approx.3 Earth masses. The broad survey will include stars spanning a wide range of ages, spectral types, metallicity, and other important parameters. Within this larger context, the Young Stars and Planets Key Project will study approx. 200 stars with ages from 1 Myr to 100 Myr to understand the formation and dynamical evolution of gas giant planets. The SIM Young Stars and Planets Project will investigate both the frequency of giant planet formation and the early dynamical history of planetary systems. We will gain insight into how common the basic architecture of our solar system is compared with recently discovered systems with close-in giant planets by examining 200 of the nearest (less than 150 pc) and youngest (1-100 Myr) solar-type stars for planets. The sensitivity of the survey for stars located 140 pc away is shown in the planet mass-separation plane. We expect to find anywhere from 10 (assuming that only the presently known fraction of stars. 5-7%, has planets) to 200 (all young stars have planets) planetary systems. W-e have set our sensitivity threshold to ensure the detection of Jupiter-mass planets in the critical orbital range of 1 to 5 AU. These observations, when combined with the results of planetary searches of mature stars, will allow us to test theories of planetary formation and early solar system evolution. By searching for planets around pre-main sequence stars carefully selected to span an age range from 1 to 100 Myr, we will learn a t what epoch and with what frequency giant planets are found at the water-ice snowline where they are expected to form. This will provide insight into the physical mechanisms by which planets form and migrate from their place of birth, and about their survival rate. With these data in hand, we will provide data, for the first time, on such important questions as: What processes affect the formation and dynamical evolution of planets? When and where do planets form? What is initial mass distribution of planetary systems around young stars? How might planets be destroyed? What is the origin of the eccentricity of planetary orbits? What is the origin of the apparent dearth of companion objects between planets and brown dwarfs seen in mature stars? The observational strategy is a compromise between the desire to extend the planetary mass function as low as possible and the essential need to build up sufficient statistics on planetary occurrence. About half of the sample will be used to address the "where" and "when" of planet formation. We will study classical T Tauri stars (cTTs) which have massive accretion disks and post- accretion, weak-lined T Tauri stars (wTTs). Preliminary estimates suggest the sample will consist of approx. 30% cTTs and approx. 70% wTTs, driven in part by the difficulty of making accurate astrometric measurements toward objects with strong variability or prominent disks.

Beichman, Charles A.↗

NASA NTRS2003

Analytical Modeling of the High Strain Rate Deformation of Polymer Matrix Composites

The results presented here are part of an ongoing research program to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric matrix materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical 5 plasticity theory definitions of effective stress and effective plastic strain are modified by applying variations of the Drucker-Prager yield criterion. To verify the revised formulation, the shear and tensile deformation of a representative toughened epoxy is analyzed across a wide range of strain rates (from quasi-static to high strain rates) and the results are compared to experimentally obtained values. For the analyzed polymers, both the tensile and shear stress-strain curves computed using the analytical model correlate well with values obtained through experimental tests. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. In the micromechanics, the unit cell is divided up into a number of independently analyzed slices, and laminate theory is then applied to obtain the effective deformation of the unit cell. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite (composed using the representative polymer analyzed for the correlation of the polymer constitutive equations) for several fiber orientation angles across a variety of strain rates. The computed values compare favorably to experimentally obtained results.

Goldberg, Robert K.↗

NASA NTRS1995

The nonlinear evolution of inviscid Goertler vortices in three-dimensional boundary layers

The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is considered. We do not follow the classical approach of weakly nonlinear stability problems and consider a mode which has just become unstable. Instead we extend the method of Blackaby, Dando, and Hall (1992), which considered the closely related nonlinear development of disturbances in stratified shear flows. The Gortler modes we consider are initially fast growing and we assume, following others, that boundary-layer spreading results in them evolving in a linear fashion until they reach a stage where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando and Hall (1993) is apparent, given the range of parameters for the Gortler problem, that there are three possible nonlinear integro-differential evolution equations for the disturbance amplitude. These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism investigated in this previous paper, and a quintic. In this paper we shall concentrate on the two cubic integro-differential equations and in particular, on the one due to the novel mechanism as this will be the first to affect a disturbance. It is found that the consideration of a spatial evolution problem as opposed to temporal (as was considered in Blackaby, Dando, and Hall, 1992) causes a number of significant changes to the evolution equations.

Blackaby, Nicholas↗

NASA NTRS2003

Competing Classical and Quantum Effects in Shape Relaxation of a Metallic Nanostructure

We demonstrate for the first time that the quantum size effect (QSE) plays a competing role along side the classical thermodynamic effect in the shape relaxation of a small metallic island. Together, these effects transforms a lead(Pb) island grown on Si(111) substrate from its initially flattop faceted morphology to a peculiar ring-shape island, a process catalysed by the tip electric field of a scanning tunnelling microscope (STM). We shall show for the first time how QSE affects the relaxation process dynamically. In particular, it leads to a novel strip-flow growth and double-step growth on selective strips of a plateau inside the ring, defined by the substrate steps more than 60?0?3 below. It appears that atoms diffusing on the plateau can clearly (sub i)(deg)sense(sub i)+/- the quantized energy states inside the island and have preferentially attached to regions that further reduces the surface energy as a result of the QSE, limiting its own growth and stabilizing the ring shape. The mechanism proposed here offers a sound explanation for ring shape metal and semiconductor islands observed in other systems as well.

Chen, Dongmin↗

NASA NTRS2003

On the Hilbert-Huang Transform Data Processing System Development

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). The Fourier view of nonlinear mechanics that had existed for a long time, and the associated FFT (fairly recent development), carry strong a-priori assumptions about the source data, such as linearity and of being stationary. Natural phenomena measurements are essentially nonlinear and nonstationary. A very recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT) proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using the Empirical Mode Decomposition (EMD) followed by the Hilbert Transform of the empirical decomposition data (HT), the HHT allows spectrum analysis of nonlinear and nonstationary data by using an engineering a-posteriori data processing, based on the EMD algorithm. This results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF) that can be further analyzed for spectrum interpretation by the classical Hilbert Transform. This paper describes phase one of the development of a new engineering tool, the HHT Data Processing System (HHTDPS). The HHTDPS allows applying the "T to a data vector in a fashion similar to the heritage FFT. It is a generic, low cost, high performance personal computer (PC) based system that implements the HHT computational algorithms in a user friendly, file driven environment. This paper also presents a quantitative analysis for a complex waveform data sample, a summary of technology commercialization efforts and the lessons learned from this new technology development.

Kizhner, Semion↗

NASA NTRS2015

Unsteady Loss in the Stator Due to the Incoming Rotor Wake in a Highly-Loaded Transonic Compressor

The present paper reports an investigation of unsteady loss generation in the stator due to the incoming rotor wake in an advanced GE transonic compressor design with a high-fidelity numerical method. This advanced compressor with high reaction and high stage loading has been investigated both experimentally and analytically in the past. The measured efficiency in this advanced compressor is significantly lower than the design intention/goal. The general understanding is that the current generation of compressor design/analysis tools miss some important flow physics in this modern compressor design. To pinpoint the source of the efficiency miss, an advanced test with a detailed flow traverse was performed for the front one and a half stage at the NASA Glenn Research Center. Detailed data-match analysis by GE identified an unexpected high loss generation in the pressure side of the stator passage. Higher total temperature and lower total pressure are measured near the pressure side of the stator. Various analyses based on the RANS and URANS of the compressor stage do not calculate the measured higher total temperature and lower total pressure on the pressure side of the stator. In the present paper, a Large Eddy Simulation (LES) is applied to find the fundamental mechanism of this unsteady loss generation in the stator due to the incoming rotor wake. The results from the LES were first compared with the NASA test results and the GE interpretation of the test data. LES calculates lower total pressure and higher total temperature on the pressure side of the stator, as the measured data showed, resulting in large loss generation on the pressure side of the stator. Detailed examination of the unsteady flow field from LES shows that the rotor wake, which has higher total temperature and higher total pressure relative to the free stream, interacts quite differently with the pressure side of the blade compared to the suction side of the blade. The higher temperature in the wake remains high as the wake passes through the pressure side of the blade. On the other hand, the total temperature diffuses as it passes through near the suction surface. For the presently investigated compressor, the classical intra-stator wake transport to the pressure side of the blade by the slip velocity in the wake seems to be minor. The main causes of this phenomenon are three-dimensional unsteady vortex interactions near the blade surface. The stabilizing effect of the concave curvature on the suction side keeps the rotor wake thin. On the other hand, the destabilizing effect of the convex curvature of the pressure side makes the rotor wake thicker, which results in a higher total temperature measurement at the stator exit. Additionally, wake stretching through the stator seems to contribute to the redistribution of the total temperature and the loss generation.

Efficiency↗

NASA NTRS1994

Micromechanical Modeling of the Thermal Expansion of Graphite/copper Composites with Nonuniform Microstructure

Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.

Bednarcyk, Brett A.↗

NASA NTRS2013

Modeling Non-homologous End Joining

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several NHEJ proteins such as Ku, DNA-PKcs, XRCC4, Ligase IV and so on. Once DSBs are generated, Ku is first recruited to the DNA end, followed by other NHEJ proteins for DNA end processing and ligation. Because of the direct ligation of break ends without the need for a homologous template, NHEJ turns out to be an error-prone but efficient repair pathway. Some mechanisms have been proposed of how the efficiency of NHEJ repair is affected. The type of DNA damage is an important factor of NHEJ repair. For instance, the length of DNA fragment may determine the recruitment efficiency of NHEJ protein such as Ku [1], or the complexity of the DNA breaks [2] is accounted for the choice of NHEJ proteins and subpathway of NHEJ repair. On the other hand, the chromatin structure also plays a role of the accessibility of NHEJ protein to the DNA damage site. In this talk, some mathematical models of NHEJ, that consist of series of biochemical reactions complying with the laws of chemical reaction (e.g. mass action, etc.), will be introduced. By mathematical and numerical analysis and parameter estimation, the models are able to capture the qualitative biological features and show good agreement with experimental data. As conclusions, from the viewpoint of modeling, how the NHEJ proteins are recruited will be first discussed for connection between the classical sequential model [4] and recently proposed two-phase model [5]. Then how the NHEJ repair pathway is affected, by the length of DNA fragment [6], the complexity of DNA damage [7] and the chromatin structure [8], will be addressed

Li, Yongfeng↗