Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximation & 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Design and Analysis of Shape Memory Spring Tires for Martian and Lunar Rover Vehicles

Shape memory alloys (SMAs) have played an important role in various innovative engineering and medical applications, such as aerospace actuators, vibration damping devices, and coronary stents. In applications, shape memory alloys are commonly utilized in two fundamentally different ways: (i) making use of the superelasticity/ pseudoelasticity (SE/PE) phenomena, as in applications in biomedical engineering, and (ii) taking advantages of the shape memory effect (SME), as is used for actuators. Their ability to act in such vastly different capacities is mainly due to their unique capability to recover large amounts of deformation produced by either applied stresses or temperature changes. One recent emerging application in use of SMAs has been in the area of non-pneumatic tire designs for Martian or Lunar roving vehicles. These vehicles require tires that are capable of traversing rugged terrain while withstanding extreme temperatures and atmospheric conditions. Inspired by the flexible wire mesh tires used on three Lunar Roving Vehicle (LRV) missions to the Moon on Apollo 15, 16, and 17, a new compliant tire technology was developed by the NASA Glenn Research Center (GRC) and Goodyear Tire & Rubber, known as the Spring Tire. The Spring Tire consists of several hundred coiled springs woven into a flexible mesh and formed into the shape of a tire. Like the LRV wire mesh tires, the original Spring Tires were made from spring steel and were prone to permanent deformation when undergoing high localized loads. Later, a new iteration of this technology was invented, known as the ‘Superelastic Tire’. This new technology incorporated the use of superelastic SMA springs, which could effectively undergo approximately 30 times more reversible deformation than the steel spring. It also provided even greater durability and allowed for more flexibility in design, such as the use of other structural forms to reduce mass or increase load carrying capacity. Because of the unique nature of both the SMA material and the complex interactions between the springs, designing Spring Tires for a specific application requires extensive effort. Historically, design decisions have relied on full-scale empirical testing; however, this is very expensive and time consuming, especially when multiple iterations are needed. Therefore, developing a large-scale, robust, and predictive numerical model entailing complex spring interactions and the shape memory material behavior within a tire construct is the first essential step towards a successful design program. The current work focuses on implementation of the user-defined Shape Memory Alloy (SMA) model, otherwise known as SMA-GVIPs, in the Finite Element analysis (FEA) program ABAQUS for large-scale simulations of the GRC-developed Spring Tires made of SMA. The novelty of this work lies in the thorough, detail-oriented, and computationally efficient finite element analysis of full-scale SMA tires. A well-thought material characterization plan followed by model validation and a systemic sensitivity study on SMA tires has never been reported in the previous literature. The main objective of this study is to help the team improve and optimize the structural design of the SMA tires through in-depth numerical analysis and sensitivity studies. Various design variables (wire diameter, coil diameter, pitch, bead angle, and number of springs) were varied to study their influence on the global load-displacement response of the tire construct. A detailed investigation of the three-dimensional stress states was also carried out to enhance our understanding of the local changes as the tire goes through global deformation. It was concluded that a robust numerical model with a good predictive capability, together with a thoughtfully crafted sensitivity study can result in improved design iterations required to reach a desired tire performance while, significantly reducing manufacturing, labor and testing expenses. A summary of the Finite Element (FE) model construction will be presented together with a description of the user-defined SMA model, characterization process, experimental results, model validation, and numerical sensitivity study results.

shape memory alloys↗
NASA NTRS2004

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗
NASA NTRS2008

The Spectrum of 1ES0229 + 200 and the Cosmic Infrared Background

Aims. We reexamine the implications of the recent HESS observations of the blazar 1ES0229+200 for constraining the extragalactic mid-infrared background radiation. Methods. We examine the effect of gamma-ray absorption by the extragalactic infrared radiation on predicted intrinsic spectra for this blazar and compare our results with the observational data. Results. We find agreement with our previous results on the shape of the infrared spectral energy distribution, contrary to the recent assertion of the HESS group. Our analysis indicates that 1ES0229+200 has a very hard intrinsic spectrum with a spectral index between 1.1 +/- 0.3 and 1.5 0.3 in the energy range between 0.5 TeV and 15 TeV. Conclusions. Under the assumptions that (1) the models of Stecker et al. as derived from numerous detailed infrared observations are reasonable, and (2) spectral indexes in the range 1 < Gamma < 1.5 are obtainable from relativistic shock acceleration under the astrophysical conditions extant in blazar flares, the fits to the observations of 1ES0229+200 using our previous infrared spectral energy distributions are consistent with both the infrared and -ray observations. Our analysis presents evidence indicating that the energy spectrum of relativistic particles in 1ES0229+200 is produced by relativistic shock acceleration, producing an intrinsic gamma-ray spectrum with index 1 < Gamma < 1.5 and with no evidence of a peak in the spectral energy distribution up to energies approximately 15 TeV.

Stecker, F. W.↗
NASA NTRS1984

Numerical optimization design of advanced transonic wing configurations

A computationally efficient and versatile technique for use in the design of advanced transonic wing configurations has been developed. A reliable and fast transonic wing flow-field analysis program, TWING, has been coupled with a modified quasi-Newton method, unconstrained optimization algorithm, QNMDIF, to create a new design tool. Fully three-dimensional wing designs utilizing both specified wing pressure distributions and drag-to-lift ration minimization as design objectives are demonstrated. Because of the high computational efficiency of each of the components of the design code, in particular the vectorization of TWING and the high speed of the Cray X-MP vector computer, the computer time required for a typical wing design is reduced by approximately an order of magnitude over previous methods. In the results presented here, this computed wave drag has been used as the quantity to be optimized (minimized) with great success, yielding wing designs with nearly shock-free (zero wave drag) pressure distributions and very reasonable wing section shapes.

Cosentino, G. B.↗
NASA NTRS1985

Numerical optimization design of advanced transonic wing configurations

A computationally efficient and versatile technique for use in the design of advanced transonic wing configurations has been developed. A reliable and fast transonic wing flow-field analysis program, TWING, has been coupled with a modified quasi-Newton method, unconstrained optimization algorithm, QNMDIF, to create a new design tool. Fully three-dimensional wing designs utilizing both specified wing pressure disributions and drag-to-lift ratio minimization as design objectives are demonstrated. Because of the high computational efficiency of each of the components of the design code, in particular the vectorization of TWING and the high speed of the Cray X-MP vector computer, the computer time required for a typical wing design is reduced by approximately an order of magnitude over previous methods. In the results presented here, this computed wave drag has been used as the quantity to be optimized (minimized) with great success, yielding wing designs with nearly shock-free (zero wave drag) pressure distributions and very reasonable wing section shapes.

Cosentino, G. B.↗
NASA NTRS1961

Method for Approximating the Vacuum Motions of Spinning Symmetrical Bodies with Nonconstant Spin Rates

A method for approximating the vacuum motions of spinning rigid symmetrical bodies with varying spin rates and inertias has been completed. The analysis includes the effects of time varying thrust misalignments, mass unbalance, and jet damping. Results are given in the form of equations for space­ referenced Euler angles, flight-path angles, body­ referenced attitude rates, and earth-referenced vehicle-trajectory coordinates. The method consists of dividing the problem into intervals during which the time-dependent variables are assumed constant at their mean interval value. In order to check this procedure, solutions for various interval sizes are compared with solutions obtained with numerical methods. Although the method is somewhat lengthy for accurate hand computation in most cases, it is readily programed for machine solutions. Probably more important, the general solutions give insight into the separate effects of the variables and, in many cases, can be quickly used to determine the approxi­mate ranges of the variables required for the desired solution to a given problem. In this respect, equations for determining maximum wobble have been derived for certain input conditions. The method has been shown to compare closely with the numerical solutions of two sample problems. The sample problems also illustrated the relatively large effect of pitch and yaw jet damping on body motions.

VACUUM↗
NASA NTRS1981

Development of a linear unsteady aerodynamic analysis for finite-deflection subsonic cascades

A linear unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, is being developed to predict the aerodynamic response to blade vibrations in the fan, compressor or turbine stages of modern jet engines. In previous work numerical solutions were restricted to cascades of sharp-edged blades aligned with the mean flow. Under the present effort the solution procedure has been extended to treat blades with rounded or blunt edges. As part of this effort an analytical model-problem study has been conducted to clarify the behavior of first-order perturbation solutions in the vicinity of airfoil edges. Further, a numerical approximation using concepts from singular perturbation theory has been developed to resolve the unsteady boundary value problem for cascades of blunt leading-edged blades. Numerical results for NACA 0012 cascades, including detailed results in the vicinity of a blade leading edge, are presented and evaluated.

Verdon, J. M.↗
NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗
NASA NTRS2006

On the Rates of Coronal Mass Ejections: Remote Solar and In Situ Observations

We compare the rates of coronal mass ejections (CMEs) as inferred from remote solar observations and interplanetary CMEs (ICMEs) as inferred from in situ observations at both 1 AU and Ulyssses from 1996 through 2004. We also distinguish between those ICMEs that contain a magnetic cloud (MC) and those that do not. While the rates of CMEs and ICMEs track each other well at solar minimum, they diverge significantly in early 1998, during the ascending phase of the solar cycle, with the remote solar observations yielding approximately 20 times more events than are seen at 1 AU. This divergence persists through 2004. A similar divergence occurs between MCs and non-MC ICMEs. We argue that these divergences are due to the birth of midlatitude active regions, which are the sites of a distinct population of CMEs, only partially intercepted by Earth, and we present a simple geometric argument showing that the CME and ICME rates are consistent with one another. We also acknowledge contributions from (1) an increased rate of high-latitude CMEs and (2) focusing effects from the global solar field. While our analysis, coupled with numerical modeling results, generally supports the interpretation that whether one observes a MC within an ICME is sensitive to the trajectory of the spacecraft through the ICME (i.e., an observational selection effect), one result directly contradicts it. Specifically, we find no systematic offset between the latitudinal origin of ICMEs that contain MCs at 1 AU in the ecliptic plane and that of those that do not.

Riley, Pete↗
NASA NTRSDate not supplied

Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel

To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.

Hypersonic Flow↗
NASA NTRSDate not supplied

Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel

To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.

Hypersonic Flow↗
NASA NTRS2024

Filament Wound Composite Analysis Using the NASA Multiscale Analysis Tool (NASMAT) and Finite Element Analysis

Fiber reinforced composite materials, owing to their tailorable thermomechanical and functional properties, allow one to produce a structure that is stronger, stiffer, and lighter than its metal counterpart while performing the same function, yielding a more efficient structure. This not only allows for the improvement of current technologies like aircraft structures, but also enables new technologies like gaseous hydrogen storage for mobility applications, which are otherwise impractical when manufactured using traditional metals due to weight and space restrictions or material embrittlement. However, the use of composites imposes greater design and manufacturing challenges on an engineer, since they are heterogenous, having a distinct structure across multiple length scale, behave generally anisotropically at the structural level and require complex manufacturing and processing methods. Capturing this complex behavior requires detailed numerical simulations, including the modeling of microstructural features like undulations, voids, and fiber alignment. In this paper, multiple repeating unit cells (RUCs), representing filament wound composites, are developed (via a script provided in the Appendix) and analyzed. The refinement of these RUCs is varied, and the analyses are performed using both the Abaqus finite element software and the NASA Multiscale Analysis Tool (NASMAT). A study is undertaken to compare the predicted effective elastic properties of the wound RUC to a laminate representation of the wound RUC, which neglects the undulations. Additionally, two different sets of periodic boundary conditions (PBCs) have been examined. One approximates the real boundary conditions using a standard approach and the other represents the PBCs exactly through the use of an offset. Lastly, a comparison of the local elastic stress fields is made among the models and approaches. Since wound structures are often approximated as laminated structures, it is important to understand the degree to which this assumption is valid, namely by first comparing the elastic constants and local elastic fields. This will provide, on the one hand, information concerning the bulk mechanical behavior and, on the other hand, insights concerning local load distributions and likely damage initiation sites.

Composites↗
NASA NTRS2025

Mind the Gap: Addressing Data Gaps and Assessing Noise Mismodeling in LISA

Due to the sheer complexity of the Laser Interferometer Space Antenna (LISA) space mission, data gaps arising from instrumental irregularities and/or scheduled maintenance are unavoidable. Focusing on merger-dominated massive black hole binary signals, we test the appropriateness of the Whittle-likelihood on gapped data in a variety of cases. From first principles, we derive the likelihood valid for gapped data in both the time and frequency domains. Cheap-to-evaluate proxies to p-p plots are derived based on a Fisher-based formalism, and verified through Bayesian techniques. Our tools allow to predict the altered variance in the parameter estimates that arises from noise mismodeling, as well as the information loss represented by the broadening of the posteriors. The result of noise mismodeling with gaps is sensitive to the characteristics of the noise model, with strong low-frequency (red) noise and strong high-frequency (blue) noise giving statistically significant fluctuations in recovered parameters. We demonstrate that the introduction of a tapering window reduces statistical inconsistency errors, at the cost of less precise parameter estimates. We also show that the assumption of independence between inter-gap segments appears to be a fair approximation even if the data set is inherently coherent. However, if one instead assumes fictitious correlations in the data stream, when the data segments are actually independent, then the resultant parameter recoveries could be inconsistent with the true parameters. The theoretical and numerical practices that are presented in this work could readily be incorporated into global-fit pipelines operating on gapped data.

LISA↗
NASA NTRS1973

On the relationship between acoustic energy density flux near the jet axis and far field acoustic intensity

By measurement and analysis, the relationship between the distribution of the outflow of acoustic energy over the jet boundary and the far-field intensity is considered. The physical quantity used is the gradient of the pressure evaluated on a geometrical plane at the smallest possible radial distance from the jet axis, but outside the vortical region, in the area where the homogeneous wave equation is reasonably well satisfied. The numerical and experimental procedures involved have been checked out by using a known source. Results indicate that the acoustic power output per unit length of the jet, in the region from which the sound emanates, peaks at approximately 9 diameters downstream. The acoustic emission for a jet Strouhal number of about 0.3 exceeds the emission for all other Strouhal numbers nearly everywhere along the measurement plane. However, the far-field peak intensity distribution obtained from the contribution of each station was found to depend on the spatial extent of the region where sound emanates from the jet, which, in turn, depends more on the far-field angle than on the Strouhal number.

Maestrello, L.↗
NASA NTRS2018

Multi-Scale Modeling of Low-Density Carbon-Phenolic Ablators

Protecting a spacecraft during atmospheric entry is one of highest risk factors that needs to be mitigated during design of a space exploration mission. At entry speeds from space, air turns into high-temperature plasma, and spacecraft Thermal Protection Systems (TPS) are needed to protect the vehicle payload. Modern successful material architectures of spacecraft shields use a porous carbon fiber substrate impregnated with phenolic as an ablator material. In the lecture, efforts to build a Predictive Material Modeling framework for porous ablators from micro-scale to macro-scale will be presented. Several numerical methods and techniques will be summarized that use voxelized images to compute geometrical properties of the porous substrate. These computed properties include porosity, specific surface area and tortuosity that are otherwise indirectly measured through experimental techniques. Direct simulation Monte Carlo (DSMC), a particle-based method for approximating the Boltzmann equation, is used to compute the permeability coefficient of the porous substrate based on its digitized representation. The method computes the flow within the microstructure, where the size of the pores may approach the mean-free-path of the flow. Finally, a high-fidelity model implemented in PATO (Porous-material Analysis Toolbox) is discussed, and some examples of ablative material response are presented, including for the first time 3D simulations of the full tiled heat shield for the Mars Science Laboratory (MSL) capsule.

Ablators↗
NASA NTRS1994

Thermoviscoplastic analysis of fibrous periodic composites by the use of triangular subvolumes

The non-linear viscoplastic behavior of fibrous periodic composites is analyzed by discretizing the unit cell into triangular subvolumes. A set of these subvolumes can be configured by the analyst to construct a representation for the unit cell of a periodic composite. In each step of the loading history the total strain increment at any point is governed by an integral equation which applies to the entire composite. A Fourier series approximation allows the incremental stresses and strains to be determined within a unit cell of the periodic lattice. The non-linearity arising from the viscoplastic behavior of the constituent materials comprising the composite is treated as a fictitious body force in the governing integral equation. Specific numerical examples showing the stress distributions in the unit cell of a fibrous tungsten/copper metal-matrix composite under viscoplastic loading conditions are given. The stress distribution resulting in the unit cell when the composite material is subjected to an overall transverse stress loading history perpendicular to the fibers is found to be highly heterogeneous, and typical homogenization techniques based on treating the stress and strain distributions within the constituent phases as homogeneous result in large errors under inelastic loading conditions.

Walker, Kevin P.↗
NASA NTRS2004

Fracture Analysis of Semi-Elliptical Surface Cracks in Ductile Materials

Accurate life assessment of structural components may require advanced life prediction criteria and methodologies. Structural components often exhibit several different types of defects, among the most prevalent being surface cracks. A semi-elliptical surface crack subjected to monotonic loading will exhibit stable crack growth until the crack has reached a critical size, at which the crack loses stability and fracture ensues (Newman, 2000). The shape and geometry of the flaw are among the most influential factors. When considering simpler crack configurations, such as a through-the-thickness crack, a three-dimensional (3D) geometry may be modeled under the approximation of two-dimensional (2D) plane stress or plane strain. The more complex surface crack is typically modeled numerically with the Finite Element Method (FEM). A semi-elliptical surface crack is illustrated in Figure 1-1.

Daniewicz, S. R.↗
NASA NTRS1977

Sonic boom analysis for high-altitude flight at high Mach number

Numerical programs for the computation of the flow field from the airplane at the flight altitude to the ground are presented. They take into account the nonlinear effects of high Mach number, the entropy change across the shock, the entropy and enthalpy variations in the atmospheric layer, and the gravitational effect. Extension of the programs for the axisymmetric problems to handle nonaxisymmetric terms is described. The asymmetry can be caused by the geometry of the body and the lift, and also by the fact that the variations in the atmospheric layer are two-dimensional. Numerical results for ground level signatures of several configurations at various flight conditions are presented and compared with existing approximate theories to demonstrate the influences of these nonlinear effects.

Ferri, A.↗