Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A COMPACT model reduction methodology for articulated, multi-flexible bodies structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one need appropriate reduced-order models for all the flexible components involved. The component modes projection and assembly model reduction methodology is one way to construct these reduced-order component models. However, the generated reduced-order system models typically contain high frequency extraneous modes. The presence of high frequency extraneous modes in a reduced-order system model forces us to use a small stepsize in the numerical integration of that model, resulting in an increase in the elapsed time of the integration process. To eliminate these high-frequency modes, we first identify and discard the component modes that 'contributed' to the generations of these modes. A sensitivity analysis is then used to determine how the remaining component modes should be perturbed in order to still closely capture the selected system modes. The effectiveness of the proposed methodology has been successfully verified using a finite-element model of the Galileo spacecraft.

Lee, Allan Y.↗
NASA NTRS1987

An evaluation of higher-order model methods for calculating transient structural response

A higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method is evaluated. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force-derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗
NASA NTRS1987

An evaluation of higher-order modal methods for calculating transient structural response

The present study evaluates a higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.

Camarda, Charles J.↗
NASA NTRS1987

Decoupled dynamic analysis of combined systems by iterative determination of interface accelerations

A dynamic analysis technique is presented that can be used to determine the response of a discrete model of a large linear structural system composed of multiple substructures. The technique circumvents the costly computation of the modal characteristics of the combined system. This is accomplished by relying on a predictor-corrector scheme to converge iteratively to the interface accelerations of the combined system, while the equations of motions of the individual structures are integrated separately. In this regard, the temporal slopes of the interface accelerations (jerks) are computed at each time point of integration to predict the interface accelerations at the next time point. The proposed technique is exemplified by conducting a Space Shuttle landing loads analysis; the obtained numerical data demonstrate its reliability and efficiency.

Spanos, P. D.↗
NASA NTRS1991

Real-time in situ sensors and control integration for life support systems

The limitations of the state-of-the-art for in situ sensors are discussed and a program of adaptation and enhancement of off-the-shelf sensor technologies and of innovation and research to develop more appropriate sensor technologies for life support systems is offered. By critically assessing the state-of-the-art in multifunctional sensors and smart sensors, research and development requirements for life support systems can be defined. Consideration is given to the desirable characteristics of smart sensors for life support applications, and some preliminary concepts for hierarchical integration of in situ sensors and control elements are presented.

Voecks, G. E.↗
NASA NTRS1995

Building a better leapfrog

In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

Hut, Piet↗
NASA NTRS2010

High-Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

It has been known for some time that Taylor series (TS) integration is among the most efficient and accurate numerical methods in solving differential equations. However, the full benefit of the method has yet to be realized in calculating spacecraft trajectories, for two main reasons. First, most applications of Taylor series to trajectory propagation have focused on relatively simple problems of orbital motion or on specific problems and have not provided general applicability. Second, applications that have been more general have required use of a preprocessor, which inevitably imposes constraints on computational efficiency. The latter approach includes the work of Berryman et al., who solved the planetary n-body problem with relativistic effects. Their work specifically noted the computational inefficiencies arising from use of a preprocessor and pointed out the potential benefit of manually coding derivative routines. In this Engineering Note, we report on a systematic effort to directly implement Taylor series integration in an operational trajectory propagation code: the Spacecraft N-Body Analysis Program (SNAP). The present Taylor series implementation is unique in that it applies to spacecraft virtually anywhere in the solar system and can be used interchangeably with another integration method. SNAP is a high-fidelity trajectory propagator that includes force models for central body gravitation with N X N harmonics, other body gravitation with N X N harmonics, solar radiation pressure, atmospheric drag (for Earth orbits), and spacecraft thrusting (including shadowing). The governing equations are solved using an eighth-order Runge-Kutta Fehlberg (RKF) single-step method with variable step size control. In the present effort, TS is implemented by way of highly integrated subroutines that can be used interchangeably with RKF. This makes it possible to turn TS on or off during various phases of a mission. Current TS force models include central body gravitation with the J2 spherical harmonic, other body gravitation, thrust, constant atmospheric drag from Earth's atmosphere, and solar radiation pressure for a sphere under constant illumination. The purpose of this Engineering Note is to demonstrate the performance of TS integration in an operational trajectory analysis code and to compare it with a standard method, eighth-order RKF. Results show that TS is 16.6 times faster on average and is more accurate in 87.5% of the cases presented.

Scott, James R.↗
NASA NTRSDate not supplied

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗
NASA NTRS2017

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗
NASA NTRS2007

Analytic Formulation and Numerical Implementation of an Acoustic Pressure Gradient Prediction

The scattering of rotor noise is an area that has received little attention over the years, yet the limited work that has been done has shown that both the directivity and intensity of the acoustic field may be significantly modified by the presence of scattering bodies. One of the inputs needed to compute the scattered acoustic field is the acoustic pressure gradient on a scattering surface. Two new analytical formulations of the acoustic pressure gradient have been developed and implemented in the PSU-WOPWOP rotor noise prediction code. These formulations are presented in this paper. The first formulation is derived by taking the gradient of Farassat's retarded-time Formulation 1A. Although this formulation is relatively simple, it requires numerical time differentiation of the acoustic integrals. In the second formulation, the time differentiation is taken inside the integrals analytically. The acoustic pressure gradient predicted by these new formulations is validated through comparison with the acoustic pressure gradient determined by a purely numerical approach for two model rotors. The agreement between analytic formulations and numerical method is excellent for both stationary and moving observers case.

Lee, Seongkyu↗
NASA NTRS2005

Development of a Real-Time GPS/Seismic Displacement Meter: Applications to Civilian Infrastructure in Orange and Western Riverside Counties, California

We propose a three-year applications project that will develop an Integrated Real-Time GPS/Seismic System and deploy it in Orange and Western Riverside Counties, spanning three major strike-slip faults in southern California (San Andreas, San Jacinto, and Elsinore) and significant populations and civilian infrastructure. The system relying on existing GPS and seismic networks will collect and analyze GPS and seismic data for the purpose of estimating and disseminating real-time positions and total ground displacements (dynamic, as well as static) during all phases of the seismic cycle, from fractions of seconds to years. Besides its intrinsic scientific use as a real-time displacement meter (transducer), the GPS/Seismic System will be a powerful tool for local and state decision makers for risk mitigation, disaster management, and structural monitoring (dams, bridges, and buildings). Furthermore, the GPS/Seismic System will become an integral part of California's spatial referencing and positioning infrastructure, which is complicated by tectonic motion, seismic displacements, and land subsidence. Finally, the GPS/Seismic system will also be applicable to navigation in any environment (land, sea, or air) by combining precise real-time instantaneous GPS positioning with inertial navigation systems. This development will take place under the umbrella of the California Spatial Reference Center, in partnership with local (Counties, Riverside County Flood and Water Conservation District, Metropolitan Water District), state (Caltrans), and Federal agencies (NGS, NASA, USGS), the geophysics community (SCIGN/SCEC2), and the private sector (RBF Consulting). The project will leverage considerable funding, resources, and R&D from SCIGN, CSRC and two NSF-funded IT projects at UCSD and SDSU: RoadNet (Real-Time Observatories, Applications and Data Management Network) and the High Performance Wireless Research and Education Network (HPWREN). These two projects are funded to develop both the wireless networks and the integrated, seamless, and transparent information management system that will deliver seismic, geodetic, oceanographic, hydrological, ecological, and physical data to a variety of end users in real-time in the San Diego region. CSRC is interested in providing users access to real-time, accurate GPS data for a wide variety of applications including RTK surveying/GIS and positioning of moving platforms such as aircraft and emergency vehicles. SCIGN is interested in upgrading sites to high-frequency real-time operations for rapid earthquake response and GPS seismology. The successful outcome of the project will allow the implementation of similar systems elsewhere, particularly in plate boundary zones with significant populations and civilian infrastructure. CSRC would like to deploy the GPS/Seismic System in other parts of California, in particular San Diego, Los Angeles County and the San Francisco Bay Area.

Bock, Yehuda↗
NASA NTRS1991

Capturing, using, and managing quality assurance knowledge for shuttle post-MECO flight design

Ascent initialization values used by the Shuttle's onboard computer for nominal and abort mission scenarios are verified by a six degrees of freedom computer simulation. The procedure that the Ascent Post Main Engine Cutoff (Post-MECO) group uses to perform quality assurance (QA) of the simulation is time consuming. Also, the QA data, checklists and associated rationale, though known by the group members, is not sufficiently documented, hindering transfer of knowledge and problem resolution. A new QA procedure which retains the current high level of integrity while reducing the time required to perform QA is needed to support the increasing Shuttle flight rate. Documenting the knowledge is also needed to increase its availability for training and problem resolution. To meet these needs, a knowledge capture process, embedded into the group activities, was initiated to verify the existing QA checks, define new ones, and document all rationale. The resulting checks were automated in a conventional software program to achieve the desired standardization, integrity, and time reduction. A prototype electronic knowledge base was developed with Macintosh's HyperCard to serve as a knowledge capture tool and data storage.

Peters, H. L.↗
NASA NTRS2011

Can CO2 Turbulent Flux Be Measured by Lidar? A Preliminary Study

The vertical profiling ofCO2 turbulent fluxes in the atmospheric boundary layer (ABL) is investigated using a coherent differential absorption lidar (CDIAL) operated nearby a tall tower in Wisconsin during June 2007. A CDIAL can perform simultaneous range-resolved CO2 DIAL and velocity measurements. The lidar eddy covariance technique is presented. The aims of the study are (i) an assessment of performance and current limitation of available CDIAL for CO2 turbulent fluxes and (ii) the derivation of instrument specifications to build a future CDIAL to perform accurate range-resolved CO2 fluxes. Experimental lidar CO2 mixing ratio and vertical velocity profiles are successfully compared with in situ sensors measurements. Time and space integral scales of turbulence in the ABL are addressed that result in limitation for time averaging and range accumulation. A first attempt to infer CO2 fluxes using an eddy covariance technique with currently available 2-mm CDIAL dataset is reported.

Gilbert, Fabien↗
NASA NTRS1969

Path integrals and Liapunov functionals.

Liapunov functionals for time delay systems by path integrals in state space, using convolution equations involving distributions with compact support

Gruber, M.↗
NASA NTRS2021

Weather Intelligent Navigation Data and Models for Aviation Planning (WINDMAP)

WINDMAP addresses the emerging needs in the aviation community of providing real-time weather forecasting to improve the safety of low altitude aircraft operations. This is accomplished through the integration of real-time observations from autonomous systems, such as drones and urban air taxis, with numerical weather prediction models and flight management and safety systems. To solve this problem, several technical challenges have been identified. These include (1) developing autonomous UAS capable of conducting observations accurately and reliably; (2) determining the number and frequency of required observations and the sensitivity of these observations in data sparse regions of the lower atmosphere;(3) assimilating dense observational data into models in real-time with sufficient resolution and accuracy; (4) developing novel physics-based reduced order models capable of incorporating diverse data sets; and (5)integrating real-time forecasting into UTM and DAA (detect-and-avoid) architectures for path planning and navigation. The goal of this proposed effort is to demonstrate the value of using small UAS to collect measurements of the dynamic and thermodynamic properties of the lower atmosphere at scales that match or exceed the spatio-temporal resolution of today’s best numerical weather prediction models

Koushik Datta↗
NASA NTRS1990

A design algorithm using Z-plane closed loop pole placement

It is argued that the design of a single input single output (SISO) control system is simplified if the designer has direct control over the closed loop poles. A space-laser-communications design example utilizing a two-transform algorithm that constructs a generic discrete compensator is presented. The plant in the example has an anti-aliasing filter, a double integrator, a torquer time constant, and a sample-and-hold as continuous elements, and a discrete proportional plus integral module. This fifth-order system is representative of many instrument pointing control systems. Poles resulting from the compensator are placed near the origin of the Z-plane so that the response is dominated by the closed loop poles from the plant. Root locus, time response, and frequency response data are given. Normalization by the sample interval is used to provide a dimensionless example.

Gatlin, James A.↗
NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗