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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A method for modelling peak signal statistics on a mobile satellite transponder

A simulation method is proposed. The simulation was developed to model the peak duration and energy content of signal peaks in a mobile communication satellite operating in a Frequency Division Multiple Access (FDMA) mode and presents an estimate of those power peaks for a system where the channels are modeled as band limited Gaussian noise, which is taken as a reasonable representation for Amplitude Commanded Single Sideband (ACSSB), Minimum Shift Keying (MSK), or Phase Shift Keying (PSK) modulated signals. The simulation results show that, under this hypothesis, the level of the signal power peaks for 10 percent, 1 percent, and 0.1 percent of the time are well described by a Rayleigh law and that their duration is extremely short and inversely proportional to the total FDM system bandwidth.

Bilodeau, Andre↗
NASA NTRS1984

Simulating a Three-Dimensional Flow in Pipes

Numerical method simulates in unsteady, incompressible flow in axisymmetric pipe in three dimensions. Method relies on velocity-vector expansion technique in which each vector function in expansion set is divergence-free and satisfies boundary conditions for viscous flow.

Leonard, A.↗
NASA NTRS1980

Vortex methods for flow simulation

Recent progress in the development of vortex methods and their applications to the numerical simulation of incompressible fluid flows are reviewed. Emphasis is on recent results concerning the accuracy of these methods, improvements in computational efficiency, and the development of three-dimensional methods. Simulations of several example flows which display some of the strengths and weaknesses of vortex methods are presented.

Leonard, A.↗
DOE OSTI2015

An Analysis Platform for Multiscale Hydrogeologic Modeling with Emphasis on Hybrid Multiscale Methods

One of the most significant challenges facing hydrogeologic modelers is the disparity between those spatial and temporal scales at which fundamental flow, transport and reaction processes can best be understood and quantified (e.g., microscopic to pore scales, seconds to days) and those at which practical model predictions are needed (e.g., plume to aquifer scales, years to centuries). While the multiscale nature of hydrogeologic problems is widely recognized, technological limitations in computational and characterization restrict most practical modeling efforts to fairly coarse representations of heterogeneous properties and processes. For some modern problems, the necessary level of simplification is such that model parameters may lose physical meaning and model predictive ability is questionable for any conditions other than those to which the model was calibrated. Recently, there has been broad interest across a wide range of scientific and engineering disciplines in simulation approaches that more rigorously account for the multiscale nature of systems of interest. In this paper, we review a number of such approaches and propose a classification scheme for defining different types of multiscale simulation methods and those classes of problems to which they are most applicable. Our classification scheme is presented in terms of a flow chart (Multiscale Analysis Platform or MAP), and defines several different motifs of multiscale simulation. Within each motif, the member methods are reviewed and example applications are discussed. We focus attention on hybrid multiscale methods, in which two or more models with different physics described at fundamentally different scales are directly coupled within a single simulation. Very recently these methods have begun to be applied to groundwater flow and transport simulations, and we discuss these applications in the context of our classification scheme. As computational and characterization capabilities continue to improve, we envision that hybrid multiscale modeling will become more common and may become a viable alternative to conventional single-scale models in the near future.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

Modeling of proton-induced radioactivation background in hard X-ray telescopes: Geant4-based simulation and its demonstration by Hitomi ’s measurement in a low Earth orbit

Hard X-ray astronomical observatories in orbit suffer from a significant amount of background due to radioactivation induced by cosmic-ray protons and/or geomagnetically trapped protons. Within the framework of a full Monte Carlo simulation, we present modeling of in-orbit instrumental background which is dominated by radioactivation. To reduce the computation time required by straightforward simulations of delayed emissions from activated isotopes, we insert a semi-analytical calculation that converts production probabilities of radioactive isotopes by interaction of the primary protons into decay rates at measurement time of all secondary isotopes. Therefore, our simulation method is separated into three steps: (1) simulation of isotope production, (2) semi-analytical conversion to decay rates, and (3) simulation of decays of the isotopes at measurement time. This method is verified by a simple setup that has a CdTe semiconductor detector, and shows a 100-fold improvement in efficiency over the straightforward simulation. To demonstrate its experimental performance, the simulation framework was tested against data measured with a CdTe sensor in the Hard X-ray Imager onboard the Hitomi X-ray Astronomy Satellite, which was put into a low Earth orbit with an altitude of 570 km and an inclination of 31°, and thus experienced a large amount of irradiation from geomagnetically trapped protons during its passages through the South Atlantic Anomaly. The simulation is able to treat full histories of the proton irradiation and multiple measurement windows. As a result, the simulation results agree very well with the measured data, showing that the measured background is well described by the combination of proton-induced radioactivation of the CdTe detector itself and thick Bi 4 Ge 3 O 12 scintillator shields, leakage of cosmic X-ray background and albedo gamma-ray radiation, and emissions from naturally contaminated isotopes in the detector system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Comparison of Optical Modelling Tools for Sunshape and Surface Slope Error

The sunshape and mirror surface slope error are two main factors that influence the characteristics of incoming and reflected solar beams in a solar concentrator. Different optical modelling tools, however, address these two factors with different methods and simplifications. In this study, simulation methods are reviewed in five optical modelling tools, namely SolTrace, Tonatiuh, Tracer, Heliosim and SolarPILOT. A series of simple test cases are simulated in each tool to compare and verify the implementation of sunshape and slope error models separately. This strategy has successfully identified discrepancies among the tools and contributed to a better understanding of their reasons. It has also uncovered few glitches in some of those tools and contributed to improve them. It is our hope that this study will help to resolve common misconceptions, to build confidence amongst CSP researches on the accuracy of tools with which they have less experience, and to ensure better agreement in optical modelling work in future.

14 SOLAR ENERGY↗
DOE OSTI2021

Electrochemical ion insertion from the atomic to the device scale

Electrochemical ion insertion involves coupled ion–electron transfer reactions, transport of guest species and redox of the host. The hosts are typically anisotropic solids with 2D conduction planes but can also be materials with 1D or isotropic transport pathways. These insertion compounds have traditionally been studied in the context of energy storage but also find extensive applications in electrocatalysis, optoelectronics and computing. Recent developments in operando, ultrafast and high-resolution characterization methods, as well as accurate theoretical simulation methods, have led to a renaissance in the understanding of ion-insertion compounds. In this Review, we present a unified framework for understanding insertion compounds across timescales and length scales ranging from atomic to device levels. Using graphite, transition metal dichalcogenides, layered oxides, oxyhydroxides and olivines as examples, we explore commonalities in these materials in terms of point defects, interfacial reactions and phase transformations. We illustrate similarities in the operating principles of various ion-insertion devices, ranging from batteries and electrocatalysts to electrochromics and thermal transistors, with the goal of unifying research across disciplinary boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Ion manipulations in structures for lossless ion manipulations (SLIM): computational evaluation of a 90° turn and a switch

The process of redirecting ions through 90° turns and ‘tee’ switches utilizing Structures for Lossless Ion Manipulations (SLIM) was evaluated using theoretical and simulation methods at 4 Torr pressure. SIMION simulations were used to optimize and evaluate conditions for performing turns without loss of signal intensity or ion mobility resolving power. Fundamental considerations indicated that the “race track” effect during ion turns may incur only small losses to the ion mobility resolving power at 4 Torr pressure for the typical plume widths predicted in an optimized SLIM ‘tee’ switch design. The dynamic switching of ions into orthogonal channels was also evaluated using SIMION ion trajectory simulations, and achieved similar performance. Simulation results were in close agreement with experimental results and were used to refine SLIM designs and applied potentials for their use.

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

Analysis of transmittance properties of surface plasmon modes on periodic solid/outline bowtie nanoantenna arrays

We numerically analyze the transmittance properties of the periodic outline bowtie nanoantenna arrays and examine the nanoantenna resonance conditions compared to their solid counterparts by means of finite element method. Simulation results show that the structural parameters may affect the nanoantenna resonance conditions, such as peak resonance wavelengths, propagation properties, and charge densities at spectral points of interest. Besides, the essences of transmittance spectra of the periodic solid/outline bowtie nanoantenna arrays corresponding to their bonding mode and anti-bonding mode are investigated as well.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

First-principles calculations of thermoelectric properties of TiN/MgO superlattices: The route for an enhancement of thermoelectric effects in artificial nanostructures

We present the thermoelectric properties of TiN/MgO superlattices employing first-principles calculation techniques. The Seebeck coefficients, the electrical conductances, the thermal conductances, and the figure of merit are investigated employing electrical and thermal transport calculations based on density functional theory combined with the nonequilibrium Green's function and nonequilibrium molecular dynamics simulation methods. The TiN/MgO superlattices with a small lattice mismatch at the interfaces are ideal systems to study the way for an enhancement of thermoelectric properties in artificial nanostructures. We find that the interfacial scattering between the two materials in the metal/insulator superlattices causes the electrical conductance to change rapidly, which enhances the Seebeck coefficient significantly. We show that the figure of merit for the artificial superlattice nanostructures has a much larger value compared with that of the bulk material and changes drastically with the superlattice configurations at the atomistic level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Investigations into the feasibility of optical-CT 3D dosimetry with minimal use of refractively matched fluids

Purpose: In optical-CT, the use of a refractively matched polyurethane solid-tank in place of a fluid bath has the potential to greatly increase practical convenience, reduce cost, and possibly improve the efficacy of flood corrections. This work investigates the feasibility of solid-tank optical-CT imaging for 3D dosimetry through computer simulation. Methods: A MATLAB ray-tracing simulation platform, ScanSim, was used to model a parallel-source telecentric optical-CT imaging system through a polyurethane solid-tank containing a central cylindrical hollow into which PRESAGE radiochromic dosimeters can be placed. A small amount of fluid fills the 1–5 mm gap between the dosimeter and the walls of the tank. The use of the solid-tank reduces the required amount of fluid by approximately 97%. To characterize the efficacy of solid-tank, optical-CT scanning simulations investigated sensitivity to refractive index (RI) mismatches between dosimeter, solid-tank, and fluid, for a variety of dosimeter (RI = 1.5–1.47) and fluid (RI = 1.55–1.0) combinations. Efficacy was evaluated through the usable radius (r{sub u}) metric, defined as the fraction of the radius of the dosimeter where measured dose is predicted to be within 2% of the ground truth entered into the simulation. Additional simulations examined the effect of increasing gap size (1–5 mm) between the dosimeter and solid-tank well. The effects of changing the lens tolerance (0.5°–5.0°) were also investigated. Results: As the RI mismatch between the dosimeter and solid-tank increased from 0 to 0.02, the usable radius decreased from 97.6% to 50.2%. The optimal fluid RI decreased nonlinearly from 1.5 to 1.34 as the mismatch increased and was up to 9% lower than the tank. Media mismatches between the dosimeter and solid-tank also exacerbate the effects of changing the gap size, with no easily quantifiable relationship with usable radius. Generally, the optimal fluid RI value increases as gap size increases and is closely matched to the dosimeter at large gap sizes (>3 mm). Increasing the telecentric lens tolerance increases the usable radius for all refractive media combinations and improves the maximum usable radius of mismatched media to that of perfectly matched media for tolerances >5.0°. The maximum usable radius can be improved up to a factor of 2 when lens tolerances are small (<1.0°). Conclusions: Dry solid-tank optical-CT imaging in a telecentric system is feasible if the dosimeter RI is a close match with the solid-tank (<0.01 difference), providing accurate dose measurements within ±2% of true dose to over 80% of the dosimeter volume. In order to achieve accurate measurements over 96% of the dosimeter volume (representing out to 2 mm from the dosimeter edge), the dosimeter-tank RI mismatch must be less than 0.005. Optimal results occur when the RI of the dosimeter and tank is the same, in which case the fluid will have the same RI. If mismatches between the tank and dosimeter RI occur, the RI of the matching fluid needs to be fine tuned to achieve the highest usable radius.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1998

PCEMCAN - Probabilistic Ceramic Matrix Composites Analyzer: User's Guide, Version 1.0

PCEMCAN (Probabalistic CEramic Matrix Composites ANalyzer) is an integrated computer code developed at NASA Lewis Research Center that simulates uncertainties associated with the constituent properties, manufacturing process, and geometric parameters of fiber reinforced ceramic matrix composites and quantifies their random thermomechanical behavior. The PCEMCAN code can perform the deterministic as well as probabilistic analyses to predict thermomechanical properties. This User's guide details the step-by-step procedure to create input file and update/modify the material properties database required to run PCEMCAN computer code. An overview of the geometric conventions, micromechanical unit cell, nonlinear constitutive relationship and probabilistic simulation methodology is also provided in the manual. Fast probability integration as well as Monte-Carlo simulation methods are available for the uncertainty simulation. Various options available in the code to simulate probabilistic material properties and quantify sensitivity of the primitive random variables have been described. The description of deterministic as well as probabilistic results have been described using demonstration problems. For detailed theoretical description of deterministic and probabilistic analyses, the user is referred to the companion documents "Computational Simulation of Continuous Fiber-Reinforced Ceramic Matrix Composite Behavior," NASA TP-3602, 1996 and "Probabilistic Micromechanics and Macromechanics for Ceramic Matrix Composites", NASA TM 4766, June 1997.

Shah, Ashwin R.↗
DOE OSTI2013

Downscaling Solar Power Output to 4-Seconds for Use in Integration Studies (Presentation)

High penetration renewable integration studies require solar power data with high spatial and temporal accuracy to quantify the impact of high frequency solar power ramps on the operation of the system. Our previous work concentrated on downscaling solar power from one hour to one minute by simulation. This method used clearness classifications to categorize temporal and spatial variability, and iterative methods to simulate intra-hour clearness variability. We determined that solar power ramp correlations between sites decrease with distance and the duration of the ramp, starting at around 0.6 for 30-minute ramps between sites that are less than 20 km apart. The sub-hour irradiance algorithm we developed has a noise floor that causes the correlations to approach ~0.005. Below one minute, the majority of the correlations of solar power ramps between sites less than 20 km apart are zero, and thus a new method to simulate intra-minute variability is needed. These intra-minute solar power ramps can be simulated using several methods, three of which we evaluate: a cubic spline fit to the one-minute solar power data; projection of the power spectral density toward the higher frequency domain; and average high frequency power spectral density from measured data. Each of these methods either under- or over-estimates the variability of intra-minute solar power ramps. We show that an optimized weighted linear sum of methods, dependent on the classification of temporal variability of the segment of one-minute solar power data, yields time series and ramp distributions similar to measured high-resolution solar irradiance data.

14 SOLAR ENERGY↗
DOE OSTI2013

Downscaling Solar Power Output to 4-Seconds for Use in Integration Studies: Preprint

High penetration renewable integration studies require solar power data with high spatial and temporal accuracy to quantify the impact of high frequency solar power ramps on the operation of the system. Our previous work concentrated on downscaling solar power from one hour to one minute by simulation. This method used clearness classifications to categorize temporal and spatial variability, and iterative methods to simulate intra-hour clearness variability. We determined that solar power ramp correlations between sites decrease with distance and the duration of the ramp, starting at around 0.6 for 30-minute ramps between sites that are less than 20 km apart. The sub-hour irradiance algorithm we developed has a noise floor that causes the correlations to approach ~0.005. Below one minute, the majority of the correlations of solar power ramps between sites less than 20 km apart are zero, and thus a new method to simulate intra-minute variability is needed. These intra-minute solar power ramps can be simulated using several methods, three of which we evaluate: a cubic spline fit to the one-minute solar power data; projection of the power spectral density toward the higher frequency domain; and average high frequency power spectral density from measured data. Each of these methods either under- or over-estimates the variability of intra-minute solar power ramps. We show that an optimized weighted linear sum of methods, dependent on the classification of temporal variability of the segment of one-minute solar power data, yields time series and ramp distributions similar to measured high-resolution solar irradiance data.

14 SOLAR ENERGY↗
DOE OSTI2017

Semiclassical approach to dynamics of interacting fermions

Understanding the behaviour of interacting fermions is of fundamental interest in many fields ranging from condensed matter to high energy physics. Developing numerically efficient and accurate simulation methods is an indispensable part of this. Already in equilibrium, fermions are notoriously hard to handle due to the sign problem. Out of equilibrium, an important outstanding problem is the efficient numerical simulation of the dynamics of these systems. In this work we develop a new semiclassical phase-space approach (a.k.a. the truncated Wigner approximation) for simulating the dynamics of interacting fermions in arbitrary dimensions. As fermions are essentially non-classical objects, a phase-space is constructed out of all fermionic bilinears. Classical phase-space is thus comprised of highly non-local (hidden) variables representing these bilinears, and the cost of the method is that it scales quadratic rather than linear with system size. We demonstrate the strength of the method by comparing the results to the exact quantum dynamics of fermion expansion in the Hubbard model and quantum thermalization in the Sachdev–Ye–Kitaev (SYK) model for small systems, where the semiclassics nearly perfectly reproduces correct results. We furthermore analyse fermion expansion in a larger, intractable by exact methods, 2D Hubbard model, which is directly relevant to recent cold atom experiments. - Highlights: • We propose a semiclassical method to simulate dynamics of interacting fermions. • A phase-space is constructed out of all fermionic bilinears. • We study quantum thermalization in the Sachdev–Ye–Kitaev model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Whole cell simulator

A whole cell model (WCM) is a comprehensive multi-scale computational model representing all the known biochemical processes in a cell. It relies on a variety of intracellular pathway models and omics data. This computational framework enables seamless integration of diverse simulation methods used in WCM such as stochastic simulation algorithm (SSA), ordinary differential equations (ODE), flux balance analysis (FBA), and logic-based approaches. These methods run simultaneously, not only for whole pathways but also for subsets of reactions. Furthermore, we aim to enable dynamic switching between the methods when beneficial. Is

GEORGAKOUDIS, GIORGIS↗
NASA NTRS1989

Perception of numerical methods in rarefied gasdynamics

The relationships between various numerical methods applied to problems in rarefied gasdynamics are discussed, with emphasis on conflicting viewpoints and computational requirements associated with physical simulation versus the numerical solution of the Boltzmann equation. The basic differences between the molecular dynamics and direct simulation methods are shown to affect their applicability to dense and rarefied flows. Methods for the probabilistic selection of representative collision in the direct simulation Monte Carlo method are reviewed. A method combining the most desirable features of the earlier methods is presented.

Bird, G. A.↗
NASA NTRS1985

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗