Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear solid 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.

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

69 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Boundary layer roll circulations during FIRE

The probable mechanism underlying the development of boundary layer roll circulations are studied using wind and temperature profiles measured by the National Center for Atmospheric Research (NCAR) Electra during the stratocumulus phase of the First ISCCP Regional Experiment (FIRE). The expected, or preferred, roll orientations, horizontal wavelengths, and propagation periods are determined by finding the minimum values of the dynamic and thermodynamic forcing parameters, which here are the eddy Reynolds number (Re) and moist Rayleigh number (Ra sub m). These minimum values depend on the height z sub T of the capping temperature inversion and on the values of the Fourier coefficients of the background height-dependent vector wind profile. As input to our nonlinear spectral model, descent and ascent runs by the Electra provide for initial estimates of the inversion height and the wind profiles. In the first phase of the investigation presented here, a mechanism is said to be a probable contributor to the development of roll circulations within the stratocumulus-topped boundary layer if the modeled roll orientation and wavelengths agree with their observed values. Preliminary results using the 14-coefficient model of Haack-Hirschberg (1988) are discussed for the 7 July 1987 Electra Mission 188-A (Flight 5). This mission was flown across a sharp cloud boundary that was within a LANDSAT/SPOT scene. The stratocumulus deck was relatively solid in the eastern part of the scene, while there was a rapid decrease in cloud cover to scattered cumulus clouds aligned in streets to the west. These cloud streets were oriented nearly parallel to the mean wind direction in the layer, which was approximately 340 degrees. The hypothesis that roll circulations occurred in both the relatively clear and the cloudy regions is investigated using as model input a descent profile obtained in the relatively clear air and an ascent profile obtained in the cloudy air. Initial results for the clear air case are that the pure inflection point mode is not possible and the pure thermal mode was oriented 35 degrees to the right of the mean wind direction. The origin of this unacceptably large discrepancy between the observed and modeled results will be investigated further and the conclusions reported at the next FIRE workshop.

Shirer, Hampton N.↗
NASA NTRS1994

Propagation of whistler waves driven by fine structured ion beams in the magnetotail

In a previous paper, which examined the propagation of low-frequency whistler waves generated by ion beams in the Earth's plasma sheet boundary layer (PSBL), it was found that whistler waves driven in the PSBL are focused toward the central plasma sheet due to the global magnetotail inhomogeneities; this finding may help explain the observations of magnetic noise bursts in the tail (Burinskaya et al., 1993). In this paper the same phenomenon is examined, but this time a much more realistic model is used for the ion beam in the PSBL. While the PSBL has been modeled as a solid, homogeneous ion beams with a width of one Earth radius, observations and theoretical considerations have shown that PSBL ion beams actually have a decreasing velocity profile toward the plasma sheet and that the density of the beams within the PSBL can vary locally. We consider again the propagation and generation of electromagnetic waves but in the presence of fine structured ion beams in the PSBL. Our results show that whistler waves, generated quasi-parallel to the background magnetic field, can be trapped locally within small spatial regions where the ion beam density is enhanced compared to the density of the adjacent PSBL region. Wave spectra and nonlinear saturation mechanisms are discussed.

Burinskaya, T.↗
NASA NTRS2016

Growth Mechanism of Nanowires: Binary and Ternary Chalcogenides

Semiconductor nanowires exhibit very exciting optical and electrical properties including high transparency and a several order of magnitude better photocurrent than thin film and bulk materials. We present here the mechanism of nanowire growth from the melt-liquid-vapor medium. We describe preliminary results of binary and ternary selenide materials in light of recent theories. Experiments were performed with lead selenide and thallium arsenic selenide systems which are multifunctional material and have been used for detectors, acousto-optical, nonlinear and radiation detection applications. We observed that small units of nanocubes and elongated nanoparticles arrange and rearrange at moderate melt undercooling to form the building block of a nanowire. Since we avoided the catalyst, we observed self-nucleation and uncontrolled growth of wires from different places. Growth of lead selenide nanowires was performed by physical vapor transport method and thallium arsenic selenide nanowire by vapor-liquid-solid (VLS) method. In some cases very long wires (>mm) are formed. To achieve this goal experiments were performed to create situation where nanowires grew on the surface of solid thallium arsenic selenide itself.

Singh, N. B.↗
NASA NTRSDate not supplied

Flight Performance and Stability of Space Launch System Core Stage Thrust Vector Control

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. The Core Stage TVC shares vehicle control authority with the SLS 5-segment Solid Rocket Boosters (SRBs) during boost phase flight, and is the sole means of vehicle flight control during in exoatmospheric flight following SRB separation. TVC responses during Green Run Hot Fire (GRHF) testing revealed that the TVC did not meet its performance specifications. Step and frequency responses exhibited unexpected departures from prior laboratory data and modeled behavior. Post-test analysis determined that the characteristics of the structure and gimbal friction are significantly influenced by the thrust-loaded conditions, and the command avionics exhibited a small but important gain nonlinearity. Using the available test data, the design team augmented the flight control TVC models to bound the observed results and include the additional fidelity needed for vehicle flight control analysis so as to build sufficient rationale for flight certification. Prior to the Green Run tests, “simplex” linear models typically used for flight control analysis did not include gimbal friction and other nonlinearities owing to long-standing assumptions that these effects were negligible in the Shuttle Orbiter TVC system. Following the Green Run findings, simulation analysis of the flight dynamics in the time and frequency domain revealed the propensity for a flight control limit cycle oscillation (LCO) if friction and structural compliances fell near the edges of test-predicted bounds. While the “most probable” models did not predict an in-flight LCO, the SLS Program conservatively proceeded with a system-wide evaluation and ultimate acceptance of the possibility for a small amplitude, low-frequency TVC LCO in flight. A final validation of the extensive test and modeling effort occurred when the first flight of SLS successfully demonstrated the fully integrated performance of the vehicle’s TVC system This paper is the final installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, the development of flight rationale in light of the TVC responses observed in Green Run is discussed, along with a review of the flight telemetry illustrating the correlation of the preflight predictions with the observed performance.

John H. Wall↗
NASA NTRS1972

The effects of various cure cycles upon the viability of Bacillus subtillis var. niger spores within solid propellant.

Saturethane solid propellant was used in all tests. The spore inoculum was evenly distributed in the propellant. Samples weighing approximately 5 g were aseptically removed, placed into curing ovens, and exposed to cure temperatures. Initial tests were conducted at 82 and 93 C. Analysis of the obtained data indicated that the survivor curves were not linear. Exposure of the inoculated propellant samples to 93 C reduced the initial population to less than 0.01% in about 20 hours. At 82 C, approximately 168 hours were required for a similar reduction. Tests involving curing temperatures of 105 and 115 C were also conducted. It is pointed out that changes in the mechanism of spore inactivation due to chemical and physical changes in the propellant could account for the nonlinear survivor curves obtained.

Brewer, W. A.↗
NASA NTRS1986

Anisotropic constitutive model for nickel base single crystal alloys: Development and finite element implementation

A tool for the mechanical analysis of nickel base single crystal superalloys, specifically Rene N4, used in gas turbine engine components is developed. This is achieved by a rate dependent anisotropic constitutive model implemented in a nonlinear three dimensional finite element code. The constitutive model is developed from metallurigical concepts utilizing a crystallographic approach. A non Schmid's law formulation is used to model the tension/compression asymmetry and orientation dependence in octahedral slip. Schmid's law is a good approximation to the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response, and strain rate sensitivity of these alloys. Methods for deriving the material constants from standard tests are presented. The finite element implementation utilizes an initial strain method and twenty noded isoparametric solid elements. The ability to model piecewise linear load histories is included in the finite element code. The constitutive equations are accurately and economically integrated using a second order Adams-Moulton predictor-corrector method with a dynamic time incrementing procedure. Computed results from the finite element code are compared with experimental data for tensile, creep and cyclic tests at 760 deg C. The strain rate sensitivity and stress relaxation capabilities of the model are evaluated.

Dame, L. T.↗
NASA NTRS2002

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Freed, Alan D.↗
NASA NTRS1987

Variational approach to probabilistic finite elements

Probabilistic finite element method (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties, and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗
NASA NTRS1991

Variational approach to probabilistic finite elements

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗
NASA NTRS1989

Assessment of SPAR elements and formulation of some basic 2-D and 3-D elements for use with testbed generic element processor

The initial Computational Structural Mechanics (CSM) Testbed was based on Level 13 of the SPAR finite element computer program. Until recently, the element library of the Testbed has been limited to those elements in Level 13 of SPAR. The development of a generic element processor has enabled element researchers to develop, implement and assess element formulations with relative ease. An assessment of new elements as well as the existing SPAR Level 13 elements has revealed some definite shortcomings with the SPAR Level 13 2-D and 3-D elements. The SPAR S81 solid element does not pass the patch test problem proposed by MacNeal-Harder. These deficiencies are identified here. The 2-D elements, however, seem to perform well taking into account the limitations imposed by the theory used to formulate them, (i.e., thin plates only). Common deficiencies of the 2-D and 3-D elements in SPAR have to do with their adaptability to the nonlinear analysis utilities developed by Lockheed Palo Alto Research Lab. Also, the EFIL format of the SPAR element data does not conform to the standard format of the Testbed.

Aminpour, Mohammad A.↗
NASA NTRS1992

Venus tectonic styles and crustal differentiation

Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.

Kaula, W. M.↗
NASA NTRS1995

Venus' Free Obliquity

The predicted orientation of Venus' rotation axis relative to its orbit ran be uniquely determined given knowledge of its J2 gravity coefficient and polar moment of inertia C if its free obliquity is fully damped. This assumption seems warranted given the dominant damping mechanism: turbulent fluid friction at a core mantle boundary (CMB). This skin friction results from differential obliquity of mantle and core spin axes, and the associated damping rate could be as short as 1/10(exp 6) year. However, the observed pole orientation indicates a free obliquity amplitude epsilon approx. = 2.1 deg compared with a nominal forced amplitude of 0.5 deg. There are two plausible explanations. The most likely is that the observed obliquity is a tidally evolved end state in which core friction, modulated by CMB ellipticity and core obliquity amplitude, counterbalances solid and atmospheric tidal torques. This concept is similar to the explanation for the retrograde spin omega as an end state in which solid and atmospheric thermal tidal torques balance at the present spin rate because of the omega(sup -1) dependence of the axial thermal torque. Weaker effects such as solid and thermal tides can then compete with core friction and for plausible models, their sum tends to increase free obliquity. The obliquity balance is controlled by the nonlinear (and nearly quadratic) dependence of the CMB turbulent "skin friction" torque on obliquity.Core ellipticity also has a profound effect on tidal evolution of Venus obliquity.

Yoder, Charles F.↗
NASA NTRS1993

Every Good Virtue You Ever Wanted in a Q-switched Solid-state Laser and More: Monolithic, Diode-pumped, Self-q-switched, Highly Reproducible, Diffraction-limited Nd:yag Laser

The applications of Q-switched lasers are well known, for example, laser radar, laser remote sensing, satellite orbit determination, Moon orbit and 'moon quake' determination, satellite laser communication, and many nonlinear optics applications. Most of the applications require additional properties of the Q-switched lasers, such as single-axial and/or single-transverse mode, high repetition rate, stable pulse shape and pulse width, or ultra compact and rugged oscillators. Furthermore, space based and airborne lasers for lidar and laser communication applications require efficient, compact, lightweight, long-lived, and stable-pulsed laser sources. Diode-pumped solid-state lasers (DPSSL) have recently shown the potential for satisfying all of these requirements. We will report on the operating characteristics of a diode-pumped, monolithic, self-Q-switched Cr,Nd:YAG laser where the chromium ions act as a saturable absorber for the laser emission at 1064 nm. The pulse duration is 3.5 ns and the output is highly polarized with an extinction ratio of 700:1. It is further shown that the output is single-longitudinal-mode with transform-limited spectral line width without pulse-to-pulse mode competition. Consequently, the pulse-to-pulse intensity fluctuation is less than the instrument resolution of 0.25 percent. This self-stabilization mechanism is because the lasing mode bleaches the distributed absorber and establishes a gain-loss grating similar to that used in the distributed feedback semiconductor lasers. A repetition rate above 5 KHz has also been demonstrated. For higher power, this laser can be used for injection seeding an amplifier (or amplifier chain) or injection locking of a power oscillator pumped by diode lasers. We will discuss some research directions on the master oscillator for higher output energy per pulse as well as how to scale the output power of the diode-pumped amplifier(s) to multi-kilowatt average power.

Chen, Y. C.↗
NASA NTRS2017

Fully Coupled Aero-Thermochemical-Elastic Simulations of an Eroding Graphite Nozzle

A multiphysics simulation capability has been developed that incorporates mutual interactions between aerodynamics, structural response from aero/thermal loading, ablation/pyrolysis, heating, and surface-to-surface radiation to perform high-fidelity, fully coupled aerothermoelastic ablation simulations, which to date had been unattainable. The multiphysics framework couples CHAR (a 3-D implicit charring ablator solver), Loci/CHEM (a computational fluid dynamics solver for high-speed chemically reacting flows), and Abaqus (a nonlinear structural dynamics solver) to create a fully coupled aerothermoelastic charring ablative solver. The solvers are tightly coupled in a fully integrated fashion to resolve the effects of the ablation pyrolysis and charring process and chemistry products upon the flow field, the changes in surface geometry due to recession upon the flow field, and thermal-structural analysis of the body from the induced aerodynamic heating from the flow field. The multiphysics framework was successfully demonstrated on a solid rocket motor graphite nozzle erosion application. Comparisons were made with available experimental data that measured the throat erosion during the motor firing. The erosion data is well characterized, as the test rig was equipped with a windowed nozzle section for real-time X-ray radiography diagnostics of the instantaneous throat variations for deducing the instantaneous erosion rates. The nozzle initially undergoes a nozzle contraction due to thermal expansion before ablation effects are able to widen the throat. A series of parameters studies were conducted using the coupled simulation capability to determine the sensitivity of the nozzle erosion to different parameters. The parameter studies included the shape of the nozzle throat (flat versus rounded), the material properties, the effect of the choice of turbulence model, and the inclusion or exclusion of the mechanical thermal expansion. Overall, the predicted results match the experiment very well, and the predictions were able to bound the data within acceptable limits.

Blades, E. L.↗
NASA NTRS1999

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.

Non-NASA Center↗