Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer simulation model 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.

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

At least 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

The Digital Astronaut Project Computational Bone Remodeling Model (Beta Version) Bone Summit Summary Report

Under the conditions of microgravity, astronauts lose bone mass at a rate of 1% to 2% a month, particularly in the lower extremities such as the proximal femur [1-3]. The most commonly used countermeasure against bone loss in microgravity has been prescribed exercise [4]. However, data has shown that existing exercise countermeasures are not as effective as desired for preventing bone loss in long duration, 4 to 6 months, spaceflight [1,3,5,6]. This spaceflight related bone loss may cause early onset of osteoporosis to place the astronauts at greater risk of fracture later in their lives. Consequently, NASA seeks to have improved understanding of the mechanisms of bone demineralization in microgravity in order to appropriately quantify this risk, and to establish appropriate countermeasures [7]. In this light, NASA's Digital Astronaut Project (DAP) is working with the NASA Bone Discipline Lead to implement well-validated computational models to help predict and assess bone loss during spaceflight, and enhance exercise countermeasure development. More specifically, computational modeling is proposed as a way to augment bone research and exercise countermeasure development to target weight-bearing skeletal sites that are most susceptible to bone loss in microgravity, and thus at higher risk for fracture. Given that hip fractures can be debilitating, the initial model development focused on the femoral neck. Future efforts will focus on including other key load bearing bone sites such as the greater trochanter, lower lumbar, proximal femur and calcaneus. The DAP has currently established an initial model (Beta Version) of bone loss due to skeletal unloading in femoral neck region. The model calculates changes in mineralized volume fraction of bone in this segment and relates it to changes in bone mineral density (vBMD) measured by Quantitative Computed Tomography (QCT). The model is governed by equations describing changes in bone volume fraction (BVF), and rates of changes in bone cell populations that remove and replace bone in packets within the bone region. The DAP bone model is unique in several respects. In particular in takes former models of volume fraction changes one step higher in fidelity and separates BVF into separate equations for mineralized and osteoid volume fractions governed by a mineralization rate. This more closely follows the physiology of the remodeling unit cycles where bone is first resorbed and then followed by the action of osteoblasts to lay down collagen matrix which eventually becomes mineralized. In another respect, the modules allow the functional description of the time rate of change of other parameters and variables in the model during a computational simulation. More detailed description of the model, preliminary validation results, current limitation and caveats, and planned advancements are provided in sections 2 through 5. The DAP bone model is being developed primarily as a research tool, and not as a clinical tool like QCT. Even if it transitions to a clinical tool, it is not intended to replace QCT or any other clinical tool. Moreover, the DAP bone model does not predict bone fracture. Its purpose is to provide valuable additional data via "forward prediction" simulations for during and after spaceflight missions to gain insight on, (1) mechanisms of bone demineralization in microgravity, and (2) the volumetric changes at the various bone sites in response to in-flight and post-flight exercise countermeasures. This data can then be used as input to the Keyak [8] (or equivalent) FE analysis method to gain insight on how bone strength may change during and after flight. This information can also be useful to help optimize exercise countermeasure protocols to minimize changes in bone strength during flight, and improve regain of bone strength post-flight. To achieve this goal, the bone model will be integrated with DAP's exercise countermeasure models to simulate the effect of exercise prescriptions on preserving bone. More specifically, the model will accept loading history due to muscle and joint force on bone and produce quantified remodeling within the bone region under influence of the applied stress. Furthermore, because they tend to respond differently, the bone remodeling model includes both trabecular bone and cortical bone.

Pennline, James↗
NASA NTRSDate not supplied

Implementation of Active Sites in DSMC to Capture Pitting of Oxidizing Carbon Materials

In this work we demonstrate a newly developed capability to capture pitting of carbon fibers in DSMC simulations, specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. State-of-the-art reactive surface models in DSMC compute collision dependent carbon consumption rates (usually through desorption of CO) based on a set of surface reactions that has been derived from molecular beam experiments. The reactivity on each carbon surface element is constant in those models, such that the carbon surface recedes uniformly as a result of ablation. However, it is well known that in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. These defective sites have a much higher reactivity than the average sites (2-3 orders of magnitude) and are first to react during ablation leading to its removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these ”active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. Within this framework, we have implemented the capability of a single surface having multiple site sets with different reactivities. Using this feature, we can simulate the presence of active sites on carbon surfaces, whose reactivity is much greater than an average site as a result of defects. We have implemented the active site fraction as a property of surface elements within SPARTA, which is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface, and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we are able to capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
NASA NTRS1999

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗
NASA NTRS1972

Utilization of a CRT display light pen in the design of feedback control systems

A hierarchical structure of the interlinked programs was developed to provide a flexible computer-aided design tool. A graphical input technique and a data structure are considered which provide the capability of entering the control system model description into the computer in block diagram form. An information storage and retrieval system was developed to keep track of the system description, and analysis and simulation results, and to provide them to the correct routines for further manipulation or display. Error analysis and diagnostic capabilities are discussed, and a technique was developed to reduce a transfer function to a set of nested integrals suitable for digital simulation. A general, automated block diagram reduction procedure was set up to prepare the system description for the analysis routines.

Thompson, J. G.↗
NASA NTRS1977

Dynamics and control simulation of the Spacelab Experiment Pointing Mount

Computer simulations were developed to evaluate the performance of four Experiment Pointing Mounts (EPM) being considered for Spacelab experiments in the 1980-1990 time frame. The system modeled compromises a multibody system consisting of the shuttle, a mechanical isolation device, the EPM, celestial and inertial sensors, bearings, gimbal torque motors and associated nonlinearities, the experiment payload, and control and estimator algorithms. Each mount was subjected to a common disturbance (shuttle vernier thruster firing and man push off) and command (stellar pointing or solar raster scan) input. The fundamental limitation common to all mounts was found to be sensor noise. System dynamics and hardware nonlinearities have secondary effects on pointing performance for sufficiently high bandwidth.

Marsh, E. L.↗
NASA NTRS1977

Atmospheric density models

An atmospheric model developed by Jacchia, quite accurate but requiring a large amount of computer storage and execution time, was found to be ill-suited for the space shuttle onboard program. The development of a simple atmospheric density model to simulate the Jacchia model was studied. Required characteristics including variation with solar activity, diurnal variation, variation with geomagnetic activity, semiannual variation, and variation with height were met by the new atmospheric density model.

Mueller, A. C.↗
NASA NTRS1978

A digital computer program for the dynamic interaction simulation of controls and structure (DISCOS), volume 1

A theoretical development and associated digital computer program system for the dynamic simulation and stability analysis of passive and actively controlled spacecraft are presented. The dynamic system (spacecraft) is modeled as an assembly of rigid and/or flexible bodies not necessarily in a topological tree configuration. The computer program system is used to investigate total system dynamic characteristics, including interaction effects between rigid and/or flexible bodies, control systems, and a wide range of environmental loadings. In addition, the program system is used for designing attitude control systems and for evaluating total dynamic system performance, including time domain response and frequency domain stability analyses.

Bodley, C. S.↗
NASA NTRS1980

Handling Qualities of Large Flexible Aircraft

The effects on handling qualities of elastic modes interaction with the rigid body dynamics of a large flexible aircraft are studied by a mathematical computer simulation. An analytical method to predict the pilot ratings when there is a severe modes interactions is developed. This is done by extending the optimal control model of the human pilot response to include the mode decomposition mechanism into the model. The handling qualities are determined for a longitudinal tracking task using a large flexible aircraft with parametric variations in the undamped natural frequencies of the two lowest frequency, symmetric elastic modes made to induce varying amounts of mode interaction.

Poopaka, S.↗
NASA NTRS1984

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.↗
NASA NTRS1987

Research relative to an advanced rod control system for quadrupole mass spectrometry applications

The design of a suitable amplifier output stage using available transistors and passive components is summarized. All of the analysis and calculation confirm that it is feasible to design the amplifier and quadrupole coupling circuit needed for the Advanced Rod Control System. The progress obtained so far concerning the three frequency tank circuits to be used in the oscillator for the mass spectrometer of the Cometary Rendezvous Asteroid Flyby (CRAF) project is presented. Results from this study look promising. However, it is not known what minimum impedance levels are required to make it possible for the oscillator to work properly. Therefore, it is necessary to construct a prototype circuit in the laboratory which can be measured and tested in an oscillator circuit. Continued attempts will be made to develop a useful inductor motor with better characteristics than the one being used at the moment. It is important that such a model be found if computer simulation is to reflect reality more closely.

Carignan, George R.↗
NASA NTRS1987

Spacecraft attitude determination using a second-order nonlinear filter

The stringent attitude determination accuracy and faster slew maneuver requirements demanded by present-day spacecraft control systems motivate the development of recursive nonlinear filters for attitude estimation. This paper presents the second-order filter development for the estimation of attitude quaternion using three-axis gyro and star tracker measurement data. Performance comparisons have been made by computer simulation of system models and filter mechanization. It is shown that the second-order filter consistently performs better than the extended Kalman filter when the performance index of the root sum square estimation error of the quaternion vector is compared. The second-order filter identifies the gyro drift rates faster than the extended Kalman filter. The uniqueness of this algorithm is the online generation of the time-varying process and measurement noise covariance matrices, derived as a function or the process and measurement nonlinearity, respectively.

Vathsal, S.↗
NASA NTRS1989

Why is CDMA the solution for mobile satellite communication

It is demonstrated that spread spectrum Code Division Multiple Access (CDMA) systems provide an economically superior solution to satellite mobile communications by increasing the system maximum capacity with respect to single channel per carrier Frequency Division Multiple Access (FDMA) systems. Following the comparative analysis of CDMA and FDMA systems, the design of a model that was developed to test the feasibility of the approach and the performance of a spread spectrum system in a mobile environment. Results of extensive computer simulations as well as laboratory and field tests results are presented.

Gilhousen, Klein S.↗
NASA NTRS2000

Acceleration of the Fast Solar Wind by Solitary Waves in Coronal Holes

The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear, time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that solitary-like waves are generated in coronal holes nonlinearly by torsional Alfven waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature, and magnetic field geometry) that will become available from the recently launched SOHO spacecraft.

Ofman, Leon↗
NASA NTRS2001

Acceleration of the Fast Solar Wind by Solitary Waves in Coronal Holes

The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear. time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that nonlinear waves are generated in coronal holes by torsional Alfv\'{e}n waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the, fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature. and magnetic field geometry,) that will become available from the recently launched SOHO spacecraft.

Ofman, Leon↗
NASA NTRS2003

LS-DYNA Implementation of Polymer Matrix Composite Model Under High Strain Rate Impact

A recently developed constitutive model is implemented into LS-DYNA as a user defined material model (UMAT) to characterize the nonlinear strain rate dependent behavior of polymers. By utilizing this model within a micromechanics technique based on a laminate analogy, an algorithm to analyze the strain rate dependent, nonlinear deformation of a fiber reinforced polymer matrix composite is then developed as a UMAT to simulate the response of these composites under high strain rate impact. The models are designed for shell elements in order to ensure computational efficiency. Experimental and numerical stress-strain curves are compared for two representative polymers and a representative polymer matrix composite, with the analytical model predicting the experimental response reasonably well.

Zheng, Xia-Hua↗
NASA NTRS2009

Creating a Lunar EVA Work Envelope

A work envelope has been defined for weightless Extravehicular Activity (EVA) based on the Space Shuttle Extravehicular Mobility Unit (EMU), but there is no equivalent for planetary operations. The weightless work envelope is essential for planning all EVA tasks because it determines the location of removable parts, making sure they are within reach and visibility of the suited crew member. In addition, using the envelope positions the structural hard points for foot restraints that allow placing both hands on the job and provides a load path for reacting forces. EVA operations are always constrained by time. Tasks are carefully planned to ensure the crew has enough breathing oxygen, cooling water, and battery power. Planning first involves computers using a virtual work envelope to model tasks, next suited crew members in a simulated environment refine the tasks. For weightless operations, this process is well developed, but planetary EVA is different and no work envelope has been defined. The primary difference between weightless and planetary work envelopes is gravity. It influences anthropometry, horizontal and vertical mobility, and reaction load paths and introduces effort into doing "overhead" work. Additionally, the use of spacesuits other than the EMU, and their impacts on range of motion, must be taken into account. This paper presents the analysis leading to a concept for a planetary EVA work envelope with emphasis on lunar operations. There is some urgency in creating this concept because NASA has begun building and testing development hardware for the lunar surface, including rovers, habitats and cargo off-loading equipment. Just as with microgravity operations, a lunar EVA work envelope is needed to guide designers in the formative stages of the program with the objective of avoiding difficult and costly rework.

Griffin, Brand N.↗
NASA NTRS2018

High-Fidelity Computational Methodology for Stitched Composite Aerospace Structures

Due to the high demands for energy efficient commercial transportation, the aviation industry has taken a leading role in the integration of composite structures. Among the leading concepts to develop lighter, more fuel-efficient commercial transport is the Pultruded Rod Stitched Efficient Unitized Structure (PRSEUS) concept, an enabling technology for hybrid wing bodies. Many proof-of-concept tests have been performed to demonstrate that the use of PRSEUS has improved the residual strength of damaged structures compared to conventional composite structures, but efficient computational tools must be developed before the concept can be commercially certified and implemented. In an attempt to address the need for efficient computational tools, a comprehensive modeling approach is developed and applied to investigate applications of PRSEUS at multiple scales. Therefore, a computational methodology has been progressively developed based on physically realistic concepts. The focus of the work described herein is to define the modeling characteristics required to accurately simulate the damage progression and failure of PRSEUS at the coupon scale. The work herein is focused on the development and analysis of a PRSEUS stringer, the methodology for which may be extended to other PRSEUS coupons and components.

Brandon Horton↗
NASA NTRS2023

CESO 22-3: Simulated Excavation Environment for Lunar Operations Follow-on

As more lunar In-Situ Resource Utilization (ISRU) missions are planned and developed, a need was identified for a simulation environment with accurate regolith interaction mechanics that is less computationally expensive than commercially available multiphysics tools. The purpose of the Simulated Excavation Environment for Lunar Operations (SEELO) project was to develop an interactive simulation environment that efficiently models lunar ISRU excavation operations, including soil deformation, while running in faster-than-real-time.

Michael DuPuis↗