Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Modeling”

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,603 records · Page 89

The encasulation effect on thermal performance of micro-encapsulated phase change materials during energy absorption

Understanding the performance of the thermal storage materials is required to successfully meet the needs of next generation concentrating solar power projects. This study focuses on the thermal characteristics of the micro-encapsulated phase change material (MEPCM) with metallic core under high temperature. A numerical model was developed and verified to study the transient heat absorption process of MEPCMs made of tin, aluminum and copper core. Convective heat transfer rate, total energy absorption and volumetric energy density were presented. The effect of encapsulation on theses variables was analyzed under various boundary and initial conditions. Its impact on characteristic time variables was discussed. The method presented can be adopted to analyzing both manufacturing and application processes for MEPCMs of various compositions, during which melting and solidification happen. Given certain operating conditions, the results from the model can be referenced to obtain optimal thermal performance of these materials.

Benner, Jingru↗

MRAC Revisited: Guaranteed Performance with Reference Model Modification

This paper presents modification of the conventional model reference adaptive control (MRAC) architecture in order to achieve guaranteed transient performance both in the output and input signals of an uncertain system. The proposed modification is based on the tracking error feedback to the reference model. It is shown that approach guarantees tracking of a given command and the ideal control signal (one that would be designed if the system were known) not only asymptotically but also in transient by a proper selection of the error feedback gain. The method prevents generation of high frequency oscillations that are unavoidable in conventional MRAC systems for large adaptation rates. The provided design guideline makes it possible to track a reference command of any magnitude form any initial position without re-tuning. The benefits of the method are demonstrated in simulations.

Stepanyan, Vahram↗

On multidomain multiscale modeling and simulation of a novel partial pressure and temperature swing adsorptive reactor ($\mathrm{PPTSAR}$) with application to the water gas shift reaction

Here in this work, a dynamic, multidomain, multiscale model of a novel Partial Pressure and Temperature Swing Adsorptive Reactor (PPTSAR) Process is developed and simulated for a Water Gas Shift Reaction (WGSR) application. The proposed model couples the reactor wall, catalyst/adsorbent pellet, and fluid mixture domains at their respective scales to comprehensively model mass, momentum, and energy transport phenomena occurring within each domain. The resulting model elucidates the PPTSAR’s start-up, transient, and long-term, periodic, dynamic behavior, during both adsorption/reaction and desorption operating modes. In the WGSR application considered, the use of two PPTSARs is shown to allow continuous operation that delivers CO conversion greater than 95% with simultaneous complete CO 2 capture.

42 ENGINEERING↗

Secondary Volcanically-Induced Lunar Atmosphere and Lunar Volatiles: 3-D Modeling and Analysis

The fact that the Moon could have a transient secondary atmosphere due to volcanic outgassing has been known for some time, though typically such an atmosphere was believed to be extremely thin (~10-8 bar) [1]. But recent research by Needham and Kring (NK) [2] suggests that during the peak of volcanic activity ~3.5 Ga such a volcanically-outgassed atmosphere could reach ~10-2 bar of surface pressure. In similar research Wilson et al. [3] proposed a more conservative estimate, arguing that the thickness of such an atmosphere would depend on the intervals between major eruptions and may not exceed microbar densities. In either case a collisional atmosphere could be present, which would control transport of outgassed volatiles (such as H2O) and their deposition in polar regions, where they could be preserved until modern day frozen in permanently shadowed regions (PSR) or buried beneath the regolith. Here we study such a hypothetical atmosphere to investigate its stability, meteorological properties and the effect on transport of volatiles. We use the ROCKE-3D planetary 3-D General Circulation Model (GCM)[4]. The insolation and orbital parameters were set to conditions 3.5 Ga. The atmospheric composition, based on the list of outgassed species presented by NK in combination with our estimates for atmospheric escape, condensation and the results from our 1-D chemistry model, was chosen to be either CO-dominated or CO2-dominated (depending on atmospheric temperature). In this study we restricted ourselves to relatively "thick" lunar atmospheres of 1-10 mb, though we believe that our results will scale to thinner atmospheres as well. We present the results for ground and atmospheric temperature for modeled atmospheres over a wide parameter space. In particular we consider different atmospheric compositions (CO or CO2 dominated), a set of atmospheric pressures from 1 mb to 10 mb and a set of obliquities from 0o to 40o. We also present an experiment of a single major eruption [5] and show that in just 3 years ~80% of the outgassed water is deposited in polar regions. This demonstrates the efficiency of such an atmosphere in delivering volatiles. We argue that a secondary lunar atmosphere could play a significant role in forming volatile deposits currently observed in the polar regions of the Moon.

Lunar atmospheres↗

An investigation into the causes of stratospheric ozone loss in the southern Australasian region

Measurements of total ozone at Macquarie Island (55 deg S, 159 deg E) reveal statistically significant reductions of approximately twelve percent during July to September when comparing the mean levels for 1987-90 with those in the seventies. In order to investigate the possibility that these ozone changes may not be a result of dynamic variability of the stratosphere, a simple linear model of ozone was created from statistical analysis of tropopause height and isentropic transient eddy heat flux, which were assumed representative of the dominant dynamic influences. Comparison of measured and modeled ozone indicates that the recent downward trend in ozone at Macquarie Island is not related to stratospheric dynamic variability and therefore suggests another mechanism, possibly changes in photochemical destruction of ozone.

Lehmann, P.↗

Pulsed Plasma Acceleration Modeling in Detonation and Deflagration Modes

Gas-fed electromagnetic pulsed plasma accelerators discharge electrical energy into a gas to ionize and electromagnetically accelerate propellant in the domain. Efforts to model pulsed accelerators have either assumed that the discharge is short and completely transient, accelerating the gas by entraining it in a moving current sheet, or that the discharge is relatively long, establishing a stable quasi-steady current distribution pattern that accelerates a plasma flowing through it. There have been pulsed plasma accelerator tests that appear to fit somewhere between these two bounds, exhibiting some properties that are associated with the purely transient devices while also showing others that are associated with quasi-steady-state plasma acceleration. A model is presented based upon the premise that all pulsed plasma accelerators first form an accelerating current sheet (detonation mode accelerator). Depending upon the pulse length and the gas conditions in the dis- charge channel, the plasma sheet may reach the end of the accelerator before the discharge has completed a full half-cycle, wherein the proposed model transitions to a quasi-steady description of the acceleration process (deflagration mode accelerator). A review of the entire model is presented, highlighting improvements and upgrades implemented to aid in the stability of the numerical scheme used to model the gas flow in the channel and an improved treatment of current sheet mass shedding, which adds gas to the wake of the sheet. The assumptions employed to model the transition from detonation to deflagration mode are presented and used to generate solutions to the governing equations. The modeling of the deflagration-mode under the present assumptions results in very low deflagration impulse bits relative to those obtained at the end of the detonation mode, implying that the present assumptions and modeling of the deflagration mode may not represent a fruitful approach.

Polzin, Kurt A.↗

Dynamic Model Agnostic Reliability Evaluation of Machine-Learning Models Integrated in Instrumentation & Control Systems

In recent years, the field of data-driven neural network-based machine learning (ML) algorithms has grown significantly and spurred research in its applicability to instrumentation and control systems. While they are promising in operational contexts, the trustworthiness of such algorithms is not adequately assessed. Failures of ML-integrated systems are poorly understood; the lack of comprehensive risk modeling can degrade the trustworthiness of these systems. In recent reports by the National Institute for Standards and Technology, trustworthiness in ML is a critical barrier to adoption and will play a vital role in intelligent systems' safe and accountable operation. Thus, in this work, we demonstrate a real-time model-agnostic method to evaluate the relative reliability of ML predictions by incorporating out-of-distribution detection on the training dataset. It is well documented that ML algorithms excel at interpolation (or near-interpolation) tasks but significantly degrade at extrapolation. This occurs when new samples are "far" from training samples. The method, referred to as the Laplacian distributed decay for reliability (LADDR), determines the difference between the operational and training datasets, which is used to calculate a prediction's relative reliability. LADDR is demonstrated on a feedforward neural network-based model used to predict safety significant factors during different loss-of-flow transients. LADDR is intended as a "data supervisor" and determines the appropriateness of well-trained ML models in the context of operational conditions. Ultimately, LADDR illustrates how training data can be used as evidence to support the trustworthiness of ML predictions when utilized for conventional interpolation tasks.

97 MATHEMATICS AND COMPUTING↗

Impeller leakage flow modeling for mechanical vibration control

HPOTP and HPFTP vibration test results have exhibited transient and steady characteristics which may be due to impeller leakage path (ILP) related forces. For example, an axial shift in the rotor could suddenly change the ILP clearances and lengths yielding dynamic coefficient and subsequent vibration changes. ILP models are more complicated than conventional-single component-annular seal models due to their radial flow component (coriolis and centrifugal acceleration), complex geometry (axial/radial clearance coupling), internal boundary (transition) flow conditions between mechanical components along the ILP and longer length, requiring moment as well as force coefficients. Flow coupling between mechanical components results from mass and energy conservation applied at their interfaces. Typical components along the ILP include an inlet seal, curved shroud, and an exit seal, which may be a stepped labyrinth type. Von Pragenau (MSFC) has modeled labyrinth seals as a series of plain annular seals for leakage and dynamic coefficient prediction. These multi-tooth components increase the total number of 'flow coupled' components in the ILP. Childs developed an analysis for an ILP consisting of a single, constant clearance shroud with an exit seal represented by a lumped flow-loss coefficient. This same geometry was later extended to include compressible flow. The objective of the current work is to: supply ILP leakage-force impedance-dynamic coefficient modeling software to MSFC engineers, base on incompressible/compressible bulk flow theory; design the software to model a generic geometry ILP described by a series of components lying along an arbitrarily directed path; validate the software by comparison to available test data, CFD and bulk models; and develop a hybrid CFD-bulk flow model of an ILP to improve modeling accuracy within practical run time constraints.

Palazzolo, Alan B.↗

Gas hydrate deposit formation in transient flowloop tests and mitigation with a surface treatment

Thermodynamic and thermal mitigation strategies are traditionally utilized to prevent gas hydrate formation but may be rendered ineffective for certain unplanned, extended shut-in scenarios. Unfortunately, the current knowledge regarding transient gas hydrate formation has been extremely limited to only a few studies and unable to address this complex issue adequately. This work, conducted in a lab-scale flowloop, aims to add to the understanding of this area by delving into the implications of transient conditions on gas hydrate formation mechanisms, while also exploring a possible mitigation strategy of a surface treatment. Specifically, the effects of phase separation on gas hydrate formation locations, transportability, and plugging are presented for an oil-dominated system. An omniphobic surface treatment is shown to affect initial gas hydrate deposition and its wall adherence strength, as well as the overall plugging outcome. Overall, this work combines observational information with process data measurements to provide a better understanding of phase separated gas hydrate systems and passive deposition mitigation strategies. This can further aid advanced gas hydrate modeling efforts, which currently do not adequately address transient flow effects on phase separation/dispersion. Here, it offers representative gas hydrate testing for coatings, which have exclusively been limited to non-flowing and/or bench-scale analysis, to give confidence in their potential application as advanced deposition prevention tools.

03 NATURAL GAS↗

Contact-impact analysis of deployable space systems

This paper presents simulation results of a flexible, deployable space system with multiple contact-impact events, using three different modeling approaches. In the first two approaches, the equations of motion of the system are written in terms of generalized coordinates of each component, with modal coordinates used to account for flexibility in the components. Two different methods are used to model contact-impact events in the systems. In the momentum balance method, as contact-impact occurs, momentum equations of the systems are solved for the jump discontinuity in velocity of generalized coordinates. Integration of the system is restarted with the new initial conditions from solutions of momentum equations. In the logical function method, a virtual spring connecting the two contact points is assumed. The virtual spring is active only when impact occurs. No special treatment is needed in the integration algorithm using this method. In the third approach, instead of using modal coordinates, a finite element based formulation for the transient dynamics of mechanical system is used to model flexibility of the components. The logical function method is used to model contact-impact in the third approach.

Wu, Shih-Chin↗

Effect of Hydrogen on Steady-State and Transient Combustion Instability Characteristics

This paper examines the effects of steady-state and transient hydrogen enrichment on thermoacoustic instability in a model gas turbine combustor. Combustion instability, a feedback loop between flame heat release rate oscillations and combustor acoustics, is characterized in a swirl-stabilized flame operated at a range of hydrogen–natural gas fuel blends and heat rates. Measurements of combustor chamber pressure fluctuations and CH* chemiluminescence imaging are used to characterize instability at a range of operating conditions. Steady-state tests show that both mixture heat rate and hydrogen content affect system stability. At a given heat rate, higher levels of hydrogen result in unstable combustion. As heat rate increases, instability occurs at lower concentrations of hydrogen in the fuel. Transient operation was tested in two directions—instability onset and decay—and two hydrogen-addition times—a short time of 1 ms and a longer time of 4 s. Results show that instability onset processes, through the transient addition of hydrogen, are highly repeatable regardless of the timescale of hydrogen addition. Certain instability decay processes are less repeatable, resulting in cases that do not fully transition from unstable to stable combustion despite similar changes in hydrogen fuel flow rate. Flame behavior before, during, and after the transient is characterized using high-speed CH* chemiluminescence imaging. Analysis of the high-speed images shows changes in flame stabilization and dynamics during the onset and decay processes. Finally, the results of this study can have implications for systems that experience variations in fuel composition, particularly in light of growing interest in hydrogen as a renewable fuel.

42 ENGINEERING↗

Numerical simulation of small-scale thermal convection in the atmosphere

A Boussinesq system is integrated numerically in three dimensions and time in a study of nonhydrostatic convection in the atmosphere. Simulation of cloud convection is achieved by the inclusion of parametrized effects of latent heat and small-scale turbulence. The results are compared with the cell structure observed in Rayleigh-Benard laboratory conversion experiments in air. At a Rayleigh number of 4000, the numerical model adequately simulates the experimentally observed evolution, including some prominent transients of a flow from a randomly perturbed initial conductive state into the final state of steady large-amplitude two-dimensional rolls. At Rayleigh number 9000, the model reproduces the experimentally observed unsteady equilibrium of vertically coherent oscillatory waves superimposed on rolls.

Somerville, R. C. J.↗

Simulation of the Space Station strut-out condition

A method is presented for reanalyzing a truss structure when one of the truss elements (struts) has failed. The method uses a modal model of the nominal structure coupled with a residual flexibility term to predict the effect of the failed strut without resolving the finite element model. By implementing the method as part of the transient simulation, it is feasible to consider a large number of potential strut failures with a minimum amount of extra effort. Preliminary application of the method to the Space Station indicates excellent agreement with results based on modifying and resolving the finite element model.

Blelloch, Paul A.↗

Dynamic simulation of compressor and gas turbine performance

Dynamic simulation of compressor and gas turbine performance is useful in a variety of design, analysis, and test applications. Mathematical models employing estimated or experimentally-derived component characteristics have been widely used, especially for near-design-point studies of acceleration and deceleration transients and control strategies. With recent progress in simulation methods and computational power, it has become possible to build models with detail at the stage level, and with fundamental fluid mechanics input. Such models are more useful because of the ability to study the details of inter-component flow property behavior during machine transients, even including stall and surge. Dynamic performance simulation models are discussed with emphasis on the fundamental principles of the models and the methods used to represent component and stage flow characteristics. Results of several simulations of the dynamic behavior of multistage compressors are shown with comparisons to experimental data. Possibilities for advanced computational techniques for near-real-time simulations of compressors and gas turbines are reviewed.

Walter F O'Brien↗

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗

Enhancement of the no-vent fill process

This paper presents an analytical and experimental evaluation of an enhanced techniques for no-vent fill. The method entails injecting liquid through the top of the receiver vessel, thereby increasing surface area and agitation of the ullage/liquid interface. Both of these factors promote condensation induced ullage collapse, and reduce compressive impedance to the incoming liquid. The enhanced process was analyzed by modifying the surface area algorithm of an existing tank thermodynamic code to model a downward-pointing, conical jet impringing on a steadily rising liquid surface. Transient pressure and temperature measurements from several tests with Freon-114 were input into the revised model to calculate condensation rate as a function of fill level. By expressing these rates in dimensionless form (i.e., in terms of Stanton number and Prandtl number), an empirical correlation similar to the submerged jet model of Brown and Sonin (1989) was derived. This provided a basis for developing an expression which relates top fill to bottom fill performance.

Vaughan, David A.↗

A Martian global groundwater model

A global groundwater flow model was constructed for Mars to study hydrologic response under a variety of scenarios, improving and extending earlier simple cross sectional models. The model is capable of treating both steady state and transient flow as well as permeability that is anisotropic in the horizontal dimensions. A single near surface confining layer may be included (representing in these simulations a coherent permafrost layer). Furthermore, in unconfined flow, locations of complete saturation and seepage are determined. The flow model assumes that groundwater gradients are sufficiently low that DuPuit conditions are satisfied and the flow component perpendicular to the ground surface is negligible. The flow equations were solved using a finite difference method employing 10 deg spacing of latitude and longitude.

Howard, Alan D.↗

Characterization of an Optically Accessible RP/GOx Rotating Detonation Engine

The accurate characterization of the dynamic response of liquid jets to transient detonation wave passage is crucial for optimizing and modeling liquid-fueled rotating detonation combustors. To address this need, we have developed an optically accessible rocket RDE (RP-RDE), which utilizes liquid RP-1 and gaseous oxygen as propellants and operates at elevated pressures. The experimental setup features optical access from the oxidizer injection plenum to the combustor exit, allowing for detailed characterization of injection, mixing, and combustion processes and their impact on the detonation wave structure. To design and develop the liquid injection scheme for the RP-RDE test article, we employed a combined experimental and numerical approach that included cold-flow simulations. Our aim was to optimize the liquid injection process and ensure that the resulting fuel-oxidizer mixture is homogeneous and stable. To this end, high-speed optical imaging will be employed to provide qualitative visualization of the liquid and gaseous phases of the propellant as well as the local heat release during detonation wave propagation through the non-premixed medium. This work will provide critical insights into the behavior of the transient detonation waves in more relevant rocket conditions, which is essential for the optimization and modeling of rotating detonation combustors.

Propulsion↗