Search NASA⌕ Search

SEARCH · Search NASA

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 793 records · Page 44

Constant Density Approximations for the Flow Behind Axisymmetric Shock Waves

The incompressible rotational flow equations are used to obtain solutions for the flow behind axisymmetric shock waves with conic longitudinal sections. The nonlinear part of the term due to rotation is retained in the analysis. Numerical results for standoff distance and stagnation point velocity gradient are presented for the case in which the shock wave is a paraboloid, a sphere, or an oblate or prolate ellipsoid. A similarity parameter is proposed which correlates approximately the flow behind geometrically similar shock waves at different free-stream conditions.

Munson, Albert G.↗

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 t.he 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/square root Hz(0.01 Hz/square root Hz) at 10 mHz (100 Hz) at room temperature. This level coincides with the world-highest level stabilization results.

Numata, Kenji↗

Effects of training on muscle O2 transport at VO2max

To quantify the relative contributions of convective and peripheral diffusive components of O2 transport to the increase in leg O2 uptake (VO2leg) at maximum O2 uptake (VO2max) after 9 wk of endurance training, 12 sedentary subjects (age 21.8 +/- 3.4 yr, VO2max 36.9 +/- 5.9 ml.min-1.kg-1) were studied. VO2max, leg blood flow (Qleg), and arterial and femoral venous PO2, and thus VO2leg, were measured while the subjects breathed room air, 15% O2, and 12% O2. The sequence of the three inspirates was balanced. After training, VO2max and VO2leg increased at each inspired O2 concentration [FIO2; mean over the 3 FIO2 values 25.2 +/- 17.8 and 36.5 +/- 33% (SD), respectively]. Before training, VO2leg and mean capillary PO2 were linearly related through the origin during hypoxia but not during room air breathing, suggesting that, at 21% O2, VO2max was not limited by O2 supply. After training, VO2leg and mean capillary PO2 at each FIO2 fell along a straight line with zero intercept, just as in athletes (Roca et al. J. Appl. Physiol. 67: 291-299, 1989). Calculated muscle O2 diffusing capacity (DO2) rose 34% while Qleg increased 19%. The relatively greater rise in DO2 increased the DO2/Qleg, which led to 9.9% greater O2 extraction. By numerical analysis, the increase in Qleg alone (constant DO2) would have raised VO2leg by 35 ml/min (mean), but that of DO2 (constant Qleg) would have increased VO2leg by 85 ml/min, more than twice as much. The sum of these individual effects (120 ml/min) was less (P = 0.013) than the observed rise of 164 ml/min (mean). This synergism (explained by the increase in DO2/Qleg) seems to be an important contribution to increases in VO2max with training.

NASA Discipline Musculoskeletal↗

Distributed Actuation and Sensing on an Uninhabited Aerial Vehicle

An array of effectors and sensors has been designed, tested and implemented on a Blended Wing Body Uninhabited Aerial Vehicle (UAV). The UAV is modified to serve as a flying, controls research, testbed. This effector/sensor array provides for the dynamic vehicle testing of controller designs and the study of decentralized control techniques. Each wing of the UAV is equipped with 12 distributed effectors that comprise a segmented array of independently actuated, contoured control surfaces. A single pressure sensor is installed near the base of each effector to provide a measure of deflections of the effectors. The UAV wings were tested in the North Carolina State University Subsonic Wind Tunnel and the pressure distribution that result from the deflections of the effectors are characterized. The results of the experiments are used to develop a simple, but accurate, prediction method, such that for any arrangement of the effector array the corresponding pressure distribution can be determined. Numerical analysis using the panel code CMARC verifies this prediction method.

Barnwell, William Garrard↗

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↗

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

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↗

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↗

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↗

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

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

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↗

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↗

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

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↗

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↗