Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Behavior modeling”

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 883 records · Page 49

Biogeochemical modeling at mass extinction boundaries

The causes of major mass extinctions is a subject of considerable interest to those concerned with the history and evolution of life on earth. The primary objectives of the proposed plan of research are: (1) to develop quantitative time-dependent biogeochemical cycle models, coupled with an ocean atmosphere in order to improve the understanding of global scale physical, chemical, and biological processes that control the distribution of elements important for life at times of mass extinctions; and (2) to develop a comprehensive data base of the best available geochemical, isotopic, and other relevant geologic data from sections across mass extinction boundaries. These data will be used to constrain and test the biogeochemical model. These modeling experiments should prove useful in: (1) determining the possible cause(s) of the environmental changes seen at bio-event boundaries; (2) identifying and quantifying little-known feedbacks among the oceans, atmosphere, and biosphere; and (3) providing additional insights into the possible responses of the earth system to perturbations of various timescales. One of the best known mass extinction events marks the Cretaceous/Tertiary (K/T) boundary (66 Myr ago). Data from the K/T boundary are used here to constrain a newly developed time-dependent biogeochemical cycle model that is designed to study transient behavior of the earth system. Model results predict significant fluctuations in ocean alkalinity, atmospheric CO2, and global temperatures caused by extinction of calcareous plankton and reduction in the sedimentation rates of pelagic carbonates and organic carbon. Oxygen-isotome and other paleoclimatic data from K/T time provide some evidence that such climatic fluctuations may have occurred, but stabilizing feedbacks may have acted to reduce the ocean alkalinity and carbon dioxide fluctuations.

Rampino, M. R.↗

Information models of software productivity - Limits on productivity growth

Research into generalized information-metric models of software process productivity establishes quantifiable behavior and theoretical bounds. The models establish a fundamental mathematical relationship between software productivity and the human capacity for information traffic, the software product yield (system size), information efficiency, and tool and process efficiencies. An upper bound is derived that quantifies average software productivity and the maximum rate at which it may grow. This bound reveals that ultimately, when tools, methodologies, and automated assistants have reached their maximum effective state, further improvement in productivity can only be achieved through increasing software reuse. The reuse advantage is shown not to increase faster than logarithmically in the number of reusable features available. The reuse bound is further shown to be somewhat dependent on the reuse policy: a general 'reuse everything' policy can lead to a somewhat slower productivity growth than a specialized reuse policy.

Tausworthe, Robert C.↗

Estimation Filter for Alignment of the Spitzer Space Telescope

A document presents a summary of an onboard estimation algorithm now being used to calibrate the alignment of the Spitzer Space Telescope (formerly known as the Space Infrared Telescope Facility). The algorithm, denoted the S2P calibration filter, recursively generates estimates of the alignment angles between a telescope reference frame and a star-tracker reference frame. At several discrete times during the day, the filter accepts, as input, attitude estimates from the star tracker and observations taken by the Pointing Control Reference Sensor (a sensor in the field of view of the telescope). The output of the filter is a calibrated quaternion that represents the best current mean-square estimate of the alignment angles between the telescope and the star tracker. The S2P calibration filter incorporates a Kalman filter that tracks six states - two for each of three orthogonal coordinate axes. Although, in principle, one state per axis is sufficient, the use of two states per axis makes it possible to model both short- and long-term behaviors. Specifically, the filter properly models transient learning, characteristic times and bounds of thermomechanical drift, and long-term steady-state statistics, whether calibration measurements are taken frequently or infrequently. These properties ensure that the S2P filter performance is optimal over a broad range of flight conditions, and can be confidently run autonomously over several years of in-flight operation without human intervention.

Bayard, David↗

Thermal crosstalk measurements and simulations for an X-ray microcalorimeter array

Arrays of high-density microcalorimeters require careful heat sinking in order to minimize the thermal crosstalk between nearby pixels. For the array of microcalorimeters developed for the Athena X-ray Integral Field Unit instrument, which has more than 3000 pixels on a 275 µm pitch, it is essential to address this problem in order to meet the energy-resolution requirements. The instrument’s energy-resolution budget requires that the impact of the thermal crosstalk on the energy resolution be a contribution that, added in quadrature to other energy-resolution contributions, is less than 0.2 eV. This value results in a derived requirement that the ratio between the amplitude of the crosstalk signal to an X-ray pulse (for example at 6 keV) is less than 1 × 10(exp −3) (for the first neighbor), less than 4 × 10(exp −4) (for the diagonal neighbor) and less than 8 × 10(exp −5) (for the second nearest neighbor). We have measured the thermal crosstalk levels between pixels in various geometries and configurations. The results show a crosstalk ratio which is at least a factor of 4 lower than the derived requirement. We also developed a finite element (FEM) 2D thermal model to predict the thermal behavior of large-scale arrays. This model successfully simulates the measured data in terms of pulse amplitude and time constants.

Antoine R. Miniussi↗

A potential-energy scaling model to simulate the initial stages of thin-film growth

A solid on solid (SOS) Monte Carlo computer simulation employing a potential energy scaling technique was used to model the initial stages of thin film growth. The model monitors variations in the vertical interaction potential that occur due to the arrival or departure of selected adatoms or impurities at all sites in the 400 sq. ft. array. Boltzmann ordered statistics are used to simulate fluctuations in vibrational energy at each site in the array, and the resulting site energy is compared with threshold levels of possible atomic events. In addition to adsorption, desorption, and surface migration, adatom incorporation and diffusion of a substrate atom to the surface are also included. The lateral interaction of nearest, second nearest, and third nearest neighbors is also considered. A series of computer experiments are conducted to illustrate the behavior of the model.

Heinbockel, J. H.↗

Faulty behavior of asynchronous storage elements

It is often assumed that the faults in storage elements (SE's) can be modeled as output/input stuck-at-faults of the element. They are implicitly considered equivalent to the stuck-at faults in the combinational logic surrounding the SE cells. A more accurate higher level fault model for elementary SE's used in asynchronous circuits is presented. This model offers better representation of the physical failures. It is shown that the stuck-at model may be adequate if only modest fault coverage is desired. The enhanced model includes some common fault behaviors of SE's that are not covered by the stuck-at model. These include data-feed-through behaviors that cause the SE to be combinational. Fault models for complex SE cells can be obtained without a significant loss of information about the structure of the circuit.

Al-Assadi, Waleed K.↗

Neural Network Models of Simple Mechanical Systems Illustrating the Feasibility of Accelerated Life Testing

A complete evaluation of the tribological characteristics of a given material/mechanical system is a time-consuming operation since the friction and wear process is extremely systems sensitive. As a result, experimental designs (i.e., Latin Square, Taguchi) have been implemented in an attempt to not only reduce the total number of experimental combinations needed to fully characterize a material/mechanical system, but also to acquire life data for a system without having to perform an actual life test. Unfortunately, these experimental designs still require a great deal of experimental testing and the output does not always produce meaningful information. In order to further reduce the amount of experimental testing required, this study employs a computer neural network model to investigate different material/mechanical systems. The work focuses on the modeling of the wear behavior, while showing the feasibility of using neural networks to predict life data. The model is capable of defining which input variables will influence the tribological behavior of the particular material/mechanical system being studied based on the specifications of the overall system.

Fusaro, Robert L.↗

New Developments in the CompDam Fiber Kinking Model for the Interaction of Kinking and Splitting Cracks with Application to Open Hole Compression Specimens

This paper describes two new additions to the NASA fiber kinking model: 1) A failure index for fiber kinking is developed to facilitate interaction between kinking and matrix cracking and 2) A criterion for fiber rupture is introduced to account for fiber breakage before or after fiber kinking. The model formulation is described with the aid of graphical Considered constructions to elucidate the behavior. Verification studies are conducted with single element models to show the behavior of the model in light of the new developments and to evaluate the effect of numerical parameters of the model. The model is used to analyze open hole compression specimens with soft, quasi-isotropic, and hard layups. Analysis results are shown for a conventional continuum damage mechanics model as well as the proposed model to highlight the differences between the two. Comparison with test data highlights the capabilities and remaining gaps in the quantitative prediction of damage evolution in open hole compression. Overall, the results show that the proposed model is capable of representing the kinematics of fiber kinking, and specifically the role of fiber kinking in open hole compression.

Andrew Cole Bergan↗

Rayleigh Light Scattering for Concentration Measurements in Turbulent Flows

Despite intensive research over a number of years, an understanding of scalar mixing in turbulent flows remains elusive. An understanding is required because turbulent mixing has a pivotal role in a wide variety of natural and technologically important processes. As an example, the mixing and transport of pollutants in the atmosphere and in bodies of water are often dependent on turbulent mixing processes. Turbulent mixing is also central to turbulent combustion which underlies most hydrocarbon energy use in modern societies as well as in unwanted fire behavior. Development of models for combusting flows is therefore crucial, however, an understanding of scalar mixing is required before useful models of turbulent mixing and, ultimately, turbulent combustion can be developed. An important subset of turbulent flows is axisymmetric turbulent jets and plumes because they are relatively simple to generate, and because the provide an appropriate test bed for the development of general theories of turbulent mixing which can be applied to more complex geometries and flows. This paper focuses on a number of experimental techniques which have been developed at the National Institute of Standards and Development for measuring concentration in binary axisymmetric turbulent jets. In order to demonstrate the value of these diagnostics, some of the more important results from earlier and on-going investigations are summarized. Topics addressed include the similarity behavior of variable density axisymmetric jets, the behavior of absolutely unstable axisymmetric helium jets, and the role of large scale structures and scalar dissipation in these flows.

Pitts, William M.↗

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗

Sea surface topography fields of the tropical Pacific from data assimilation

Time series of maps of monthly tropical Pacific dynamic topography anomalies from 1979 through 1985 were constructed by means of assimilation of the tide gauge and expendable bathythermograph (XBT) data into a linear model driven by observed winds. Estimates of error statistics were calculated and compared to actual differences between hindcasts and observations. Four experiements were performed as follows: one with no assimilation, one with assimiation of sea level anomaly data from eight selected island tide gauge stations, one with assimilation of dynamic height anomalies derived from XBT data, and one with both XBT and tide gauge data assimilated. Data from seven additional tide gauge stations were withheld from the assimilation process and used for verification in all four experiements. Statistical objective maps based on data alone were also constructed for comparison purposes. The dynamic response of the model without assimilation was, in general, weaker than the observed response. Assimilation resulted in enhanced signal amplitude in all three assimilation experiments. With few exceptions the error estimated generated by the Kalman filter appeared quite reasonable. Since the error processes cannot be assumed to be white or stationary, we could find no straightforward way to test the formal statistical hypothesis that the time series of differences between the filter ouput and the actual observations were drawn from a population with statistics given by the Kalman filter estimaes. We found that the best filter we could devise was still short of the goal of producing a white innovation sequence. Extensive changes in the assumed error statistics make only marginal differences. The same is true for long time and space scale behavior of different models with richer physics and finer resolution. Better data assimilation results will probably require relaxation of the assumptions of stationarity and serial independence of the errors. Formation of such detailed noise models will require longer time series, with the attendant problems of matching very different data sets.

Miller, Robert N.↗

Stochastic differential equations and turbulent dispersion

Aspects of the theory of continuous stochastic processes that seem to contribute to an understanding of turbulent dispersion are introduced and the theory and philosophy of modelling turbulent transport is emphasized. Examples of eddy diffusion examined include shear dispersion, the surface layer, and channel flow. Modeling dispersion with finite-time scale is considered including the Langevin model for homogeneous turbulence, dispersion in nonhomogeneous turbulence, and the asymptotic behavior of the Langevin model for nonhomogeneous turbulence.

Durbin, P. A.↗

Modeling the Benchmark Active Control Technology Wind-Tunnel Model for Application to Flutter Suppression

This paper describes the formulation of a model of the dynamic behavior of the Benchmark Active Controls Technology (BACT) wind-tunnel model for application to design and analysis of flutter suppression controllers. The model is formed by combining the equations of motion for the BACT wind-tunnel model with actuator models and a model of wind-tunnel turbulence. The primary focus of this paper is the development of the equations of motion from first principles using Lagrange's equations and the principle of virtual work. A numerical form of the model is generated using values for parameters obtained from both experiment and analysis. A unique aspect of the BACT wind-tunnel model is that it has upper- and lower-surface spoilers for active control. Comparisons with experimental frequency responses and other data show excellent agreement and suggest that simple coefficient-based aerodynamics are sufficient to accurately characterize the aeroelastic response of the BACT wind-tunnel model. The equations of motion developed herein have been used to assist the design and analysis of a number of flutter suppression controllers that have been successfully implemented.

Waszak, Martin R.↗

Design for Validation

The use of computer hardware and software in life-critical applications, such as for civil air transports, demands the use of rigorous formal mathematical validation procedures. However, formal specification and verification will only be tractable if the system is designed in a manner that lends itself to formal methods. Likewise, accurate reliability analysis will only be tractable if the number of interacting components that must be individually included in a single reliability model is kept to a low number and if their failure behavior interactions can be modeled simply. Also, the system must be designed such that the system reliability does not directly depend on system parameters that cannot be accurately determined. This paper presents a design methodology based on the concept of designing a system in such a manner that it can be rigorously validated, or "design for validation."

Johnson, Sally C.↗

MISTE flight experiment status

The MISTE flight experiment has made significant scientific and technical progress in preparing for a future microgravity flight. We are collaborating with several theoretical modeling groups that have developed crossover (equation-of-state) models for predicting thermodynamic behavior near the liquid-gas critical point. Several of these models have already been used to test experimental measurements of the heat capacity at constant volume, isothermal susceptibility, and coexistence curve in the crossover region near the 3He liquid-gas critical point. A brief description of these models and a representative fit to experimental data will be presented. In collaboration with Mission Research Corporation, MISTE has been testing a new small pneumatic valve for use at low temperatures. The results of recent successful low temperature actuation tests will also be discussed.

critical phenomena liquid-gas critical point↗

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗

Stardust-NExT, Deep Impact, and the Accelerating Spin of 9P/Tempel One

The evolution of the spin rate of comet 9P/Tempel 1 through two perihelion passages (CYs 2000 and 2005) is determined from 1922 Earth-based observations taken over a period of 13y as part of a World-Wide observing campaign and 2888 observations taken over a period of 50d from the Deep Impact spacecraft. We determine the following sidereal spin rates (periods): 209.023 +/- 0.025 degrees /day (41.335 +/- 0.005 h) prior to the 2000 perihelion passage, 210.448 +/- 0.016 degrees/day (41.055 +/- 0.003 h) for the interval between the 2000 and 2005 perihelion passages, 211.856 +/- 0.030 degrees/day (40.783 +/- 0.006 h) from Deep Impact photometry just prior to the 2005 perihelion passage, and 211.625 +/- 0.012 degrees /day (40.827 +/- 0.002 h) in the interval 2006-2010 following the 2005 perihelion passage. The period decreased by 16.8 +/- 0.3 min during the 2000 passage and by 13.7 +/- 0.2 min during the 2005 passage suggesting a secular decrease in the net torque. The change in spin rate is asymmetric with respect to perihelion with the maximum net torque being applied on approach to perihelion. The Deep Impact data alone show that the spin rate was increasing at a rate of 0.024 +/- 0.003 degree/d/d at JD2453530.60510 (i.e., 25.134 d before impact) and provides independent confirmation of the change seen in the Earth-based observations. The rotational phase of the nucleus at times before and after each perihelion and at the Deep Impact encounter is estimated based on the Thomas et al. pole and longitude system. The possibility of a 180 degree error in the rotational phase is assessed and found to be significant. Analytical and physical modeling of the behavior of the spin rate through of each perihelion is presented and used as a basis to predict the rotational state of the nucleus at the time of the nominal (i.e., prior to February 2010) Stardust-NExT encounter on 2011 February 14 20:42. We find that a net torque in the range of 0.3 - 2.5 x 10(exp 7) kg.square m2/square s acts on the nucleus during perihelion passage. The spin rate initially slows down on approach to perihelion and then passes through a minimum. It then accelerates rapidly as it passes through perihelion eventually reaching a maximum post-perihelion. It then decreases to a stable value as the nucleus moves away from the sun. We find that the pole direction is unlikely to precess by more than approximately 1 degree/perihelion passage. The trend of the period with time and the fact that the modeled peak torque that occurs before perihelion is in agreement with published accounts of trends in water production rate and suggests that widespread H2O out-gassing from the surface is largely responsible for the observed spin-up.

Belton, Michael J. S.↗

Stardust-NExT, Deep Impact, and the Accelerating Spin of 9P/Tempel 1

The evolution of the spin rate of Comet 9P/Tempel 1 through two perihelion passages (in 2000 and 2005) is determined from 1922 Earth-based observations taken over a period of 13 year as part of a World-Wide observing campaign and from 2888 observations taken over a period of 50 days from the Deep Impact spacecraft. We determine the following sidereal spin rates (periods): 209.023 +/- 0.025deg/dy (41.335 0.005 h) prior to the 2000 perihelion passage, 210.448 +/- 0.016deg/dy (41.055 +/- 0.003 h) for the interval between the 2000 and 2005 perihelion passages, 211.856 +/- 0.030deg/dy (40.783 +/- 0.006 h) from Deep Impact photometry just prior to the 2005 perihelion passage, and 211.625 +/- 0.012deg/dy (40.827 +/- 0.002 h) in the interval 2006-2010 following the 2005 perihelion passage. The period decreased by 16.8 +/- 0.3 min during the 2000 passage and by 13.7 +/- 0.2 min during the 2005 passage suggesting a secular decrease in the net torque. The change in spin rate is asymmetric with respect to perihelion with the maximum net torque being applied on approach to perihelion. The Deep Impact data alone show that the spin rate was increasing at a rate of 0.024 +/- 0.003deg/dy/dy at JD2453530.60510 (i.e., 25.134 dy before impact), which provides independent confirmation of the change seen in the Earth-based observations. The rotational phase of the nucleus at times before and after each perihelion and at the Deep Impact encounter is estimated based on the Thomas et al. (Thomas et al. [2007]. Icarus 187, 4-15) pole and longitude system. The possibility of a 180deg error in the rotational phase is assessed and found to be significant. Analytical and physical modeling of the behavior of the spin rate through of each perihelion is presented and used as a basis to predict the rotational state of the nucleus at the time of the nominal (i.e., prior to February 2010) Stardust-NExT encounter on 2011 February 14 at 20:42. We find that a net torque in the range of 0.3-2.5 x 10(exp 7) kg/sq m/sq s acts on the nucleus during perihelion passage. The spin rate initially slows down on approach to perihelion and then passes through a minimum. It then accelerates rapidly as it passes through perihelion eventually reaching a maximum post-perihelion. It then decreases to a stable value as the nucleus moves away from the Sun. We find that the pole direction is unlikely to precess by more than approx. 1deg per perihelion passage. The trend of the period with time and the fact that the modeled peak torque occurs before perihelion are in agreement with published accounts of trends in water production rate and suggests that widespread H2O out-gassing from the surface is largely responsible for the observed spin-up.

Belton, Michael J. S.↗