Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time discretization”

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 595 records · Page 33

Phase transformation kinetics in finite inhomogeneously nucleated systems

Phase transformation kinetics that occur by a nucleation and growth process are investigated. A simple discrete space and time model is used for the dynamics and analytical results are obtained for the volume fraction of the material transformed for both finite systems and a special example of an inhomogeneously nucleated system. The theory is developed for two cases, initial nucleation, and continuous nucleation. The results are compared with simulations of the model.

Weinberg, Michael↗

CFD validation for aerodynamic flows - Challenge for the '90's

The process of computational fluid dynamics (CFD) validation is described from two perspectives, numerical error validation and physical model validation. Errors associated with time and space discretization, grid refinement, numerical dissipation, and level of solution convergence, are all possible sources of numerical error. The second aspect of CFD code validation is associated with the governing equations and the physical models, e.g., chemistry and turbulence models. This type of validation requires comparisons with carefully conducted experiments and is the primary aspect of the validation process discussed in this paper. Examples of validation are shown for selected Reynolds-averaged Navier-Stokes codes for conditions from transonic through hypersonic. Since code applications are becoming more complex, it no longer suffices to use data from surface or integral quantities alone to provide the required validation. Flow field surveys and experimental boundary condition measurements are emerging as critical data that must be obtained for CFD code validation. Progress in the overall level of validation is shown to be improving. However, future challenges remain, some of which are enumerated.

Marvin, Joseph G.↗

A method for solution of the Euler-Bernoulli beam equation in flexible-link robotic systems

An efficient numerical method for solving the partial differential equation (PDE) governing the flexible manipulator control dynamics is presented. A finite-dimensional model of the equation is obtained through discretization in both time and space coordinates by using finite-difference approximations to the PDE. An expert program written in the Macsyma symbolic language is utilized in order to embed the boundary conditions into the program, accounting for a mass carried at the tip of the manipulator. The advantages of the proposed algorithm are many, including the ability to (1) include any distributed actuation term in the partial differential equation, (2) provide distributed sensing of the beam displacement, (3) easily modify the boundary conditions through an expert program, and (4) modify the structure for running under a multiprocessor environment.

Tzes, Anthony P.↗

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗

Dynamics of spacecraft with deploying flexible appendages

A computational formulation for the dynamic analysis of spacecraft with deploying appendages is presented. The appendage model is based on a geometrically nonlinear beam formulation which accurately accounts for large rotational and large deformation motions. A moving finite element reference grid is incorporated within the nonlinear beam formulation to model the deployment motion. Hamilton's Law is used to formulate the general equations of motion, and a transient integration solution procedure is derived from a space-time finite element discretization of the Hamiltonian variational statement. Computational results of the methodology are presented for a classical gravity gradient stabilized satellite configuration.

Downer, Janice D.↗

Hidden Markov models for fault detection in dynamic systems

The invention is a system failure monitoring method and apparatus which learns the symptom-fault mapping directly from training data. The invention first estimates the state of the system at discrete intervals in time. A feature vector x of dimension k is estimated from sets of successive windows of sensor data. A pattern recognition component then models the instantaneous estimate of the posterior class probability given the features, p(w(sub i) perpendicular to x), 1 less than or equal to i is less than or equal to m. Finally, a hidden Markov model is used to take advantage of temporal context and estimate class probabilities conditioned on recent past history. In this hierarchical pattern of information flow, the time series data is transformed and mapped into a categorical representation (the fault classes) and integrated over time to enable robust decision-making.

Smyth, Padhraic J.↗

A finite element study of incompressible flows past oscillating cylinders and aerofoils

We present our numerical results for certain unsteady flows past oscillating cylinders and aerofoils. The computations are based on the stabilized space-time finite element formulation. The implicit equation systems resulting from the space-time finite element discretizations are solved using iterative solution techniques. One of the problems studied is flow past a cylinder which is forced to oscillate in the horizontal direction. In this case we observe a change from an unsymmetric mode of vortex shedding to a symmetric one. An extensive study was carried out for the case in which a cylinder is mounted on lightly damped springs and allowed to oscillate in the vertical direction. In this case the motion of the cylinder needs to be determined as part of the solution, and under certain conditions this motion changes the vortex-shedding pattern of the flow field significantly. This nonlinear fluid-structure interaction exhibits certain interesting behavior such as 'lock-in' and 'hysteresis', which are in good agreement with the laboratory experiments carried out by other researchers in the past. Preliminary results for flow past a pitching aerofoil are also presented.

Mittal, S.↗

Experiment and simulation for CSI: What are the missing links?

Viewgraphs on experiment and simulation for control structure interaction (CSI) are presented. Topics covered include: control structure interaction; typical control/structure interaction system; CSI problem classification; actuator/sensor models; modeling uncertainty; noise models; real-time computations; and discrete versus continuous.

Belvin, W. Keith↗

Reliability of Autonomic Responses and Malaise Across Multiple Motion Sickness Stimulation Tests

There is general agreement that a high degree of variability exists between subjects in their autonomic nervous system responses to motion sickness stimulation. Additionally, a paucity of data exists that examines the variability within an individual across repeated motion sickness tests. Investigators have also examined the relationship of autonomic responses to motion sickness development. These investigations have used analyses at discrete points in time to describe this relationship. This approach fails to address the time course of autonomic responses and malaise development throughout the motion sickness test. Our objectives were to examine the reliability of autonomic responses and malaise using the final minute of the motion sickness test across five testing occasions, to examine the reliability of the change in autonomic responses and the change in malaise across five testing occasions, and to examine the relationship between changes in autonomic responses and changes in malaise level across the entire motion sickness test. Our results indicate that, based on the final minute of testing, the autonomic responses of heart rate, blood volume pulse, and respiration rate are moderately stable across multiple tests. Changes in heart rate, blood volume pulse, respiration rate, and malaise throughout the test duration were less stable across the tests. We attribute this instability to variations in individual susceptibility and the error associated with estimating a measure of autonomic gain.

Stout, Cynthia S.↗

Global changes in biogeochemical cycles in response to human activities

The main objective of our research was to characterize biogeochemical cycles at continental and global scales in both terrestrial and aquatic ecosystems. This characterization applied to both natural ecosystems and those disturbed by human activity. The primary elements of interest were carbon and nitrogen and the analysis sought to quantify standing stocks and dynamic cycling processes. The translocation of major nutrients from the terrestrial landscape to the atmosphere (via trace gases) and to fluvial systems (via leaching, erosional losses, and point source pollution) were of particular importance to this study. Our aim was to develop the first generation of Earth System Models. Our research was organized around the construction and testing of component biogeochemical models which treated terrestrial ecosystem processes, aquatic nutrient transport through drainage basins, and trace gas exchanges at the continental and global scale. A suite of three complementary models were defined within this construct. The models were organized to operate at a 1/2 degree latitude by longitude level of spatial resolution and to execute at a monthly time step. This discretization afforded us the opportunity to understand the dynamics of the biosphere down to subregional scales, while simultaneously placing these dynamics into a global context.

Moore, Berrien, III↗

An efficient solution procedure for the thermoelastic analysis of truss space structures

A solution procedure is proposed for the thermal and thermoelastic analysis of truss space structures in periodic motion. In this method, the spatial domain is first descretized using a consistent finite element formulation. Then the resulting semi-discrete equations in time are solved analytically by using Fourier decomposition. Geometrical symmetry is taken advantage of completely. An algorithm is presented for the calculation of heat flux distribution. The method is demonstrated via a numerical example of a cylindrically shaped space structure.

Givoli, D.↗

Hidden Markov models for fault detection in dynamic systems

The invention is a system failure monitoring method and apparatus which learns the symptom-fault mapping directly from training data. The invention first estimates the state of the system at discrete intervals in time. A feature vector x of dimension k is estimated from sets of successive windows of sensor data. A pattern recognition component then models the instantaneous estimate of the posterior class probability given the features, p(w(sub i) (vertical bar)/x), 1 less than or equal to i isless than or equal to m. Finally, a hidden Markov model is used to take advantage of temporal context and estimate class probabilities conditioned on recent past history. In this hierarchical pattern of information flow, the time series data is transformed and mapped into a categorical representation (the fault classes) and integrated over time to enable robust decision-making.

Smyth, Padhraic J.↗

Optimal Approximation of Quadratic Interval Functions

Measurements are never absolutely accurate, as a result, after each measurement, we do not get the exact value of the measured quantity; at best, we get an interval of its possible values, For dynamically changing quantities x, the additional problem is that we cannot measure them continuously; we can only measure them at certain discrete moments of time t(sub 1), t(sub 2), ... If we know that the value x(t(sub j)) at a moment t(sub j) of the last measurement was in the interval [x-(t(sub j)), x + (t(sub j))], and if we know the upper bound D on the rate with which x changes, then, for any given moment of time t, we can conclude that x(t) belongs to the interval [x-(t(sub j)) - D (t - t(sub j)), x + (t(sub j)) + D (t - t(sub j))]. This interval changes linearly with time, an is, therefore, called a linear interval function. When we process these intervals, we get an expression that is quadratic and higher order w.r.t. time t, Such "quadratic" intervals are difficult to process and therefore, it is necessary to approximate them by linear ones. In this paper, we describe an algorithm that gives the optimal approximation of quadratic interval functions by linear ones.

Koshelev, Misha↗

A Discontinuous Galerkin Method for Parabolic Problems with Modified hp-Finite Element Approximation Technique

A recent paper is generalized to a case where the spatial region is taken in R(sup 3). The region is assumed to be a thin body, such as a panel on the wing or fuselage of an aerospace vehicle. The traditional h- as well as hp-finite element methods are applied to the surface defined in the x - y variables, while, through the thickness, the technique of the p-element is employed. Time and spatial discretization scheme based upon an assumption of certain weak singularity of double vertical line u(sub t) double vertical line 2, is used to derive an optimal a priori error estimate for the current method.

Kaneko, Hideaki↗

Pressure Oscillations and Structural Vibrations in Space Shuttle RSRM and ETM-3 Motors

The complex interactions between internal motor pressure oscillations resulting from vortex shedding, the motor's internal acoustic modes, and the motor's structural vibration modes were assessed for the Space Shuttle four-segment booster Reusable Solid Rocket Motor and for the five-segment engineering test motor ETM-3. Two approaches were applied 1) a predictive procedure based on numerically solving modal representations of a solid rocket motor s acoustic equations of motion and 2) a computational fluid dynamics two-dimensional axi-symmetric large eddy simulation at discrete motor burn times.

Mason, D. R.↗

Continuum Damage Mechanics Used to Predict the Creep Life of Monolithic Ceramics

Significant improvements in propulsion and power generation for the next century will require revolutionary advances in high-temperature materials and structural design. Advanced ceramics are candidate materials for these elevated temperature applications. High-temperature and long-duration applications of monolithic ceramics can place their failure mode in the creep rupture regime. An analytical methodology in the form of the integrated design program-Ceramics Analysis and Reliability Evaluation of Structures/Creep (CARES/Creep) has been developed by the NASA Lewis Research Center to predict the life of ceramic structural components subjected to creep rupture conditions. This program utilizes commercially available finite element packages and takes into account the transient state of stress and creep strain distributions (stress relaxation as well as the asymmetric response to tension and compression). The creep life of a component is discretized into short time steps, during which the stress distribution is assumed constant. Then, the damage is calculated for each time step on the basis of a modified Monkman-Grant (MMG) creep rupture criterion. The cumulative damage is subsequently calculated as time elapses in a manner similar to Miner's rule for cyclic fatigue loading. Failure is assumed to occur when the normalized cumulative damage at any point in the component reaches unity. The corresponding time is the creep rupture life for that component.

Powers, Lynn M.↗

Statistical and Probabilistic Extensions to Ground Operations' Discrete Event Simulation Modeling

NASA's human exploration initiatives will invest in technologies, public/private partnerships, and infrastructure, paving the way for the expansion of human civilization into the solar system and beyond. As it is has been for the past half century, the Kennedy Space Center will be the embarkation point for humankind's journey into the cosmos. Functioning as a next generation space launch complex, Kennedy's launch pads, integration facilities, processing areas, launch and recovery ranges will bustle with the activities of the world's space transportation providers. In developing this complex, KSC teams work through the potential operational scenarios: conducting trade studies, planning and budgeting for expensive and limited resources, and simulating alternative operational schemes. Numerous tools, among them discrete event simulation (DES), were matured during the Constellation Program to conduct such analyses with the purpose of optimizing the launch complex for maximum efficiency, safety, and flexibility while minimizing life cycle costs. Discrete event simulation is a computer-based modeling technique for complex and dynamic systems where the state of the system changes at discrete points in time and whose inputs may include random variables. DES is used to assess timelines and throughput, and to support operability studies and contingency analyses. It is applicable to any space launch campaign and informs decision-makers of the effects of varying numbers of expensive resources and the impact of off nominal scenarios on measures of performance. In order to develop representative DES models, methods were adopted, exploited, or created to extend traditional uses of DES. The Delphi method was adopted and utilized for task duration estimation. DES software was exploited for probabilistic event variation. A roll-up process was used, which was developed to reuse models and model elements in other less - detailed models. The DES team continues to innovate and expand DES capabilities to address KSC's planning needs.

Trocine, Linda↗

The Evolving Polarized Jet of Black Hole Candidate Swift J1745-26

Swift J1745−26 is an X-ray binary towards the Galactic Centre that was detected when it went into outburst in September 2012. This source is thought to be one of a growing number of sources that display "failed outbursts", in which the self-absorbed radio jets of the transient source are never fully quenched and the thermal emission from the geometrically-thin inner accretion disk never fully dominates the X-ray flux. We present multifrequency data from the Very Large Array, Australia Telescope Compact Array and Karoo Array Telescope (KAT- 7) radio arrays, spanning the entire period of the outburst. Our rich data set exposes radio emission that displays a high level of large scale variability compared to the X-ray emission and deviations from the standard radio-X-ray correlation that are indicative of an unstable jet and confirm the outburst's transition from the canonical hard state to an intermediate state. We also observe steepening of the spectral index and an increase of the linear polarization to a large fraction (is approx. equal to 50%) of the total flux, as well as a rotation of the electric vector position angle. These are consistent with a transformation from a self-absorbed compact jet to optically-thin ejecta - the first time such a discrete ejection has been observed in a failed outburst - and may imply a complex magnetic field geometry.

black↗