Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPUTER 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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Computational Fluid Dynamics Simulations of the Transonic Dynamics Tunnel Airstream Oscillator System

This paper presents computational fluid dynamics simulations of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) with the airstream oscillator system (AOS) in operation. Test section flow angle can be modulated using sinusoidally oscillating vanes which are located in pairs within the contraction section of the tunnel circuit. A compilation of previously published experimental data for validation has been reviewed that contains hot-wire anemometer and fast-response probe data acquired in an empty test section with the vanes oscillating in-phase and out-of-phase. These experimental data have not previously been compared with computational results. The present results provide an opportunity for validation of the simulations of the sinusoidal AOS flow field in the TDT test section. Computed tunnel centerline, lateral and cross-section contour data of flow angle and vorticity show that resolving the vane tip vortex is critical to accurately predicting the test section center flow angle. The computed results show, as do the experimental data, what appears to be a resonance at higher frequencies. It is currently not known what is the cause of that resonance. Further refinement of the mesh is expected to improve the comparisons with experimental data

Robert E Bartels↗
DOE OSTI2025

A Computational Framework for Simulations of Dissipative Nonadiabatic Dynamics on Hybrid Oscillator-Qubit Quantum Devices

Here, we introduce a computational framework for simulating nonadiabatic vibronic dynamics on circuit quantum electrodynamics (cQED) platforms. Our approach leverages hybrid oscillator-qubit quantum hardware with midcircuit measurements and resets, enabling the incorporation of environmental effects such as dissipation and dephasing. To demonstrate its capabilities, we simulate energy transfer dynamics in a triad model of photosynthetic chromophores inspired by natural antenna systems. We specifically investigate the role of dissipation during the relaxation dynamics following photoexcitation, where electronic transitions are coupled to the evolution of quantum vibrational modes. Our results indicate that hybrid oscillator-qubit devices, operating with noise levels below the intrinsic dissipation rates of typical molecular antenna systems, can achieve the simulation fidelity required for practical computations on near-term and early fault-tolerant quantum computing platforms.

Hamiltonians↗
DOE OSTI2021

Discussion on molecular dynamics (MD) simulations of the asphalt materials

The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.

36 MATERIALS SCIENCE↗
DOE OSTI2016

WE-H-207A-01: Computational Evaluation of High-Resolution 18F Positron Imaging Using Radioluminescence Microscopy with Lu2O3: Eu Thin-Film Scintillator

Purpose: Radioluminescence microscopy, an emerging and powerful tool for high resolution beta imaging, has been applied to molecular imaging of cellular metabolism to understand tumor biology. A novel thin-film (10 µm thickness) scintillator made of Lu{sub 2}O{sub 3}: Eu has been developed to enhance the system performance. However the advances of radioluminescence imaging with Lu{sub 2}O{sub 3}scintillator compared with that using conventional scintillator have not been explored theoretically to date. To validate the advantages of the thin-film scintillator, this study uses a novel computational simulation framework to evaluate the performance of radioluminescence microscopy using both conventional and thin-film scintillators. Methods: Numerical models for different stages of positron imaging are established. Positron from {sup 18}F passing through the scintillator and its neighbor structures are modeled by Monte-Carlo simulation using Geant4. The propagation and focus of photons by the microscope are modeled by convolution with a depth-varying point spread function generated by the Gibson-Lanni model. Photons focused on the detector plane are then captured and converted into electronic signals by an electron multiplication (EM) CCD camera, which is described by a photosensor model considering various noises and charge amplification. Results: The performance metrics of radioluminescence imaging with a thin-film Lu{sub 2}O{sub 3} and conventional CdWO{sub 4} scintillator are compared, including spatial resolution, sensitivity, positron track area and intensity. The spatial resolution of Lu{sub 2}O{sub 3} system can achieve 10 µm maximally, a 12 µm enhancement from that obtained from CdWO{sub 4} system. Meanwhile, the system with Lu{sub 2}O{sub 3} scintillator can provide a higher mean sensitivity: 40% compared with that (21.5%) obtained from CdWO{sub 4} system. Moreover, the simulation results are in good agreement with previous experimental measurements. Conclusion: This study provides a new theoretical understanding of our imaging system and has the potential to promote the development of radioluminescence microscopy for more reliable and robust application on the functional imaging of delicate biological structures. The authors acknowledge funding from NIH grant R01CA186275 and SBIR grant 1R43GM110888-01.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

(Project 18-15502) Reducing Uncertainty in Radionuclide Transport Prediction Using Multiple Environmental Tracers (Final Report)

In order to successfully site and design a nuclear waste disposal facility, DOE scientists are required to show safe containment of the radioactive waste for up to 1 million years. A significant hurdle to accurately and convincingly demonstrating disposal safety is predicting the fate of radioactive elements once they enter the groundwater system surrounding the repository. These predictions are often made with computer models, which contain accurate physics of groundwater movement and chemical reactions that occur during groundwater flow. As computer power increases, the physics and chemistry of these computer simulators can become more and more realistic and the physically based error decreases. However, the ability of these computer models to provide accurate predictions in a specific place, over long time periods, requires the scientists and engineers to know the subsurface properties of the Earth that control groundwater movement. In particular, groundwater scientists and engineers need to know the groundwater fluid velocity, which can change over time and strongly vary with location within the groundwater system. Because the Earth’s subsurface cannot be directly seen, and can only be sampled at drilling locations, the properties of the groundwater system are never known completely, and computer models of groundwater transport will always have some amount of uncertainty. This project’s principal goal was to use chemicals and isotope “tracers”, which have been introduced to the groundwater system by natural processes over long time periods, to help inform computer models of the groundwater velocity and subsurface properties. The goal was to calculate how much better predictions of groundwater transport were when these “tracers” were used to inform the computer models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Novel Relativistic Electronic Structure Theories for Actinide-Containing Compounds

Actinides of importance to basic energy sciences contain electrons moving at speed comparable to the speed of light. Reliable computational simulation of these electrons and hence actinide chemistry requires accurate description of relativistic effects. The present project advances computational actinide chemistry with development of new methodologies, algorithms, and computer programs in relativistic quantum chemistry, as well as applications to actinide chemistry and spectroscopy. A new “electrons-only” exact two-component approach has been developed to provide efficient treatments of relativistic effects, while maintaining chemical accuracy. New computational algorithms developed here extend the applicability of relativistic electron-correlation methods to larger molecules. The method-development work in this project also features the first implementation of analytic gradient technique for relativistic electron-correlation methods, which provides significantly enhanced ability to compute properties for molecules containing actinides. The applicability and usefulness of these new methods and computer programs have been demonstrated in calculations of actinide-containing molecules to facilitate understanding of actinide chemistry and spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1973

Pilot performance during a simulated standard instrument procedure turn with and without a predictor display

A simulator study was conducted to measure the effectiveness of predictor information incorporated into a CRT display of a computer simulated aircraft's horizontal and vertical situation. Professional pilots served as subjects for the task of executing a standard instrument procedure turn at constant altitude in constant crosswinds with and without their predicted ground track displayed. The results showed that the display with the predicted ground track was markedly and significantly superior to the display without this information and that the subjects were generally satisfied with this type of information. Mean rms lateral path error was independent of the crosswind velocity with the predictor information, and increased without it and with increasing wind velocity. Rms stick activity decreased with the predictor display which also uncoupled aileron and elevator activity.

Kreifeldt, J. G.↗
NASA NTRS1973

Calculated performance map of a 4 1/2-stage 15.0 centimeter (5.9 inch) mean diameter turbine designed for a turbofan simulator

The overall performance of an existing high-ratio turbine is calculated analytically over a range of speed and pressure ratio in order to determine its capability for other applications. The analytical performance covers a speed range from 50 to 120 percent of design and a pressure-ratio range from 5.0 to 35.0. The turbine was designed for a 50.8 centimeter (20.0 in.) tip diameter turbofan simulator. Computed results are compared with the experimental turbine data obtained from testing three fan configurations with the turbofan simulator in air. The comparison indicates good agreement over the range of speeds and pressure ratios covered by the experimental data.

Wasserbauer, C. A.↗
NASA NTRS1973

Pilot Performance During A Simulated Standard Instrument Procedure Turn With and Without A Predictor Display

A simulator study was conducted to measure the effectiveness of predictor information incorporated into a CRT display of a computer simulated aircraft's horizontal and vertical situation. Professional pilots served as subjects for the task of executing a standard instrument procedure turn at constant altitude in constant crosswinds with and without their predicted ground track displayed. The results show that the display with the predicted ground track was markedly and significantly superior to the display without this information and that the subjects were generally satisfied with this type of information. Mean rms lateral path error was independent of the crosswind velocity with the predictor information, and increased without it with increasing wind velocity. Rms stick activity decreased with the predictor display which also uncoupled aileron and elevator activity.

John G Kreifeldt↗
NASA NTRS1996

Coherent Lidar Design and Performance Verification

This final report summarizes the investigative results from the 3 complete years of funding and corresponding publications are listed. The first year saw the verification of beam alignment for coherent Doppler lidar in space by using the surface return. The second year saw the analysis and computerized simulation of using heterodyne efficiency as an absolute measure of performance of coherent Doppler lidar. A new method was proposed to determine the estimation error for Doppler lidar wind measurements without the need for an independent wind measurement. Coherent Doppler lidar signal covariance, including wind shear and turbulence, was derived and calculated for typical atmospheric conditions. The effects of wind turbulence defined by Kolmogorov spatial statistics were investigated theoretically and with simulations. The third year saw the performance of coherent Doppler lidar in the weak signal regime determined by computer simulations using the best velocity estimators. Improved algorithms for extracting the performance of velocity estimators with wind turbulence included were also produced.

Frehlich, Rod↗
NASA NTRS1997

Initialization and Simulation of Three-Dimensional Aircraft Wake Vortices

This paper studies the effects of axial velocity profiles on vortex decay, in order to properly initialize and simulate three-dimensional wake vortex flow. Analytical relationships are obtained based on a single vortex model and computational simulations are performed for a rather practical vortex wake, which show that the single vortex analytical relations can still be applicable at certain streamwise sections of three-dimensional wake vortices.

Ash, Robert L.↗
NASA NTRS2012

Observing Planets and Small Bodies in Sputtered High Energy Atom (SHEA) Fluxes

The evolution of the surfaces of bodies unprotected by either strong magnetic fields or thick atmospheres in the Solar System is caused by various processes, induced by photons, energetic ions and micrometeoroids. Among these processes, the continuous bombardment of the solar wind or energetic magnetospheric ions onto the bodies may significantly affect their surfaces, with implications for their evolution. Ion precipitation produces neutral atom releases into the exosphere through ion sputtering, with velocity distribution extending well above the particle escape limits. We refer to this component of the surface ejecta as sputtered high-energy atoms (SHEA). The use of ion sputtering emission for studying the interaction of exposed bodies (EB) with ion environments is described here. Remote sensing in SHEA in the vicinity of EB can provide mapping of the bodies exposed to ion sputtering action with temporal and mass resolution. This paper speculates on the possibility of performing remote sensing of exposed bodies using SHEA The evolution of the surfaces of bodies unprotected by either strong magnetic fields or thick atmospheres in the Solar System is caused by various processes, induced by photons, energetic ions and micrometeoroids. Among these processes, the continuous bombardment of the solar wind or energetic magnetospheric ions onto the bodies may significantly affect their surfaces, with implications for their evolution. Ion precipitation produces neutral atom releases into the exosphere through ion sputtering, with velocity distribution extending well above the particle escape limits. We refer to this component of the surface ejecta as sputtered high-energy atoms (SHEA). The use of ion sputtering emission for studying the interaction of exposed bodies (EB) with ion environments is described here. Remote sensing in SHEA in the vicinity of EB can provide mapping of the bodies exposed to ion sputtering action with temporal and mass resolution. This paper speculates on the possibility of performing remote sensing of exposed bodies using SHEA and suggests the need for quantitative results from laboratory simulations and molecular physic modeling in order to understand SHEA data from planetary missions. In the Appendix, referenced computer simulations using existing sputtering data are reviewed.

Milillo, A.↗
DOE OSTI2015

CAMELOT: A machine learning approach for coarse-grained simulations of aggregation of block-copolymeric protein sequences

We report the development and deployment of a coarse-graining method that is well suited for computer simulations of aggregation and phase separation of protein sequences with block-copolymeric architectures. Our algorithm, named CAMELOT for Coarse-grained simulations Aided by MachinE Learning Optimization and Training, leverages information from converged all atom simulations that is used to determine a suitable resolution and parameterize the coarse-grained model. To parameterize a system-specific coarse-grained model, we use a combination of Boltzmann inversion, non-linear regression, and a Gaussian process Bayesian optimization approach. The accuracy of the coarse-grained model is demonstrated through direct comparisons to results from all atom simulations. We demonstrate the utility of our coarse-graining approach using the block-copolymeric sequence from the exon 1 encoded sequence of the huntingtin protein. This sequence comprises of 17 residues from the N-terminal end of huntingtin (N17) followed by a polyglutamine (polyQ) tract. Simulations based on the CAMELOT approach are used to show that the adsorption and unfolding of the wild type N17 and its sequence variants on the surface of polyQ tracts engender a patchy colloid like architecture that promotes the formation of linear aggregates. These results provide a plausible explanation for experimental observations, which show that N17 accelerates the formation of linear aggregates in block-copolymeric N17-polyQ sequences. The CAMELOT approach is versatile and is generalizable for simulating the aggregation and phase behavior of a range of block-copolymeric protein sequences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Building complex simulations rapidly using MATRIX(x): The Space Station redesign

MSFC's quick response to the Space Station redesign effort last year required the development of a computer simulation to model the attitude and station-keeping dynamics of a complex body with rotating solar arrays in orbit around the Earth. The simulation was written using a rapid-prototyping graphical simulation and design tool called MATRIX(x) and provided the capability to quickly remodel complex configuration changes by icon manipulation using a mouse. The simulation determines time-dependent inertia properties, and models forces and torques from gravity-gradient, solar radiation, and aerodynamic disturbances. Surface models are easily built from a selection of beams, plates, tetrahedrons, and cylinders. An optimization scheme was written to determine the torque equilibrium attitudes that balance gravity-gradient and aerodynamic torques over an orbit, and propellant-usage estimates were determined. The simulation has been adapted to model the attitude dynamics for small spacecraft.

Carrington, C. K.↗
DOE OSTI2014

Model-based feasibility assessment and evaluation of prostate hyperthermia with a commercial MR-guided endorectal HIFU ablation array

Purpose: Feasibility of targeted and volumetric hyperthermia (40–45 °C) delivery to the prostate with a commercial MR-guided endorectal ultrasound phased array system, designed specifically for thermal ablation and approved for ablation trials (ExAblate 2100, Insightec Ltd.), was assessed through computer simulations and tissue-equivalent phantom experiments with the intention of fast clinical translation for targeted hyperthermia in conjunction with radiotherapy and chemotherapy. Methods: The simulations included a 3D finite element method based biothermal model, and acoustic field calculations for the ExAblate ERUS phased array (2.3 MHz, 2.3 × 4.0 cm{sup 2}, ∼1000 channels) using the rectangular radiator method. Array beamforming strategies were investigated to deliver protracted, continuous-wave hyperthermia to focal prostate cancer targets identified from representative patient cases. Constraints on power densities, sonication durations and switching speeds imposed by ExAblate hardware and software were incorporated in the models. Preliminary experiments included beamformed sonications in tissue mimicking phantoms under MR temperature monitoring at 3 T (GE Discovery MR750W). Results: Acoustic intensities considered during simulation were limited to ensure mild hyperthermia (T{sub max} < 45 °C) and fail-safe operation of the ExAblate array (spatial and time averaged acoustic intensity I{sub SATA} < 3.4 W/cm{sup 2}). Tissue volumes with therapeutic temperature levels (T > 41 °C) were estimated. Numerical simulations indicated that T > 41 °C was calculated in 13–23 cm{sup 3} volumes for sonications with planar or diverging beam patterns at 0.9–1.2 W/cm{sup 2}, in 4.5–5.8 cm{sup 3} volumes for simultaneous multipoint focus beam patterns at ∼0.7 W/cm{sup 2}, and in ∼6.0 cm{sup 3} for curvilinear (cylindrical) beam patterns at 0.75 W/cm{sup 2}. Focused heating patterns may be practical for treating focal disease in a single posterior quadrant of the prostate and diffused heating patterns may be useful for heating quadrants, hemigland volumes or even bilateral targets. Treatable volumes may be limited by pubic bone heating. Therapeutic temperatures were estimated for a range of physiological parameters, sonication duty cycles and rectal cooling. Hyperthermia specific phasing patterns were implemented on the ExAblate prostate array and continuous-wave sonications (∼0.88 W/cm{sup 2}, 15 min) were performed in tissue-mimicking material with real-time MR-based temperature imaging (PRFS imaging at 3.0 T). Shapes of heating patterns observed during experiments were consistent with simulations. Conclusions: The ExAblate 2100, designed specifically for thermal ablation, can be controlled for delivering continuous hyperthermia in prostate while working within operational constraints.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Spanning the Core Design Space with Binary Perturbations - A Mathematical Investigation

A reactor core design engineer, even one with many years of experience or one using advanced fuel pattern optimization software, will likely fine tune the pattern before finalizing the design. Fine tuning can further minimize local enthalpy peaks, optimize reactivity to extend burnup, and add additional margin to limits on reactivity coefficients or heat flux. For reactors that have reached an equilibrium fuel cycle, the number of feeds each cycle is typically the similar to previous cycles as is their placement in the core. It is then the engineer's responsibly to efficiently select the placement of burned assemblies to ensure optimal fuel utilization, satisfy safety limits, and meet cycle energy requirements. One of the engineer's most commonly used tools in this fine tuning investigation is to swap the positions of two assemblies, most likely two previously irradiated ('burned') assemblies since their reactivity is lower than fresh fuel. Such an operation is referred to in this work as a binary perturbation or a 'swap'. If quarter core symmetry is maintained, the perturbation translates to the other three quadrants. Likewise if eighth core symmetry is maintained, the perturbation translates to other seven octants. The question is then, how much of the core design space can binary perturbations possibly span? This work seeks to quantitatively define the space spanned by m number of binary perturbations as a function of variable size through combinatorial principles and computer simulation. The results of this study are useful in understanding the in-core fuel management problem faced by both professionals who employ hand manipulations and those who attempt stochastic solutions of the problem.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2008

On-Orbit Multi-Field Wavefront Control with a Kalman Filter

A document describes a multi-field wavefront control (WFC) procedure for the James Webb Space Telescope (JWST) on-orbit optical telescope element (OTE) fine-phasing using wavefront measurements at the NIRCam pupil. The control is applied to JWST primary mirror (PM) segments and secondary mirror (SM) simultaneously with a carefully selected ordering. Through computer simulations, the multi-field WFC procedure shows that it can reduce the initial system wavefront error (WFE), as caused by random initial system misalignments within the JWST fine-phasing error budget, from a few dozen micrometers to below 50 nm across the entire NIRCam Field of View, and the WFC procedure is also computationally stable as the Monte-Carlo simulations indicate. With the incorporation of a Kalman Filter (KF) as an optical state estimator into the WFC process, the robustness of the JWST OTE alignment process can be further improved. In the presence of some large optical misalignments, the Kalman state estimator can provide a reasonable estimate of the optical state, especially for those degrees of freedom that have a significant impact on the system WFE. The state estimate allows for a few corrections to the optical state to push the system towards its nominal state, and the result is that a large part of the WFE can be eliminated in this step. When the multi-field WFC procedure is applied after Kalman state estimate and correction, the stability of fine-phasing control is much more certain. Kalman Filter has been successfully applied to diverse applications as a robust and optimal state estimator. In the context of space-based optical system alignment based on wavefront measurements, a KF state estimator can combine all available wavefront measurements, past and present, as well as measurement and actuation error statistics to generate a Maximum-Likelihood optimal state estimator. The strength and flexibility of the KF algorithm make it attractive for use in real-time optical system alignment when WFC alone cannot effectively align the system.

Lou, John↗