Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 91 records · Page 5

Impeller shroud to casing leakage flow simulations in the Space Shuttle Main Engine high pressure fuel pump

Quasi-three-dimensional Navier-Stokes calculations were carried out for the Space Shuttle Main Engine high-pressure fuel pump to simulate the impeller shroud to casing leakage flow. This flow geometry was modeled as an axisymmetric cavity flow with a stationary surface representing the casing, and a rotating surface denoting the impeller. A 63 x 81-node mesh provided sufficient resolution in the regions of greatest flow variations and reduced the effects of numerical diffusion. The turbulence field was closed with the high Reynolds number form of the k-epsilon model supplemented with wall functions in the vicinity of the walls. Finally, a parametric study quantified the effects of through mass flow changes on this leakage flow.

Sindir, Munir M.↗

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Water Content↗

Physical and Electronic Isolation of Carbon Nanotube Conductors

Multi-walled nanotubes are proposed as a method to electrically and physically isolate nanoscale conductors from their surroundings. We use tight binding (TB) and density functional theory (DFT) to simulate the effects of an external electric field on multi-wall nanotubes. Two categories of multi-wall nanotube are investigated, those with metallic and semiconducting outer shells. In the metallic case, simulations show that the outer wall effectively screens the inner core from an applied electric field. This offers the ability to reduce crosstalk between nanotube conductors. A semiconducting outer shell is found not to perturb an electric field incident on the inner core, thereby providing physical isolation while allowing the tube to remain electrically coupled to its surroundings.

OKeeffe, James↗

Impact of the Transformer Magnetizing Inductance on the Performance of the Dual-Active Bridge Converter

The magnetizing inductance of the high-frequency transformer in the dual-active bridge (DAB) converter is usually assumed to be much larger than its leakage inductance, thereby the magnetizing current can be negligible. However, ignoring the magnetizing current will lead to an error in the DAB converter design process. This paper analyzes the operation and performance of the DAB converter for a general case in which the transformer magnetizing inductance is taken into consideration. The analysis shows that the magnetizing inductance impacts the soft-switching region, the power transfer capability, and the leakage inductor current. A simulation case study is provided to verify the analysis. A DAB converter prototype is built to experimentally validate the analysis as well.

Yin, Shiyuan↗

Developing parallel GeoFEST(P) using the PYRAMID AMR library

The PYRAMID parallel unstructured adaptive mesh refinement (AMR) library has been coupled with the GeoFEST geophysical finite element simulation tool to support parallel active tectonics simulations. Specifically, we have demonstrated modeling of coseismic and postseismic surface displacement due to a simulated Earthquake for the Landers system of interacting faults in Southern California. The new software demonstrated a 25-times resolution improvement and a 4-times reduction in time to solution over the sequential baseline milestone case. Simulations on workstations using a few tens of thousands of stress displacement finite elements can now be expanded to multiple millions of elements with greater than 98% scaled efficiency on various parallel platforms over many hundreds of processors. Our most recent work has demonstrated that we can dynamically adapt the computational grid as stress grows on a fault. In this paper, we will describe the major issues and challenges associated with coupling these two programs to create GeoFEST(P). Performance and visualization results will also be described.

adaptive mesh refinement (AMR)↗

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Applications of active adaptive noise control to jet engines

During phase 2 research on the application of active noise control to jet engines, the development of multiple-input/multiple-output (MIMO) active adaptive noise control algorithms and acoustic/controls models for turbofan engines were considered. Specific goals for this research phase included: (1) implementation of a MIMO adaptive minimum variance active noise controller; and (2) turbofan engine model development. A minimum variance control law for adaptive active noise control has been developed, simulated, and implemented for single-input/single-output (SISO) systems. Since acoustic systems tend to be distributed, multiple sensors, and actuators are more appropriate. As such, the SISO minimum variance controller was extended to the MIMO case. Simulation and experimental results are presented. A state-space model of a simplified gas turbine engine is developed using the bond graph technique. The model retains important system behavior, yet is of low enough order to be useful for controller design. Expansion of the model to include multiple stages and spools is also discussed.

Shoureshi, Rahmat↗

Numerical Simulations of A Conceptual MSR-EES Shoulder Recession

The present study demonstrates our in-house material response solver, Icarus's capability to simulate the ablative conditions, including pyrolysis effects due to the interaction between hypersonic boundary layers and the thermal protection system (TPS). A conceptual aeroshell shoulder design, which undergoes mission-relevant flow and material conditions, is selected for the demonstration purpose. LAURA, a structured flow solver, is used to solve flow around the shoulder at several trajectory points of a flight path. The aerothermal dataset obtained from LAURA is used to enforce boundary conditions on the aeroshell wall to simulate ablation processes. The two-layered material system is stacked with HEEET (TPS) and Aluminum (actual material). A set of base parameters, such as the angle of material orientation with respect to the flow direction, convective heat-transfer co-efficient, and aerothermal boundary conditions on the rear end of the HEEET material layer, is selected to conduct a base-case simulation. A small amount of recession was observed, indicating that the design should be fine to go through the selected flow and material conditions. Furthermore, the parameters mentioned above are individually altered and are observed to affect the shoulder design's recession compared to the base case. Verification of our Icarus setup was also carried out using a one-dimensional grid and FIAT, a one-dimensional material-response solver, to build credibility for our results. As the current work is an uncoupled fluid-material-response simulation, a local sharp mesh deformation in the recessed surface near the shoulder corner is observed.

Prakash Shrestha↗

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulation studies for surfaces and materials strength

Computer simulation studies were conducted to investigate materials properties. All the calculations were carried out using atomic level simulation techniques which are based on semiempirical or model functions. During this project, functions with varying degrees of complexity were derived and employed in simulations. Investigations include covalently bonded materials as well as systems involving metal atoms. For small clusters calculated results provided information on various energy- and structure-related properties along with vibrational characteristics. Also, energy barriers for configurational transitions were calculated for selected cases. Simulation calculations for surfaces produced new results in areas related to surface energies, reconstructions and relaxations, surface defects, surface stresses as well as adsorption and nucleation processes. Simulations were also performed on bulk materials. Calculated results provided an atomic level understanding on energetics and structures of point defects, crystal stability, elastic properties, and materials strength for various systems. Calculations involving polymeric materials include studies of polar polymer melts and polymer/solid interfaces. Simulations employing specially developed codes provided significant information about energetics and conformational characteristics of different polymeric chains.

Halicioglu, Timur↗

Tracking performance and cycle slipping in the all-digital symbol synchronizer loop of the block 5 receiver

Computer simulated noise performance of the symbol synchronizer loop (SSL) in the Block 5 receiver is compared with the theoretical noise performance. Good agreement is seen at the higher loop SNR's (SNR(sub L)'s), with gradual degradation as the SNR(sub L) is decreased. For the different cases simulated, cycle slipping is observed (within the simulation time of 10(exp 4) seconds) at SNR(sub L)'s below different thresholds, ranging from 6 to 8.5 dB, comparable to that of a classical phase-locked loop. An important point, however, is that to achieve the desired loop SNR above the seemingly low threshold to avoid cycle slipping, a large data-to-loop-noise power ratio, P(sub D)/(N(sub 0)B(sub L)), is necessary (at least 13 dB larger than the desired SNR(sub L) in the optimum case and larger otherwise). This is due to the large squaring loss (greater than or equal to 13 dB) inherent in the SSL. For the special case of symbol rates approximately equaling the loop update rate, a more accurate equivalent model accounting for an extra loop update period delay (characteristic of the SSL phase detector design) is derived. This model results in a more accurate estimation of the noise-equivalent bandwidth of the loop.

Aung, M.↗

Forward Flight Rotor Performance at Martian Atmospheric Densities and Sensitivity to Low Reynolds Numbers

Much effort has been made to enhance exploration on Mars. In addition to a rover and Mars-orbiting satellites, a Mars Helicopter Technology Demonstrator was proposed by the NASA Jet Propulsion Laboratory (JPL) to augment planetary research for the Mars 2020 Mission. Understanding rotor performance is vital for operations at Martian atmospheric conditions. The work presented is a study investigating rotor performance at Martian atmospheric conditions. Forward flight rotor tests were conducted in the Planetary Aeolian Laboratory (PAL) at NASA Ames Research Center, which has the capability to evacuate the air in the chamber to reach Martian atmospheric densities. A 1-meter-diameter rotor, roughly approximating the Mars Helicopter Technology Demonstrator, was tested at multiple atmospheric densities, including that of Mars. Rotor rotational speed, thrust, torque, power, and airspeed measurements were collected during the test. These results were then correlated with simulated cases using a mid-fidelity computational fluid dynamics software, Rotorcraft CFD (RotCFD). C81Generator (C81Gen) was used to generate airfoil aerodynamic coefficient for the spanwise locations along the rotor. To observe the differences between the C81Gen flow type modes at low Reynolds number, the simulations at Martian atmospheric densities were run under the “fully turbulent”, and the fully laminar flow type. In addition, Reynolds number effects (within 2x104 to 9x104) on experimental thrust coefficient, power coefficient, and figure of merit were analyzed. Within this chord- based Reynolds number range, CT and FM decreased around 26% and 36%, respectively, while CP remained fairly constant, exhibiting variations of no more than 5.5%. Despite the challenges involved in testing at a large difference of atmospheric ensities between Earth and Mars, repeatable data was obtained in all the measurements at Martian atmospheric conditions.

Perez Perez, Brenda NAtalia↗

Surrogate Model of Flexible Research Platform EnergyPlus Models to Enable Sensitivity Analysis

This letter report describes the surrogate models developed from the EnergyPlus model of Oak Ridge National Laboratory’s Flexible Research Platform. Two data-driven black-box models were developed, and the outputs of the surrogate models were compared with the EnergyPlus model. The two models developed are a multilayer perceptron deep learning model, and a long short-term memory (LSTM) neural network model. The three factors for selecting the black-box models are scalability, computation time, and accuracy. A total of 107 input variables were the dominant variables in determining the outputs of building energy consumptions and thermal comfort. A total of 54 output variables were identified as the prediction targets, including the system- and zone-level outputs. The large set of the simulation cases were generated by integrating sensor errors into an emulator based on EnergyPlus and Python EMS, which includes advanced control sequences from ASHRAE Guideline 36-2018: High-Performance Sequences of Operation. The surrogate models were developed based on a set of large-scale simulation runs (i.e., 4,000 runs) on a cloud platform. The comparison analysis shows that the two black-box models had good accuracy for predicting new outputs for sensitivity analysis using the root mean square error metric. As a next step, the developed surrogate models will be used to perform sensitivity analysis for different sensor impacts (e.g., sensor types, sensor locations).

42 ENGINEERING↗

Flexible Service Contracting for Risk Management within Integrated Transmission and Distribution Systems

The general objective of our project has been to investigate the ability of Independent Distribution System Operators (IDSOs), functioning as linkage agents for Integrated Transmission and Distribution (ITD) systems, to facilitate the flexible availability and usage of reserve in support of ITD system operations. This general objective is in accordance with Order 2222 of the U.S. Federal Energy Regulatory Commission (FERC), titled “Participation of Distributed Energy Resource Aggregations in Markets Operated by Regional Transmission Organizations and Independent System Operators”. The Final Rule for FERC Order 2222 was issued on September 17, 2020. The primary contribution of our project is that we have formulated an innovative energy management Transactive Energy Design (TES) approach for ITD systems that provides promising support for our project objective. Specifically, we have developed a new type of IDSO-managed TES design for distribution system operations, as well as new types of swing contracts permitting IDSOs to participate in transmission system operations as providers of reserve harnessed from distribution system resources in return for appropriate compensation. Together, these design elements constitute a scalable market-based approach facilitating efficient reserve procurement for ITD system operations from a fuller range of power resources. The efficacy of our approach has been demonstrated by means of detailed conceptual analyses as well as test case simulations. To implement the latter, we have developed the ITD TES Platform V2.0, a computational platform that permits the modeling and software implementation of ITD system operations over time. The Electric Reliability Council of Texas (ERCOT) energy region has been used as the empirical anchor for the development of this platform. Our conceptual and test-case work has been reported in a Wiley/IEEE Press book, two refereed book chapters, and seven refereed journal articles. The key components comprising the ITD TES Platform V2.0 have been released as documented open-source software at online GitHub repositories.

24 POWER TRANSMISSION AND DISTRIBUTION↗

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance Characterization and Simulation of Amine-Based Vacuum Swing Adsorption Units for Spacesuit Carbon Dioxide and Humidity Control

Controlling carbon dioxide (CO2) and water (H2O) concentrations in the vapor phase of a space suit is critical to ensuring an astronauts safety, comfortability, and capability to perform extra-vehicular activity (EVA) tasks. Historically, this has been accomplished using lithium hydroxide (LiOH) and metal oxides (MetOx). Lithium hydroxide is a consumable material and requires priming with water before it becomes effective at removing carbon dioxide. MetOx is regenerable through a power-intensive thermal cycle but is significantly heavier on a volume basis than LiOH. As an alternative, amine-based vacuum swing beds are under aggressive development for EVA applications which control atmospheric concentrations of both CO2 and H2O through a fully-regenerative process. The current concept, referred to as the rapid cycle amine (RCA), has resulted in numerous laboratory prototypes. Performance of these prototypes have been assessed and documented from experimental and theoretical perspectives. To support developmental efforts, a first principles model has also been established for the vacuum swing adsorption technology. The efforts documented herein summarize performance characterization and simulation results for several variable metabolic profiles subjected to the RCA. Furthermore, a variety of control methods are explored including timed swing cycles, instantaneous CO2 feedback control, and time-averaged CO2 feedback control. A variety of off-nominal tests are also explored including high/low suit temperatures, increasingly high humidity cases, and dynamic pressure cases simulating the suit pre-breathe protocol. Consequently, this work builds on efforts previous efforts to fully bound the performance of the rapid cycle amine under a variety of nominal and off-nominal conditions.

Swickrath, Michael J.↗

Dynamic mode decomposition of nonequilibrium electron-phonon dynamics: accelerating the first-principles real-time Boltzmann equation

Abstract Nonequilibrium dynamics governed by electron–phonon ( e -ph) interactions plays a key role in electronic devices and spectroscopies and is central to understanding electronic excitations in materials. The real-time Boltzmann transport equation (rt-BTE) with collision processes computed from first principles can describe the coupled dynamics of electrons and atomic vibrations (phonons). Yet, a bottleneck of these simulations is the calculation of e –ph scattering integrals on dense momentum grids at each time step. Here we show a data-driven approach based on dynamic mode decomposition (DMD) that can accelerate the time propagation of the rt-BTE and identify dominant electronic processes. We apply this approach to two case studies, high-field charge transport and ultrafast excited electron relaxation. In both cases, simulating only a short time window of ~10% of the dynamics suffices to predict the dynamics from initial excitation to steady state using DMD extrapolation. Analysis of the momentum-space modes extracted from DMD sheds light on the microscopic mechanisms governing electron relaxation to a steady state or equilibrium. The combination of accuracy and efficiency makes our DMD-based method a valuable tool for investigating ultrafast dynamics in a wide range of materials.

36 MATERIALS SCIENCE↗