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 1,153 records · Page 64

Wind-US Results for the AIAA 2nd Propulsion Aerodynamics Workshop

This presentation contains Wind-US results presented at the 2nd Propulsion Aerodynamics Workshop. The workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion Systems Integration Technical Committee with the purpose of assessing the accuracy of computational fluid dynamics for air breathing propulsion applications. Attendees included representatives from government, industry, academia, and commercial software companies. Participants were encouraged to explore and discuss all aspects of the simulation process including the effects of mesh type and refinement, solver numerical schemes, and turbulence modeling. The first set of challenge cases involved computing the thrust and discharge coefficients for a 25deg conical nozzle for a range of nozzle pressure ratios between 1.4 and 7.0. Participants were also asked to simulate two cases in which the 25deg conical nozzle was bifurcated by a solid plate, resulting in vortex shedding (NPR=1.6) and shifted plume shock (NPR=4.0). A second set of nozzle cases involved computing the discharge and thrust coefficients for a convergent dual stream nozzle for a range of subsonic nozzle pressure ratios. The workshop committee also compared the plume mixing of these cases across various codes and models. The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 and an offset of 1.34 times the inlet diameter. Boundary layer profiles, wall static pressure, and total pressure at downstream rake locations were examined.

Wind-US↗

Apollo experience report: A use of network simulation techniques in the design of the Apollo lunar surface experiments package support system

A case study of data-communications network modeling and simulation is presented. The applicability of simulation techniques in early system design phases is demonstrated, and the ease with which model parameters can be changed and comprehensive statistics gathered is shown. The discussion of the model design and application also yields an insight into the design and implementation of the Apollo lunar surface experiments package ground-support system.

Gustafson, R. A.↗

Development of an Inertial Sensor-based Methodology for Spacesuited Geology Task Assessments during Simulated Lunar Extravehicular Activities

Lunar surface exploration during Artemis missions will require the specific skill set of geology sampling. Apollo astronauts had extensive training and used specialized tools to collect lunar rocks, core samples, pebbles, sand, and dust. The inflexibility of the pressurized Apollo spacesuits forced sampling to be taken at a standstill posture. However, new exploration spacesuits are expected to incorporate advanced materials and joint bearings, allowing for greater mobility and a wider range of functional postures. Thus, science and exploration during Artemis missions will likely involve a variety of standing, squatting, and kneeling postures. In preparation for future lunar exploration missions, NASA provides geologic training to astronauts and other mission personnel. This professional training with a spacesuit in simulated lunar environments will enhance performance and reduce risk of injury to astronauts on the lunar surface. However, anecdotally, untrained or newly trained people wearing prototype planetary spacesuits have been observed to performing motions differently than a trained geologist would when conducting the same geology sampling tasks. Therefore, a tool for evaluating geology postures at extravehicular activity (EVA) training facilities becomes required. In this paper, we introduce a novel inertial measurement unit (IMU)-based method of geology task assessments in spacesuited conditions during simulated lunar EVAs. As a case study, two subjects (one geologist and one non-geologist) participated and donned the Mark III prototype planetary spacesuit during offloading with the spreader bar gimbal in NASA’s Active Response Gravity Offload System (ARGOS). For automated geology task assessments, the spacesuit was instrumented with three wireless IMUs (APDM Opal, OR, USA): one on the chest and one each on the left and right ankle bearings. Then subjects performed geology tasks using various tools (rake, trench, hammer chisel, scoop, and drive tube) for 45 minutes each. The chest IMU measured the torso tilt angle in the sagittal plane. We used an ensemble learning method with the ankle IMUs to discriminate between standing and kneeling activities. IMU data were processed using custom MATLAB (Mathworks, MA, USA) software. In our case study, the developed method was able to discriminate differences in standing and kneeling activity levels between subjects who were all highly experienced with spacesuited testing. Our preliminary data showed one subject maintained the constant and lower range of the upper body tilt angle while both standing and kneeling, while the other subject showed more variation of the upper body tilt angle and preferred bending the upper body rather than changing from standing to kneeling posture and vice versa. While geology experience may be a factor, these results need further investigation as suit sizing and ARGOS offloading configurations have been proven to have a significant influence on suited ARGOS tasks. Also, more subjects will be needed to complete these tasks for validation. IMU-based geology task assessments can provide useful information for geology training programs. Additionally, our IMU-based posture analysis can provide new insights into how to evaluate spacesuited geology task characteristics of astronauts during simulated lunar EVAs.

Kyoung Jae Kim↗

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗

A Domain-Decomposed Multi-Level Method for Adaptively Refined Cartesian Grids with Embedded Boundaries

The work presents a new method for on-the-fly domain decomposition technique for mapping grids and solution algorithms to parallel machines, and is applicable to both shared-memory and message-passing architectures. It will be demonstrated on the Cray T3E, HP Exemplar, and SGI Origin 2000. Computing time has been secured on all these platforms. The decomposition technique is an outgrowth of techniques used in computational physics for simulations of N-body problems and the event horizons of black holes, and has not been previously used by the CFD community. Since the technique offers on-the-fly partitioning, it offers a substantial increase in flexibility for computing in heterogeneous environments, where the number of available processors may not be known at the time of job submission. In addition, since it is dynamic it permits the job to be repartitioned without global communication in cases where additional processors become available after the simulation has begun, or in cases where dynamic mesh adaptation changes the mesh size during the course of a simulation. The platform for this partitioning strategy is a completely new Cartesian Euler solver tarcreted at parallel machines which may be used in conjunction with Ames' "Cart3D" arbitrary geometry simulation package.

Aftosmis, M. J.↗

Validation of the imaging neutral particle analyzer in nearly MHD quiescent plasmas using injected beam ions on DIII-D

The imaging neutral particle analyzer (INPA) is a scintillator-based diagnostic that provides energy and radially resolved measurements of confined fast ions. To verify operation, understanding, and modeling of the diagnostic performance, distributions of deuterium beam ions in low density, nearly MHD-quiescent plasmas are measured on the DIII-D tokamak using the INPA. Here, the distribution is altered by changing the pitch-angle scattering rate with electron cyclotron heating (ECH). Experimental images are mostly reproduced with simulation of core fast ions charge exchanging with injected beam and edge cold neutrals in classical cases. Image comparisons between the experiment and simulation show agreement with less than 25% discrepancies near the beam injection energy of 50 keV. Other strong passive signals, observed in a localized region on the image, are suggested to come from charge exchange between cold neutrals and fast ions around the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Effect of Boundary Layer Character on Stochastic Rotor Blade Vortex Shedding Noise

This work illustrates the effect of a rotor blade’s boundary layer on the broadband LBL-VS self-noise emitted from an optimum hovering rotor through experimental and multi-fidelity computational studies. Blade surface roughness effects associated with different manufacturing techniques and the effect of adding a boundary layer trip were shown to decrease LBL-VS noise by upwards of 30 dB at the frequency of maximum emission with a slight penalty in aerodynamic performance when compared with smooth rotor blades. Low-fidelity 2-D viscous flow analysis verified the presence of laminar separation bubbles on the rotor blades which are responsible for LBL-VS noise. Three high-fidelity lattice-Boltzmann simulations were conducted with different wall-functions to predict the boundary layer character correspondent to their experimental counterpart and the resultant presence or absence of LBL-VS noise. Excellent aerodynamic and aeroacoustic agreement was seen between the lattice-Boltzmann simulations and the experimental data for the cases with surface roughness and the boundary layer trip, with the broadband noise results from the simulation with fully turbulent wall-functions diverging from the experimental results above 5 kHz. The transitional wall-function simulation which emulated the smooth experimental blades underpredicted thrust by 14% and broadband noise by a minimum of 10 dB with an accurately predicted broadband noise trend.

Christopher S Thurman↗

Neutrino heating in 1D, 2D, and 3D core-collapse supernovae: characterizing the explosion of high-compactness stars

Massive stars can end their lives with a successful supernova explosion (leaving behind a neutron star or, more rarely, a black hole), or a failed explosion that leaves behind a black hole. The density structure of the pre-collapse progenitor star already encodes much of the information regarding the outcome and properties of the explosion. However, the complexity of the collapse and subsequent shock expansion phases prevents drawing a straightforward connection between the pre-collapse and post-explosion properties. In order to derive such a connection several explodability studies have been performed in recent years. However, different studies can predict different explosion outcomes. In this article, we show how compactness, which is related to the average density of the star’s core, has an important role in determining the efficiency of neutrino heating, and therefore the outcome of the explosion. Commonly, high-compactness progenitors are assumed to yield failed explosions, due to their large mass accretion rates, preventing the shock from expanding. We show by analysing ~150 2D flash and F ornax simulations and 20 3D F ornax simulations that this is not the case. Instead, due to the rapid increase of neutrino heating with compactness, high-compactness progenitors lead to successful shock revival. We also show that 1D+ simulations that include v-driven convection using a mixing-length theory approach correctly reproduce this trend. Finally, we compare 1D+ models, which we show can reproduce some aspects of multi-D simulations with reasonable accuracy, with other widely used 1D models in the literature.

neutrinos↗

Simulation of a spatial, servo-hydraulic test facility for space structures

In this paper different control concepts for servo-hydraulic test facilities are derived using exact linearization techniques. Based on different linear and nonlinear models of the test table and the actuator dynamics several nonlinear controllers of different complexity are derived. The closed loop system performance of the controlled servo-hydraulic test facility is tested in various computer simulations using both, standard test signals and large test signals as system inputs. The simulation results turn out, that in case of standard input signals the test facility controller must include a linear test table mechanics model and a nonlinear servo-hydraulic actuator model. Additional simulations demonstrate the robustness of the control concept selected for standard test signals with respect to variations of plant parameters.

Leimbach, K.-D.↗

Global MHD Simulation of Mesoscale Structures at the Magnetospheric Boundary

We have continued our modeling of the solar wind/magnetosphere interaction for steady solar Wind conditions by carrying out a series of 3 D global magneto-hydrodynamics (MHD) simulations using a set of predetermined solar wind parameters as input to the code. With the results from the simulations that we ran last year for two other dipole tilts (0 and 35 deg), we have 36 data sets to study the dynamics of mesoscale structures at the magnetospheric boundary. We have started the analysis of these runs to investigate geometrical properties and the topology of the magnetic and electric fields for the different solar wind regimes and orientations considered. Preliminary results indicate that merging sites are consistent with patterns proposed for antiparallel merging at the dayside magnetopause. Another goal of this investigation is to establish the displacement of the cusp region as a function of the solar wind dynamic pressure, IMF direction and magnetic field dipole tilt. One of the difficulties of that study is to locate precisely the cusp in the simulation results. We have used several case studies to establish a series of criteria in to permit a routine identification of that location. We need now to process the rest of the simulation results using these criteria to obtain statistical results.

Berchem, Jean↗

Mixed Phase Modeling in GlennICE with Application to Engine Icing

A capability for modeling ice crystals and mixed phase icing has been added to GlennICE. Modifications have been made to the particle trajectory algorithm and energy balance to model this behavior. This capability has been added as part of a larger effort to model ice crystal ingestion in aircraft engines. Comparisons have been made to four mixed phase ice accretions performed in the Cox icing tunnel in order to calibrate an ice erosion model. A sample ice ingestion case was performed using the Energy Efficient Engine (E3) model in order to illustrate current capabilities. Engine performance characteristics were supplied using the Numerical Propulsion System Simulation (NPSS) model for this test case.

Wright, William B.↗

Measurement And Simulation of Signal Fluctuations Caused By Propagation Through Trees

We present measured fluctuations of UHF, L-band, and C-band signals that were transmitted from the ground through a forest canopy to an airborne radar. By use of a transmitter in the open, we are able to derive both magnitude and phase fluctuations. We then propose a simple Monte Carlo simulation for amplitude and phase fluctuations cased by propagation of a signal through a forest canopy. The measured RMS magnitude and phase fluctuations are similar to those predicted by the Monte Carlo simulation.

Signal↗

Quantum Computing for High-Energy Physics: State of the Art and Challenges

Quantum computers offer an intriguing path for a paradigmatic change of computing in the natural sciences and beyond, with the potential for achieving a so-called quantum advantage—namely, a significant (in some cases exponential) speedup of numerical simulations. The rapid development of hardware devices with various realizations of qubits enables the execution of small-scale but representative applications on quantum computers. In particular, the high-energy physics community plays a pivotal role in accessing the power of quantum computing, since the field is a driving source for challenging computational problems. This concerns, on the theoretical side, the exploration of models that are very hard or even impossible to address with classical techniques and, on the experimental side, the enormous data challenge of newly emerging experiments, such as the upgrade of the Large Hadron Collider. In this Roadmap paper, led by CERN, DESY, and IBM, we provide the status of high-energy physics quantum computations and give examples of theoretical and experimental target benchmark applications, which can be addressed in the near future. Having in mind hardware with about 100 qubits capable of executing several thousand two-qubit gates, where possible, we also provide resource estimates for the examples given using error-mitigated quantum computing. The ultimate declared goal of this task force is therefore to trigger further research in the high-energy physics community to develop interesting use cases for demonstrations on near-term quantum computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large-Eddy Simulations of Convection Initiation over Heterogeneous, Low Terrain

Abstract Large-eddy simulations are conducted to investigate and physically interpret the impacts of heterogeneous, low terrain on deep-convection initiation (CI). The simulations are based on a case of shallow-to-deep convective transition over the Amazon River basin, and use idealized terrains with varying levels of ruggedness. The terrain is designed by specifying its power-spectral shape in wavenumber space, inverting to physical space assuming random phases for all wave modes, and scaling the terrain to have a peak height of 200 m. For the case in question, these modest terrain fields expedite CI by up to 2–3 h, largely due to the impacts of the terrain on the size of, and subcloud support for, incipient cumuli. Terrain-induced circulations enhance subcloud kinetic energy on the mesoscale, which is realized as wider and longer-lived subcloud circulations. When the updraft branches of these circulations breach the level of free convection, they initiate wider and more persistent cumuli that subsequently undergo less entrainment-induced cloud dilution and detrainment-induced mass loss. As a result, the clouds become more vigorous and penetrate deeper into the troposphere. Larger-scale terrains are more effective than smaller-scale terrains in promoting CI because they induce larger enhancements in both the width and the persistence of subcloud updrafts.

54 ENVIRONMENTAL SCIENCES↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Irrigation Requirement Estimation using MODIS Vegetation Indices and Inverse Biophysical Modeling; A Case Study for Oran, Algeria

Human demand for food influences the water cycle through diversion and extraction of fresh water needed to support agriculture. Future population growth and economic development alone will substantially increase water demand and much of it for agricultural uses. For many semi-arid lands, socio-economic shifts are likely to exacerbate changes in climate as a driver of future water supply and demand. For these areas in particular, where the balance between water supply and demand is fragile, variations in regional climate can have potentially predictable effect on agricultural production. Satellite data and biophysically-based models provide a powerful method to quantify the interactions between local climate, plant growth and water resource requirements. In irrigated agricultural lands, satellite observations indicate high vegetation density while the precipitation amount indicates otherwise. This inconsistency between the observed precipitation and the observed canopy leaf density triggers the possibility that the observed high leaf density is due to an alternate source of water, irrigation. We explore an inverse process approach using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS), climatological data, and the NASA's Simple Biosphere model, SiB2, to quantitatively assess water demand in a semi-arid agricultural land by constraining the carbon and water cycles modeled under both equilibrium (balance between vegetation and prevailing local climate) and nonequilibrium (water added through irrigation) conditions. We postulate that the degree to which irrigated lands vary from equilibrium conditions is related to the amount of irrigation water used. We added water using two distribution methods: The first method adds water on top of the canopy and is a proxy for the traditional spray irrigation. The second method allows water to be applied directly into the soil layer and serves as proxy for drip irrigation. Our approach indicates that over the study site, for the month of July, spray irrigation resulted in an irrigation amount of about 1.4 mm per occurrence with an average frequency of occurrence of 24.6 hours. The simulated total monthly irrigation for July was 34.85 mm. In contrast, the drip irrigation resulted in less frequent irrigation events with an average water requirement about 57% less than that simulated during the spray irrigation case. The efficiency of the drip irrigation method rests on its reduction of the canopy interception loss compared to the spray irrigation method. When compared to a country-wide average estimate of irrigation water use, our numbers are quite low. We would have to revise the reported country level estimates downward to 17% or less

Bounoua, L.↗

An attempt to simulate laser-induced all-optical spin switching in a crystalline ferrimagnet

Interest in all-optical spin switching (AOS) is growing rapidly. The recent discovery of AOS in Mn 2 RuGa provides a much needed clean case of crystalline ferrimagnets for theoretical simulations. Here, we attempt to simulate it using the state-of-the-art first-principles method combined with the Heisenberg exchange model. We first compute the spin moments at two inequivalent manganese sites and then feed them into our model Hamiltonian. We employ an ultrafast laser pulse to switch the spins. We find that there is a similar optimal laser field amplitude to switch spins. However, we find that the exchange interaction has a significant effect on the system switchability. Weakening the exchange interaction could make the system unswitchable. This provides a crucial insight into the switching mechanism in ferrimagnets.

36 MATERIALS SCIENCE↗