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

Ellerman bombs—evidence for magnetic reconnection in the lower solar atmosphere

The presence of photospheric magnetic reconnection has long been thought to give rise to short and impulsive events, such as Ellerman bombs (EBs) and Type II spicules. In this article, we combine high-resolution, high-cadence observations from the Interferometric BIdimensional Spectrometer and Rapid Oscillations in the Solar Atmosphere instruments at the Dunn Solar Telescope, National Solar Observatory, New Mexico, with co-aligned Solar Dynamics Observatory Atmospheric Imaging Assembly and Hinode Solar Optical Telescope (SOT) data to observe small-scale events situated within an active region. These data are then compared with state-of-the-art numerical simulations of the lower atmosphere made using the MURaM code. It is found that brightenings, in both the observations and the simulations, of the wings of the Hα line profile, interpreted as EBs, are often spatially correlated with increases in the intensity of the Fe I λ6302.5 line core. Bipolar regions inferred from Hinode/SOT magnetic field data show evidence of flux cancellation associated, co-spatially, with these EBs, suggesting that magnetic reconnection could be a driver of these high-energy events. Through the analysis of similar events in the simulated lower atmosphere, we are able to infer that line profiles analogous to the observations occur co-spatially with regions of strong opposite-polarity magnetic flux. These observed events and their simulated counterparts are interpreted as evidence of photospheric magnetic reconnection at scales observable using current observational instrumentation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Magnetic field amplification and evolution in turbulent collisionless magnetohydrodynamics: An application to the intracluster medium

The amplification of magnetic fields (MFs) in the intracluster medium (ICM) is attributed to turbulent dynamo (TD) action, which is generally derived in the collisional-MHD framework. However, this assumption is poorly justified a priori, since in the ICM the ion mean free path between collisions is of the order of the dynamical scales, thus requiring a collisionless MHD description. The present study uses an anisotropic plasma pressure that brings the plasma within a parametric space where collisionless instabilities take place. In this model, a relaxation term of the pressure anisotropy simulates the feedback of the mirror and firehose instabilities, in consistency with empirical studies. Our three-dimensional numerical simulations of forced transonic turbulence, aiming the modeling of the turbulent ICM, were performed for different initial values of the MF intensity and different relaxation rates of the pressure anisotropy. We found that in the high-β plasma regime corresponding to the ICM conditions, a fast anisotropy relaxation rate gives results that are similar to the collisional-MHD model, as far as the statistical properties of the turbulence are concerned. Also, the TD amplification of seed MFs was found to be similar to the collisional-MHD model. The simulations that do not employ the anisotropy relaxation deviate significantly from the collisional-MHD results and show more power at the small-scale fluctuations of both density and velocity as a result of the action of the instabilities. For these simulations, the large-scale fluctuations in the MF are mostly suppressed and the TD fails in amplifying seed MFs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Detection of coherent structures in photospheric turbulent flows

We study coherent structures in solar photospheric flows in a plage in the vicinity of the active region AR 10930 using the horizontal velocity data derived from Hinode/Solar Optical Telescope magnetograms. Eulerian and Lagrangian coherent structures (LCSs) are detected by computing the Q-criterion and the finite-time Lyapunov exponents of the velocity field, respectively. Our analysis indicates that, on average, the deformation Eulerian coherent structures dominate over the vortical Eulerian coherent structures in the plage region. We demonstrate the correspondence of the network of high magnetic flux concentration to the attracting Lagrangian coherent structures (aLCSs) in the photospheric velocity based on both observations and numerical simulations. In addition, the computation of aLCS provides a measure of the local rate of contraction/expansion of the flow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

NUMERICAL STUDY ON IN SITU PROMINENCE FORMATION BY RADIATIVE CONDENSATION IN THE SOLAR CORONA

We propose an in situ formation model for inverse-polarity solar prominences and demonstrate it using self-consistent 2.5 dimensional MHD simulations, including thermal conduction along magnetic fields and optically thin radiative cooling. The model enables us to form cool dense plasma clouds inside a flux rope by radiative condensation, which is regarded as an inverse-polarity prominence. Radiative condensation is triggered by changes in the magnetic topology, i.e., formation of the flux rope from the sheared arcade field, and by thermal imbalance due to the dense plasma trapped inside the flux rope. The flux rope is created by imposing converging and shearing motion on the arcade field. Either when the footpoint motion is in the anti-shearing direction or when heating is proportional to local density, the thermal state inside the flux rope becomes cooling-dominant, leading to radiative condensation. By controlling the temperature of condensation, we investigate the relationship between the temperature and density of prominences and derive a scaling formula for this relationship. This formula suggests that the proposed model reproduces the observed density of prominences, which is 10–100 times larger than the coronal density. Moreover, the time evolution of the extreme ultraviolet emission synthesized by combining our simulation results with the response function of the Solar Dynamics Observatory Atmospheric Imaging Assembly filters agrees with the observed temporal and spatial intensity shift among multi-wavelength extreme ultraviolet emission during in situ condensation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SMALL-SCALE MAGNETIC ISLANDS IN THE SOLAR WIND AND THEIR ROLE IN PARTICLE ACCELERATION. I. DYNAMICS OF MAGNETIC ISLANDS NEAR THE HELIOSPHERIC CURRENT SHEET

Increases of ion fluxes in the keV–MeV range are sometimes observed near the heliospheric current sheet (HCS) during periods when other sources are absent. These resemble solar energetic particle events, but the events are weaker and apparently local. Conventional explanations based on either shock acceleration of charged particles or particle acceleration due to magnetic reconnection at interplanetary current sheets (CSs) are not persuasive. We suggest instead that recurrent magnetic reconnection occurs at the HCS and smaller CSs in the solar wind, a consequence of which is particle energization by the dynamically evolving secondary CSs and magnetic islands. The effectiveness of the trapping and acceleration process associated with magnetic islands depends in part on the topology of the HCS. We show that the HCS possesses ripples superimposed on the large-scale flat or wavy structure. We conjecture that the ripples can efficiently confine plasma and provide tokamak-like conditions that are favorable for the appearance of small-scale magnetic islands that merge and/or contract. Particles trapped in the vicinity of merging islands and experiencing multiple small-scale reconnection events are accelerated by the induced electric field and experience first-order Fermi acceleration in contracting magnetic islands according to the transport theory of Zank et al. We present multi-spacecraft observations of magnetic island merging and particle energization in the absence of other sources, providing support for theory and simulations that show particle energization by reconnection related processes of magnetic island merging and contraction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Distribution function approach to redshift space distortions. Part V: perturbation theory applied to dark matter halos

Numerical simulations show that redshift space distortions (RSD) introduce strong scale dependence in the power spectra of halos, with ten percent deviations relative to linear theory predictions even on relatively large scales (k < 0.1h/Mpc) and even in the absence of satellites (which induce Fingers-of-God, FoG, effects). If unmodeled these effects prevent one from extracting cosmological information from RSD surveys. In this paper we use Eulerian perturbation theory (PT) and Eulerian halo biasing model and apply it to the distribution function approach to RSD, in which RSD is decomposed into several correlators of density weighted velocity moments. We model each of these correlators using PT and compare the results to simulations over a wide range of halo masses and redshifts. We find that with an introduction of a physically motivated halo biasing, and using dark matter power spectra from simulations, we can reproduce the simulation results at a percent level on scales up to k ∼ 0.15h/Mpc at z = 0, without the need to have free FoG parameters in the model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Scaling relations of halo cores for self-interacting dark matter

Using a simple analytic formalism, we demonstrate that significant dark matter self-interactions produce halo cores that obey scaling relations nearly independent of the underlying particle physics parameters such as the annihilation cross section and the mass of the dark matter particle. For dwarf galaxies, we predict that the core density ρ{sub c} and the core radius r{sub c} should obey ρ{sub c} r{sub c} ≈ 41 M{sub ⊙} pc{sup −2} with a weak mass dependence ∼ M{sup 0.2}. Remarkably, such a scaling relation has recently been empirically inferred. Scaling relations involving core mass, core radius, and core velocity dispersion are predicted and agree well with observational data. By calibrating against numerical simulations, we predict the scatter in these relations and find them to be in excellent agreement with existing data. Future observations can test our predictions for different halo masses and redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

CAN AGN FEEDBACK BREAK THE SELF-SIMILARITY OF GALAXIES, GROUPS, AND CLUSTERS?

It is commonly thought that active galactic nucleus (AGN) feedback can break the self-similar scaling relations of galaxies, groups, and clusters. Using high-resolution three-dimensional hydrodynamic simulations, we isolate the impact of AGN feedback on the L {sub x}-T {sub x} relation, testing the two archetypal and common regimes, self-regulated mechanical feedback and a quasar thermal blast. We find that AGN feedback has severe difficulty in breaking the relation in a consistent way. The similarity breaking is directly linked to the gas evacuation within R {sub 500}, while the central cooling times are inversely proportional to the core density. Breaking self-similarity thus implies breaking the cool core, morphing all systems to non-cool-core objects, which is in clear contradiction with the observed data populated by several cool-core systems. Self-regulated feedback, which quenches cooling flows and preserves cool cores, prevents dramatic evacuation and similarity breaking at any scale; the relation scatter is also limited. The impulsive thermal blast can break the core-included L {sub x}-T {sub x} at T {sub 500} ≲ 1 keV, but substantially empties and overheats the halo, generating a perennial non-cool-core group, as experienced by cosmological simulations. Even with partial evacuation, massive systems remain overheated. We show that the action of purely AGN feedback is to lower the luminosity and heat the gas, perpendicular to the fit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

RANDOM WALKS AND EFFECTIVE OPTICAL DEPTH IN RELATIVISTIC FLOW

We investigate the random walk process in relativistic flow. In the relativistic flow, photon propagation is concentrated in the direction of the flow velocity due to the relativistic beaming effect. We show that in the pure scattering case, the number of scatterings is proportional to the size parameter ξ ≡ L/l {sub 0} if the flow velocity β ≡ v/c satisfies β/Γ >> ξ{sup –1}, while it is proportional to ξ{sup 2} if β/Γ << ξ{sup –1}, where L and l {sub 0} are the size of the system in the observer frame and the mean free path in the comoving frame, respectively. We also examine the photon propagation in the scattering and absorptive medium. We find that if the optical depth for absorption τ{sub a} is considerably smaller than the optical depth for scattering τ{sub s} (τ{sub a}/τ{sub s} << 1) and the flow velocity satisfies β≫√(2τ{sub a}/τ{sub s}), then the effective optical depth is approximated by τ{sub *} ≅ τ{sub a}(1 + β)/β. Furthermore, we perform Monte Carlo simulations of radiative transfer and compare the results with the analytic expression for the number of scatterings. The analytic expression is consistent with the results of the numerical simulations. The expression derived in this study can be used to estimate the photon production site in relativistic phenomena, e.g., gamma-ray burst and active galactic nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

MOSFIRE ABSORPTION LINE SPECTROSCOPY OF z > 2 QUIESCENT GALAXIES: PROBING A PERIOD OF RAPID SIZE GROWTH

Using the MOSFIRE near-infrared multi-slit spectrograph on the Keck 1 Telescope, we have secured high signal-to-noise ratio absorption line spectra for six massive galaxies with redshift 2 < z < 2.5. Five of these galaxies lie on the red sequence and show signatures of passive stellar populations in their rest-frame optical spectra. By fitting broadened spectral templates we have determined stellar velocity dispersions and, with broad-band Hubble Space Telescope and Spitzer photometry and imaging, stellar masses and effective radii. Using this enlarged sample of galaxies, we confirm earlier suggestions that quiescent galaxies at z > 2 have small sizes and large velocity dispersions compared to local galaxies of similar stellar mass. The dynamical masses are in very good agreement with stellar masses (log M {sub *}/M {sub dyn} = –0.02 ± 0.03), although the average stellar-to-dynamical mass ratio is larger than that found at lower redshift (–0.23 ± 0.05). By assuming evolution at fixed velocity dispersion, not only do we confirm a surprisingly rapid rate of size growth but we also consider the necessary evolutionary track on the mass-size plane and find a slope α = dlog R{sub e} /dlog M {sub *} ≳ 2 inconsistent with most numerical simulations of minor mergers. Both results suggest an additional mechanism may be required to explain the size growth of early galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Geometrical splitting technique to improve the computational efficiency in Monte Carlo calculations for proton therapy

Purpose: To present the implementation and validation of a geometrical based variance reduction technique for the calculation of phase space data for proton therapy dose calculation. Methods: The treatment heads at the Francis H Burr Proton Therapy Center were modeled with a new Monte Carlo tool (TOPAS based on Geant4). For variance reduction purposes, two particle-splitting planes were implemented. First, the particles were split upstream of the second scatterer or at the second ionization chamber. Then, particles reaching another plane immediately upstream of the field specific aperture were split again. In each case, particles were split by a factor of 8. At the second ionization chamber and at the latter plane, the cylindrical symmetry of the proton beam was exploited to position the split particles at randomly spaced locations rotated around the beam axis. Phase space data in IAEA format were recorded at the treatment head exit and the computational efficiency was calculated. Depth-dose curves and beam profiles were analyzed. Dose distributions were compared for a voxelized water phantom for different treatment fields for both the reference and optimized simulations. In addition, dose in two patients was simulated with and without particle splitting to compare the efficiency and accuracy of the technique. Results: A normalized computational efficiency gain of a factor of 10-20.3 was reached for phase space calculations for the different treatment head options simulated. Depth-dose curves and beam profiles were in reasonable agreement with the simulation done without splitting: within 1% for depth-dose with an average difference of (0.2 {+-} 0.4)%, 1 standard deviation, and a 0.3% statistical uncertainty of the simulations in the high dose region; 1.6% for planar fluence with an average difference of (0.4 {+-} 0.5)% and a statistical uncertainty of 0.3% in the high fluence region. The percentage differences between dose distributions in water for simulations done with and without particle splitting were within the accepted clinical tolerance of 2%, with a 0.4% statistical uncertainty. For the two patient geometries considered, head and prostate, the efficiency gain was 20.9 and 14.7, respectively, with the percentages of voxels with gamma indices lower than unity 98.9% and 99.7%, respectively, using 2% and 2 mm criteria. Conclusions: The authors have implemented an efficient variance reduction technique with significant speed improvements for proton Monte Carlo simulations. The method can be transferred to other codes and other treatment heads.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-C-204-06: Monte Carlo Dose Calculation for Kilovoltage X-Ray-Psoralen Activated Cancer Therapy (X-PACT): Preliminary Results

Purpose: X-PACT is an experimental cancer therapy where kV x-rays are used to photo-activate anti-cancer therapeutics through phosphor intermediaries (phosphors that absorb x-rays and re-radiate as UV light). Clinical trials in pet dogs are currently underway (NC State College of Veterinary Medicine) and an essential component is the ability to model the kV dose in these dogs. Here we report the commissioning and characterization of a Monte Carlo (MC) treatment planning simulation tool to calculate X-PACT radiation doses in canine trials. Methods: FLUKA multi-particle MC simulation package was used to simulate a standard X-PACT radiation treatment beam of 80kVp with the Varian OBI x-ray source geometry. The beam quality was verified by comparing measured and simulated attenuation of the beam by various thicknesses of aluminum (2–4.6 mm) under narrow beam conditions (HVL). The beam parameters at commissioning were then corroborated using MC, characterized and verified with empirically collected commissioning data, including: percent depth dose curves (PDD), back-scatter factors (BSF), collimator scatter factor(s), and heel effect, etc. All simulations were conducted for N=30M histories at M=100 iterations. Results: HVL and PDD simulation data agreed with an average percent error of 2.42%±0.33 and 6.03%±1.58, respectively. The mean square error (MSE) values for HVL and PDD (0.07% and 0.50%) were low, as expected; however, longer simulations are required to validate convergence to the expected values. Qualitatively, pre- and post-filtration source spectra matched well with 80kVp references generated via SPEKTR software. Further validation of commissioning data simulation is underway in preparation for first-time 3D dose calculations with canine CBCT data. Conclusion: We have prepared a Monte Carlo simulation capable of accurate dose calculation for use with ongoing X-PACT canine clinical trials. Preliminary results show good agreement with measured data and hold promise for accurate quantification of dose for this novel psoralen X-ray therapy. Funding Support, Disclosures, & Conflict of Interest: The Monte Carlo simulation work was not funded; Drs. Adamson & Oldham have received funding from Immunolight LLC for X-PACT research.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Peculiarities of Forming a Multiband Transmission Function Based on Multifrequency Acousto-Optic Diffraction

An analytical solution for the transmission function of noncoherent wideband radiation is obtained under acousto-optic (AO) filtering using a discrete set of monochromatic AO waves with a small spectral overlap. We studied characteristics of the AO transformation of a continuous spectrum of noncoherent radiation into a given set of discrete narrow bands of spectral transmission by excitation of a discrete set of sound frequencies. We carried out the analysis of transmission functions of individual channels taking into account a partial overlap of their spectra and possible intermodulation distortions. It is shown that a stationary value of the root-mean-square light power is found at the electronic output due to the photoelectric transformation and detecting diffracted light. Based on this, a necessary stationary, multiband, and nearly equidistant transmission function of a device can be formed by using a relevant spectrum of acoustic excitation. Peculiarities of this way of forming the multiband transmission function are revealed: the limitation of diffraction efficiency for an individual channel, the possibility of decoupling side lobes of adjacent channels, etc. A multiband acousto-optic filter (MAOF) was simulated that was based on a paratellurite monocrystal (TeO{sub 2}), which was previously used for experimental optical encoding. The theoretical and experimental results are in gratifying agreement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Gyrotron generation of broadband chaotic radiation under overlapping of high- and low-frequency resonances

Significant (up to 10%) broadening of the band of noise-like radiation in gyrotrons is possible due to specific tuning of gyrofrequency relative to the critical frequency of the working mode when the high- and low-frequency cyclotron resonances are substantially separated at the given translation and transverse velocities of electrons. The generation bands at the resonances are overlapped under operation above the threshold. The analysis with allowance for finiteness of the electron transit time in the interaction space and, hence, different slopes of the dispersion characteristics of electron beam and wave is critically important for correct investigation of the generation regimes in the framework of the average approach. The results are verified using direct 3D particle-in-cell simulation of the chaotic generation regimes of the gyrotron generation in the 8-mm-wavelength range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Light Water Reactor Sustainability Program: Analysis of Pressurized Water Reactor Station Blackout Caused by External Flooding Using the RISMC Toolkit

The existing fleet of nuclear power plants is in the process of extending its lifetime and increasing the power generated from these plants via power uprates. In order to evaluate the impact of these factors on the safety of the plant, the Risk Informed Safety Margin Characterization (RISMC) project aims to provide insight to decision makers through a series of simulations of the plant dynamics for different initial conditions (e.g., probabilistic analysis and uncertainty quantification). This report focuses, in particular, on the application of a RISMC detailed demonstration case study for an emergent issue using the RAVEN and RELAP-7 tools. This case study looks at the impact of a couple of challenges to a hypothetical pressurized water reactor, including: (1) a power uprate, (2) a potential loss of off-site power followed by the possible loss of all diesel generators (i.e., a station black-out event), (3) and earthquake induces station-blackout, and (4) a potential earthquake induced tsunami flood. The analysis is performed by using a set of codes: a thermal-hydraulic code (RELAP-7), a flooding simulation tool (NEUTRINO) and a stochastic analysis tool (RAVEN) – these are currently under development at the Idaho National Laboratory.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

CRAB FLARES DUE TO TURBULENT DISSIPATION OF THE PULSAR STRIPED WIND

We interpret γ-ray flares from the Crab Nebula as the signature of turbulence in the pulsar’s electromagnetic outflow. Turbulence is triggered upstream by dynamical instability of the wind’s oscillating magnetic field and accelerates non-thermal particles. On impacting the wind-termination shock, these particles emit a distinct synchrotron component F{sub ν,flare}, which is constantly modulated by intermittency of the upstream plasma flow. Flares are observed when the high-energy cutoff of F{sub ν,flare} emerges above the fast-declining nebular emission around 0.1–1 GeV. Simulations carried out in the force-free electrodynamics approximation predict the striped wind to become fully turbulent well ahead of the wind-termination shock, provided its terminal Lorentz factor is ≲10{sup 4}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

USING LEAKED POWER TO MEASURE INTRINSIC AGN POWER SPECTRA OF RED-NOISE TIME SERIES

Fluxes emitted at different wavebands from active galactic nuclei (AGNs) fluctuate at both long and short timescales. The variation can typically be characterized by a broadband power spectrum, which exhibits a red-noise process at high frequencies. The standard method of estimating the power spectral density (PSD) of AGN variability is easily affected by systematic biases such as red-noise leakage and aliasing, in particular when the observation spans a relatively short period and is gapped. Focusing on the high-frequency PSD that is strongly distorted due to red-noise leakage and usually not significantly affected by aliasing, we develop a novel and observable normalized leakage spectrum (NLS), which sensitively describes the effects of leaked red-noise power on the PSD at different temporal frequencies. Using Monte Carlo simulations, we demonstrate how an AGN underlying PSD sensitively determines the NLS when there is severe red-noise leakage, and thereby how the NLS can be used to effectively constrain the underlying PSD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SEMI-ANALYTIC GALAXY EVOLUTION (SAGE): MODEL CALIBRATION AND BASIC RESULTS

This paper describes a new publicly available codebase for modeling galaxy formation in a cosmological context, the “Semi-Analytic Galaxy Evolution” model, or sage for short.{sup 5} sage is a significant update to the 2006 model of Croton et al. and has been rebuilt to be modular and customizable. The model will run on any N-body simulation whose trees are organized in a supported format and contain a minimum set of basic halo properties. In this work, we present the baryonic prescriptions implemented in sage to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include the following: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling–radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population.

79 ASTRONOMY AND ASTROPHYSICS↗