Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Procedures”

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 55 records · Page 3

Evaluation of constraint stabilization procedures for multibody dynamical systems

Comparative numerical studies of four constraint treatment techniques for the simulation of general multibody dynamic systems are presented, and results are presented for the example of a classical crank mechanism and for a simplified version of the seven-link manipulator deployment problem. The staggered stabilization technique (Park, 1986) is found to yield improved accuracy and robustness over Baumgarte's (1972) technique, the singular decomposition technique (Walton and Steeves, 1969), and the penalty technique (Lotstedt, 1979). Furthermore, the staggered stabilization technique offers software modularity, and the only data each solution module needs to exchange with the other is a set of vectors plus a common module to generate the gradient matrix of the constraints, B.

Park, K. C.↗

A computational procedure for large rotational motions in multibody dynamics

A computational procedure suitable for the solution of equations of motion for multibody systems is presented. The present procedure adopts a differential partitioning of the translational motions and the rotational motions. The translational equations of motion are then treated by either a conventional explicit or implicit direct integration method. A principal feature of this procedure is a nonlinearly implicit algorithm for updating rotations via the Euler four-parameter representation. This procedure is applied to the rolling of a sphere through a specific trajectory, which shows that it yields robust solutions.

Park, K. C.↗

Stabilization of computational procedures for constrained dynamical systems

A new stabilization method of treating constraints in multibody dynamical systems is presented. By tailoring a penalty form of the constraint equations, the method achieves stabilization without artificial damping and yields a companion matrix differential equation for the constraint forces; hence, the constraint forces are obtained by integrating the companion differential equation for the constraint forces in time. A principal feature of the method is that the errors committed in each constraint condition decay with its corresponding characteristic time scale associated with its constraint force. Numerical experiments indicate that the method yields a marked improvement over existing techniques.

Park, K. C.↗

A computational procedure for multibody systems including flexible beam dynamics

A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. The flexible beams are modeled using a fully nonlinear theory which accounts for both finite rotations and large deformations. The present formulation incorporates physical measures of conjugate Cauchy stress and covariant strain increments. As a consequence, the beam model can easily be interfaced with real-time strain measurements and feedback control systems. A distinct feature of the present work is the computational preservation of total energy for undamped systems; this is obtained via an objective strain increment/stress update procedure combined with an energy-conserving time integration algorithm which contains an accurate update of angular orientations. The procedure is demonstrated via several example problems.

Downer, J. D.↗

Developments and applications of dynamic models for large eddy simulation of complex flows

The dynamic modeling procedure for large eddy simulation of turbulent flows is reviewed and recent developments in the theoretical aspects and applications are described. Methods for inclusion of backscatter of energy from small to large scale motions are presented. New formulations of the dynamic procedure are proposed which are optimized based on the subgrid scale flux vector or the energy dissipation rate instead of the subgrid scale stress tensor. Recent results from application of the model to forced isotropic turbulence with an inertial subrange, flow over a backward facing step at Reynolds number of 28000, and flow over a concave curved surface are presented.

Moin, P.↗

Verification of an improved computational design procedure for TWT-dynamic refocuser-MDC systems with secondary electron emission losses

A computational procedure for the design of TWT-refocuser-MDC systems was used to design a short 'dynamic' refocusing system and highly efficient four-stage depressed collector for a 200-W 8-18-GHz TWT. The computations were carried out with advanced multidimensional computer programs which model the electron beam as a series of disks of charge and follow their trajectories from the RF input of the TWT, through the slow-wave structure and refocusing section, to their points of impact in the depressed collector. Secondary emission losses in the MDC were treated semiquantitatively by injecting a representative beam of secondary electrons into the MDC analysis at the point of impact of each primary beam. A comparison of computed and measured TWT and MDC performance showed very good agreement. The electrodes of the MDC were fabricated from a particular form of isotropic graphite that was selected for its low secondary electron yield, thermal expansion characteristics, ease of machinability and vacuum properties. This MDC was tested at CW for more than 1000 h with negligible degradation in TWT and MDC performances.

Ramins, P.↗

Advanced superposition methods for high speed turbopump vibration analysis

The small, high pressure Mark 48 liquid hydrogen turbopump was analyzed and dynamically tested to determine the cause of high speed vibration at an operating speed of 92,400 rpm. This approaches the design point operating speed of 95,000 rpm. The initial dynamic analysis in the design stage and subsequent further analysis of the rotor only dynamics failed to predict the vibration characteristics found during testing. An advanced procedure for dynamics analysis was used in this investigation. The procedure involves developing accurate dynamic models of the rotor assembly and casing assembly by finite element analysis. The dynamically instrumented assemblies are independently rap tested to verify the analytical models. The verified models are then combined by modal superposition techniques to develop a completed turbopump model where dynamic characteristics are determined. The results of the dynamic testing and analysis obtained are presented and methods of moving the high speed vibration characteristics to speeds above the operating range are recommended. Recommendations for use of these advanced dynamic analysis procedures during initial design phases are given.

Nielson, C. E.↗

Constraint treatment techniques and parallel algorithms for multibody dynamic analysis

Computational procedures for kinematic and dynamic analysis of three-dimensional multibody dynamic (MBD) systems are developed from the differential-algebraic equations (DAE's) viewpoint. Constraint violations during the time integration process are minimized and penalty constraint stabilization techniques and partitioning schemes are developed. The governing equations of motion, a two-stage staggered explicit-implicit numerical algorithm, are treated which takes advantage of a partitioned solution procedure. A robust and parallelizable integration algorithm is developed. This algorithm uses a two-stage staggered central difference algorithm to integrate the translational coordinates and the angular velocities. The angular orientations of bodies in MBD systems are then obtained by using an implicit algorithm via the kinematic relationship between Euler parameters and angular velocities. It is shown that the combination of the present solution procedures yields a computationally more accurate solution. To speed up the computational procedures, parallel implementation of the present constraint treatment techniques, the two-stage staggered explicit-implicit numerical algorithm was efficiently carried out. The DAE's and the constraint treatment techniques were transformed into arrowhead matrices to which Schur complement form was derived. By fully exploiting the sparse matrix structural analysis techniques, a parallel preconditioned conjugate gradient numerical algorithm is used to solve the systems equations written in Schur complement form. A software testbed was designed and implemented in both sequential and parallel computers. This testbed was used to demonstrate the robustness and efficiency of the constraint treatment techniques, the accuracy of the two-stage staggered explicit-implicit numerical algorithm, and the speed up of the Schur-complement-based parallel preconditioned conjugate gradient algorithm on a parallel computer.

Chiou, Jin-Chern↗

A simplified signal analysis technique for obtaining optimal estimates of system dynamics

A procedure is presented which provides an algorithm for the special case where the unknown system has n inputs, and n finite length realizations for each of the n input signals are available. Given an n input system and n experimental records for each input, it is shown that the optimum Wiener filter is easily computed using only the experimental data. No assumptions concerning the analytical structure of the spectral matrix are necessary.

Whitbeck, R. F.↗

An efficient procedure for computing partial ionization of a gas in numerical fluid dynamics

A numerical procedure is described that allows efficient solution for the three nonlinear and implicit equations of state relating the equilibrium properties of a partially (and singly) ionized neutral gas. In the computation, the values for the intensive variables, that is pressure, temperature and ionization coefficient, are iteratively obtained with a given accuracy starting from assigned values for energy and volume per unit mass. Iterative computation is performed only when the value for the ionization coefficient lies within the lower and the upper bounds fixed by the required number of significant digits; otherwise the two explicit equations of state for nonionized or for fully ionized gases are employed.

Garribba, S.↗

Diabatic dynamic initialization

A generalized dynamical adjustment procedure was applied to a diabatic model to produce balanced initial conditions. Namely, backward adiabatic model integration is followed by forward diabatic model integration, with a high frequency (low pass) filter in the form of the Euler-backward time differencing scheme being applied throughout the whole integration. As a result of the application of such a diabatic dynamic initialization procedure within the Goddard Laboratory for Atmospheres (GLA) 4-D data assimilation system, the following properties of forecasts from initialized fields are achieved right from the beginning of the usual forecast integration: (1) the forecast tendencies (and fields) are free of any noise due to imbalance in initial conditions; (2) the shocks related to an initial imbalance between model physics and dynamics, and especially the substantial initial imbalance of precipitation and evaporation fields, or the initial spin-up effect, are practically removed. Diabatic dynamic initialization has been compared with implicit nonlinear normal mode initialization, and found to be superior in removing the initial spin-up effect, and also in improving the tropical structure. The diabatic dynamic initialization procedure was successfully tested for the GLA system with the use of all conventional data and the GLA satellite data retrievals. It allows a smooth data insertion without any shocks or imbalances, which is highly desirable for efficient functioning of 4-D data assimilation systems. The developed initialization procedure is computationally efficient and in principle easily applicable to different forecast models.

Fox-Rabinovitz, Michael S.↗

eVTOL Vehicle-Agnostic Instrument Flight Procedures Test Plan

NASA Advanced Air Mobility National Campaign is researching the utility of electric vertical takeoff and land (eVTOL) advanced air mobility (AAM) instrument flight procedures. The result will be dynamic and tailored procedures that align to the following modus operandi: maximize safety, optimize efficiency, support passenger comfort and minimize acoustics. This is achieved through dynamic airspace procedure design, which is a modular approach to create an airspace construct that customizes procedures to vehicle design and configuration, operation and environmental conditions. The test plan supports different eVTOL platforms and envisioned operations for flight test or simulation and may be leveraged by AAM aircraft manufacturers and operators for any given aircraft, location and operation.

eVTOL test plan↗

A Framework for Dynamic Constraint Reasoning Using Procedural Constraints

Many complex real-world decision and control problems contain an underlying constraint reasoning problem. This is particularly evident in a recently developed approach to planning, where almost all planning decisions are represented by constrained variables. This translates a significant part of the planning problem into a constraint network whose consistency determines the validity of the plan candidate. Since higher-level choices about control actions can add or remove variables and constraints, the underlying constraint network is invariably highly dynamic. Arbitrary domain-dependent constraints may be added to the constraint network and the constraint reasoning mechanism must be able to handle such constraints effectively. Additionally, real problems often require handling constraints over continuous variables. These requirements present a number of significant challenges for a constraint reasoning mechanism. In this paper, we introduce a general framework for handling dynamic constraint networks with real-valued variables, by using procedures to represent and effectively reason about general constraints. The framework is based on a sound theoretical foundation, and can be proven to be sound and complete under well-defined conditions. Furthermore, the framework provides hybrid reasoning capabilities, as alternative solution methods like mathematical programming can be incorporated into the framework, in the form of procedures.

Jonsson, Ari K.↗

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory-Helioseismic and Magnetic Imager Time-Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a time - distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time - distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference traveltime perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.

HMI↗

Galerkin/least-squares procedures in computational fluid dynamics

The application of Galerkin/least-squares methodology to computational fluid dynamics is reviewed. Applications to advective-diffusive systems of equations and incompressible and compressible Euler and Navier-Stokes equations are described. Space-time finite-element methodology for time-dependent problems is emphasized and the status and role of so-called 'discontinuity-capturing operators' is reviewed.

Hughes, Thomas J. R.↗