SEARCH · Search NASA

Results for “force method”

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

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1960

An Approximate Analytical Method for Studying Atmosphere Entry of Vehicles with Modulated Aerodynamic Forces

The dimensionless, transformed, nonlinear differential equation developed in NASA TR R-11 for describing the approximate motion and heating during entry into planetary atmospheres for constant aerodynamic coefficients and vehicle shape has been modified to include entries during which the aerodynamic coefficients and the vehicle shape are varied. The generality of the application of the original equation to vehicles of arbitrary weight, size, and shape and to arbitrary atmospheres is retained. A closed-form solution for the motion, heating, and the variation of drag loading parameter m/C(D)A has been obtained for the case of constant maximum resultant deceleration during nonlifting entries. This solution requires certain simplifying assumptions which do not compromise the accuracy of the results. The closed-form solution has been used to determine the variation of m/C(D)A required to reduce peak decelerations and to broaden the corridor for nonlifting entry into the earth's atmosphere at escape velocity. The attendant heating penalty is also studied.

Levy, Lionel L., Jr.↗

NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

NASA NTRS1995

Progress in HTS Trapped Field Magnets: J(sub c), Area, and Applications

Progress in trapped field magnets is reported. Single YBCO grains with diameters of 2 cm are made in production quantities, while 3 cm, 4 1/2 cm and 6 cm diameters are being explored. For single grain tiles: J(sub c) - 10,000 A/sq cm for melt textured grains; J(sub c) - 40,000 A/sq cm for light ion irradiation; and J(sub c) - 85,000 A/J(sub c) for heavy ion irradiation. Using 2 cm diameter tiles bombarded by light ions, we have fabricated a mini-magnet which trapped 2.25 Tesla at 77K, and 5.3 Tesla at 65K. A previous generation of tiles, 1 cm x 1 cm, was used to trap 7.0 Tesla at 55K. Unirradiated 2.0 cm tiles were used to provide 8 magnets for an axial gap generator, in a collaborative experiment with Emerson Electric Co. This generator delivered 100 Watts to a resistive load, at 2265 rpm. In this experiment, activation of the TFMs was accomplished by a current pulse of 15 ms duration. Tiles have also been studied for application as a bumper-tether system for the soft docking of spacecraft. A method for optimizing tether forces, and mechanisms of energy dissipation are discussed. A bus bar was constructed by welding three crystals while melt-texturing, such that their a,b planes were parallel and interleaved. The bus bar, of area approx. 2 sq cm, carried a transport current of 1000 amps, the limit of the testing equipment available.

Weinstein, Roy↗

NASA NTRS1999

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

NASA NTRS1999

Delamination Modeling of Composites for Improved Crash Analysis

Fleming, David C.↗

NASA NTRS2003

Application of a Split-Fiber Probe to Velocity Measurement in the NASA Research Compressor

A split-fiber probe was used to acquire unsteady data in a research compressor. The probe has two thin films deposited on a quartz cylinder 200 microns in diameter. A split-fiber probe allows simultaneous measurement of velocity magnitude and direction in a plane that is perpendicular to the sensing cylinder, because it has its circumference divided into two independent parts. Local heat transfer considerations indicated that the probe direction characteristic is linear in the range of flow incidence angles of +/- 35. Calibration tests confirmed this assumption. Of course, the velocity characteristic is nonlinear as is typical in thermal anemometry. The probe was used extensively in the NASA Glenn Research Center (GRC) low-speed, multistage axial compressor, and worked reliably during a test program of several months duration. The velocity and direction characteristics of the probe showed only minute changes during the entire test program. An algorithm was developed to decompose the probe signals into velocity magnitude and velocity direction. The averaged unsteady data were compared with data acquired by pneumatic probes. An overall excellent agreement between the averaged data acquired by a split-fiber probe and a pneumatic probe boosts confidence in the reliability of the unsteady content of the split-fiber probe data. To investigate the features of unsteady data, two methods were used: ensemble averaging and frequency analysis. The velocity distribution in a rotor blade passage was retrieved using the ensemble averaging method. Frequencies of excitation forces that may contribute to high cycle fatigue problems were identified by applying a fast Fourier transform to the absolute velocity data.

Lepicovsky, Jan↗

NASA NTRS2003

Inflow/Outflow Conditions for Unsteady Aerodynamics and Aeroacoustics in Nonuniform Flow

The effect of a nonuniform mean flow on the normal modes; the inflow/outflow nonreflecting boundary conditions; and the sound power are studied. The normal modes in an annular duct are computed using a spectral method in combination with a shooting method. The swirl causes force imbalance which couples the acoustic and vortical modes. The acoustic modes are distinguished from the vortical modes by their large pressure and small vorticity content. The mean swirl also produces a Doppler shift in frequency. This results in more counter-spinning modes cut-on at a given frequency than modes spinning with the swirl. Nonreflecting boundary conditions are formulated using the normal mode solutions. The inflow/outflow boundary conditions are implemented in a linearized Euler scheme and validated by computing the propagation of acoustic and vortical waves in a duct for a variety of swirling mean flows. Numerical results show that the evolution of the vortical disturbances is sensitive to the inflow conditions and the details of the wake excitations. All three components of the wake velocity must be considered to correctly compute the wake evolution and the blade upwash. For high frequencies, the acoustic-vortical mode coupling is weak and a conservation equation for the acoustic energy can be derived. Sound power calculations show significant mean flow swirl effects, but mode interference effects are small.

Atassi, Oliver V.↗

NASA NTRS2004

Ocean Carbon Flux, Transport, and Burial Within the Western and Eastern US Coastal Zones

This project has been to develop and apply a regional. eddy-resolving circulation and biogeochemistry model of both the western and eastern U.S. coastal regions, capable of simulating the processes that control the carbon cycle. Validation has been by statistical comparison with analyses from various satellite measurements, including those from EOS sensors, as well as from in situ measurements. Sensitivity studies were carried out to investigate how the coastal ecosystem and biogeochemical cycles respond to changes in climate, large-scale eutrophication from indus- trial pollution, and other anthropogenic induced changes. The research has been conducted in collaboration with research groups at UCLA. NASA/GSFC (Wallops), Rutgers, and SIO. Overall. the project was focused on several key modeling issues, each of which tie back into completing the primary task of developing a coastal carbon model for both the eastern and western US. coasts. Individual groups within the entire program are still collaborating to address these specific tasks. These include: implementation of the coupled circulation/biogeochemical model within the U.S. West Coast. including high-resolution, embedded subdomains for the Southern California Bight and Monterey Bay region; development of a biogeochemical model with resolved carbon, nitrogen and oxygen cycles; development of data assimilation techniques for use of satellite data sets; reconfiguration of the model domain to U.S. East Coast; development of coastal forcing fields: development of methods to compare the model against remotely sensed data; and, the test of model sensitivity to environmental conditions. Below, we present a summary of the progress made toward achieving these soak. Because this has been a multi-institutional, collaborative effort, we note the groups involved with particular activities.

McWilliams, James C.↗

NASA NTRS2000

Calcaneal loading during walking and running

PURPOSE: This study of the foot uses experimentally measured kinematic and kinetic data with a numerical model to evaluate in vivo calcaneal stresses during walking and running. METHODS: External ground reaction forces (GRF) and kinematic data were measured during walking and running using cineradiography and force plate measurements. A contact-coupled finite element model of the foot was developed to assess the forces acting on the calcaneus during gait. RESULTS: We found that the calculated force-time profiles of the joint contact, ligament, and Achilles tendon forces varied with the time-history curve of the moment about the ankle joint. The model predicted peak talocalcaneal and calcaneocuboid joint loads of 5.4 and 4.2 body weights (BW) during walking and 11.1 and 7.9 BW during running. The maximum predicted Achilles tendon forces were 3.9 and 7.7 BW for walking and running. CONCLUSIONS: Large magnitude forces and calcaneal stresses are generated late in the stance phase, with maximum loads occurring at approximately 70% of the stance phase during walking and at approximately 60% of the stance phase during running, for the gait velocities analyzed. The trajectories of the principal stresses, during both walking and running, corresponded to each other and qualitatively to the calcaneal trabecular architecture.

NASA Center ARC↗

NASA NTRS2004

Anisotropy of MHD Turbulence at Low Magnetic Reynolds Number

Turbulent fluctuations in MHD flows are known to become dimensionally anisotropic under the action of a sufficiently strong magnetic field. We consider the technologically relevant case of low magnetic Reynolds number and apply the method of DNS of forced flow in a periodic box to generate velocity fields. The analysis based on different anisotropy characteristics shows that the dimensional anisotropy is virtually scale-independent. We also find that, except for the case of very strong magnetic field, the flow is componentally isotropic. Its kinetic energy is practically uniformly distributed among the velocity components.

Zikanov, O.↗

NASA NTRS2008

A Dynamic Analysis of Hydrodynamic Wave Journal Bearings

The purpose of this paper is to study the dynamic behavior of a three-wave journal bearing using a transient approach. The transient analysis permits the determination of the rotor behavior after the fractional frequency whirl appears. The journal trajectory is determined by solving a set of nonlinear equations of motion using the Runge-Katta method. The fluid film forces are computed by integrating the transient Reynolds equation at each time step location of the shaft with respect to the bearing. Because of the large values of the rotational speeds, turbulent effects were included in the computations. The influence of the temperature on the viscosity was also considered. Numerical results were compared to experimenta1 results obtained at the NASA Glenn Research Center. Comparisons of the theoretical results with experimental data were found to be in good agreement. The numerical and experimental results showed that the fluid film of a three-wave journal bearing having a diameter of 30 mm, a length of 27 mm, and a wave amplitude ratio greater than 0.15 is stable even at rotational speeds of 60,000 RPM. For lower wave amplitude ratios, the threshold speed at which the fluid film becomes unstable depends on the wave amplitude and on the supply pocket pressure. Even if the fluid film is unstable, the wave bearing maintains the whirl orbit inside the bearing clearance.

Ene, Nicoleta M.↗

NASA NTRS2014

Data Mining of Historical Human Data to Assess the Risk of Injury due to Dynamic Loads

The NASA Occupant Protection Group is charged with ensuring crewmembers are protected during all dynamic phases of spaceflight. Previous work with outside experts has led to the development of a definition of acceptable risk (DAR) for space capsule vehicles. The DAR defines allowable probability rates for various categories of injuries. An important question is how to validate these probabilities for a given vehicle. One approach is to impact test human volunteers under projected nominal landing loads. The main drawback is the large number of subject tests required to attain a reasonable level of confidence that the injury probability rates would meet those outlined in the DAR. An alternative is to mine existing databases containing human responses to impact. Testing an anthropomorphic test device (ATD) at the same human‐exposure levels could yield a range of ATD responses that would meet DAR. As one aspect of future vehicle validation, the ATD could be tested in the vehicle's seat and suit configuration at nominal landing loads and compared with the ATD responses supported by the human data set. This approach could reduce the number of human‐volunteer tests NASA would need to conduct to validate that a vehicle meets occupant protection standards. METHODS: The U.S. Air Force has recorded hundreds of human responses to frontal, lateral, and spinal impacts at many acceleration levels and pulse durations. All of this data are stored on the Collaborative Biomechanics Data Network (CBDN), which is maintained by the Wright Patterson Air Force Base (WPAFB). The test device for human occupant restraint (THOR) ATD was impact tested on WPAFB's horizontal impulse accelerator (HIA) matching human‐volunteer exposures on the HIA to 5 frontal and 3 spinal loading conditions. No human injuries occurred as a result of these impact conditions. Peak THOR response variables for neck axial tension and compression, and thoracic‐spine axial compression were collected. Maximal chest deflection was determined from motion capture video of the impact test. HIC‐ 15 and BRIC were calculated from head acceleration responses. Given the number of human subjects for each test condition a confidence interval of injury probability will be obtained. RESULTS: Results will be discussed in terms of injury‐risk probability estimates based on the human data set evaluated. Also, gaps in the data set will be identified. These gaps could be one of two types. One is areas where additional THOR testing would increase the comparable human data set, thereby improving confidence in the injury probability rate. The other is where additional human testing would assist in obtaining information on other acceleration levels or directions. DISCUSSION: The historical human data showed validity of the THOR ATD for supplemental testing. The historical human data are limited in scope, however. Further data are needed to characterize the effects of sex, age, anthropometry, and deconditioning due to spaceflight on risk of injury

Wells, Jesica↗

NASA NTRS2017

Aeroservoelastic Modeling of Body Freedom Flutter for Control System Design

One of the most severe forms of coupling between aeroelasticity and flight dynamics is an instability called freedom flutter. The existing tools often assume relatively weak coupling, and are therefore unable to accurately model body freedom flutter. Because the existing tools were developed from traditional flutter analysis models, inconsistencies in the final models are not compatible with control system design tools. To resolve these issues, a number of small, but significant changes have been made to the existing approaches. A frequency domain transformation is used with the unsteady aerodynamics to ensure a more physically consistent stability axis rational function approximation of the unsteady aerodynamic model. The aerodynamic model is augmented with additional terms to account for limitations of the baseline unsteady aerodynamic model and to account for the gravity forces. An assumed modes method is used for the structural model to ensure a consistent definition of the aircraft states across the flight envelope. The X-56A stiff wing flight-test data were used to validate the current modeling approach. The flight-test data does not show body-freedom flutter, but does show coupling between the flight dynamics and the aeroelastic dynamics and the effects of the fuel weight.

aeroservoelastic↗

NASA NTRS2018

High Fidelity Collision Probabilities Estimated Using Brute Force Monte Carlo Simulations

The NASA Conjunction Assessment Risk Analysis team has implemented new software to estimate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm employs a brute force Monte Carlo (BFMC) method which differs from most other methods because it uses orbital states and covariances propagated from their orbit determination epoch times using the full set of Astrodynamics Support Workstation higher order theory models, including the High Accuracy Satellite Drag Model. This paper de-scribes the BFMC algorithm, presents comparisons of BFMC P (sub c) estimates to those calculated using other methods, and discusses the implications for conjunction risk assessment.

Probabilities↗

NASA NTRSDate not supplied

Predicting Char Yield of High-Temperature Resins

A simulation protocol is developed to predict the char yield of organic resins during high-temperature processing. Such in silico methods can help screen promising new formulations for advanced materials, but previously no chemistry-sensitive technique existed to predict the important experimental value of char yield. The method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to processing temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small byproduct molecules to mimic outgassing, and the application of high pressure to eliminate porosity and encourage graphitization. More than ten different resin chemistries are investigated, including arylacetylenes, cyanate esters, phthalonitriles, and polyimides, representing a diverse group of precursors with respect to initial cyclic content, heteroatoms and reactive groups. The protocol correctly predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. The properties of the resins during processing are compared, including outgassing products, morphology of the final chemical configurations, cyclic content and mechanical properties.

polymer simulations, molecular dynamics↗

NASA NTRSDate not supplied

Predicting Char Yield of High-Temperature Resins

A simulation technique has been developed for predicting the char yield of organic resins during high-temperature processing. In silico methods can aid in the screening of new advanced materials for a number of important properties, but no chemistry-sensitive protocol currently exists for predicting the important experimental value of char yield. The proposed method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to pyrolysis temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small by-product molecules to mimic outgassing, and the application of high pressure to encourage the formation of a dense, glassy network. Six different resin chemistries were investigated: an ethynyl, a phenylethynyl, a cyanate ester, a phthalonitrile, acrylonitrile, and adamantane. These candidates represent a diverse group of precursors with respect to initial cyclic content, presence of heteroatoms, and types of reactive groups. The protocol developed accurately predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. Several simulated properties of the carbonized structures are compared with experimental results, including outgassing products, morphology of the final chemical configurations, cyclic content, and mechanical properties.

polymer simulations, molecular dynamics↗

NASA NTRSDate not supplied

Predicting Char Yield of High-Temperature Resins

A simulation protocol is developed to predict the char yield of organic resins during high-temperature processing. Simulation protocol in silico methods can help screen promising new formulations for advanced materials, but previously no chemistry-sensitive technique existed to predict the important experimental value of char yield. The method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to processing temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small byproduct molecules to mimic outgassing, and the application of high pressure to eliminate porosity and encourage graphitization. More than ten different resin chemistries are investigated, including arylacetylenes, cyanate esters, phthalonitriles, and polyimides, representing a diverse group of precursors with respect to initial cyclic content, heteroatoms and reactive groups. The protocol correctly predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. The properties of the resins during processing are compared, including outgassing products, morphology of the final chemical configurations, cyclic content and mechanical properties.

polymer simulations, molecular dynamics↗

NASA NTRS2023

Fast Coupled Loads Analysis Method: Norton-Thevenin Receptance Coupling

A new method called Norton-Thevenin Receptance Coupling (NTRC) has been developed to perform coupled loads analysis (CLA). NTRC provides a tool that payload developers can use to obtain launch loads at a fraction of the cost of a CLA any time it is required in the payload design cycle. NTRC combines the frequency domain component coupling method of Receptance Coupling with the Norton and Thevenin theory used in force limiting to derive an alternate method for performing CLA.

Coupled Loads Analysis↗