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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Magnetic field line random walk in models and simulations of reduced magnetohydrodynamic turbulence

The random walk of magnetic field lines is examined numerically and analytically in the context of reduced magnetohydrodynamic (RMHD) turbulence, which provides a useful description of plasmas dominated by a strong mean field, such as in the solar corona. A recently developed non-perturbative theory of magnetic field line diffusion is compared with the diffusion coefficients obtained by accurate numerical tracing of magnetic field lines for both synthetic models and direct numerical simulations of RMHD. Statistical analysis of an ensemble of trajectories confirms the applicability of the theory, which very closely matches the numerical field line diffusion coefficient as a function of distance z along the mean magnetic field for a wide range of the Kubo number R. This theory employs Corrsin's independence hypothesis, sometimes thought to be valid only at low R. However, the results demonstrate that it works well up to R = 10, both for a synthetic RMHD model and an RMHD simulation. The numerical results from the RMHD simulation are compared with and without phase randomization, demonstrating a clear effect of coherent structures on the field line random walk for a very low Kubo number.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Comprehensive rate coefficients for electron-collision-induced transitions in hydrogen

Energy-changing electron-hydrogen atom collisions are crucial to regulating the energy balance in astrophysical and laboratory plasmas and are relevant to the formation of stellar atmospheres, recombination in H II clouds, primordial recombination, three-body recombination, and heating in ultracold and fusion plasmas. Computational modeling of electron-hydrogen collision has been attempted through quantum mechanical scattering state-to-state calculations of transitions involving low-lying energy levels in hydrogen (with principal quantum number n < 7) and at large principal quantum numbers using classical trajectory techniques. Analytical expressions are proposed that interpolate the current quantum mechanical and classical trajectory results for electron-hydrogen scattering in the entire range of energy levels for nearly the entire temperature range of interest in astrophysical environments. An asymptotic expression for the Born cross section is interpolated with a modified expression previously derived for electron-hydrogen scattering in the Rydberg regime using classical trajectory Monte Carlo simulations. The derived formula is compared to existing numerical data for transitions involving low principal quantum numbers, and the dependence of the deviations on temperature is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Hubble space telescope/advanced camera for surveys confirmation of the dark substructure in A520

We present a weak-lensing study of the cluster A520 based on Advanced Camera for Surveys (ACS) data. The excellent data quality provides a mean source density of ∼109 arcmin{sup –2}, which improves both resolution and significance of the mass reconstruction compared to a previous study based on Wide Field Planetary Camera 2 (WFPC2) images. We take care in removing instrumental effects such as the charge trailing due to radiation damage of the detector and the position-dependent point-spread function. This new ACS analysis confirms the previous claims that a substantial amount of dark mass is present between two luminous subclusters where we observe very little light. The centroid of the dark peak in the current ACS analysis is offset to the southwest by ∼1' with respect to the centroid from the WFPC2 analysis. Interestingly, this new centroid is in better agreement with the location where the X-ray emission is strongest, and the mass-to-light ratio estimated with this centroid is much higher (813 ± 78 M {sub ☉}/L {sub R☉}) than the previous value; the aperture mass with the WFPC2 centroid provides a consistent mass. Although we cannot provide a definite explanation for the dark peak, we discuss a revised scenario, wherein dark matter with a more conventional range (σ{sub DM}/m {sub DM} < 1 cm{sup 2} g{sup –1}) of self-interacting cross-section can lead to the detection of this dark substructure. If supported by detailed numerical simulations, this hypothesis opens up the possibility that the A520 system can be used to establish a lower limit of the self-interacting cross-section of dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE LAUNCHING OF COLD CLOUDS BY GALAXY OUTFLOWS. I. HYDRODYNAMIC INTERACTIONS WITH RADIATIVE COOLING

To better understand the nature of the multiphase material found in outflowing galaxies, we study the evolution of cold clouds embedded in flows of hot and fast material. Using a suite of adaptive mesh refinement simulations that include radiative cooling, we investigate both cloud mass loss and cloud acceleration under the full range of conditions observed in galaxy outflows. The simulations are designed to track the cloud center of mass, enabling us to study the cloud evolution at long disruption times. For supersonic flows, a Mach cone forms around the cloud, which damps the Kelvin–Helmholtz instability but also establishes a streamwise pressure gradient that stretches the cloud apart. If time is expressed in units of the cloud crushing time, both the cloud lifetime and the cloud acceleration rate are independent of cloud radius, and we find simple scalings for these quantities as a function of the Mach number of the external medium. A resolution study suggests that our simulations accurately describe the evolution of cold clouds in the absence of thermal conduction and magnetic fields, physical processes whose roles will be studied in forthcoming papers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE KOZAI–LIDOV MECHANISM IN HYDRODYNAMICAL DISKS. III. EFFECTS OF DISK MASS AND SELF-GRAVITY

Previously we showed that a substantially misaligned viscous accretion disk with pressure that orbits around one component of a binary system can undergo global damped Kozai–Lidov (KL) oscillations. These oscillations produce periodic exchanges of the disk eccentricity with inclination. The disk KL mechanism is quite robust and operates over a wide range of binary and disk parameters. However, the effects of self-gravity, which are expected to suppress the KL oscillations for sufficiently massive disks, were ignored. Here, we analyze the effects of disk self-gravity by means of hydrodynamic simulations and compare the results with the expectations of analytic theory. The disk mass required for suppression in the simulations is a few percent of the mass of the central star and this roughly agrees with an analytical estimate. The conditions for suppression of the KL oscillations in the simulations are close to requiring that the disk be gravitationally unstable. We discuss some implications of our results for the dynamics of protoplanetary disks and the related planet formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A Hubble diagram for quasars

We present a new method to test the ΛCDM cosmological model and to estimate cosmological parameters based on the nonlinear relation between the ultraviolet and X-ray luminosities of quasars. We built a data set of 1138 quasars by merging several samples from the literature with X-ray measurements at 2 keV and SDSS photometry, which was used to estimate the extinction-corrected 2500 Å flux. We obtained three main results: (1) we checked the nonlinear relation between X-ray and UV luminosities in small redshift bins up to z∼6, confirming that the relation holds at all redshifts with the same slope; (2) we built a Hubble diagram for quasars up to z∼6, which is well matched to that of supernovae in the common z = 0–1.4 redshift interval and extends the test of the cosmological model up to z∼6; and (3) we showed that this nonlinear relation is a powerful tool for estimating cosmological parameters. Using the present data and assuming a ΛCDM model, we obtain Ω{sub M} = 0.22{sub −0.08}{sup +0.10} and Ω{sub Λ} = 0.92{sub −0.30}{sup +0.18} (Ω{sub M} = 0.28 ± 0.04 and Ω{sub Λ} = 0.73 ± 0.08 from a joint quasar-SNe fit). Much more precise measurements will be achieved with future surveys. A few thousand SDSS quasars already have serendipitous X-ray observations from Chandra or XMM-Newton, and at least 100,000 quasars with UV and X-ray data will be made available by the extended ROentgen Survey with an Imaging Telescope Array all-sky survey in a few years. The Euclid, Large Synoptic Survey Telescope, and Advanced Telescope for High ENergy Astrophysics surveys will further increase the sample size to at least several hundred thousand. Our simulations show that these samples will provide tight constraints on the cosmological parameters and will allow us to test for possible deviations from the standard model with higher precision than is possible today.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

GALAXY SPIN ALIGNMENT IN FILAMENTS AND SHEETS: OBSERVATIONAL EVIDENCE

The properties of galaxies are known to be affected by their environment. One important question is how their angular momentum reflects the surrounding cosmic web. We use the Sloan Digital Sky Survey to investigate the spin axes of spiral and elliptical galaxies relative to their surrounding filament/sheet orientations. To detect filaments, a marked point process with interactions (the {sup B}isous model{sup )} is used. Sheets are found by detecting 'flattened' filaments. The minor axes of ellipticals are found to be preferentially perpendicular to hosting filaments. A weak correlation is found with sheets. These findings are consistent with the notion that elliptical galaxies formed via mergers, which predominantly occurred along the filaments. The spin axis of spiral galaxies is found to align with the host filament, with no correlation between spiral spin and sheet normal. When examined as a function of distance from the filament axis, a much stronger correlation is found in the outer parts, suggesting that the alignment is driven by the laminar infall of gas from sheets to filaments. When compared with numerical simulations, our results suggest that the connection between dark matter halo and galaxy spin is not straightforward. Our results provide an important input to the understanding of how galaxies acquire their angular momentum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-T-304: Study of Secondary Neutrons From Uniform Scanning Proton Beams

Purpose: Secondary neutrons are unwanted byproducts from proton therapy and exposure from secondary radiation during treatment could increase risk of developing a secondary cancer later in a patient's lifetime. The purpose of this study is to investigate secondary neutrons from uniform scanning proton beams under various beam conditions using both measurements and Monte Carlo simulations. Methods: CR-39 Plastic Track Nuclear Detectors (PNTD) were used for the measurement. CR-39 PNTD has tissue like sensitivity to the secondary neutrons but insensitive to the therapeutic protons. In this study, we devised two experimental conditions: a) hollow-phantom; phantom is bored with a hollow cylinder along the direction of the beam so that the primary proton passes through the phantom without interacting with the phantom material, b) cylindrical-phantom; a solid cylinder of diameter close to the beam diameter is placed along the beam path. CR-39 PNTDs were placed laterally inside a 60X20X35 cm3 phantom (hollow-phantom) and in air (cylindrical-phantom) at various angles with respect to the primary beam axis. We studied for three different proton energies (78 MeV, 162 MeV and 226 MeV), using a 4 cm modulation width and 5cm diameter brass aperture for the entire experiment and simulation. A comparison of the experiment was performed using the Monte Carlo code FLUKA. Results: The measured secondary neutron dose equivalent per therapeutic primary proton dose (H/D) ranges from 2.1 ± 0.2 to 25.42 ± 2.3 mSv/Gy for the hollow phantom study, and 2.7 ± 0.3 to 46.4 ± 3.4 mSv/Gy for the cylindrical phantom study. Monte Carlo simulations predicated neutron dose equivalent from measurements within a factor of 5. Conclusion: The study suggests that the production of external neutrons is significantly higher than the production of internal neutrons.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

WE-E-17A-07: Patient-Specific Mathematical Neuro-Oncology: Biologically-Informed Radiation Therapy and Imaging Physics

Purpose: To demonstrate a method of generating patient-specific, biologically-guided radiation therapy (RT) plans and to quantify and predict response to RT in glioblastoma. We investigate the biological correlates and imaging physics driving T2-MRI based response to radiation therapy using an MRI simulator. Methods: We have integrated a patient-specific biomathematical model of glioblastoma proliferation, invasion and radiotherapy with a multiobjective evolutionary algorithm for intensity-modulated RT optimization to construct individualized, biologically-guided plans. Patient-individualized simulations of the standard-of-care and optimized plans are compared in terms of several biological metrics quantified on MRI. An extension of the PI model is used to investigate the role of angiogenesis and its correlates in glioma response to therapy with the Proliferation-Invasion-Hypoxia- Necrosis-Angiogenesis model (PIHNA). The PIHNA model is used with a brain tissue phantom to predict tumor-induced vasogenic edema, tumor and tissue density that is used in a multi-compartmental MRI signal equation for generation of simulated T2- weighted MRIs. Results: Applying a novel metric of treatment response (Days Gained) to the patient-individualized simulation results predicted that the optimized RT plans would have a significant impact on delaying tumor progression, with Days Gained increases from 21% to 105%. For the T2- MRI simulations, initial validation tests compared average simulated T2 values for white matter, tumor, and peripheral edema to values cited in the literature. Simulated results closely match the characteristic T2 value for each tissue. Conclusion: Patient-individualized simulations using the combination of a biomathematical model with an optimization algorithm for RT generated biologically-guided doses that decreased normal tissue dose and increased therapeutic ratio with the potential to improve survival outcomes for treatment of glioblastoma. Simulated T2-MRI is shown to be consistent with known physics of MRI and can be used to further investigate biological drivers of imaging-based response to RT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-242: A Method for Collision Avoidance in External Beam Radiation Therapy

Purpose: We proposed a method for collision avoidance (CA) in external beam radiation therapy (EBRT). The method encompasses the analysis of all positions of the moving components of the beam delivery system such as the treatment table and gantry, including patient specific information obtained from the CT images. This method eliminates the need for time-consuming dry-runs prior to the actual treatments. Methods: The QA procedure for EBRT requires that the collision should be checked prior to treatment. We developed a system capable of a rigorous computer simulation of all moving components including positions of the couch and gantry during the delivery, position of the patients, and imaging equipment. By running this treatment simulation it is possible to quantify and graphically represent all positions and corresponding trajectories of all points of the moving parts during the treatment delivery. The development of the workflow for implementation of the CA includes several steps: a) derivation of combined dynamic equation of motion of the EBRT delivery systems, b) developing the simulation model capable of drawing the motion trajectories of the specific points, c) developing the interface between the model and the treatment plan parameters such as couch and gantry parameters for each field. Results: The patient CT images were registered to the treatment couch so the patient dimensions were included into the simulation. The treatment field parameters were structured in the xml-file which was used as the input into the dynamic equations. The trajectories of the moving components were plotted on the same graph using the dynamic equations. If the trajectories intersect that was the signal that collision exists. Conclusion: This CA method was proved to be effective in the simulation of treatment delivery. The proper implementation of this system can potentially improve the QA program and increase the efficacy in the clinical setup.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Thermal fluctuations of pair interaction forces in liquid Yukawa systems

The dynamics of charged particles in two- and three-dimensional Yukawa systems is studied using numerical methods. Nonideal systems are simulated in a wide range of their parameters. The density of thermal fluctuations of pair interaction forces in liquid structures is obtained and compared with the internal energy density of the systems under study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Evaluation of DUSTRAN Software System for Modeling Chloride Deposition on Steel Canisters

The degradation of steel by stress corrosion cracking (SCC) when exposed to atmospheric conditions for decades is a significant challenge in the fossil fuel and nuclear industries. SCC can occur when corrosive contaminants such as chlorides are deposited on a susceptible material in a tensile stress state. The Nuclear Regulatory Commission has identified chloride-induced SCC as a potential cause for concern in stainless steel used nuclear fuel (UNF) canisters in dry storage. The modeling of contaminant deposition is the first step in predictive multiscale modeling of SCC that is essential to develop mitigation strategies, prioritize inspection, and ensure the integrity and performance of canisters, pipelines, and structural materials. A multiscale simulation approach can be developed to determine the likelihood that a canister would undergo SCC in a certain period of time. This study investigates the potential of DUSTRAN, a dust dispersion modeling system developed by Pacific Northwest National Laboratory, to model the deposition of chloride contaminants from sea salt aerosols on a steel canister. Results from DUSTRAN simulations run with historical meteorological data were compared against measured chloride data at a coastal site in Maine. DUSTRAN’s CALPUFF model tended to simulate concentrations higher than those measured; however, the closest estimations were within the same order of magnitude as the measured values. The decrease in discrepancies between measured and simulated values as the level of abstraction in wind speed decreased suggest that the model is very sensitive to wind speed. However, the influence of other parameters such as the distinction between open-ocean and surf-zone sources needs to be explored further. Deposition values predicted by the DUSTRAN system were not in agreement with concentration values and suggest that the deposition calculations may not fully represent physical processes. Overall, results indicate that with parameter refinement, DUSTRAN has the potential to simulate atmospheric chloride dispersion on steel canisters.

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

Community Project for Accelerator Science and Simulation (ComPASS)

After concluding our initial exercise (solving a simplified statistical inverse problem with unknown parameter laser intensity) of coupling Vorpal and our parallel statistical library QUESO, we shifted the application focus to DLA. Our efforts focused on developing a Gaussian process (GP) emulator within QUESO for efficient optimization of power couplers within woodpiles. The smaller simulation size (compared with LPA) allows for sufficient “training runs” to develop a reasonable GP statistical emulator for a parameter space of moderate dimension.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Evaluation of CFD Methods for Simulation of Two-Phase Boiling Flow Phenomena in a Helical Coil Steam Generator

The U.S. Department of Energy, Office of Nuclear Energy charges participants in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program with the development of advanced modeling and simulation capabilities that can be used to address design, performance and safety challenges in the development and deployment of advanced reactor technology. The NEAMS has established a high impact problem (HIP) team to demonstrate the applicability of these tools to identification and mitigation of sources of steam generator flow induced vibration (SGFIV). The SGFIV HIP team is working to evaluate vibration sources in an advanced helical coil steam generator using computational fluid dynamics (CFD) simulations of the turbulent primary coolant flow over the outside of the tubes and CFD simulations of the turbulent multiphase boiling secondary coolant flow inside the tubes integrated with high resolution finite element method assessments of the tubes and their associated structural supports. This report summarizes the demonstration of a methodology for the multiphase boiling flow analysis inside the helical coil steam generator tube. A helical coil steam generator configuration has been defined based on the experiments completed by Polytecnico di Milano in the SIET helical coil steam generator tube facility. Simulations of the defined problem have been completed using the Eulerian-Eulerian multi-fluid modeling capabilities of the commercial CFD code STAR-CCM+. Simulations suggest that the two phases will quickly stratify in the slightly inclined pipe of the helical coil steam generator. These results have been successfully benchmarked against both empirical correlations for pressure drop and simulations using an alternate CFD methodology, the dispersed phase mixture modeling capabilities of the open source CFD code Nek5000.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Variations of the ISM compactness across the main sequence of star forming galaxies: Observations and simulations

The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass (M{sub *}) plane, of the form SFR∝M{sub ∗}{sup α}, usually referred to as the star formation main sequence (MS). The physics that sets the properties of the MS is currently a subject of debate, and no consensus has been reached regarding the fundamental difference between members of the sequence and its outliers. Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present Chiburst, a Markov Chain Monte Carlo spectral energy distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, SFRs, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate, and the compactness parameter C, that parametrizes this geometry and hence the evolution of dust temperature (T{sub dust}) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of luminous infrared galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Fast radio bursts as probes of magnetic fields in the intergalactic medium

We examine the proposal that the dispersion measures (DMs) and Faraday rotation measures (RMs) of extragalactic linearly polarized fast radio bursts (FRBs) can be used to probe the intergalactic magnetic field (IGMF) in filaments of galaxies. The DM through the cosmic web is dominated by contributions from the warm-hot intergalactic medium (WHIM) in filaments and from the gas in voids. On the other hand, RM is induced mostly by the hot medium in galaxy clusters, and only a fraction of it is produced in the WHIM. We show that if one excludes FRBs whose sightlines pass through galaxy clusters, the line of sight (LOS) strength of the IGMF in filaments, B{sub ||}, is approximately C(〈1+z〉/f{sub DM})(RM/DM), where C is a known constant. Here, the redshift of the FRB is not required to be known; f {sub DM} is the fraction of total DM due to the WHIM, while 〈1+z〉 is the redshift of interevening gas weighted by the WHIM gas density, both of which can be evaluated for a given cosmology model solely from the DM of an FRB. Using data on structure formation simulations and a model IGMF, we show that C(〈1+z〉/f{sub DM})(RM/DM) closely reproduces the density-weighted LOS strength of the IGMF in filaments of the large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A NUMERICAL SIMULATION OF COSMIC RAY MODULATION NEAR THE HELIOPAUSE. II. SOME PHYSICAL INSIGHTS

Cosmic ray (CR) transport near the heliopause (HP) is studied using a hybrid transport model, with the parameters constrained by observations from the Voyager 1 spacecraft. We simulate the CR radial flux along different directions in the heliosphere. There is no well-defined thin layer between the solar wind region and the interstellar region along the tail and polar directions of the heliosphere. By analyzing the radial flux curve along the direction of Voyager 2 , together with its trajectory information, the crossing time of the HP by Voyager 2 is predicted to be in 2017.14. We simulate the CR radial flux for different energy values along the direction of Voyager 1 . We find that there is only a modest modulation region of about 10 au wide beyond the HP, so that Voyager 1 observing the Local Interstellar Spectra is justified in numerical modeling. We analyze the heliospheric exit information of pseudo-particles in our stochastic numerical (time-backward) method, conjecturing that they represent the behavior of CR particles, and we find that pseudo-particles that have been traced from the nose region exit in the tail region. This implies that many CR particles diffuse directly from the heliospheric tail region to the nose region near the HP. In addition, when pseudo-particles were traced from the Local Interstellar Medium (LISM), it is found that their exit location (entrance for real particles) from the simulation domain is along the prescribed Interstellar Magnetic Field direction. This indicates that parallel diffusion dominates CR particle transport in the LISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

AGE dependence of wind properties for solar-type stars: a 3D study

Young and rapidly rotating stars are known for intense, dynamo-generated magnetic fields. Spectropolarimetric observations of those stars in precisely aged clusters are key input for gyrochronology and magnetochronology. We use Zeeman Doppler imaging maps of several young K-type stars of similar mass and radius but with various ages and rotational periods to perform three-dimensional (3D) numerical MHD simulations of their coronae and follow the evolution of their magnetic properties with age. Those simulations yield the coronal structure as well as the instant torque exerted by the magnetized, rotating wind on the star. As stars get older, we find that the angular momentum loss decreases with Ω{sub ⋆}{sup 3}, which is the reason for the convergence on the Skumanich law. For the youngest stars of our sample, the angular momentum loss shows signs of saturation around 8Ω{sub ⊙}, which is a common value used in spin evolution models for K-type stars. We compare these results to semianalytical models and existing braking laws. We observe a complex wind-speed distribution for the youngest stars with slow, intermediate, and fast wind components, which are the result of interaction with intense and nonaxisymmetric magnetic fields. Consequently, in our simulations, the stellar wind structure in the equatorial plane of young stars varies significantly from a solar configuration, delivering insight about the past of the solar system interplanetary medium.

79 ASTRONOMY AND ASTROPHYSICS↗