Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic models”

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

Clumpy molecular clouds - A dynamic model self-consistently regulated by T Tauri star formation

A model is proposed which can account for the longevity, energetics, and dynamical structure of dark molecular clouds. Recent observations of a high space density of T Tauri stars in some dark clouds provide the basis for the assertion that high-velocity winds from these low mass pre-main-sequence stars provide a continuous dynamic input into molecular clouds. It is postulated that as clumps are driven above the Jeans mass, both by coalescence and the enhancement of ram pressure through continuing acceleration by protostellar winds, collapse is followed by the formation of low-mass stars that generate additional protostellar winds. Finally, it is found that star formation occurs on a relatively slow time scale, comparable to the cloud lifetime

Norman, C.↗

Improved Temperature Dynamic Model of Turbine Subcomponents for Facilitation of Generalized Tip Clearance Control

This paper describes efforts conducted to improve dynamic temperature estimations of a turbine tip clearance system to facilitate design of a generalized tip clearance controller. This work builds upon research previously conducted and presented in and focuses primarily on improving dynamic temperature estimations of the primary components affecting tip clearance (i.e. the rotor, blades, and casing/shroud). The temperature profiles estimated by the previous model iteration, specifically for the rotor and blades, were found to be inaccurate and, more importantly, insufficient to facilitate controller design. Some assumptions made to facilitate the previous results were not valid, and thus improvements are presented here to better match the physical reality. As will be shown, the improved temperature sub- models, match a commercially validated model and are sufficiently simplified to aid in controller design.

Kypuros, Javier A.↗

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

Modeling Dynamic Fracture of Cryogenic Pellets

This work is part of an investigation with the long-range objective of predicting the size distribution function and velocity dispersion of shattered pellet fragments after a large cryogenic pellet impacts a solid surface at high velocity. The study is vitally important for the shattered pellet injection (SPI) technique, one of the leading technologies being implemented at ORNL for the mitigation of disruption damage on current tokamaks and ITER. The report contains three parts that are somewhat interwoven. In Part I we formulated a self-similar model for the expansion dynamics and velocity dispersion of the debris cloud following pellet impact against a thick (rigid) target plate. Also presented in Part I is an analytical fracture model that predicts the nominal or mean size of the fragments in the debris cloud and agrees well with known SPI data. The aim of Part II is to gain an understanding of the pellet fracturing process when a pellet is shattered inside a miter tube with a sharp bend. Because miter tubes have a thin stainless steel (SS) wall a permanent deformation (dishing) of the wall is produced at the site of the impact. A review of the literature indicates that most projectile impact on thin plates are those for which the target is deformed and the projectile is perfectly rigid. Such impacts result in “projectile embedding” where the projectile speed is reduced to zero during the interaction so that all the kinetic energy (KE) of the projectile goes into the energy stored in plastic deformation. Much of the literature deals with perforation of the target. The problem here is quite different; the softer pellet easily undergoes complete material failure causing only a small transfer of KE to stored energy of wall deformation. For the real miter tube, we derived a strain energy function for the wall deflection using a non-linear (plastic) stress-strain relation for 304 SS. Using a dishing profile identical to the linear Kirchkoff-Love profile (for lack of a rigorously derived profile) we derived the strain energy associated with the deflection and applied a virtual work principle to find a relationship between the impact (load) pressure to the measured wall deflection depth. The inferred impact pressure was in good agreement with the expected pressure for oblique cryogenic pellet impacts where the pellet shear stress causing cleavage fracture is well above the yield stress for pure shear. The section is concluded with additional discussion on how this wall deformation data lends further support to the analytical fracture model presented in Part I. In Part III we present three different size distribution models. A summary, with a few brief suggestions for a follow on study, is provided at the end of this report.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multibody dynamics: Modeling component flexibility with fixed, free, loaded, constraint, and residual modes

The assumed-modes method in multibody dynamics allows the elastic deformation of each component in the system to be approximated by a sum of products of spatial and temporal functions commonly known as modes and modal coordinates respectively. The choice of component modes used to model articulating and non-articulating flexible multibody systems is examined. Attention is directed toward three classical Component Mode Synthesis (CMS) methods whereby component normal modes are generated by treating the component interface (I/F) as either fixed, free, or loaded with mass and stiffness contributions from the remaining components. The fixed and free I/F normal modes are augmented by static shape functions termed constraint and residual modes respectively. A mode selection procedure is outlined whereby component modes are selected from the Craig-Bampton (fixed I/F plus constraint), MacNeal-Rubin (free I/F plus residual), or Benfield-Hruda (loaded I/F) mode sets in accordance with a modal ordering scheme derived from balance realization theory. The success of the approach is judged by comparing the actuator-to-sensor frequency response of the reduced order system with that of the full order system over the frequency range of interest. A finite element model of the Galileo spacecraft serves as an example in demonstrating the effectiveness of the proposed mode selection method.

Spanos, John T.↗

Stellar and gas dynamical model for tidal disruption events in a quiescent galaxy

A detailed model of the tidal disruption events (TDEs) has been constructed using stellar dynamical and gas dynamical inputs that include black hole (BH) mass M{sub •}, specific orbital energy E and angular momentum J, star mass M{sub ⋆} and radius R{sub ⋆,} and the pericenter of the star orbit r{sub p}(E, J, M{sub •}). We solved the steady state Fokker–Planck equation using the standard loss cone theory for the galactic density profile ρ (r) ∝ r{sup −γ} and stellar mass function ξ(m) where m = M{sub ⋆}/M{sub ⊙} and obtained the feeding rate of stars to the BH integrated over the phase space as N-dot {sub t}∝M{sub ∙}{sup β}, where β=−0.3±0.01 for M{sub •} > 10{sup 7}M{sub ⊙} and ∼6.8 × 10{sup −5} Yr{sup −1} for γ = 0.7. We use this to model the in-fall rate of the disrupted debris, M-dot (E,J,m,t), and discuss the conditions for the disk formation, finding that the accretion disk is almost always formed for the fiduciary range of the physical parameters. We also find the conditions under which the disk formed from the tidal debris of a given star with a super Eddington accretion phase. We have simulated the light curve profiles in the relevant optical g band and soft X-rays for both super and sub-Eddington accretion disks as a function of M-dot (E,J,t). Using this, standard cosmological parameters, and mission instrument details, we predict the detectable TDE rates for various forthcoming surveys as a function of γ.

79 ASTRONOMY AND ASTROPHYSICS↗

Methodology for Thermodynamic Analysis Coupled with Computational Fluid Dynamics Modeling for Casting a Novel Aluminum–Cerium Alloy

In this paper, we present the results of computational fluid dynamics (CFD) analysis to assess castability, porosity, molten metal fluidity and other technological properties of a new Al-Ce alloy. A thermodynamic analysis of a new Al-Ce alloy is done to obtain the physical properties as a function of temperature across both solid and liquid phases. These properties are then used to build the CFD of model the full casting process from initial pouring through to final solidification of the part, here an example of a heavy duty mor mount is used. While several excellent commercial CFD codes exists this process shows that CFD can be used to assess the capability of the alloy to properly fill the mold, as well as give predictions where scattered porosity or large-scale defects may occur in the casting. Further, unlike the commercially available software (e.g., ProCAST, SOLIDCast , MAGMASOFT®) the complex 3D-analysis of the stress / strain fields in cast parts is not performed at this time. However, the availability of free software for assessing the required thermodynamic / thermo-physical properties of new alloys (OpenCALPHAD) and CFD codes such as OpenFOAM® makes the developed option attractive and economical, especially for the analysis of new Al-Ce alloys, for which the available data does not exist.

36 MATERIALS SCIENCE↗

Crowd Dynamics: Modeling and Control of Multiagent Systems

This review aims to present recent developments in modeling and control of multiagent systems. A particular focus is set on crowd dynamics characterized by complex interactions among agents, also called social interactions, and large-scale systems. Specifically, in a crowd each individual agent interacts with a field generated by the other agents and the environment. These systems can be modeled at the microscopic scale by ordinary differential equations, while an alternative description at the mesoscopic scale is given by a partial differential equation for the propagation of the probability density of the agents. Control actions can be applied at the individual level as well as at the level of the corresponding fields. This article presents and compares different control types, and the specific application to multilane, multiclass traffic is developed in some detail, showing the main tools at work in a hybrid setting with relevant impacts on autonomous driving.

33 ADVANCED PROPULSION SYSTEMS↗

Cluster Dynamics Modeling with Bubble Nucleation, Growth and Coalescence

The topic of this communication pertains to defect formation in irradiated solids such as plasma-facing tungsten submitted to helium implantation in fusion reactor com- ponents, and nuclear fuel (metal and oxides) submitted to volatile ssion product generation in nuclear reactors. The purpose of this progress report is to describe ef- forts towards addressing the prediction of long-time evolution of defects via continuum cluster dynamics simulation. The di culties are twofold. First, realistic, long-time dynamics in reactor conditions leads to a non-dilute di usion regime which is not accommodated by the prevailing dilute, stressless cluster dynamics theory. Second, long-time dynamics calls for a large set of species (ideally an in nite set) to capture all possible emerging defects, and this represents a computational bottleneck. Extensions beyond the dilute limit is a signi cant undertaking since no model has been advanced to extend cluster dynamics to non-dilute, deformable conditions. Here our proposed approach to model the non-dilute limit is to monitor the appearance of a spatially localized void volume fraction in the solid matrix with a bell shape pro le and insert an explicit geometrical bubble onto the support of the bell function. The newly cre- ated internal moving boundary provides the means to account for the interfacial ux of mobile species into the bubble, and the growth of bubbles allows for coalescence phenomena which captures highly non-dilute interactions. We present a preliminary interfacial kinematic model with associated interfacial di usion transport to follow the evolution of the bubble in any number of spatial dimensions and any number of bubbles, which can be further extended to include a deformation theory. Finally we comment on a computational front-tracking method to be used in conjunction with conventional cluster dynamics simulations in the non-dilute model proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A coupled vacancy diffusion-dislocation dynamics model for the climb-glide motion of jogged screw dislocations

Here, we develop a novel model to study the climb/glide motion of jogged screw dislocations within the discrete dislocation dynamics (DDD) framework. We present results for the dependence of the climb velocity on the applied stress and on the jog size and distribution statistics. We show that the model predictions are consistent with experimental data in both γ-TiAl and Zircaloy-4. The ranges of the applied stress and jog spacing that determine the dominance of one of three dislocation mechanisms are identified. These are the jog dragging, dipole dragging, and dipole bypass mechanisms, respectively. The overall dislocation motion in the jog dragging regime is composed of glide of screw segments and climb of jogs, controlling the plastic strain and the creep rate, respectively. Based on current simulations and on a detailed examination of the predicted jog heights compared to experiments, we advance the hypothesis that a combination of jog dragging and dipole bypass mechanisms is necessary to reproduce the high creep rate observed in some experiments.

36 MATERIALS SCIENCE↗

SSME structural dynamic model development

The high pressure fuel turbopump (HPFTP) is a major component of the Space Shuttle Main Engine (SSME) powerhead. The device is a three stage centrifugal pump that is directly driven by a two stage hot gas turbine. The purpose of the pump is to deliver fuel (liquid hydrogen) from the low pressure fuel turbopump (LPFTP) through the main fuel valve (MFV) to the thrust chamber coolant circuits. In doing so, the pump pressurizes the fuel from an inlet pressure of approximately 178 psi to a discharge pressure of over 6000 psi. At full power level (FPL), the pump rotates at a speed of over 37,000 rpm while generating approximately 77,000 horsepower. Obviously, a pump failure at these speeds and power levels could jeopardize the mission. Results are summarized for work in which the solutions obtained from analytical models of the fuel turbopump impellers are compared with the results obtained from dynamic tests.

Foley, Michael J.↗

A dynamic model of Venus's gravity field

Unlike Earth, long wavelength gravity anomalies and topography correlate well on Venus. Venus's admittance curve from spherical harmonic degree 2 to 18 is inconsistent with either Airy or Pratt isostasy, but is consistent with dynamic support from mantle convection. A model using whole mantle flow and a high viscosity near surface layer overlying a constant viscosity mantle reproduces this admittance curve. On Earth, the effective viscosity deduced from geoid modeling increases by a factor of 300 from the asthenosphere to the lower mantle. These viscosity estimates may be biased by the neglect of lateral variations in mantle viscosity associated with hot plumes and cold subducted slabs. The different effective viscosity profiles for Earth and Venus may reflect their convective styles, with tectonism and mantle heat transport dominated by hot plumes on Venus and by subducted slabs on Earth. Convection at degree 2 appears much stronger on Earth than on Venus. A degree 2 convective structure may be unstable on Venus, but may have been stabilized on Earth by the insulating effects of the Pangean supercontinental assemblage.

Kiefer, W. S.↗

A dynamic model of Venus's gravity field

Unlike Earth, long wavelength gravity anomalies and topography correlate well on Venus. Venus's admittance curve from spherical harmonic degree 2 to 18 is inconsistent with either Airy or Pratt isostasy, but is consistent with dynamic support from mantle convection. A model using whole mantle flow and a high viscosity near surface layer overlying a constant viscosity mantle reproduces this admittance curve. On Earth, the effective viscosity deduced from geoid modeling increases by a factor of 300 from the asthenosphere to the lower mantle. These viscosity estimates may be biased by the neglect of lateral variations in mantle viscosity associated with hot plumes and cold subducted slabs. The different effective viscosity profiles for Earth and Venus may reflect their convective styles, with tectonism and mantle heat transport dominated by hot plumes on Venus and by subducted slabs on Earth. Convection at degree 2 appears much stronger on Earth than on Venus. A degree 2 convective structure may be unstable on Venus, but may have been stabilized on Earth by the insulating effects of the Pangean supercontinental assemblage.

Kiefer, W. S.↗

Dynamical modeling of serial manipulators with flexible links and joints using the method of kinematic influence

A method of formulating the dynamical equations of a flexible, serial manipulator is presented, using the Method of Kinematic Influence. The resulting equations account for rigid body motion, structural motion due to link and joint flexibilities, and the coupling between these two motions. Nonlinear inertial loads are included in the equations. A finite order mode summation method is used to model flexibilities. The structural data may be obtained from experimental, finite element, or analytical methods. Nonlinear flexibilities may be included in the model.

Graves, Philip L.↗

Identification of helicopter rotor dynamic models

A recursive, extended Kalman-filter approach is applied to the identifiction of rotor damping levels of representative helicopter dynamic systems. The general formulation of the approach is presented in the context of a typically posed stochastic estimation problem, and the method is analytically applied to determining the damping levels of a coupled rotor-body system. The identified damping covergence characteristics are studied for sensitivity to both constant-coefficient and periodic-coefficient measurement models, process-noise covariance levels, and specified initial estimates of the rotor-system damping. A second application of the method to identifying the plant model for a highly damped, isolated flapping blade with a constant-coefficient state model (hover) and a periodic-coefficient state model (forward flight) is also investigated. The parameter-identification capability is evaluated for the effect of periodicity on the plant model coefficients and the influence of different measurement noise levels.

Molusis, J. A.↗