Search NASA⌕ Search

SEARCH · Search NASA

Results for “engine 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,729 records · Page 96

FLUKA simulation validation to measure the beam energy of a cyclotron

FLUKA is a fully integrated particle physics MonteCarlo simulation package which has many applications in high energy experimental physics and engineering. When using a cyclotron for medical radionuclide production, it is important to know that the beam energy is suitable to cause the intended nuclear reaction and minimize unwanted impurities. Here, this work focuses on the comparison of FLUKA as a tool to measure the beam energy of a cyclotron with experimental results using the method of stacked copper foils interspaced by an aluminum degrader. We compared two methods to determine the activity ratio: determining the activity in each foil using the experimentally measured cross-sections from the National Nuclear Data Center (NNDC) cross-section data and calculating the activity in the two Cu foils using the FLUKA program. These results were compared to those reported by Gagnon et al., 2011. The FLUKA results gave energy measurements that were different than the true energy by a percentage in the range of 2.4% to 13%, suggesting that FLUKA is not recommended for pre-calculating the cyclotron beam energy using these low Z materials without experimental confirmation. It also confirms that the experimental cross-section data from the NNDC database is superior for this application giving energy values that ranged from 0.1% to 1.6% from the true energy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Detailed chemistry modeling of rotating detonations with dilute n -heptane sprays and preheated air

Utilization of liquid fuels is crucial to enabling commercialization of rotating detonation engines (RDE) in the near future. In this study, Eulerian-Lagrangian simulations are conducted for rotating detonative combustion with dilute n-heptane sprays and preheated air. Two-dimensional flattened configuration is used and a skeletal chemical mechanism with 44 species and 112 elementary reactions for n-heptane combustion is adopted. The flow structure, droplet distribution, and thermochemical parameters in the refill zone are first analyzed. It is shown that the mixture in the refill zone is heterogeneous, including evaporating droplets, vapor, and air. When the total temperature is below 950 K, the average equivalence ratio increases with the total temperature. When it is higher than 950 K, the average equivalence ratio is almost constant. Subsequently, the chemical explosive mode analysis is applied to identify the controlling reactions and dominant combustion modes in the fuel refill zone and reaction fronts. Results demonstrate that the initiation reaction (R104: n-C 7 H 16 + O 2 2-C 7 H 15 + HO 2 ) and low-temperature reaction (R107: RO 2 R'O 2 H) are dominant in the upstream and downstream of the refill zone, respectively. The intermediate species from low-temperature chemistry, R'O 2 H, is found to be important for the chemical explosive mode in the undetonated mixture. The influence of species diffusion and dispersed droplets is further analyzed. Results show that vapor autoignition facilitated by droplet evaporation occurs in the refill zone. Finally, the effects of the air total temperature on the detonation propagation speed and RDE propulsion performance are investigated. Further, it is found that the detonation propagation speed and specific impulse increase with air total temperature. The total pressure ratio first increases and then decreases as the air total temperature increases. Moreover, when the total temperature of the preheated air is above 1,300 K, the effects of low-temperature chemistry are negligible.

33 ADVANCED PROPULSION SYSTEMS↗

Multi-fidelity physics-informed machine learning for probabilistic damage diagnosis

Machine learning (ML) models are gaining popularity in structural health monitoring (SHM) because of their ability to learn the complex relationship between damage and sensor data. However, the lack of sufficient experimental data for structures with different degrees of damage is a key problem in training ML models for SHM. This problem can be alleviated by using physics-based models to generate the required training data to build physics-informed ML (PIML) models for SHM. However, it takes significant computational effort to perform enough high-fidelity simulations of the diagnostic test. It is thus desirable to know whether the available computational resource budget should be expended on numerous low-fidelity physics simulations, or a small number of high-fidelity simulations, or their combination. In this paper, we investigate this aspect of generating adequate training data for PIML, by constructing multi-fidelity PIML models. We evaluate the performance of several PIML models, trained with different amounts of low-fidelity and high-fidelity data, in locating hidden cracks in concrete structures using a nonlinear dynamics-based diagnosis technique. Here, we find that high-fidelity physics simulations that do not cover the (test and damage) parameter space do not improve the performance of diagnostic PIML models built using data from many low-fidelity physics simulations.

42 ENGINEERING↗

Net-zero energy building design and life-cycle cost analysis with air-source variable refrigerant flow and distributed photovoltaic systems

Net-zero energy building (NZEB) design has been gaining widespread acceptance as a major driver to reach feasible long-term goals of energy reductions in a building sector for the coming decades. Accordingly, net-zero energy targets for commercial buildings have been specified by federal and/or local government agencies in U.S. Although moving toward NZEB goals has technically feasible long-term advantages, economic issues are one of primarily concerns by building owners and designers for the application of advanced building technologies in NZEB design of their retrofit or new building construction. However, the lack of life cycle cost (LCC) data has still hindered reliable understanding of adopting NZEB design by building owners in a cost-effective manner. In response to this gap, this study conducts life-cycle cost analysis (LCCA) of NZEB design with an energy-efficient heating, ventilation, and air conditioning (HVAC) system and a solar power generation system. Here, to enable NZEB design, an air-source VRF heat pump (HP) type system is considered for energy-efficient HVAC equipment in an office building model with a grid-tied photovoltaic (PV) system. Results reveal that the LCC of NZEBs varies primarily due to a combined effect of annual site EUIs, on-site power generation, and cost-items of each location in U.S. climates. Some hot and mild climate zones in the U.S. represent that LCC values of NZEBs are considerably lower. In contrast, some other cold zones result in unattractive LCC values for NZEBs, caused by a combined effect of relatively lower potential of on-site PV electricity generation and city cost indices. Yet, with subsidized financing by the federal government for PV installation costs, LCC values of NZEBs become significantly attractive in most climate locations.

42 ENGINEERING↗

The Magnetic Genome of Two-Dimensional van der Waals Materials

Magnetism in two-dimensional (2D) van der Waals (vdW) materials has recently emerged as one of the most promising areas in condensed matter research, with many exciting emerging properties and significant potential for applications ranging from topological magnonics to low-power spintronics, quantum computing, and optical communications. In the brief time after their discovery, 2D magnets have blossomed into a rich area for investigation, where fundamental concepts in magnetism are challenged by the behavior of spins that can develop at the single layer limit. However, much effort is still needed in multiple fronts before 2D magnets can be routinely used for practical implementations. In this comprehensive review, prominent authors with expertise in complementary fields of 2D magnetism (i.e., synthesis, device engineering, magneto-optics, imaging, transport, mechanics, spin excitations, and theory and simulations) have joined together to provide a genome of current knowledge and a guideline for future developments in 2D magnetic materials research.

36 MATERIALS SCIENCE↗

Nanoscale roughness effect on Maxwell-like boundary conditions for the Boltzmann equation

It is well known that the roughness of the wall has an effect on microscale gas flows. This effect can be shown for large Knudsen numbers by using a numerical solution of the Boltzmann equation. However, when the wall is rough at a nanometric scale, it is necessary to use a very small mesh size which is much too expansive. An alternative approach is to incorporate the roughness effect in the scattering kernel of the boundary condition, such as the Maxwell-like kernel introduced by the authors in a previous paper. Here, we explain how this boundary condition can be implemented in a discrete velocity approximation of the Boltzmann equation. Moreover, the influence of the roughness is shown by computing the structure scattering pattern of mono-energetic beams of the incident gas molecules. The effect of the angle of incidence of these molecules, of their mass, and of the morphology of the wall is investigated and discussed in a simplified two-dimensional configuration. The effect of the azimuthal angle of the incident beams is shown for a three-dimensional configuration. Finally, the case of non-elastic scattering is considered. All these results suggest that our approach is a promising way to incorporate enough physics of gas-surface interaction, at a reasonable computing cost, to improve kinetic simulations of micro- and nano-flows.

42 ENGINEERING↗

Toward digital design at the exascale: An overview of project ICECap

High performance computing has entered the Exascale Age. Capable of performing over 1018 floating point operations per second, exascale computers, such as El Capitan, the National Nuclear Security Administration's first, have the potential to revolutionize the detailed in-depth study of highly complex science and engineering systems. However, in addition to these kind of whole machine “hero” simulations, exascale systems could also enable new paradigms in digital design by making petascale hero runs routine. Currently, untenable problems in complex system design, optimization, model exploration, and scientific discovery could all become possible. Motivated by the challenge of uncovering the next generation of robust high-yield inertial confinement fusion (ICF) designs, project ICECap (Inertial Confinement on El Capitan) attempts to integrate multiple advances in machine learning (ML), scientific workflows, high performance computing, GPU-acceleration, and numerical optimization to prototype such a future. Built on a general framework, ICECap is exploring how these technologies could broadly accelerate scientific discovery on El Capitan. In addition to our requirements, system-level design, and challenges, we describe some of the key technologies in ICECap, including ML replacements for multiphysics packages, tools for human-machine teaming, and algorithms for multifidelity design optimization under uncertainty. As a test of our prototype pre-El Capitan system, we advance the state-of-the art for ICF hohlraum design by demonstrating the optimization of a 17-parameter National Ignition Facility experiment and show that our ML-assisted workflow makes design choices that are consistent with physics intuition, but in an automated, efficient, and mathematically rigorous fashion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evidence of non-isentropic release from high residual temperatures in shocked metals measured with ultrafast x-ray diffraction

Shock experiments are widely used to understand the mechanical and electronic properties of matter under extreme conditions. However, after shock loading to a Hugoniot state, a clear description of the post-shock thermal state and its impacts on materials is still lacking. We used diffraction patterns from 100-fs x-ray pulses to investigate the temperature evolution of laser-shocked Al–Zr metal film composites at time delays ranging from 5 to 75 ns driven by a 120-ps short-pulse laser. We found significant heating of both Al and Zr after shock release, which can be attributed to heat generated by inelastic deformation. A conventional hydrodynamic model that employs (i) typical descriptions of Al and Zr mechanical strength and (ii) elevated strength responses (which might be attributed to an unknown strain rate dependence) did not fully account for the measured temperature increase, which suggests that other strength-related mechanisms (such as fine-scale void growth) could play an important role in thermal responses under shock wave loading/unloading cycles. Our results suggest that a significant portion of the total shock energy delivered by lasers becomes heat due to defect-facilitated plastic work, leaving less converted to kinetic energy. This heating effect may be common in laser-shocked experiments but has not been well acknowledged. High post-shock temperatures may induce phase transformation of materials during shock release. Another implication for the study is the preservability of magnetic records from planetary surfaces that have a shock history from frequent impact events.

42 ENGINEERING↗

Leader–follower fixed-time consensus of multi-agent systems with high-order integrator dynamics

The leader-follower fixed-time consensus of high-order multi-agent systems with external disturbances is investigated in this paper. A novel sliding manifold is designed to ensure that the tracking errors converge to zero in a fixed-time during the sliding motion. Then, a distributed control law is designed based on Lyapunov technique to drive the system states to the sliding manifold in finite-time independent of initial conditions. Finally, the efficiency of the proposed method is illustrated by numerical simulations.

42 ENGINEERING↗

Effect of hollow fiber membrane packing on the performance of modules formed with fiber tows

Membrane processes are the preferred option for many gas separations. Hollow fiber membranes are commonly used in these processes. Commercial hollow fiber membrane modules are often fabricated from small groups of fibers (i.e., a tow) that are arranged into a larger cylindrical bundle. The existence of fiber tows can lead to packing non-uniformity. This nonuniformity can result in nonideal fluid distribution in the shell and be detrimental to module performance. In this work, full three-dimensional models of flow and mass transfer within fiber bundles formed from tows have been developed. These models account for the detailed fiber arrangement within modules and are solved to evaluate module separation performance in terms of stage cut and required membrane area as a function product gas purity for carbon capture applications. Results are presented illustrating the effects of varying fiber packing arrangements within a tow and varying tow packing arrangements. Further, the results suggest smaller tows are preferred over larger tows, at fixed overall fiber packing and lower intra- and inter-tow packings are preferred. To reduce the computational cost of simulating large fiber bundles, the use of an equivalent planar bundle is proposed (EPB). A simple procedure for determining the geometry of the EPB is described and validated through comparisons of calculated performance metrics with results obtained from full three-dimensional simulations.

42 ENGINEERING↗

Scoping study of lower hybrid current drive for CFETR

The paper assesses the applicability of lower hybrid current drive (LHCD) for two potential operating scenarios for the China Fusion Engineering Test Reactor (CFETR): the “hybrid” scenario in which some of the plasma current is sustained by the Ohmic transformer, and the fully non-inductive “steady state” scenario. Here, the πScope workflow engine was used to set up a large number of ray tracing/Fokker-Planck simulations (> 10 4 ) with parametric scans in the antenna poloidal position and launched parallel refractive index (n || ) for both the hybrid and steady state scenarios. Modeling predicts efficient off-axis current drive (1.3 MA for 20 MW launched power) with a peak near ρ of 0.6-0.65 for waves launched from the high field side (HFS). Waves launched from the low field side (LFS) damp at larger radius (ρ > 0.73) with similar efficiency to HFS launch. Stability analysis of the CFETR scenarios favors current drive profiles peaked near the mid-radius, suggesting that HFS launch is preferable due to the current drive location. The effect of wave scattering from density blobs in the edge/scrape-off-layer region was assessed through rotation of the perpendicular wavenumber at the ray origin. Simulations show that this effect can be quite large both in efficiency and damping location, however by adjusting the launched n|| much of the unperturbed performance can be recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling and Analysis of DC Microgrids as Stochastic Hybrid Systems

This study proposes a method of predicting the influence of random load behavior on the dynamics of dc microgrids and distribution systems. This is accomplished by combining stochastic load models and deterministic microgrid models. Together, these elements constitute a stochastic hybrid system. The resulting model enables straightforward calculation of dynamic state moments, which are used to assess the probability of desirable operating conditions. Specific consideration is given to systems based on the dual active bridge (DAB) topology. Bounds are derived for the probability of zero voltage switching (ZVS) in DAB converters. A simple example is presented to demonstrate how these bounds may be used to improve ZVS performance as an optimization problem. In conclusion, predictions of state moment dynamics and ZVS probability assessments are verified through comparisons to Monte Carlo simulations.

42 ENGINEERING↗

Global System Reduction Order Modeling for Localized Feature Inclusion

The development of reduced-order models remains an active research area, despite advances in computational resources. Here we develop a novel order-reduction approach that is designed to incorporate isolated regions that contain, for example, nonlinearitites or accumulating damage. The approach is designed to use global modes of the overall system response, which are then naturally coupled to the response in the isolated region of interest. Two examples are provided to demonstrate both the accuracy and the computational efficiency of the proposed approach. The performance of this approach is compared to the exact response corresponding to a finite element simulation for the chosen problems. In addition, the accuracy and computational efficiency are shown relative to a standard Galerkin reduction based on the linear normal modes. It is found that the proposed reduction offer computational efficiency comparable to a Galerkin reduction, but more accurately represents the response of the system when both are compared to the finite element simulation.

42 ENGINEERING↗

Phonon localization in heat conduction

Nondiffusive phonon thermal transport, extensively observed in nanostructures, has largely been attributed to classical size effects, ignoring the wave nature of phonons. We report localization behavior in phonon heat conduction due to multiple scattering and interference events of broadband phonons, by measuring the thermal conductivities of GaAs/AlAs superlattices with ErAs nanodots randomly distributed at the interfaces. With an increasing number of superlattice periods, the measured thermal conductivities near room temperature increased and eventually saturated, indicating a transition from ballistic to diffusive transport. In contrast, at cryogenic temperatures the thermal conductivities first increased but then decreased, signaling phonon wave localization, as supported by atomistic Green's function simulations. The discovery of phonon localization suggests a new path forward for engineering phonon thermal transport.

36 MATERIALS SCIENCE↗

Effects of Hot Streak and Phantom Cooling on Heat Transfer in a Cooled Turbine Stage Including Particulate Deposition

The objective of this research effort was to develop a validated computational modeling capability for the characterization of the effects of hot streaks and particulate deposition on the heat load of modern gas turbines. This was accomplished with a multi-faceted approach including analytical, experimental, and computational components. A 1-year no cost extension request was approved for this effort, so the total duration was 4 years. The research effort succeeded in its ultimate objective by leveraging extensive experimental deposition studies complemented by computational modeling. Experiments were conducted with hot streaks, vane cooling, and combinations of hot streaks with vane cooling. These studies contributed to a significant body of corporate knowledge of deposition, in combination with particle rebound and deposition studies funded by other agencies, to provide suitable conditions for the development of a new model. The model includes the following physical phenomena: elastic deformation, plastic deformation, adhesion, and shear removal. It also incorporates material property sensitivity to temperature and tangential-normal velocity rebound cross-dependencies observed in experiments. The model is well-suited for incorporation in CFD simulations of complex gas turbine flows due to its algebraic (explicit) formulation. This report contains model predictions compared to coefficient of restitution data available in the open literature as well as deposition results from two different high temperature turbine deposition facilities. While the model comparisons with experiments are in many cases promising, several key aspects of particle deposition remain elusive. The simple phenomenological nature of the model allows for parametric dependencies to be evaluated in a straightforward manner. This effort also included the first-ever full turbine stage deposition model published in the open literature. The simulations included hot streaks and simulated vane cooling. The new deposition model was implemented into the CFD model as a wall boundary condition, with various particle sizes investigated in the simulation. Simulations utilizing a steady mixing plane formulation and an unsteady sliding mesh were conducted and the flow solution of each was validated against experimental data. Results from each of these simulations, including impact and capture distributions and efficiencies, were compared and potential reasons for differences discussed in detail. The inclusion of a large range of particle sizes allowed investigation of trends with particle size, such as increased radial migration and reduced sticking efficiency at the larger particle sizes. The unsteady simulation predicted lower sticking efficiencies on the rotor blades than the mixing plane simulation for the majority of particle sizes. This is postulated to be due to the preservation of the hot streak and cool vane wake through the vane-rotor interface (which are smeared out circumferentially in the mixing-plane simulation). The results reported here represent the successful implementation of a novel deposition model into validated vane-rotor flow solutions that include a non-uniform inlet temperature profile and simulated vane cooling.

42 ENGINEERING↗