SEARCH · Search NASA

Results for “Beam code development”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2012

GERMcode: A Stochastic Model for Space Radiation Risk Assessment

A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and high charge and energy (HZE) particles that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of HZE particles in tissue and shielding materials is made with a stochastic approach that includes both particle track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections. For NSRL applications, the GERMcode evaluates a set of biophysical properties, such as the Poisson distribution of particles or delta-ray hits for a given cellular area and particle dose, the radial dose on tissue, and the frequency distribution of energy deposition in a DNA volume. By utilizing the ProE/Fishbowl ray-tracing analysis, the GERMcode will be used as a bi-directional radiation transport model for future spacecraft shielding analysis in support of Mars mission risk assessments. Recent radiobiological experiments suggest the need for new approaches to risk assessment that include time-dependent biological events due to the signaling times for activation and relaxation of biological processes in cells and tissue. Thus, the tracking of the temporal and spatial distribution of events in tissue is a major goal of the GERMcode in support of the simulation of biological processes important in GCR risk assessments. In order to validate our approach, basic radiobiological responses such as cell survival curves, mutation, chromosomal aberrations, and representative mouse tumor induction curves are implemented into the GERMcode. Extension of these descriptions to other endpoints related to non-targeted effects and biochemical pathway responses will be discussed.

Kim, Myung-Hee Y.↗

NASA NTRS1975

ISOFINEL: Isoparametric finite element code for elastic analysis of two-dimensional bodies

A formulation is presented for the development of a finite element program for the elastic analysis of two-dimensional bodies using the eight-node isoparametric quadrilateral. The program solves for both plane stress and plane strain problems. The finite element formulation based on the isoparametric displacement functions is presented. The program structure is given in the form of flow diagrams with descriptions of the numerical procedure used to obtain the element stiffness matrix, and the solution method employed to solve for nodal displacements. Three numerical examples (a plate under uniaxial tension, a plate under pure shear, and a beam under pure bending) are presented to illustrate the capability and limitations of the element implementation. The first problem is solved exactly by the element, as predicted by the form of its displacement functions. In the other two problems the accuracy of the solution is highly dependent upon the slenderness of the element, the number of elements in the map, and the numerical integration scheme used to build the element stiffness matrix.

Marino, C.↗

NASA NTRS1990

A rate equation approach to gain saturation effects in laser mode calculations

Space exploration and research require large amounts of power. Solar pumped laser systems have been shown to have the potential for meeting the performance requirements necessary for power transmission in space. The successful design of laser systems for these applications will depend, in part, on having a clear understanding of the development of the dynamical processes in the laser cavity and on the effects that changes in physical and design parameters have on laser performance. In particular, it is necessary to know the amplitude and phase distributions of the laser beam at the output aperture when steady state operation is achieved in order to determine the far-field power distribution. The output from the laser will depend on the active medium, the optical environment of the active material and on the gain distribution in the active region as laser action builds up and reaches steady state. An important component in the design process will be a realistic model of the active laser cavity. A computer model of the laser cavity, based on Huygens' principle was developed. The code calculates the amplitude and phase of an optical wave reflected back and forth between the mirrors of a laser cavity. The original code assumes a gain distribution which does not change with the buildup of oscillations in the cavity. A step in the direction of realism is the inclusion of a saturable gain medium in the cavity. The objective is to incorporate saturation effects into the existing computer model.

Roberts, Lila F.↗

NASA NTRS1994

Computer simulation of Alfven waves and double layers along auroral magnetic field lines

A plasma simulation has been developed to model interactions between inertial Alfven waves and double layers and to investigate their relative contributions at auroral particle acceleration. We use a novel one-dimensional particle-in-cell code, with periodic boundary conditions, to model the nonlinear excitation of current-driven weak double layers via free energy supplied by an inertial Alfven wave. Analysis of the simulation output shows that double layers are not the agent primarily responsible for electron acceleration. Rather, the inertial Alfven wave accelerates groups of electrons into a steepening beam as it encounters them. As the beam electrons reenter the main distribution, decelerated by anomalous resistive effects, they are replaced by electrons farther downstream. Hence, the particles do not free-stream over the length of the channel. Furthermore, this wave action persists even when the system is linearly stable to ion-acoustic modes, precluding the possibility that this behavior is brought about by the formation of ion-acoustic double layers.

Silberstein, M.↗

NASA NTRS2018

Numerical Studies of Hall Thruster Acceleration Region Electron Transport

Azimuthal direction 1-D and azimuthal and axial direction 2-D Particle in Cell (PIC) simulations of the Hall thruster acceleration and near-plume regions are presented. The axial extent of the 2-D computational domain is based on the ExB drift parameter, 𝑫𝑬𝒙𝑩, the number of times around the channel an electron ExB drifts before being scattered. With only high energy main beam ions, the simulations produced a broad potential profile in the acceleration region. The addition of a small amount of ionization in the acceleration domain at a rate consistent with Hall thruster parameters causes the potential profile to develop two distinct regions. In the upstream portion of the computational domain, the potential profile is very steep, typical of Hall thruster acceleration regions. In the downstream portion, the potential profile is almost flat, like potentials in Hall thruster near plume regions. This bifurcation happens naturally; there is no other change in code parameters between the two results. Without the addition of a small amount of ionization in the acceleration region, the axial electron transport is very low. With ionization, the calculated ratio of electron current to ion current is the same order of those measured in Hall thrusters. This suggests that in regions with lower energy electrons and ions, plasma waves are effective in allowing crossfield transport, and axial electric fields are small. In the acceleration region, the plasma waves are less effective at scattering the energetic electrons, and high axial electric fields are needed to produce the observed cross-field electron transport.

Ortega, Alejandro Lopez↗

NASA NTRS2015

Multi-Level Experimental and Analytical Evaluation of Two Composite Energy Absorbers

Two composite energy absorbers were developed and evaluated at NASA Langley Research Center through multi-level testing and simulation performed under the Transport Rotorcraft Airframe Crash Testbed (TRACT) research program. A conical-shaped energy absorber, designated the conusoid, was evaluated that consisted of four layers of hybrid carbon-Kevlar plain weave fabric oriented at [+45 deg/-45 deg/-45 deg/+45 deg] with respect to the vertical, or crush, direction. A sinusoidal-shaped energy absorber, designated the sinusoid, was developed that consisted of hybrid carbon-Kevlar plain weave fabric face sheets, two layers for each face sheet oriented at +/-45deg with respect to the vertical direction and a closed-cell ELFOAM P200 polyisocyanurate (2.0-lb/cu ft) foam core. The design goal for the energy absorbers was to achieve average floor-level accelerations of between 25- and 40-g during the full-scale crash test of a retrofitted CH-46E helicopter airframe, designated TRACT 2. Variations in both designs were assessed through dynamic crush testing of component specimens. Once the designs were finalized, subfloor beams of each configuration were fabricated and retrofitted into a barrel section of a CH-46E helicopter. A vertical drop test of the barrel section was conducted onto concrete to evaluate the performance of the energy absorbers prior to retrofit into TRACT 2. The retrofitted airframe was crash tested under combined forward and vertical velocity conditions onto soil, which is characterized as a sand/clay mixture. Finite element models were developed of all test articles and simulations were performed using LS-DYNA, a commercial nonlinear explicit transient dynamic finite element code. Test-analysis results are presented for each energy absorber as comparisons of time-history responses, as well as predicted and experimental structural deformations and progressive damage under impact loading for each evaluation level.

Jackson, Karen E.↗

NASA NTRS1996

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗

NASA NTRS1986

The LATDYN user's manual

The LATDYN User's Manual presents the capabilities and instructions for the LATDYN (Large Angle Transient DYNamics) computer program. The LATDYN program is a tool for analyzing the controlled or uncontrolled dynamic transient behavior of interconnected deformable multi-body systems which can undergo large angular motions of each body relative other bodies. The program accommodates large structural deformation as well as large rigid body rotations and is applicable, but not limited to, the following areas: (1) development of large flexible space structures; (2) slewing of large space structure components; (3) mechanisms with rigid or elastic components; and (4) robotic manipulations of beam members. Presently the program is limited to two dimensional problems, but in many cases, three dimensional problems can be exactly or approximately reduced to two dimensions. The program uses convected finite elements to affect the large angular motions involved in the analysis. General geometry is permitted. Detailed user input and output specifications are provided and discussed with example runstreams. To date, LATDYN has been configured for CDC/NOS and DEC VAX/VMS machines. All coding is in ANSII-77 FORTRAN. Detailed instructions regarding interfaces with particular computer operating systems and file structures are provided.

Housner, J. M.↗

NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

NASA NTRS1999

Determination of Nonlinear Stiffness Coefficients for Finite Element Models with Application to the Random Vibration Problem

In this paper, a method for obtaining nonlinear stiffness coefficients in modal coordinates for geometrically nonlinear finite-element models is developed. The method requires application of a finite-element program with a geometrically non- linear static capability. The MSC/NASTRAN code is employed for this purpose. The equations of motion of a MDOF system are formulated in modal coordinates. A set of linear eigenvectors is used to approximate the solution of the nonlinear problem. The random vibration problem of the MDOF nonlinear system is then considered. The solutions obtained by application of two different versions of a stochastic linearization technique are compared with linear and exact (analytical) solutions in terms of root-mean-square (RMS) displacements and strains for a beam structure.

Muravyov, Alexander A.↗

NASA NTRS1987

Modeling and control of beam-like structures

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for beam-like continua. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

NASA NTRS1992

Arcing in space structures in low Earth orbit

This report describes results of an experimental and theoretical program to investigate arcing of structures containing dielectric and conducting materials when they are biased negatively with respect to a plasma. An argon ion source generated Ar(+) ions of directed energy 20 to 40 eV and density approximately 10(exp 7) cm(exp -3) that impinged upon samples containing a dielectric material on top of a negatively biased Kovar plate. Arcing events were studied for bias voltages between -300 and -1000V with respect to the ion beam. The samples were Dow Corning 93-500 adhesive on Kovar, fused silica cover slips bonded on Kovar, and silicon solar cells mounted on Kovar. Measurements of discharge current, Kovar plate voltage, and radiation from the arc versus time were carried out. Microsecond duration exposure images and optical spectra in the 0.24 to 0.40 micron band were also acquired during arcing events. Arcing events were found to be associated with exposed adhesive and means were found to eliminate arcing altogether. The charging of a silica cover plate and the fields around the plate were calculated using a particle-in-cell code. Models were developed to explain the ignition of the arc and the physical processes occurring during the discharge.

Upschulte, B. L.↗

NASA NTRS2012

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engine turbines is to decompose a computational fluid dynamics (CFD).generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies using cyclically symmetric structural dynamic models. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non ]harmonic excitation sources that become present in complex flows. This complex content can only be captured by a CFD flow field encompassing at least an entire revolution. A substantial development effort to create a series of software programs to enable application of the 360 degree forcing function in a frequency response analysis on cyclic symmetric models has been completed (to be described in a future paper), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements and the disk modeled with plates (using the finite element code MSC/NASTRAN). The focus of this model is to be representative of response of realistic bladed disks, and so the dimensions are roughly equivalent to the new J2X rocket engine 1st stage fuel pump turbine. The simplicity of the model allows the CFD load to be able to be readily applied, along with analytical and experimental variations in both the temporal and spatial fourier components of the excitation. In addition, this model is a first step in identifying response differences between transient and frequency forced response analysis techniques. The second phase assesses this difference for a much more realistic solid model of a bladed-disk in order to evaluate the effect of the spatial variation in loading on blade dominated modes. Neither research on the accuracy of the frequency response method when used in this context or a comprehensive study of the effect of test-observed variation on blade forced response have been found in the literature, so this research is a new contribution to practical structural dynamic analysis of gas turbines. The primary excitation of the upstream nozzles interacts with the blades on fuel pump of the J2X causes the 5th Nodal diameter modes to be excited, as explained by Tyler and Sofrin1, so a modal analysis was first performed on the beam/plate model and the 5ND bladed-disk mode at 40167 hz was identified and chosen to be the one excited at resonance (see figure 1). The first forced response analysis with this model focuses on identifying differences between frequency and transient response analyses. A hypothesis going into the analysis was that perhaps the frequency response was enforcing a temporal periodicity that did not really exist, and so therefore it would overestimate the response. As high dynamic response was a considerable source of stress in the J2X, examining this concept could potentially be beneficial for the program.

Brown, Andrew M.↗

NASA NTRS1988

Modal cost analysis for simple continua

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations, it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for simple continua such as beam-like structures. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

NASA NTRS1988

Control of optical systems

Some of the current and planned activities at the Air Force Systems Command in structures and controls for optical-type systems are summarized. Many of the activities are contracted to industry; one task is an in-house program which includes a hardware test program. The objective of the in-house program, referred to as the Aluminum Beam Expander Structure (ABES), is to address issues involved in on-orbit system identification. The structure, which appears similar to the LDR backup structure, is about 35 feet tall. The activity to date has been limited to acquisition of about 250 hours of test data. About 30 hours of data per excitation force is gathered in order to obtain sufficient data for a good statistical estimate of the structural parameters. The development of an Integrated Structural Modeling (ISM) computer program is being done by Boeing Aerospace Company. The objective of the contracted effort is to develop a combined optics, structures, thermal, controls, and multibody dynamics simulation code.

Founds, D.↗

NASA NTRS1990

Solar flare gamma-ray line shapes

A computer code has been developed which is used to calculate ab initio the laboratory shapes and energy shifts of gamma-ray lines from (C-12)(p, gamma/4.438/)p-prime(C-12) and (O-16)(p, gamma/6.129/)p-prime(O-16) reactions and to calculate the expected shapes of these lines from solar flares. The sensitivity of observable solar flare gamma-ray line shapes to the directionality of the incident particles is investigated for several projectile angular distributions. Shapes of the carbon and oxygen lines are calculated assuming realistic proton energy spectra for particles in circular orbits at the mirror points of magnetic loops, for particle beams directed downward into the photosphere, and for isotropic particle distributions. Line shapes for flare sites near the center of the sun and on the limb are shown for both thin-target and thick-target interaction models.

Werntz, C.↗

NASA NTRS2001

A Bell-Curved Based Algorithm for Mixed Continuous and Discrete Structural Optimization

An evolutionary based strategy utilizing two normal distributions to generate children is developed to solve mixed integer nonlinear programming problems. This Bell-Curve Based (BCB) evolutionary algorithm is similar in spirit to (mu + mu) evolutionary strategies and evolutionary programs but with fewer parameters to adjust and no mechanism for self adaptation. First, a new version of BCB to solve purely discrete optimization problems is described and its performance tested against a tabu search code for an actuator placement problem. Next, the performance of a combined version of discrete and continuous BCB is tested on 2-dimensional shape problems and on a minimum weight hub design problem. In the latter case the discrete portion is the choice of the underlying beam shape (I, triangular, circular, rectangular, or U).

Kincaid, Rex K.↗

NASA NTRS1999

Recent Development in the Analytical Investigation of Shipboard Rotorcraft Engage and Disengage Operations

This paper presents an overview of recent developments in an effort to predict transient aeroelastic rotor response during shipboard engage and disengage sequences. The blade is modeled as an elastic beam undergoing in flap, lag, extension and torsion. The blade equations of motion are formulated using Hamilton's principle and they are spatially discretized using the finite element method. The discretized blade equations of motion are integrated for a specified rotor speed run-up or run-down profile. Blade element theory is used to calculate quasi-steady or unsteady aerodynamic loads in linear and nonlinear regimes. The analysis is capable of simulating both articulated, hingeless, and gimballed rotor systems. Validation of the rotor code is discussed, including correlation with droop stop impact tests and wind tunnel experiments. Predictions of safe engagement and disengagement envelopes, limited by excessive blade tip deflections or hub moments, are presented. Future directions of study are also discussed.

Smith, Edward C.↗