Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Simulation”

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,117 records · Page 62

A PC-based collapse factor analysis for high pressure cryogenic systems

A collapse factor simulation computer program, suitable for use on personal computers, for supercritical cryogen system pressurization is examined. The available literature is previewed, individual component models are developed, and the overall simulation program logic is examined. Extensive use is made of existing thermodynamic and transport property evaluation packages. Heat transfer is estimated using natural convection correlations, and transient wall heat conduction is considered. Examples of the simulation are demonstrated for two facilities. Partial validation is achieved by comparing the simulation output with data from an existing test facility. Preliminary results indicate that the collapse factors will seldom exceed 1.15 during supercritical operation of high-pressure cryogenic systems. Recommendations for improvements in the component models and the simulation are made.

Hodge, B. K.↗

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗

A study of ignition and simulation circuits for arcjet thrusters, part 1

A 1 kW electronic load was programmed to simulate the nonlinear i-v (volt-ampere) characteristics of an arcjet, both ignited and unignited. The simulator was tested and found to closely resemble an arcjet both for large transients and small perturbances up to about 40 kHz. No attempt was made to simulate the ignition process itself. The dynamic behavior of the arcjet (and the simulator) was shown to differ significantly from that of a resistor bank. Previous research led to the design and construction of a 1 kW arcjet power supply. A high voltage ignition circuit was added to this hardware, and tests on a 1 kW arcjet were performed at NASA-Lewis. All tests were successful and no ignition failures were observed. Circuit documentation and test results are included.

Stuart, Thomas A.↗

Design Power Control Strategies of Grid-Forming Inverters for Microgrid Application: Preprint

This paper compares two control schemes in the application control layer of a non-PLL grid-forming (GFM) inverter to gain the insights and understanding of how the two schemes affect the dynamic responses of the GFM inverters and the transition operation of microgrids. Both schemes can achieve smooth microgrid transition operation and power tracking, however, the first scheme expects to exhibit more transients due to the feedforward of the inverse dynamic model of the output filter. Simulation results show that both control schemes achieve successful transition operation with smooth transients, and the second control scheme exhibits slightly better transients in the system frequency during islanding operation and the first one has slightly superior transients in the system frequency during reconnection operation. The power tracking performance in the grid-connected mode is evaluated, which shows that the first scheme has better active power tracking and the second one has better reactive power tracking. The analysis and results are useful to develop reliable control schemes for non-PLL GFM inverters because more and more inverters will work as non-PLL grid-forming sources in the future grids due to better stability and reliability.

grid-forming inverter↗

Dynamic Altitude Simulation System Performance Modeling

The "Dynamic Altitude Simulation Prediction Program," DASSPP, is a program to predict the transient response of an engine, test cell, ejector system under engine shutdown conditions. These transients are important to know so that corrective modifications can be adapted to prevent any damage to the engine or test cell. The "Dynamic Altitude System Simulation Prediction Program," DASSPP, is a major rewrite of the existing program "RL-1000" written in BASIC. The RL-1000 program was written to analyze the transients of the RL-10 system only. The new program is written to run in Excell 97 and utilizes the Visual BASIC language in Excell. The program has many added features not included in the original "Rl-1000" program. The program utilizes the ejector models developed during the summer of 1997. The new program is very user friendly and utilizes a dialog box for data input.

LaFrance, Leo J.↗

Crash Simulation of a Boeing 737 Fuselage Section Vertical Drop Test

A 30-ft/s vertical drop test of a fuselage section of a Boeing 737 aircraft was conducted in October of 1999 at the FAA Technical Center in Atlantic City, NJ. This test was performed to evaluate the structural integrity of a conformable auxiliary fuel tank mounted beneath the floor and to determine its effect on the impact response of the airframe structure and the occupants. The test data were used to compare with a finite element simulation of the fuselage structure and to gain a better understanding of the impact physics through analytical/experimental correlation. To perform this simulation, a full-scale 3-dimensional finite element model of the fuselage section was developed using the explicit, nonlinear transient-dynamic finite element code, MSC.Dytran. The emphasis of the simulation was to predict the structural deformation and floor-level acceleration responses obtained from the drop test of the B737 fuselage section with the auxiliary fuel tank.

Fasanella, Edwin L.↗

Imaging Sensor Flight and Test Equipment Software

The Lightning Imaging Sensor (LIS) is one of the components onboard the Tropical Rainfall Measuring Mission (TRMM) satellite, and was designed to detect and locate lightning over the tropics. The LIS flight code was developed to run on a single onboard digital signal processor, and has operated the LIS instrument since 1997 when the TRMM satellite was launched. The software provides controller functions to the LIS Real-Time Event Processor (RTEP) and onboard heaters, collects the lightning event data from the RTEP, compresses and formats the data for downlink to the satellite, collects housekeeping data and formats the data for downlink to the satellite, provides command processing and interface to the spacecraft communications and data bus, and provides watchdog functions for error detection. The Special Test Equipment (STE) software was designed to operate specific test equipment used to support the LIS hardware through development, calibration, qualification, and integration with the TRMM spacecraft. The STE software provides the capability to control instrument activation, commanding (including both data formatting and user interfacing), data collection, decompression, and display and image simulation. The LIS STE code was developed for the DOS operating system in the C programming language. Because of the many unique data formats implemented by the flight instrument, the STE software was required to comprehend the same formats, and translate them for the test operator. The hardware interfaces to the LIS instrument using both commercial and custom computer boards, requiring that the STE code integrate this variety into a working system. In addition, the requirement to provide RTEP test capability dictated the need to provide simulations of background image data with short-duration lightning transients superimposed. This led to the development of unique code used to control the location, intensity, and variation above background for simulated lightning strikes at user-selected locations.

Freestone, Kathleen↗

Use of Heated Helium to Simulate Surface Pressure Fluctuations on the Launch Abort Vehicle During Abort Motor Firing

The solid-rocket plumes from the Abort motor of the Multi-Purpose Crew Vehicle (MPCV, also know as Orion) were simulated using hot, high pressure, Helium gas to determine the surface pressure fluctuations on the vehicle in the event of an abort. About 80 different abort situations over a wide Mach number range, (0.3< or =M< or =1.2) and vehicle attitudes (+/-15deg) were simulated inside the NASA Ames Unitary Plan, 11-Foot Transonic Wind Tunnel. For each abort case, typically two different Helium plume and wind tunnel conditions were used to bracket different flow matching critera. This unique, yet cost-effective test used a custom-built hot Helium delivery system, and a 6% scale model of a part of the MPCV, known as the Launch Abort Vehicle. The test confirmed the very high level of pressure fluctuations on the surface of the vehicle expected during an abort. In general, the fluctuations were found to be dominated by the very near-field hydrodynamic fluctuations present in the plume shear-layer. The plumes were found to grow in size for aborts occurring at higher flight Mach number and altitude conditions. This led to an increase in the extent of impingement on the vehicle surfaces; however, unlike some initial expectations, the general trend was a decrease in the level of pressure fluctuations with increasing impingement. In general, the highest levels of fluctuations were found when the outer edges of the plume shear layers grazed the vehicle surface. At non-zero vehicle attitudes the surface pressure distributions were found to become very asymmetric. The data from these wind-tunnel simulations were compared against data collected from the recent Pad Abort 1 flight test. In spite of various differences between the transient flight situation and the steady-state wind tunnel simulations, the hot-Helium data were found to replicate the PA1 data fairly reasonably. The data gathered from this one-of-a-kind wind-tunnel test fills a gap in the manned-space programs, and will be used to establish the acoustic environment for vibro-acoustic qualification testing of the MPCV.

Panda, Jayanta↗

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Modeling Technique and Representation of Failure in the Analysis of Triaxial Braided Carbon Fiber Composites

Quasi-static tests have been performed on triaxially braided carbon fiber composite materials with large unit cell sizes. The effects of different fibers and matrix materials on the failure mode were investigated. Simulations of the tests have been performed using the transient dynamic finite element code, LS-DYNA. However, the wide range of failure modes observed for the triaxial braided carbon fiber composites during tests could not be simulated using composite material models currently available within LS-DYNA. A macroscopic approach has been developed that provides better simulation of the material response in these materials. This approach uses full-field optical measurement techniques to measure local failures during quasi-static testing. Information from these experiments is then used along with the current material models available in LS-DYNA to simulate the influence of the braided architecture on the failure process. This method uses two-dimensional shell elements with integration points through the thickness of the elements to represent the different layers of braid along with a new analytical method for the import of material stiffness and failure data directly. The present method is being used to examine the effect of material properties on the failure process. The experimental approaches used to obtain the required data will be described, and preliminary results of the numerical analysis will be presented.

Littell, Justin D.↗

A Modeling Technique and Representation of Failure in the Analysis of Triaxial Braided Carbon Fiber Composites

Quasi-static tests have been performed on triaxially braided carbon fiber composite materials with large unit cell sizes. The effects of different fibers and matrix materials on the failure mode were investigated. Simulations of the tests have been performed using the transient dynamic finite element code, LS-DYNA. However, the wide range of failure modes observed for the triaxial braided carbon fiber composites during tests could not be simulated using composite material models currently available within LS-DYNA. A macroscopic approach has been developed that provides better simulation of the material response in these materials. This approach uses full-field optical measurement techniques to measure local failures during quasi-static testing. Information from these experiments is then used along with the current material models available in LS-DYNA to simulate the influence of the braided architecture on the failure process. This method uses two-dimensional shell elements with integration points through the thickness of the elements to represent the different layers of braid along with a new analytical method for the import of material stiffness and failure data directly. The present method is being used to examine the effect of material properties on the failure process. The experimental approaches used to obtain the required data will be described, and preliminary results of the numerical analysis will be presented.

Littell, Justin D.↗

Bridging microscopic dynamics and rheology in the yielding of charged colloidal suspensions

The yielding of soft materials is critical to many natural and industrial processes, yet experimental insights into microscopic aspects of yielding are limited. This study combines angle X-ray scattering, X-ray photon correlation spectroscopy, and in situ rheology (Rheo-SAXS-XPCS) with fast lubrication dynamics simulations to examine how interparticle interactions influence yielding in charged colloidal suspensions. By tuning attraction through salt addition, we compare repulsive and attractive systems under deformation. Repulsive suspensions yield uniformly with Andrade-like creep and minimal structural change. In contrast, attractive suspensions show complex behaviors, including shear banding, delayed yielding, and resolidification, governed by transient dynamics at shear band interfaces. These results directly link microscopic particle dynamics to macroscopic flow and demonstrate how interaction potentials control rheological behavior. This work offers a framework for designing soft materials with tailored properties for applications in coatings, food processing, drug delivery, and other technologies requiring precise mechanical control.

Molecular Dynamics Simulation↗

Three-dimensional numerical simulation of interplanetary magnetic field changes at 1 AU as a consequence of simulated solar flares

In order to study the relationship between different forms of activity and transient variations of the north-south component, B(sub z), of the Interplanetary Magnetic Field (IMF), at 1 AU, a three dimensional numerical simulation code was employed to study several aspects of this problem. One form of solar activity, the flare, is simulated by using a pressure pulse at different locations near the solar surface and observing the simulated IMF evolution of B(sub theta)(-B(sub z)) at 1 AU. It was found, for a given pressure pulse, that the orientation of the corresponding transient variation of B(sub z) has a strong relationship with the location of the pressure pulse and the initial condition of IMF.

Wu, S. T.↗

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

ABSTRACT Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell’s equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly rotating magnetized neutron stars, taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. Rs ∼ 0.5 [R*], we show that the new terms lead to a per cent order correction in the magnetic field orientation and strength compared to the case, with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable, when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

Torres, R. (ORCID:0000000291820228)↗

Determination of Thermal State of Charge in Solar Heat Receivers

The research project at Howard University seeks to develop analytical and numerical capabilities to study heat transfer and fluid flow characteristics, and the prediction of the performance of solar heat receivers for space applications. Specifically, the study seeks to elucidate the effects of internal and external thermal radiation, geometrical and applicable dimensionless parameters on the overall heat transfer in space solar heat receivers. Over the last year, a procedure for the characterization of the state-of-charge (SOC) in solar heat receivers for space applications has been developed. By identifying the various factors that affect the SOC, a dimensional analysis is performed resulting in a number of dimensionless groups of parameters. Although not accomplished during the first phase of the research, data generated from a thermal simulation program can be used to determine values of the dimensionless parameters and the state-of-charge and thereby obtain a correlation for the SOC. The simulation program selected for the purpose is HOTTube, a thermal numerical computer code based on a transient time-explicit, axisymmetric model of the total solar heat receiver. Simulation results obtained with the computer program are presented the minimum and maximum insolation orbits. In the absence of any validation of the code with experimental data, results from HOTTube appear reasonable qualitatively in representing the physical situations modeled.

Glakpe, E. K.↗

The Sensitivity of African Easterly Waves to Eastern Tropical Atlantic Sea-Surface Temperatures

The results of two regional atmospheric model simulations are compared to assess the influence of the eastern tropical Atlantic sea-surface temperature maximum on local precipitation, transient easterly waves and the West African summer monsoon. Both model simulations were initialized with reanalysis 2 data (US National Center for Environmental Prediction and Department of Energy) on 15 May 2006 and extended through 6 October 2006, forced by synchronous reanalysis 2 lateral boundary conditions introduced four times daily. One simulation uses 2006 reanalysis 2 sea-surface temperatures, also updated four times daily, while the second simulation considers ocean forcing absent the sea-surface temperature maximum, achieved here by subtracting 3 K at every ocean grid point between 0 and 15 N during the entire simulation. The simulation with 2006 sea-surface temperature forcing produces a realistic distribution of June-September mean precipitation and realistic westward propagating swaths of maximum rainfall, based on validation against Tropical Rainfall Measuring Mission (TRMM) estimates. The simulation without the sea-surface temperature maximum produces only 57% of the control June-September total precipitation over the eastern tropical Atlantic and about 83% of the Sahel precipitation. The simulation with warmer ocean temperatures generates generally stronger circulation, which in turn enhances precipitation by increasing moisture convergence. Some local precipitation enhancement is also attributed to lower vertical thermal stability above the warm water. The study shows that the eastern tropical Atlantic sea-surface temperature maximum enhances the strength of transient easterly waves and broadens the spatial extent of associated precipitation. However, large-scale circulation and its interaction with the African continent, and not sea-surface temperatures, control the timing and trajectories of the waves.

Druyan, Leonard M.↗

An Investigation Regarding the Effects of Volume Dynamics on Gas Turbine Engine Performance and Operability

This work explores the impact of adding volume dynamics to an electrified turbofan engine model. The zero-dimensional single-input single-output volume elements capture the physics of fluid compressibility of flow in the engine core gas path turbomachinery during transient operation with greater fidelity. Some background regarding the motivation for this effort is provided. A model of the volume element is developed from first principles. The baseline engine model is described and an explanation is given of the modifications made to incorporate the volume elements into the new model. Details of the simulated inputs, designed to elicit exaggerated dynamic responses and to facilitate comparison between the baseline and modified engine models, are explicated. The results of these simulations are presented and discussed, demonstrating that the transient performance and operability of the modified engine model, measured by parameters such as shaft speeds, thrust, and high-pressure compressor (HPC) stall margin, are impacted by the additional dynamics. Commentary is provided regarding potential explanations for the differences in behavior exhibited by the two models. Finally, suggestions are offered for future work that might improve the volume models developed in this effort and leverage the improved fidelity enabled by the modeling of turbomachinery volume dynamics.

Volume dynamics↗