Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized 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
DOE OSTI2013

A NEW MODEL FOR MIXING BY DOUBLE-DIFFUSIVE CONVECTION (SEMI-CONVECTION). II. THE TRANSPORT OF HEAT AND COMPOSITION THROUGH LAYERS

Regions of stellar and planetary interiors that are unstable according to the Schwarzschild criterion, but stable according to the Ledoux criterion, are subject to a form of oscillatory double-diffusive (ODD) convection often called ''semi-convection''. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and composition by ODD convection, and ultimately propose a new 1D prescription that can be used in stellar and planetary structure and evolution models. The first paper in this series demonstrated that under certain conditions ODD convection spontaneously transitions from an initial homogeneous state of weak wave-breaking turbulence into a staircase of fully convective layers, which results in a substantial increase in the transport of heat and composition. Here, we present simulations of ODD convection in this layered regime, we describe the dynamical behavior of the layers, and we derive empirical scaling laws for the transport through layered convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A NUMERICAL METHOD FOR STUDYING SUPER-EDDINGTON MASS TRANSFER IN DOUBLE WHITE DWARF BINARIES

We present a numerical method for the study of double white dwarf (DWD) binary systems at the onset of super-Eddington mass transfer. We incorporate the physics of ideal inviscid hydrodynamical flow, Newtonian self-gravity, and radiation transport on a three-dimensional uniformly rotating cylindrical Eulerian grid. Care has been taken to conserve the key physical quantities such as angular momentum and energy. Our new method conserves total energy to a higher degree of accuracy than other codes that are presently being used to model mass transfer in DWD systems. We present the results of verification tests and simulate the first 20 + orbits of a binary system of mass ratio q 0.7 at the onset of dynamically unstable direct impact mass transfer. The mass transfer rate quickly exceeds the critical Eddington limit by many orders of magnitude, and thus we are unable to model a trans-Eddington phase. It appears that radiation pressure does not significantly affect the accretion flow in the highly super-Eddington regime. An optically thick common envelope forms around the binary within a few orbits. Although this envelope quickly exceeds the spatial domain of the computational grid, the fraction of the common envelope that exceeds zero gravitational binding energy is extremely small, suggesting that radiation-driven mass loss is insignificant in this regime. It remains to be seen whether simulations that capture the trans-Eddington phase of such flows will lead to the same conclusion or show that substantial material gets expelled.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

On the soft limit of the large scale structure power spectrum: UV dependence

We derive a non-perturbative equation for the large scale structure power spectrum of long-wavelength modes. Thereby, we use an operator product expansion together with relations between the three-point function and power spectrum in the soft limit. The resulting equation encodes the coupling to ultraviolet (UV) modes in two time-dependent coefficients, which may be obtained from response functions to (anisotropic) parameters, such as spatial curvature, in a modified cosmology. We argue that both depend weakly on fluctuations deep in the UV. As a byproduct, this implies that the renormalized leading order coefficient(s) in the effective field theory (EFT) of large scale structures receive most of their contribution from modes close to the non-linear scale. Consequently, the UV dependence found in explicit computations within standard perturbation theory stems mostly from counter-term(s). We confront a simplified version of our non-perturbative equation against existent numerical simulations, and find good agreement within the expected uncertainties. Our approach can in principle be used to precisely infer the relevance of the leading order EFT coefficient(s) using small volume simulations in an 'anisotropic separate universe' framework. Our results suggest that the importance of these coefficient(s) is a ∼ 10% effect, and plausibly smaller.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Nonequilibrium and local detection of the normal fraction of a trapped two-dimensional Bose gas

We propose a method to measure the normal fraction of a two-dimensional Bose gas, a quantity that generally differs from the noncondensed fraction. The idea is based on applying a spatially oscillating artificial gauge field to the atoms. The response of the atoms to the gauge field can be read out either mechanically from the deposited energy into the cloud or optically from the macroscopic optical properties of the atomic gas. The local nature of the proposed scheme allows one to reconstruct the spatial profile of the superfluid component; furthermore, the proposed method does not require having established thermal equilibrium in the gas in the presence of the gauge field. The theoretical description of the system is based on a generalization of the Dum-Olshanii theory of artificial gauge fields to the interacting many-body context. The efficiency of the proposed measurement scheme is assessed by means of classical field numerical simulations. An explicit atomic level scheme minimizing disturbing effects such as spontaneous emission and light shifts is proposed for {sup 87}Rb atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2011

Magnetic-field generation and electron-collimation analysis for propagating fast electron beams in overdense plasmas

An analytical fluid model is proposed for artificially collimating fast electron beams produced in the interaction of ultraintense laser pulses with specially engineered low-density-core-high-density-cladding structure targets. Since this theory clearly predicts the characteristics of the spontaneously generated magnetic field and its dependence on the plasma parameters of the targets transporting fast electrons, it is of substantial relevance to the target design for fast ignition. The theory also reveals that the rapid changing of the flow velocity of the background electrons in a transverse direction (perpendicular to the flow velocity) caused by the density jump dominates the generation of a spontaneous interface magnetic field for these kinds of targets. It is found that the spontaneously generated magnetic field reaches as high as 100 MG, which is large enough to collimate fast electron transport in overdense plasmas. This theory is also supported by numerical simulations performed using a two-dimensional particle-in-cell code. It is found that the simulation results agree well with the theoretical analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Optimized generation of high resolution breast anthropomorphic software phantoms

Purpose: The authors present an efficient method for generating anthropomorphic software breast phantoms with high spatial resolution. Employing the same region growing principles as in their previous algorithm for breast anatomy simulation, the present method has been optimized for computational complexity to allow for fast generation of the large number of phantoms required in virtual clinical trials of breast imaging. Methods: The new breast anatomy simulation method performs a direct calculation of the Cooper's ligaments (i.e., the borders between simulated adipose compartments). The calculation corresponds to quadratic decision boundaries of a maximum a posteriori classifier. The method is multiscale due to the use of octree-based recursive partitioning of the phantom volume. The method also provides user-control of the thickness of the simulated Cooper's ligaments and skin. Results: Using the proposed method, the authors have generated phantoms with voxel size in the range of (25-1000 {mu}m){sup 3}/voxel. The power regression of the simulation time as a function of the reciprocal voxel size yielded a log-log slope of 1.95 (compared to a slope of 4.53 of our previous region growing algorithm). Conclusions: A new algorithm for computer simulation of breast anatomy has been proposed that allows for fast generation of high resolution anthropomorphic software phantoms.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Rapid Monte Carlo simulation of detector DQE(f)

Purpose: Performance optimization of indirect x-ray detectors requires proper characterization of both ionizing (gamma) and optical photon transport in a heterogeneous medium. As the tool of choice for modeling detector physics, Monte Carlo methods have failed to gain traction as a design utility, due mostly to excessive simulation times and a lack of convenient simulation packages. The most important figure-of-merit in assessing detector performance is the detective quantum efficiency (DQE), for which most of the computational burden has traditionally been associated with the determination of the noise power spectrum (NPS) from an ensemble of flood images, each conventionally having 10{sup 7} − 10{sup 9} detected gamma photons. In this work, the authors show that the idealized conditions inherent in a numerical simulation allow for a dramatic reduction in the number of gamma and optical photons required to accurately predict the NPS. Methods: The authors derived an expression for the mean squared error (MSE) of a simulated NPS when computed using the International Electrotechnical Commission-recommended technique based on taking the 2D Fourier transform of flood images. It is shown that the MSE is inversely proportional to the number of flood images, and is independent of the input fluence provided that the input fluence is above a minimal value that avoids biasing the estimate. The authors then propose to further lower the input fluence so that each event creates a point-spread function rather than a flood field. The authors use this finding as the foundation for a novel algorithm in which the characteristic MTF(f), NPS(f), and DQE(f) curves are simultaneously generated from the results of a single run. The authors also investigate lowering the number of optical photons used in a scintillator simulation to further increase efficiency. Simulation results are compared with measurements performed on a Varian AS1000 portal imager, and with a previously published simulation performed using clinical fluence levels. Results: On the order of only 10–100 gamma photons per flood image were required to be detected to avoid biasing the NPS estimate. This allowed for a factor of 10{sup 7} reduction in fluence compared to clinical levels with no loss of accuracy. An optimal signal-to-noise ratio (SNR) was achieved by increasing the number of flood images from a typical value of 100 up to 500, thereby illustrating the importance of flood image quantity over the number of gammas per flood. For the point-spread ensemble technique, an additional 2× reduction in the number of incident gammas was realized. As a result, when modeling gamma transport in a thick pixelated array, the simulation time was reduced from 2.5 × 10{sup 6} CPU min if using clinical fluence levels to 3.1 CPU min if using optimized fluence levels while also producing a higher SNR. The AS1000 DQE(f) simulation entailing both optical and radiative transport matched experimental results to within 11%, and required 14.5 min to complete on a single CPU. Conclusions: The authors demonstrate the feasibility of accurately modeling x-ray detector DQE(f) with completion times on the order of several minutes using a single CPU. Convenience of simulation can be achieved using GEANT4 which offers both gamma and optical photon transport capabilities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-J-38: Comparison of 6MV Photon Dose in a Perpendicular and Parallel Magnetic Field

Purpose: Integrating a linac with an MRI system would allow for real time tumour tracking however the patient will be irradiated in the presence of a magnetic field. The present study experimentally investigates the magnetic field effects on entrance, exit, and interface dose for both transverse and parallel magnetic fields. Methods: Polystyrene was used to construct a set of phantoms for Gafchromic film measurements. One phantom had an adjustable air gap and four other phantoms had one surface at various angles. The linac-MR prototype consisting of a biplanar permanent magnet coupled to a linac was used for the transverse magnetic field measurements. A couple of solenoid electromagnets, stacked on top of each other and irradiated along their bore, were used for the parallel field measurements. Results: All doses are relative to no magnetic field. The transverse magnetic field reduced the entrance dose for all surface angles by strongly deflecting the contaminant electrons. The exit dose in a transverse magnetic field was found to be significantly higher. The entrance dose with a parallel magnetic field present is higher due to the contaminant electrons being concentrated within the beam area. The air gap phantom measurements, done in a transverse magnetic field, show a significant increase of the dose at the proximal side of the air gap and a decrease at the distal side. The measurements, done in the parallel magnetic field, show the concentration of secondary electrons in the air gap. Conclusion: The radiation dose measurements of a 6MV beam in a parallel and transverse magnetic field presented here are currently being replicated using Monte Carlo simulations. This verified Monte Carlo system could provide the dose calculation basis for future linac-MR systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Experimental depth dose curves of a 67.5 MeV proton beam for benchmarking and validation of Monte Carlo simulation

Purpose: To measure depth dose curves for a 67.5 ± 0.1 MeV proton beam for benchmarking and validation of Monte Carlo simulation. Methods: Depth dose curves were measured in 2 beam lines. Protons in the raw beam line traversed a Ta scattering foil, 0.1016 or 0.381 mm thick, a secondary emission monitor comprised of thin Al foils, and a thin Kapton exit window. The beam energy and peak width and the composition and density of material traversed by the beam were known with sufficient accuracy to permit benchmark quality measurements. Diodes for charged particle dosimetry from two different manufacturers were used to scan the depth dose curves with 0.003 mm depth reproducibility in a water tank placed 300 mm from the exit window. Depth in water was determined with an uncertainty of 0.15 mm, including the uncertainty in the water equivalent depth of the sensitive volume of the detector. Parallel-plate chambers were used to verify the accuracy of the shape of the Bragg peak and the peak-to-plateau ratio measured with the diodes. The uncertainty in the measured peak-to-plateau ratio was 4%. Depth dose curves were also measured with a diode for a Bragg curve and treatment beam spread out Bragg peak (SOBP) on the beam line used for eye treatment. The measurements were compared to Monte Carlo simulation done with GEANT4 using TOPAS. Results: The 80% dose at the distal side of the Bragg peak for the thinner foil was at 37.47 ± 0.11 mm (average of measurement with diodes from two different manufacturers), compared to the simulated value of 37.20 mm. The 80% dose for the thicker foil was at 35.08 ± 0.15 mm, compared to the simulated value of 34.90 mm. The measured peak-to-plateau ratio was within one standard deviation experimental uncertainty of the simulated result for the thinnest foil and two standard deviations for the thickest foil. It was necessary to include the collimation in the simulation, which had a more pronounced effect on the peak-to-plateau ratio for the thicker foil. The treatment beam, being unfocussed, had a broader Bragg peak than the raw beam. A 1.3 ± 0.1 MeV FWHM peak width in the energy distribution was used in the simulation to match the Bragg peak width. An additional 1.3–2.24 mm of water in the water column was required over the nominal values to match the measured depth penetration. Conclusions: The proton Bragg curve measured for the 0.1016 mm thick Ta foil provided the most accurate benchmark, having a low contribution of proton scatter from upstream of the water tank. The accuracy was 0.15% in measured beam energy and 0.3% in measured depth penetration at the Bragg peak. The depth of the distal edge of the Bragg peak in the simulation fell short of measurement, suggesting that the mean ionization potential of water is 2–5 eV higher than the 78 eV used in the stopping power calculation for the simulation. The eye treatment beam line depth dose curves provide validation of Monte Carlo simulation of a Bragg curve and SOBP with 4%/2 mm accuracy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-H-BRA-04: Biological Geometries for the Monte Carlo Simulation Toolkit TOPASNBio

Purpose: New advances in radiation therapy are most likely to come from the complex interface of physics, chemistry and biology. Computational simulations offer a powerful tool for quantitatively investigating radiation interactions with biological tissue and can thus help bridge the gap between physics and biology. The aim of TOPAS-nBio is to provide a comprehensive tool to generate advanced radiobiology simulations. Methods: TOPAS wraps and extends the Geant4 Monte Carlo (MC) simulation toolkit. TOPAS-nBio is an extension to TOPAS which utilizes the physics processes in Geant4-DNA to model biological damage from very low energy secondary electrons. Specialized cell, organelle and molecular geometries were designed for the toolkit. Results: TOPAS-nBio gives the user the capability of simulating biological geometries, ranging from the micron-scale (e.g. cells and organelles) to complex nano-scale geometries (e.g. DNA and proteins). The user interacts with TOPAS-nBio through easy-to-use input parameter files. For example, in a simple cell simulation the user can specify the cell type and size as well as the type, number and size of included organelles. For more detailed nuclear simulations, the user can specify chromosome territories containing chromatin fiber loops, the later comprised of nucleosomes on a double helix. The chromatin fibers can be arranged in simple rigid geometries or within factual globules, mimicking realistic chromosome territories. TOPAS-nBio also provides users with the capability of reading protein data bank 3D structural files to simulate radiation damage to proteins or nucleic acids e.g. histones or RNA. TOPAS-nBio has been validated by comparing results to other track structure simulation software and published experimental measurements. Conclusion: TOPAS-nBio provides users with a comprehensive MC simulation tool for radiobiological simulations, giving users without advanced programming skills the ability to design and run complex simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Ring-shaped spatial pattern of exciton luminescence formed due to the hot carrier transport in a locally photoexcited electron-hole bilayer

A consistent explanation of the formation of a ring-shaped pattern of exciton luminescence in GaAs/AlGaAs double quantum wells is suggested. The pattern consists of two concentric rings around the laser excitation spot. It is shown that the luminescence rings appear due to the in-layer transport of hot charge carriers at high photoexcitation intensity. Interestingly, one of two causes of this transport might involve self-organized criticality (SOC) that would be the first case of the SOC observation in semiconductor physics. We test this cause in a many-body numerical model by performing extensive molecular dynamics simulations. The results show good agreement with experiments. Moreover, the simulations have enabled us to identify the particular kinetic processes underlying the formation of each of these two luminescence rings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

X-ray optics simulation using Gaussian superposition technique

We present an efficient method to perform x-ray optics simulation with high or partially coherent x-ray sources using Gaussian superposition technique. In a previous paper, we have demonstrated that full characterization of optical systems, diffractive and geometric, is possible by using the Fresnel Gaussian Shape Invariant (FGSI) previously reported in the literature. The complex amplitude distribution in the object plane is represented by a linear superposition of complex Gaussians wavelets and then propagated through the optical system by means of the referred Gaussian invariant. This allows ray tracing through the optical system and at the same time allows calculating with high precision the complex wave-amplitude distribution at any plane of observation. This technique can be applied in a wide spectral range where the Fresnel diffraction integral applies including visible, x-rays, acoustic waves, etc. We describe the technique and include some computer simulations as illustrative examples for x-ray optical component. We show also that this method can be used to study partial or total coherence illumination problem.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Dense Plasma Focus Modeling

The mechanisms for pinch formation in Dense Plasma Focus (DPF) devices, with the generation of high-energy ions beams and subsequent neutron production over a relatively short distance, are not fully understood. Here we report on high-fidelity 2D and 3D numerical magnetohydrodynamic (MHD) simulations using the LA-COMPASS code to study the pinch formation dynamics and its associated instabilities and neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

User Guidelines and Best Practices for CASL VUQ Analysis Using Dakota.

Sandia's Dakota software (available at http://dakota.sandia.gov) supports science and engineering transformation through advanced exploration of simulations. Specifically, it manages and analyzes ensembles of simulations to provide broader and deeper perspective for analysts and decision makers. This enables them to enhance understanding of risk, improve products, and assess simulation credibility. This manual offers Consortium for Advanced Simulation of Light Water Reactors (LWRs) (CASL) partners a guide to conducting Dakota-based VUQ studies for CASL problems. It motivates various classes of Dakota methods and includes examples of their use on representative application problems. On reading, a CASL analyst should understand why and how to apply Dakota to a simulation problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Nanoindentation of ZrN Coatings on Silicon and Copper Substrates

An effective mathematical model is proposed to describe experiments on the indentation of samples with layered or functionally graded coatings. The model is based on the elastic contact problem corresponding to the penetration of a punch in a coated elastic half-space. A Berkovich indenter is considered. The results of the simulation and the experiments are compared for the indentation of two zirconium nitride (ZrN) coatings on substrates of monocrystalline (100) silicon and annealed copper. The agreement of the results is good.

42 ENGINEERING↗
DOE OSTI2017

Charged Particle Diffusion in Isotropic Random Magnetic Fields

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider the diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of Larmor radius to the characteristic outer length scale of turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared to those from computer simulations, showing good agreement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Evolution of Mass and Velocity Field in the Cosmic Web: Comparison Between Baryonic and Dark Matter

We investigate the evolution of the cosmic web since z = 5 in grid-based cosmological hydrodynamical simulations, focusing on the mass and velocity fields of both baryonic and cold dark matter. The tidal tensor of density is used as the main method for web identification, with λ {sub th} = 0.2–1.2. The evolution trends in baryonic and dark matter are similar, although moderate differences are observed. Sheets appear early, and their large-scale pattern may have been set up by z = 3. In terms of mass, filaments supersede sheets as the primary collapsing structures from z ∼ 2–3. Tenuous filaments assembled with each other to form prominent ones at z < 2. In accordance with the construction of the frame of the sheets, the cosmic divergence velocity, v {sub div}, was already well-developed above 2–3 Mpc by z = 3. Afterwards, the curl velocity, v {sub curl}, grew dramatically along with the rising of filaments, becoming comparable to v {sub div}, for <2–3 Mpc at z = 0. The scaling of v {sub curl} can be described by the hierarchical turbulence model. The alignment between the vorticity and the eigenvectors of the shear tensor in the baryonic matter field resembles that in the dark matter field, and is even moderately stronger between ω and e{sub 1}, and ω and e{sub 3}. Compared with dark matter, there is slightly less baryonic matter found residing in filaments and clusters, and its vorticity developed more significantly below 2–3 Mpc. These differences may be underestimated because of the limited resolution and lack of star formation in our simulation. The impact of the change of dominant structures in overdense regions at z ∼ 2–3 on galaxy formation and evolution is shortly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗