Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient dynamic analysis”

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 235 records · Page 13

Hybrid methods for rotordynamic analysis

Effective procedures are presented for the response analysis of the Space Shuttle Main Engine turbopumps under transient loading conditions. Of particular concern is the determination of the nonlinear response of the systems to rotor imbalance in presence of bearing clearances. The proposed procedures take advantage of the nonlinearities involved being localized at only a few rotor/housing coupling joints. The methods include those based on integral formulations for the incremental solutions involving the transition matrices of the rotor and housing. Alternatively, a convolutional representation of the housing displacements at the coupling points is proposed which would allow performing the transient analysis on a reduced model of the housing. The integral approach is applied to small dynamical models to demonstrate the efficiency of the approach. For purposes of assessing the numerical integration results for the nonlinear rotor/housing systems, a numerical harmonic balance procedure is developed to enable determining all possible harmonic, subharmonic, and nonperiodic solutions of the systems. A brief account of the Fourier approach is presented as applied to a two degree of freedon rotor-support system.

Noah, Sherif T.↗

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗

On the use of first order rotor dynamics in multiblade coordinates

This paper is directed to the question of how to represent most efficiently rotor/body coupling in a linear flight dynamics analysis. Rigid body pitch, roll and vertical motions are considered for the rotor/body coupling studies. Flapping stability limits, eigenvalues, transient responses to control step inputs, to step gusts and to random gusts are determined for a hypothetical hingeless compound helicopter operating up to .8 advance ratio. Data are obtained for the basic helicopter and for the craft with two simple control feedback systems. While complete periodic system modeling is necessary for determining flapping stability limits and vibrations, constant system modeling using first order dynamics in each of the multiblade rotor coordinates was found to be adequate for rotor-craft stability and response computations.

Hohenemser, K. H.↗

Treatment of Transient Pressure Events in Space Flight Pressurized Systems

The general physics of pressure transients, sources of pressure transients, and major factors influencing pressure transients are addressed. Mitigation strategies that aid in reducing the magnitude of pressure transients are discussed. Methods to characterize pressure transients and their effects on the structure are presented. Case studies illustrate the structural analysis methodology used to predict the structural dynamic response. The structural verification process involves formulating acceptance and qualification programs, which requires knowledge of the stress state. The maximum expected operating pressure (MEOP) is established such that the maximum stress produced by static pressure is equivalent to the maximum stress at the same critical location produced by the combined effect of steady state pressure and the magnitude of the pressure transient. A roadmap on how to treat transient pressures in the structural verification process of spaceflight pressurized systems (e.g., valves, lines, pressure vessels, pressurized structures) is presented.

Transient Pressure Events↗

Electron–Ion Covariance Reveals Geometry-Resolved Inner-Shell Spectra in CS 2 Photodissociation

The chemical shifts of inner-shell atomic orbitals are highly sensitive to chemical bonding and molecular structure. In evolving systems, however, the spectra of distinct molecular species and geometries overlap, obscuring the underlying chemical dynamics. Here we demonstrate the use of electron–ion covariance analysis to combine the structural sensitivity of Coulomb explosion imaging with inner-shell spectroscopy, yielding geometry-specific spectra of transient and product species. We apply this approach to the excited state dynamics of CS 2 probed by ionization above the S 2p edge. Electron–ion covariance with time- and momentum-selected S + and S 2+ ions isolates distinct S 2p photoelectron spectra for ground-state CS 2 , bent photoexcited CS 2 , the CS photoproduct, and bare atomic sulfur ─ species whose spectra overlap strongly in the channel-averaged measurement. Clear chemical shifts are observed in the covariance photoelectron spectrum for each of these species, all of which are consistent with high-level calculations. Here, by extracting the atomic S contribution to the photoelectron spectrum in a finely time-resolved manner we can disentangle this contribution to the overall time-resolved photoelectron spectrum as the photodissociation proceeds. These results demonstrate the promise of electron–ion covariance as a general approach to geometry-resolved inner-shell spectroscopy, opening a route to tracking structural evolution through chemical shifts in complex photoexcited molecules.

Ionization↗

A hierarchical gray-box dynamic modeling methodology for direct-expansion cooling systems to support control stability analysis; Méthodologie de modélisation dynamique hiérarchique de type boîte grise pour les systèmes de refroidissement à détente directe afin de soutenir l'analyse de stabilité de la commande

In this paper, a gray-box dynamic modeling approach for direct-expansion cooling systems is presented. The overall approach incorporates a multi-stage training procedure that consists of 1) identification of component sub-models from quasi-steady-state performance data, 2) system model integration with estimation of refrigerant charge and 3) fine tuning of thermal capacitances of the evaporator and condenser to capture the system dynamic responses. Compared to traditional physics-based models, the proposed modeling approach has advantages including reduced engineering efforts in the model development phase, improved computational efficiency and enhanced prediction accuracy. The modeling method was validated using a 3-ton variable-speed heat pump and proved to be capable of accurately predicting the system transient behaviors over a wide range of operating conditions. The established dynamic model was then applied for control stability analysis, with a specific goal of determining a proper control execution time step. The case study results showed that the stable control execution time step could change significantly, from 3 sec to 19 sec, as the operating conditions and control settings vary, and a proper selection of the execution time step is critical to ensure stable and reliable operations.

42 ENGINEERING↗

A Flexible Multibody Approach to Space Launch System Liftoff Pad Separation with Umbilical Disconnect

A flexible multibody dynamic framework is applied to the Space Launch System(SLS)liftoff Coupled Loads Analysis(CLA), enabling computational efficiencies and systematic inclusion of component nonlinearities. This work simulates the SLS liftoff transient pad separation event inclusive of umbilical disconnects and preload strain energy release (twang). For this system level nonlinear dynamic simulation, the flexible multibody framework provides a systematic approach for the addition of component specific algorithms to enforce the subject nonlinearities. As such, the complexities associated with the nonlinear CLA are greatly reduced, facilitating simplified bookkeeping as well as accelerating the execution of the nonlinear time-domain simulations.

Joel W Sills↗

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗

High-Performance Computing Based EMT Simulation of Large PV or Hybrid PV Plants

Faults in the transmission grid have led to reduced power generation from power electronics resources that are typically not connected to the faulted transmission line. In many of the cases, partial loss of power is observed within the power electronics resources like large photovoltaic (PV) power plants. This phenomena is not captured in existing simulation models and/or simulators. High-fidelity switched system electromagnetic transient (EMT) dynamic models of PV power plants can improve the fidelity of models available for accurate analysis of the impact on PV plants during simulation of faults. However, these models are extremely computationally expensive and take a long time to simulate. Long simulation times limit the ability to use these models as larger regions are studied in EMT simulations with more power electronics resources. In this paper, numerical simulation algorithms are combined with high-performance computing techniques and applied to the high-fidelity switched system EMT model of PV plants. Using these techniques, a speed-up of up to 58x is obtained, while preserving the accuracy of the simulation at greater than 98%.

Debnath, Suman↗

Stabilized Hyperfoam Modeling of the General Plastics EF4003 (3 PCF) Flexible Foam

Constitutive model parameterizations for the General Plastics EF4003 low density 3 pound per cubic foot are needed for design and qualification purposes in normal and abnormal mechanical simulations. The material is expected to be deformed in two ways: first during preloading, and second under impact conditions of the system (transient dynamic). All analyses are to be performed at room temperature. The goal is to provide the analysis community a robust constitutive model parameterization to represent the compression behavior of the EF4003 foam from small deformations up to massive compressive deformations when the foam is densifying. It is worth noting the EF4003 exhibits anisotropy in its stress-strain behavior between the rise and transverse directions (See figure 2.8c-d) as well as plateau behavior that is very likely to cause material stability issues, due to the buckling transition, (and has historically done so) when using Sandia’s current workhorse models for flexible foams, Hyperfoam and Flex Foam. A Stability-informed Hyperfoam parameterization procedure is developed and executed to calibrate a hyperfoam model for the EF4003 room temperature, transversely loaded data. A rise orientation parameterization was not attempted due to localization in the experiments.

36 MATERIALS SCIENCE↗

Thermal Analysis and Design of an Advanced Space Suit

The thermal dynamics and design of an Advanced Space Suit are considered. A transient model of the Advanced Space Suit has been developed and implemented using MATLAB/Simulink to help with sizing, with design evaluation, and with the development of an automatic thermal comfort control strategy. The model is described and the thermal characteristics of the Advanced Space suit are investigated including various parametric design studies. The steady state performance envelope for the Advanced Space Suit is defined in terms of the thermal environment and human metabolic rate and the transient response of the human-suit-MPLSS system is analyzed.

Lin, Chin H.↗

Three-Dimensional, Multi-Phase Unsteady Flow of Water Undergoing Multiple Flow Regimes (Free Surface, Droplet, and Fluid Film)

Upgrades to the Mobile Launcher Rainbird nozzles were implemented for the Artemis II launch to increase water flow rates and improve the mobile launcher deck water coverage while minimizing the water throw distance to avoid impacting the Space Launch System (SLS) rocket nozzles. This analysis describes the use of a three-dimensional, multi-phase, free surface transient computational fluid dynamics (CFD) model to simulate the water flowing out through the Rainbird nozzle and ensure compliance with current Artemis requirements. Flow rates, pressures, velocities, nozzle throw distances and coverage areas were analyzed to ensure that system requirements are met. The completed results were used to make necessary modifications to the Ignition Overpressure Protection and Sound Suppression (IOP/SS) system Rainbird design by reshaping the nozzle used in the system.

CFD↗

Dynamic Modeling and Testing of MSRR-1 for Use in Microgravity Environments Analysis

Delicate microgravity science is unlikely to succeed on the International Space Station if vibratory and transient disturbers corrupt the environment. An analytical approach to compute the on-orbit acceleration environment at science experiment locations within a standard payload rack resulting from these disturbers is presented. This approach has been grounded by correlation and comparison to test verified transfer functions. The method combines the results of finite element and statistical energy analysis using tested damping and modal characteristics to provide a reasonable approximation of the total root-mean-square (RMS) acceleration spectra at the interface to microgravity science experiment hardware.

Gattis, Christy↗

DAC-3 Pointing Stability Analysis Results for SAGE 3 and Other Users of the International Space Station (ISS) Payload Attachment Sites (PAS)

The purpose of this paper is to provide final results of a pointing stability analysis for external payload attachment sites (PAS) on the International Space Station (ISS). As a specific example, the pointing stability requirement of the SAGE III atmospheric science instrument was examined in this paper. The instrument requires 10 arcsec stability over 2 second periods. SAGE 3 will be mounted on the ISS starboard side at the lower, outboard FIAS. In this engineering analysis, an open-loop DAC-3 finite element model of ISS was used by the Microgravity Group at Johnson Space Flight Center to generate transient responses at PAS to a limited number of disturbances. The model included dynamics up to 50 Hz. Disturbance models considered included operation of the solar array rotary joints, thermal radiator rotary joints, and control moment gyros. Responses were filtered to model the anticipated vibration attenuation effects of active control systems on the solar and thermal radiator rotary joints. A pointing stability analysis was conducted by double integrating acceleration transient over a 2 second period. Results of the analysis are tabulated for ISS X, Y, and Z Axis rotations. These results indicate that the largest excursions in rotation during pointing occurred due to rapid slewing of the thermal radiator. Even without attenuation at the rotary joints, the resulting pointing error was limited to less than 1.6 arcsec. With vibration control at the joints, to a maximum 0.5 arcsec over a 2 second period. Based on this current level of model definition, it was concluded that between 0 - 50 Hz, the pointing stability requirement for SAGE 3 will not be exceeded by the disturbances evaluated in this study.

Woods-Vedeler, Jessica A.↗

Method Improvements in Thermal Analysis of Mach 10 Leading Edges

Several improvements have recently been made in the thermal analysis methods for leading edges of a hypersonic vehicle. The leading edges of this vehicle undergo exceptionally high heat loads that incorporate extreme spatial gradients as well as severe transients. Due to the varying flight conditions, complex geometry, and need for thermal loads at many points along the trajectory, full computational fluid dynamics (CFD) analysis of the aeroheating; loads is not feasible. Thus, engineering methods must be used to determine the aeroheating on the vehicle surfaces, and that must be utilized in the thermal analysis. Over the last year, the thermal analysis of a hypersonic vehicle has been enhanced in several ways. Two different engineering codes are used to predict aeroheating loads: one over the curve near the stagnation point, and the other on flat surfaces downstream of the leading edge. These two are matched together at the intersection point using a method that allows closer approximation of CFD results. User-developed FORTRAN, which is part of the thermal solver PATRAN Thermal, is used to accomplish this. The customizable FORTRAN code also allows use of many different time- and space-dependent factors, interpolation of the heat load in time and space, and inclusion of both highly swept and unswept grid structures. This FORTRAN is available to other PATRAN users who may want to accomplish a similar objective in analysis. Flux, rather than convective coefficient, is used to define heat loads, which allows more accurate analysis as well as better application of margins. Improvements have also been made in more efficient utilization of imported CAD geometry, by creating faces on solids to facilitate load application.

Amundsen, Ruth M.↗

Simulation and generalized Langevin equation study of lipid subdiffusion in biomembrane phases

In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.

Malik, Sheeba [Oak Ridge National Laboratory (ORNL↗

Real-World 20-Hz IBR Subsynchronous Oscillations: Signatures and Mechanism Analysis

This article provides analyses of subsynchronous oscillations associated with solar photovoltaics (PVs) observed in real world. All have oscillation frequencies at approximately 20 Hz in voltage and current magnitudes. On the other hand, in phase voltage and current, the oscillations show very different signatures. Through constructing simplified linear models and conducting sensitivity analysis, we are able to associate dynamic phenomena with the underlying mechanism and speculate the possible causes of the 20-Hz oscillations. The analysis results have all been validated via simulation results from electromagnetic transient (EMT) test beds constructed in MATLAB/Simscape.

inverter-based resources↗