Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATION METHODS”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Hardware-in-the-loop grid simulator system and method

A hardware-in-the-loop (HIL) electrical grid simulation system and method that combines a reactive divider with a variable frequency converter to better mimic and control expected and unexpected parameters in an electrical grid. The invention provides grid simulation in a manner to allow improved testing of variable power generators, such as wind turbines, and their operation once interconnected with an electrical grid in multiple countries. The system further comprises an improved variable fault reactance (reactive divider) capable of providing a variable fault reactance power output to control a voltage profile, therein creating an arbitrary recovery voltage. The system further comprises an improved isolation transformer designed to isolate zero-sequence current from either a primary or secondary winding in a transformer or pass the zero-sequence current from a primary to a secondary winding.

17 WIND ENERGY↗
DOE OSTI2016

First-principle optimal local pseudopotentials construction via optimized effective potential method

The local pseudopotential (LPP) is an important component of orbital-free density functional theory, a promising large-scale simulation method that can maintain information on a material’s electron state. The LPP is usually extracted from solid-state density functional theory calculations, thereby it is difficult to assess its transferability to cases involving very different chemical environments. Here, we reveal a fundamental relation between the first-principles norm-conserving pseudopotential (NCPP) and the LPP. On the basis of this relationship, we demonstrate that the LPP can be constructed optimally from the NCPP for a large number of elements using the optimized effective potential method. Specially, our method provides a unified scheme for constructing and assessing the LPP within the framework of first-principles pseudopotentials. Our practice reveals that the existence of a valid LPP with high transferability may strongly depend on the element.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Analysis of nanoscale two-phase flow of argon using molecular dynamics

Two phase flows through micro and nanochannels have attracted a lot of attention because of their immense applicability to many advanced fields such as MEMS/NEMS, electronic cooling, bioengineering etc. In this work, a molecular dynamics simulation method is employed to study the condensation process of superheated argon vapor force driven flow through a nanochannel combining fluid flow and heat transfer. A simple and effective particle insertion method is proposed to model phase change of argon based on non-periodic boundary conditions in the simulation domain. Starting from a crystalline solid wall of channel, the condensation process evolves from a transient unsteady state where we study the influence of different wall temperatures and fluid wall interactions on interfacial and heat transport properties of two phase flows. Subsequently, we analyzed transient temperature, density and velocity fields across the channel and their dependency on varying wall temperature and fluid wall interaction, after a dynamic equilibrium is achieved in phase transition. Quasi-steady nonequilibrium temperature profile, heat flux and interfacial thermal resistance were analyzed. The results demonstrate that the molecular dynamics method, with the proposed particle insertion method, effectively solves unsteady nonequilibrium two phase flows at nanoscale resolutions whose interphase between liquid and vapor phase is typically of the order of a few molecular diameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Imaging the photochemistry of cyclobutanone using ultrafast electron diffraction: Experimental results

We investigated the ultrafast structural dynamics of cyclobutanone following photoexcitation at λ = 200 nm using gas-phase megaelectronvolt ultrafast electron diffraction. Our investigation complements the simulation studies of the same process within this special issue. It provides information about both electronic state population and structural dynamics through well-separable inelastic and elastic electron scattering signatures. We observe the depopulation of the photoexcited S 2 state of cyclobutanone with n3s Rydberg character through its inelastic electron scattering signature with a time constant of (0.29 ± 0.2) ps toward the S 1 state. The S 1 state population undergoes ring-opening via a Norrish Type-I reaction, likely while passing through a conical intersection with S 0 . The corresponding structural changes can be tracked by elastic electron scattering signatures. These changes appear with a delay of (0.14 ± 0.05) ps with respect to the initial photoexcitation, which is less than the S 2 depopulation time constant. This behavior provides evidence for the ballistic nature of the ring-opening once the S 1 state is reached. The resulting biradical species react further within (1.2 ± 0.2) ps via two rival fragmentation channels yielding ketene and ethylene, or propene and carbon monoxide. Furthermore, our study showcases the value of both gas-phase ultrafast diffraction studies as an experimental benchmark for nonadiabatic dynamics simulation methods and the limits in the interpretation of such experimental data without comparison with such simulations.

Carbon monoxide↗
DOE OSTI2017

Monte Carlo simulations of secondary electron emission due to ion beam milling

We present a Monte Carlo simulation study of secondary electron emission resulting from focused ion beam milling of a copper target. The basis of this study is a simulation code which simulates ion induced excitation and emission of secondary electrons, in addition to simulating focused ion beam sputtering and milling. This combination of features permits the simulation of the interaction between secondary electron emission, and the evolving target geometry as the ion beam sputters material. Previous ion induced SE Monte Carlo simulation methods have been restricted to predefined target geometries, while the dynamic target in the presented simulations makes this study relevant to image formation in ion microscopy, and chemically assisted ion beam etching, where the relationship between sputtering, and its effects on secondary electron emission, is important. We focus on a copper target, and validate our simulation against experimental data for a range of: noble gas ions, ion energies, ion/substrate angles and the energy distribution of the secondary electrons. We then provide a detailed account of the emission of secondary electrons resulting from ion beam milling; we quantify both the evolution of the yield as high aspect ratio valleys are milled, as well as the emission of electrons within these valleys that do not escape the target, but which are important to the secondary electron contribution to chemically assisted ion induced etching.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Phoenix v. 1.0-SNAPSHOT

Phoenix is a Java Virtual Machine (JVM) based library for performing mathematical and astrodynamics calculations. It consists of two primary sub-modules, phoenix-math and phoenix-astrodynamics. The mathematics package has a variety of mathematical classes for performing 3D transformations, geometric reasoning, and numerical analysis. The astrodynamics package has various classes and methods for computing locations, attitudes, accesses, and other values useful for general satellite modeling and simulation. Methods for computing celestial locations, such as the location of the Sun and Moon, are also included. Phoenix is meant to be used as a library within the context of a larger application. For example, it could be used for a web service, desktop client, or to compute simple values in a scripting environment.

Bastian, Mark↗
NASA NTRS2010

Assessment of Spatial Navigation and Docking Performance During Simulated Rover Tasks

INTRODUCTION: Following long-duration exploration transits, pressurized rovers will enhance surface mobility to explore multiple sites across Mars and other planetary bodies. Multiple rovers with docking capabilities are envisioned to expand the range of exploration. However, adaptive changes in sensorimotor and cognitive function may impair the crew s ability to safely navigate and perform docking tasks shortly after transition to the new gravitoinertial environment. The primary goal of this investigation is to quantify post-flight decrements in spatial navigation and docking performance during a rover simulation. METHODS: Eight crewmembers returning from the International Space Station will be tested on a motion simulator during four pre-flight and three post-flight sessions over the first 8 days following landing. The rover simulation consists of a serial presentation of discrete tasks to be completed within a scheduled 10 min block. The tasks are based on navigating around a Martian outpost spread over a 970 sq m terrain. Each task is subdivided into three components to be performed as quickly and accurately as possible: (1) Perspective taking: Subjects use a joystick to indicate direction of target after presentation of a map detailing current orientation and location of the rover with the task to be performed. (2) Navigation: Subjects drive the rover to the desired location while avoiding obstacles. (3) Docking: Fine positioning of the rover is required to dock with another object or align a camera view. Overall operator proficiency will be based on how many tasks the crewmember can complete during the 10 min time block. EXPECTED RESULTS: Functionally relevant testing early post-flight will develop evidence regarding the limitations to early surface operations and what countermeasures are needed. This approach can be easily adapted to a wide variety of simulated vehicle designs to provide sensorimotor assessments for other operational and civilian populations.

Wood, S. J.↗
NASA NTRS2013

NanoTHOR: Low-Cost Launch of Nanosatellites to Deep Space

The rapid development of high-performance nanosatellite platforms is enabling NASA and commercial ventures to consider performing missions to the asteroids, the Moon, and Mars at lower cost and on shorter timelines than traditional large spacecraft platforms. Currently, however, opportunities to launch secondary payloads to Earth escape are rare, and using chemical rockets to propel secondary payloads from LEO rideshares to escape is problematic due to the risks posed to primary payloads. The NanoTHOR effort has explored the technical feasibility and value proposition for using a simple momentum-exchange tether system to scavenge orbital energy from an upper stage in geostationary transfer orbit in order to boost nanosatellites to Earth escape. A NanoTHOR module will accomplish rapid transfer of a nanosatellite to an escape trajectory by deploying the nanosat at the end of a long, slender, high-strength tether and then using winching in the Earth's gravity!gradient to convert orbital angular momentum into rotational angular momentum. In the Phase I effort, we developed and simulated methods for controlling tether deployment and retraction to spin up a tether system, and these simulations demonstrated the feasibility of providing delta-Vs on the order of 800 m/s with a simple, low-mass tether system. Moreover, the NanoTHOR tether can act as a reusable in-space upper stage, boosting multiple nanosatellites on a single launch and doing so with a mass requirement lower than that of conventional rocket technologies. Serving as an escape-injection stage, NanoTHOR can enable a 6U CubeSat to deliver small payloads to Mars orbit, lunar orbit, and rendezvous with at least 110 of the known near-Earth asteroids. Evaluation of the technology readiness of the component technologies required for NanoTHOR indicate that the hardware required is all mid-TRL, and the lower-TRL controls and integration components can be advanced to mid-TRL with modest investment. By scavenging orbital energy from upper stages without any stored energy or propellant requirements, NanoTHOR permits deep-space nanosat missions to launch on rideshare opportunities, enabling NASA and commercial ventures to conduct affordable and frequent missions to explore deep space destinations.

Hoyt, Robert↗
DOE OSTI2021

NUMERICAL NOISE STUDY IN EIC BEAM-BEAM SIMULATIONS

In the Electron-Ion Collider (EIC) design, a flat beam collision scheme is adopted to achieve 1034 cm -2 s -1 luminosity. We found that the vertical growth of the proton beam is much larger than of the round beam. In this article we present the numerical noise study about the number of macroparticles, the electron slice number, and the electron bunch length. Both weak-strong and strong-strong simulation methods are used. It turns out the proton emittance growth in the strong-strong simulation mainly comes from the numerical noise. This study helps us to perform beam-beam simulation correctly for the EIC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Modelling and simulation of a porous core photonic crystal fibre for terahertz wave propagation

A porous core photonic crystal fibre based on conventional hexagonal lattice cladding is proposed for propagating terahertz radiation. The structure is designed and theoretically investigated using full vectorial finite element method. Simulation results show that at 300 µm core diameter, with a high porosity of 85%, and an operating frequency of 1.3 THz, the proposed fibre reduces the bulk absorption loss of cyclic olefin copolymer (TOPAS) by about 81%, which corresponds to an ultra-low effective material loss value of 0.039 cm{sup −1}. Furthermore, the proposed fibre shows near zero dispersion coefficient of 0.47 ps/THz/cm with an extremely small variation of 0.05 over a broad 1.3 THz bandwidth; with confinement and bending losses investigated and found to be negligibly low. It is anticipated that the proposed waveguide can potentially be used for short range transmission of terahertz radiation in the communication window.

42 ENGINEERING↗
DOE OSTI2015

GPU-based prompt gamma ray imaging from boron neutron capture therapy

Purpose: The purpose of this research is to perform the fast reconstruction of a prompt gamma ray image using a graphics processing unit (GPU) computation from boron neutron capture therapy (BNCT) simulations. Methods: To evaluate the accuracy of the reconstructed image, a phantom including four boron uptake regions (BURs) was used in the simulation. After the Monte Carlo simulation of the BNCT, the modified ordered subset expectation maximization reconstruction algorithm using the GPU computation was used to reconstruct the images with fewer projections. The computation times for image reconstruction were compared between the GPU and the central processing unit (CPU). Also, the accuracy of the reconstructed image was evaluated by a receiver operating characteristic (ROC) curve analysis. Results: The image reconstruction time using the GPU was 196 times faster than the conventional reconstruction time using the CPU. For the four BURs, the area under curve values from the ROC curve were 0.6726 (A-region), 0.6890 (B-region), 0.7384 (C-region), and 0.8009 (D-region). Conclusions: The tomographic image using the prompt gamma ray event from the BNCT simulation was acquired using the GPU computation in order to perform a fast reconstruction during treatment. The authors verified the feasibility of the prompt gamma ray image reconstruction using the GPU computation for BNCT simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-FG-BRB-07: GPU-Based Prompt Gamma Ray Imaging From Boron Neutron Capture Therapy

Purpose: The purpose of this research is to perform the fast reconstruction of a prompt gamma ray image using a graphics processing unit (GPU) computation from boron neutron capture therapy (BNCT) simulations. Methods: To evaluate the accuracy of the reconstructed image, a phantom including four boron uptake regions (BURs) was used in the simulation. After the Monte Carlo simulation of the BNCT, the modified ordered subset expectation maximization reconstruction algorithm using the GPU computation was used to reconstruct the images with fewer projections. The computation times for image reconstruction were compared between the GPU and the central processing unit (CPU). Also, the accuracy of the reconstructed image was evaluated by a receiver operating characteristic (ROC) curve analysis. Results: The image reconstruction time using the GPU was 196 times faster than the conventional reconstruction time using the CPU. For the four BURs, the area under curve values from the ROC curve were 0.6726 (A-region), 0.6890 (B-region), 0.7384 (C-region), and 0.8009 (D-region). Conclusion: The tomographic image using the prompt gamma ray event from the BNCT simulation was acquired using the GPU computation in order to perform a fast reconstruction during treatment. The authors verified the feasibility of the prompt gamma ray reconstruction using the GPU computation for BNCT simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Electrolytes between dielectric charged surfaces: Simulations and theory

We present a simulation method to study electrolyte solutions in a dielectric slab geometry using a modified 3D Ewald summation. The method is fast and easy to implement, allowing us to rapidly resum an infinite series of image charges. In the weak coupling limit, we also develop a mean-field theory which allows us to predict the ionic distribution between the dielectric charged plates. The agreement between both approaches, theoretical and simulational, is very good, validating both methods. Examples of ionic density profiles in the strong electrostatic coupling limit are also presented. Finally, we explore the confinement of charge asymmetric electrolytes between neutral surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Design-Space Data: Informing Common Design Decisions with Pre-Simulated Data

Design Space Exploration (DSE) analysis techniques represent a data-centric approach to integrating performance analysis in early design phases when there is the greatest potential to cheaply improve the energy efficiency of a building. We focus on a novel extension of DSE called Universal Design Space Exploration (UDSE), which leverages massive databases of pre-simulated analysis that represent all possible outcomes of common analysis workflows. These databases, called Design Spaces, become “universal” when a single pre-simulated design space can be re-applied to future unknown projects. Unlike current simulation methods, which require a design to exist before it can be analyzed and often take minutes or hours to simulate, UDSE leverages pre-simulation to deliver rapid and relevant insight as new designs are conceptualized. The data underpinning UDSE enables advanced statistical and Artificial Intelligence methods, allowing UDSE to deliver a greater understanding of the larger problem being explored, rather than simply delivering analysis of several pre-conceived design options. We believe that UDSE can provide instantaneous, relevant analysis for all building design projects at negligible cost. This paper has two main goals, to develop a relevant Universal Design Space that showcases the potential of UDSE and to release this data freely to industry and academia; thereby lowering the barrier to entry to digital literacy in statistics, ML and AI within the architecture, engineering, construction (AEC) industry.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Formulation of an experimental substructure model using a Craig-Bampton based transmission simulator

An experimental–analytical substructuring is attractive when there is motivation to replace one or more system subcomponents with an experimental model. This experimentally derived substructure can then be coupled to finite element models of the rest of the structure to predict the system response. The transmission simulator method couples a fixture to the component of interest during a vibration test in order to improve the experimental model for the component. The transmission simulator is then subtracted from the tested system to produce the experimental component. This method reduces ill-conditioning by imposing a least squares fit of constraints between substructure modal coordinates to connect substructures, instead of directly connecting physical interface degrees of freedom. This paper presents an alternative means of deriving the experimental substructure model, in which a Craig–Bampton representation of the transmission simulator is created and subtracted from the experimental measurements. The corresponding modal basis of the transmission simulator is described by the fixed-interface modes, rather than free modes that were used in the original approach. Moreover, these modes do a better job of representing the shape of the transmission simulator as it responds within the experimental system, leading to more accurate results using fewer modes. The new approach is demonstrated using a simple finite element model based example with a redundant interface.

42 ENGINEERING↗
NASA NTRS2018

Recent Advances to Estimation of Fixed-Interface Modal Models Using Dynamic Substructuring

In 2010, Allen & Mayes proposed to estimate the fixed-interface modes of a structure by measuring the modes of the structure bolted to a fixture and then applying constraints to the fixture using the transmission simulator method. While the method proved useful, and has indeed been used in studies since that point, a few peculiarities were noted. First, in some cases the estimated fixed-base natural frequencies were observed to converge very slowly to the true values (in simulated experiments) as the number of constraints was increased. To formulate these constraints, prior studies used only the free-interface modes of the fixture or the measured modes of the assembly. This work extends that to consider other sets of constraints, showing improved results. Furthermore, in some prior studies it has been observed that there were errors of more than 10% in the natural frequencies even when the fixture motion was hundreds of times smaller than the motion of the structure of interest (and so it had presumably been removed). This work explores this phenomenon, seeking to use the strain energy in the fixture, to the extent that it can be estimated using a test-analysis model for the fixture, as a metric to predict frequency error. The proposed methods are explored by applying them to simulated measurements from a beam and from the NASA Space Launch System coupled to the Mobile Launcher.

Allen, Mathew S.↗
DOE OSTI2024

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗