Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical analysis”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Thermal Noise Limit in Frequency Stabilization of Lasers with Rigid Cavities

We evaluated thermal noise (Brownian motion) in a rigid reference cavity Used for frequency stabilization of lasers, based on the mechanical loss of cavity materials and the numerical analysis of the mirror-spacer mechanics with the direct application of the fluctuation dissipation theorem. This noise sets a fundamental limit for the frequency stability achieved with a rigid frequency-reference cavity of order 1 Hz/rtHz at 10mHz at room temperature. This level coincides with the world-highest level stabilization results.

Numata, Kenji↗
NASA NTRS2005

System-Level Design of a Shape Memory Alloy Actuator for Active Clearance Control in the High-Pressure Turbine

This paper describes results of a numerical analysis evaluating the feasibility of high-temperature shape memory alloys (HTSMA) for active clearance control actuation in the high-pressure turbine section of a modern turbofan engine. The prototype actuator concept considered here consists of parallel HTSMA wires attached to the shroud that is located on the exterior of the turbine case. A transient model of an HTSMA actuator was used to evaluate active clearance control at various operating points in a test bed aircraft engine simulation. For the engine under consideration, each actuator must be designed to counteract loads from 380 to 2000 lbf and displace at least 0.033 inches. Design results show that an actuator comprised of 10 wires 2 inches in length is adequate for control at critical engine operating points and still exhibits acceptable failsafe operability and cycle life. A proportional-integral-derivative (PID) controller with integrator windup protection was implemented to control clearance amidst engine transients during a normal mission. Simulation results show that the control system exhibits minimal variability in clearance control performance across the operating envelope. The final actuator design is sufficiently small to fit within the limited space outside the high-pressure turbine case and is shown to consume only small amounts of bleed air to adequately regulate temperature.

DeCastro, Jonathan A.↗
NASA NTRS2005

Torsion of Noncircular Composite Cylinders

The paper presents a brief overview of the predicted deformation and failure characteristics of noncircular composite cylinders subjected to torsion. Using a numerical analysis, elliptical cylinders with a minor-to-major diameter ratio of 0.7 are considered. Counterpart circular cylinders with the same circumference as the elliptical cylinders are included for comparison. The cylinders are constructed of a medium-modulus graphite-epoxy material in a quasi-isotropic lay-up. Imperfections generated from the buckling mode shapes are included in the initial cross-sectional geometry of the cylinders. Deformations until first fiber failure, as predicted using the maximum stress failure criterion and a material degradation scheme, are presented. For increasing levels of torsion, the deformations of the elliptical cylinders, in the form of wrinkling of the cylinder wall, occur primarily in the flatter regions of the cross section. By comparison the wrinkling deformations of the circular cylinders are more uniformly distributed around the circumference. Differences in the initial failure and damage progression and the overall torque vs. twist relationship between the elliptical and circular cylinders are presented. Despite differences in the response as the cylinders are being loaded, at first fiber failure the torque and twist for the elliptical and circular cylinders nearly coincide.

Rouse, Marshall↗
NASA NTRS2006

Progress on Shape Memory Alloy Actuator Development for Active Clearance Control

Results of a numerical analysis evaluating the feasibility of high-temperature shape memory alloys (HTSMA) for active clearance control actuation in the high-pressure turbine section of a modern turbofan engine has been conducted. The prototype actuator concept considered here consists of parallel HTSMA wires attached to the shroud that is located on the exterior of the turbine case. A transient model of an HTSMA actuator was used to evaluate active clearance control at various operating points in a test bed aircraft engine simulation. For the engine under consideration, each actuator must be designed to counteract loads from 380 to 2000 lbf and displace at least 0.033 in. Design results show that an actuator comprised of 10 wires 2 in. in length is adequate for control at critical engine operating points and still exhibit acceptable failsafe operability and cycle life. A proportional-integral-derivative (PID) controller with integrator windup protection was implemented to control clearance amidst engine transients during a normal mission. Simulation results show that the control system exhibits minimal variability in clearance control performance across the operating envelope. The final actuator design is sufficiently small to fit within the limited space outside the high-pressure turbine case and is shown to consume only small amounts of bleed air to adequately regulate temperature.

DeCastro, Jonathan↗
NASA NTRS2006

Percentage Contributions from Atmospheric and Surface Features to Computed Brightness Temperatures

Over the past few years, there has become an increasing interest in the use of millimeter-wave (mm-wave) and sub-millimeter-wave (submm-wave) radiometer observations to investigate the properties of ice particles in clouds. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. One methodology used since the 1960's to discern the relationship between the physical state observed and the brightness temperature (TB) is through the temperature weighting function profile. In this research, the temperature weighting function concept is exploited to analyze the sensitivity of various characteristics of the cloud profile, such as relative humidity, ice water path, liquid water path, and surface emissivity. In our numerical analysis, we compute the contribution (in Kelvin) from each of these cloud and surface characteristics, so that the sum of these various parts equals the computed TB. Furthermore, the percentage contribution from each of these characteristics is assessed. There is some intermingling/contamination of the contributions from various components due to the integrated nature of passive observations and the absorption and scattering between the vertical layers, but all in all the knowledge gained is useful. This investigation probes the sensitivity over several cloud classifications, such as cirrus, blizzards, light snow, anvil clouds, and heavy rain. The focus is on mm-wave and submm-wave frequencies, however discussions of the effects of cloud variations to frequencies as low as 10 GHz and up to 874 GHz will also be presented. The results show that nearly 60% of the TB value at 89 GHz comes from the earth's surface for even the heaviest blizzard snow rates. On the other hand, a significant percentage of the TB value comes from the snow in the cloud for 166, and 183 plus or minus 7 GHz for the heavy and medium snow rates. For submm-wave channels, there is no contribution from the surface because these channels cannot probe through clouds, nor normal water vapor amounts in clear air regions. This work is extremely valuable in physically-based retrieval algorithm development research.

Gail Skofronick-Jackson↗
NASA NTRS2007

Investigation of a Parabolic Iterative Solver for Three-dimensional Configurations

A parabolic iterative solution procedure is investigated that seeks to extend the parabolic approximation used within the internal propagation module of the duct noise propagation and radiation code CDUCT-LaRC. The governing convected Helmholtz equation is split into a set of coupled equations governing propagation in the positive and negative directions. The proposed method utilizes an iterative procedure to solve the coupled equations in an attempt to account for possible reflections from internal bifurcations, impedance discontinuities, and duct terminations. A geometry consistent with the NASA Langley Curved Duct Test Rig is considered and the effects of acoustic treatment and non-anechoic termination are included. Two numerical implementations are studied and preliminary results indicate that improved accuracy in predicted amplitude and phase can be obtained for modes at a cut-off ratio of 1.7. Further predictions for modes at a cut-off ratio of 1.1 show improvement in predicted phase at the expense of increased amplitude error. Possible methods of improvement are suggested based on analytic and numerical analysis. It is hoped that coupling the parabolic iterative approach with less efficient, high fidelity finite element approaches will ultimately provide the capability to perform efficient, higher fidelity acoustic calculations within complex 3-D geometries for impedance eduction and noise propagation and radiation predictions.

Nark, Douglas M.↗
NASA NTRS2007

Local Elastic Constants for Epoxy-Nanotube Composites from Molecular Dynamics Simulation

A method from molecular dynamics simulation is developed for determining local elastic constants of an epoxy/nanotube composite. The local values of C11, C33, K12, and K13 elastic constants are calculated for an epoxy/nanotube composite as a function of radial distance from the nanotube. While the results possess a significant amount of statistical uncertainty resulting from both the numerical analysis and the molecular fluctuations during the simulation, the following observations can be made. If the size of the region around the nanotube is increased from shells of 1 to 6 in thickness, then the scatter in the data reduces enough to observe trends. All the elastic constants determined are at a minimum 20 from the center of the nanotube. The C11, C33, and K12 follow similar trends as a function of radial distance from the nanotube. The K13 decreases greater distances from the nanotube and becomes negative which may be a symptom of the statistical averaging.

Frankland, S. J. V.↗
NASA NTRS2007

Automated Detection of Events of Scientific Interest

A report presents a slightly different perspective of the subject matter of Fusing Symbolic and Numerical Diagnostic Computations (NPO-42512), which appears elsewhere in this issue of NASA Tech Briefs. Briefly, the subject matter is the X-2000 Anomaly Detection Language, which is a developmental computing language for fusing two diagnostic computer programs one implementing a numerical analysis method, the other implementing a symbolic analysis method into a unified event-based decision analysis software system for real-time detection of events. In the case of the cited companion NASA Tech Briefs article, the contemplated events that one seeks to detect would be primarily failures or other changes that could adversely affect the safety or success of a spacecraft mission. In the case of the instant report, the events to be detected could also include natural phenomena that could be of scientific interest. Hence, the use of X- 2000 Anomaly Detection Language could contribute to a capability for automated, coordinated use of multiple sensors and sensor-output-data-processing hardware and software to effect opportunistic collection and analysis of scientific data.

James, Mark↗
NASA NTRS2011

Steady and Unsteady Flow Effects of Circumferential Grooves Casing Treatment in a Transonic Compressor Rotor

The current paper reports on an investigation of steady and unsteady flow effects of circumferential grooves casing treatment in a transonic compressor rotor. Circumferential grooves casing treatment is used mainly to increase stall margin in axial compressors with a relatively small decrease in aerodynamic efficiency. It is widely believed that flow mechanisms of circumferential grooves casing treatment near stall conditions are not yet well understood even though this treatment has been used widely in real engines. Numerical analysis based on steady Reynolds-averaged Navier-Stokes (RANS) has been the primary tool used to understand flow mechanism for circumferential grooves casing treatment. Although steady RANS explains some flow effects of circumferential grooves casing treatment, it does not calculate all the measured changes in the compressor characteristics. Therefore, design optimization of circumferential grooves with steady RANS has not been very successful. As a compressor operates toward the stall condition, the flow field becomes transient. Major sources of self-generated flow unsteadiness are shock oscillation and interaction between the passage shock and the tip leakage vortex. In the present paper, an unsteady Reynolds-averaged Navier-Stokes (URANS) approach is applied to study the effects of circumferential grooves in a transonic compressor. The results from URANS are compared with the results from RANS and measured data. The current investigation shows that there are significant unsteady flow effects on the performance of the circumferential grooves casing treatment. For the currently investigated rotor, the unsteady effects are of the same magnitude as the steady effects in terms of extending the compressor stall margin.

Hah, Chunill↗
NASA NTRS2012

An Experimental and Numerical Study of Roughness-Induced Instabilities in a Mach 3.5 Boundary Layer

Progress on a joint experimental and numerical study of laminar-to-turbulent transition induced by an isolated roughness element in a high-speed laminar boundary layer is reported in this paper. The numerical analysis suggests that transition is driven by the instability of high- and low-speed streaks embedded in the wake of the isolated roughness element. In addition, spatial stability analysis revealed that the wake flow supports multiple modes (even and odd) of convective instabilities that experience strong enough growth to cause transition. The experimental measurements, which included hot-wire and pitot-probe surveys, confirmed the existence of embedded high- and low-speed streaks in the roughness wake. Furthermore, the measurements indicate the presence of both even and odd modes of instability, although their relative magnitude depends on the specifics of the roughness geometry and flow conditions (e.g., the value of Re(sub kk) or k/delta. For the two test cases considered in the measurements (Re(sub kk) values of 462 and 319), the even mode and the odd mode were respectively dominant and appear to play a primary role in the transition process.

Kegerise, Michael A.↗
NASA NTRS2013

Alternate Methods in Refining the SLS Nozzle Plug Loads

Numerical analysis has shown that the SLS nozzle environmental barrier (nozzle plug) design is inadequate for the prelaunch condition, which consists of two dominant loads: 1) the main engines startup pressure and 2) an environmentally induced pressure. Efforts to reduce load conservatisms included a dynamic analysis which showed a 31% higher safety factor compared to the standard static analysis. The environmental load is typically approached with a deterministic method using the worst possible combinations of pressures and temperatures. An alternate probabilistic approach, utilizing the distributions of pressures and temperatures, resulted in a 54% reduction in the environmental pressure load. A Monte Carlo simulation of environmental load that used five years of historical pressure and temperature data supported the results of the probabilistic analysis, indicating the probabilistic load is reflective of a 3-sigma condition (1 in 370 probability). Utilizing the probabilistic load analysis eliminated excessive conservatisms and will prevent a future overdesign of the nozzle plug. Employing a similar probabilistic approach to other design and analysis activities can result in realistic yet adequately conservative solutions.

Burbank, Scott↗
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↗
NASA NTRS2015

An Automated DAKOTA and VULCAN-CFD Framework with Application to Supersonic Facility Nozzle Flowpath Optimization

Removing human interaction from design processes by using automation may lead to gains in both productivity and design precision. This memorandum describes efforts to incorporate high fidelity numerical analysis tools into an automated framework and applying that framework to applications of practical interest. The purpose of this effort was to integrate VULCAN-CFD into an automated, DAKOTA-enabled framework with a proof-of-concept application being the optimization of supersonic test facility nozzles. It was shown that the optimization framework could be deployed on a high performance computing cluster with the flow of information handled effectively to guide the optimization process. Furthermore, the application of the framework to supersonic test facility nozzle flowpath design and optimization was demonstrated using multiple optimization algorithms.

Axdahl, Erik L.↗
NASA NTRS2015

X-Band CubeSat Communication System Demonstration

Today's CubeSats mostly operate their communications at UHF- and S-band frequencies. UHF band is presently crowded, thus downlink communications are at lower data rates due to bandwidth limitations and are unreliable due to interference. This research presents an end-to-end robust, innovative, compact, efficient and low cost S-band uplink and X-band downlink CubeSat communication system demonstration between a balloon and a Near Earth Network (NEN) ground system. Since communication systems serve as umbilical cords for space missions, demonstration of this X-band communication system is critical for successfully supporting current and future CubeSat communication needs. This research has three main objectives. The first objective is to design, simulate, and test a CubeSat S- and X-band communication system. Satellite Tool Kit (STK) dynamic link budget calculations and HFSS Simulations and modeling results have been used to trade the merit of various designs for small satellite applications. S- and X-band antennas have been tested in the compact antenna test range at Goddard Space Flight Center (GSFC) to gather radiation pattern data. The second objective is simulate and test a CubeSat compatible X-band communication system at 12.5Mbps including S-band antennas, X-band antennas, Laboratory for Atmospheric and Space Physics (LASP) /GSFC transmitter and an S-band receiver from TRL-5 to TRL-8 by the end of this effort. Different X-band communication system components (antennas, diplexers, etc.) from GSFC, other NASA centers, universities, and private companies have been investigated and traded, and a complete component list for the communication system baseline has been developed by performing analytical and numerical analysis. This objective also includes running simulations and performing trades between different X-band antenna systems to optimize communication system performance. The final objective is to perform an end-to-end X-band CubeSat communication system demonstration between a balloon and/or a sounding rocket and a Near Earth Network (NEN) ground system. This paper presents CubeSat communication systems simulation results, analysis of X-band and S-band antennas and RF front-end components, transceiver design, analysis and optimization of space-to-ground communication performance, subsystem development, as well as the test results for an end-to-end X-band CubeSat communication system demonstration. The outcome of this work will be used to pave the way for next generation NEN-compatible X-band CubeSat communication systems to support higher data rates with more advanced modulation and forward error correction (FEC) coding schemes, and to support and attract new science missions at lower cost. It also includes an abbreviated concept of operations for CubeSat users to utilize the NEN, starting from first contact with NASA's communication network and continuing through on-orbit operations.

NEN↗
NASA NTRS2017

Micromechanics of Composite Materials Governed by Vector Constitutive Laws

The high-fidelity generalized method of cells micromechanics theory has been extended for the prediction of the effective property tensor and the corresponding local field distributions for composites whose constituents are governed by vector constitutive laws. As shown, the shear analogy, which can predict effective transverse properties, is not valid in the general three-dimensional case. Consequently, a general derivation is presented that is applicable to both continuously and discontinuously reinforced composites with arbitrary vector constitutive laws and periodic microstructures. Results are given for thermal and electric problems, effective properties and local field distributions, ordered and random microstructures, as well as complex geometries including woven composites. Comparisons of the theory's predictions are made to test data, numerical analysis, and classical expressions from the literature. Further, classical methods cannot provide the local field distributions in the composite, and it is demonstrated that, as the percolation threshold is approached, their predictions are increasingly unreliable. XXXX It has been observed that the bonding between the fibers and matrix in composite materials can be imperfect. In the context of thermal conductivity, such imperfect interfaces have been investigated in micromechanical models by Dunn and Taya (1993), Duan and Karihaloo (2007), Nan et al. (1997) and Hashin (2001). The present HFGMC micromechanical method, derived for perfectly bonded composite materials governed by vector constitutive laws, can be easily generalized to include the effects of weak bonding between the constituents. Such generalizations, in the context of the mechanical micromechanics problem, involve introduction of a traction-separation law at the fiber/matrix interface and have been presented by Aboudi (1987), Bednarcyk and Arnold (2002), Bednarcyk et al. (2004) and Aboudi et al. (2013) and will be addressed in the future.

micromechanics↗
NASA NTRS2018

Meander-Like Shape Optimization of the Stator-Rotor Platform Cavity

Recent progress in additive manufacturing has enabled opportunities to explore novel stator rim geometries which can be implemented to improve cooling strategies in turbomachinery. This paper presents a simplified stationary geometry optimization strategy to produce enhanced stator-rotor cavity sealing and highlights main driving mechanisms.The stator and rotor rims were designed using a design strategy based on inspiration from the meandering of rivers. A minimum thickness of 2 millimeters was maintained throughout the cavity to ensure a practical implementation. The computational domain comprised of the stator outlet, hub disk leakage cavity, and rotor platform was meshed using NUMECA Int. package, Hexpress. The numerical analysis required 3D Unsteady Reynolds Average Navier-Stokes to replicate vorticial structures using Ansys Fluent. The operating conditions were representative of engine-like conditions, exploring a wide range of mass flow ratios from 1 to 3 percent. The optimization yielded designs that provide 30 percent reduction in rear platform temperature while minimizing coolant mass flow. The applicability of the design was compared against 3D sector in both stationary and in rotation.

Shape Optimization↗