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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Exploring dark matter microphysics with galaxy surveys

We use present cosmological observations and forecasts of future experiments to illustrate the power of large-scale structure (LSS) surveys in probing dark matter (DM) microphysics and unveiling potential deviations from the standard ΛCDM scenario. To quantify this statement, we focus on an extension of ΛCDM with DM-neutrino scattering, which leaves a distinctive imprint on the angular and matter power spectra. After finding that future CMB experiments (such as COrE+) will not significantly improve the constraints set by the Planck satellite, we show that the next generation of galaxy clustering surveys (such as DESI) could play a leading role in constraining alternative cosmologies and even have the potential to make a discovery. Typically we find that DESI would be an order of magnitude more sensitive to DM interactions than Planck, thus probing effects that until now have only been accessible via N-body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Dissipative dark matter and the Andromeda plane of satellites

We show that dissipative dark matter can potentially explain the large observed mass to light ratio of the dwarf satellite galaxies that have been observed in the recently identified planar structure around Andromeda, which are thought to result from tidal forces during a galaxy merger. Whereas dwarf galaxies created from ordinary disks would be dark matter poor, dark matter inside the galactic plane not only provides a source of dark matter, but one that is more readily bound due to the dark matter's lower velocity. This initial N-body study shows that with a thin disk of dark matter inside the baryonic disk, mass-to-light ratios as high as O(90) can be generated when tidal forces pull out patches of sizes similar to the scales of Toomre instabilities of the dark disk. A full simulation will be needed to confirm this result.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Gaussian variational ansatz in the problem of anomalous sea waves: Comparison with direct numerical simulation

The nonlinear dynamics of an obliquely oriented wave packet on a sea surface is analyzed analytically and numerically for various initial parameters of the packet in relation to the problem of the so-called rogue waves. Within the Gaussian variational ansatz applied to the corresponding (1+2)-dimensional hyperbolic nonlinear Schrödinger equation (NLSE), a simplified Lagrangian system of differential equations is derived that describes the evolution of the coefficients of the real and imaginary quadratic forms appearing in the Gaussian. This model provides a semi-quantitative description of the process of nonlinear spatiotemporal focusing, which is one of the most probable mechanisms of rogue wave formation in random wave fields. The system of equations is integrated in quadratures, which allows one to better understand the qualitative differences between linear and nonlinear focusing regimes of a wave packet. Predictions of the Gaussian model are compared with the results of direct numerical simulation of fully nonlinear long-crested waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Visualization of Dynamic Vortex Structures in Magnetic Films with Uniaxial Anisotropy (Micromagnetic Simulation)

Three-dimensional computer simulation of dynamic processes in a moving domain boundary separating domains in a soft magnetic uniaxial film with planar anisotropy is performed by numerical solution of Landau-Lifshitz-Gilbert equations. The developed visualization methods are used to establish the connection between the motion of surface vortices and antivortices, singular (Bloch) points, and core lines of intrafilm vortex structures. A relation between the character of magnetization dynamics and the film thickness is found. The analytical models of spatial vortex structures for imitation of topological properties of the structures observed in micromagnetic simulation are constructed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Monte Carlo Techniques for Nuclear Systems - Theory Lectures

These are lecture notes for a Monte Carlo class given at the University of New Mexico. The following topics are covered: course information; nuclear eng. review & MC; random numbers and sampling; computational geometry; collision physics; tallies and statistics; eigenvalue calculations I; eigenvalue calculations II; eigenvalue calculations III; variance reduction; parallel Monte Carlo; parameter studies; fission matrix and higher eigenmodes; doppler broadening; Monte Carlo depletion; HTGR modeling; coupled MC and T/H calculations; fission energy deposition. Solving particle transport problems with the Monte Carlo method is simple - just simulate the particle behavior. The devil is in the details, however. These lectures provide a balanced approach to the theory and practice of Monte Carlo simulation codes. The first lectures provide an overview of Monte Carlo simulation methods, covering the transport equation, random sampling, computational geometry, collision physics, and statistics. The next lectures focus on the state-of-the-art in Monte Carlo criticality simulations, covering the theory of eigenvalue calculations, convergence analysis, dominance ratio calculations, bias in Keff and tallies, bias in uncertainties, a case study of a realistic calculation, and Wielandt acceleration techniques. The remaining lectures cover advanced topics, including HTGR modeling and stochastic geometry, temperature dependence, fission energy deposition, depletion calculations, parallel calculations, and parameter studies. This portion of the class focuses on using MCNP to perform criticality calculations for reactor physics and criticality safety applications. It is an intermediate level class, intended for those with at least some familiarity with MCNP. Class examples provide hands-on experience at running the code, plotting both geometry and results, and understanding the code output. The class includes lectures & hands-on computer use for a variety of Monte Carlo calculations. Beginning MCNP users are encouraged to review LA-UR-09-00380, "Criticality Calculations with MCNP: A Primer (3nd Edition)" (available at http:// mcnp.lanl.gov under "Reference Collection") prior to the class. No Monte Carlo class can be complete without having students write their own simple Monte Carlo routines for basic random sampling, use of the random number generator, and simplified particle transport simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Flexible Production Systems with Probabilistic Perturbations

Probabilistic perturbations associated with scattering of the operating periods and equipment faults in flexible production systems are considered. The perturbations are taken into account in computer simulation of the operation of a flexible production system. The proposed modeling application is described.

42 ENGINEERING↗
DOE OSTI2016

TURBULENT TRANSPORT IN A STRONGLY STRATIFIED FORCED SHEAR LAYER WITH THERMAL DIFFUSION

This work presents numerical results on the transport of heat and chemical species by shear-induced turbulence in strongly stratified, thermally diffusive environments. The shear instabilities driven in this regime are sometimes called “secular” shear instabilities, and can take place when the Richardson number of the flow is large, provided the Péclet number is small. We have identified a set of simple criteria to determine whether these instabilities can take place or not. Generally speaking, we find that they may be relevant whenever the thermal diffusivity of the fluid is very large (typically larger than 10{sup 14} cm{sup 2} s{sup −1}), which is the case in the outer layers of high-mass stars (M ≥ 10 M{sub ⊙}), for instance. Using a simple model setup in which the shear is forced by a spatially sinusoidal, constant-amplitude body-force, we have identified several regimes ranging from effectively unstratified to very strongly stratified, each with its own set of dynamical properties. Unless the system is in one of the two extreme regimes (effectively unstratified or completely stable), however, we find that (1) only about 10% of the input power is used toward heat transport, while the remaining 90% is viscously dissipated; (2) that the effective compositional mixing coefficient is well-approximated by the model of Zahn, with D ≃ 0.02κ{sub T}/J where κ{sub T} is the thermal diffusivity and J is the Richardson number. These results need to be confirmed, however, with simulations in different model setups and at higher effective Reynolds number.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Gas dynamics in the milky way: a low pattern speed model

We present gas flow models for the Milky Way based on high-resolution grid-based hydrodynamical simulations. The basic galactic potential we use is from an N-body model constrained by the density of red clump giants in the Galactic bulge. We augment this potential with a nuclear bulge, two pairs of spiral arms, and additional mass at the bar end to represent the long bar component. With this combined model we can reproduce many features in the observed (l,v) diagram with a bar pattern speed of 33 km s{sup −1} kpc{sup −1} and a spiral pattern speed of 23 km s{sup −1} kpc{sup −1}. The shape and kinematics of the nuclear ring, Bania’s Clump 2, the Connecting arm, the Near and Far 3 kpc arms, the Molecular Ring, and the spiral arm tangent points in our simulations are comparable to those in the observations. Our results imply that a low pattern speed model for the bar in our Milky Way reproduces the observations for a suitable Galactic potential. Our best model gives a better match to the (l,v) diagram than previous high pattern speed hydrodynamical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Alignments of dark matter halos with large-scale tidal fields: mass and redshift dependence

Large-scale tidal fields estimated directly from the distribution of dark matter halos are used to investigate how halo shapes and spin vectors are aligned with the cosmic web. The major, intermediate, and minor axes of halos are aligned with the corresponding tidal axes, and halo spin axes tend to be parallel with the intermediate axes and perpendicular to the major axes of the tidal field. The strengths of these alignments generally increase with halo mass and redshift, but the dependence is only on the peak height, ν≡δ{sub c}/σ(M{sub h},z). The scaling relations of the alignment strengths with the value of ν indicate that the alignment strengths remain roughly constant when the structures within which the halos reside are still in a quasi-linear regime, but decreases as nonlinear evolution becomes more important. We also calculate the alignments in projection so that our results can be compared directly with observations. Finally, we investigate the alignments of tidal tensors on large scales, and use the results to understand alignments of halo pairs separated at various distances. Our results suggest that the coherent structure of the tidal field is the underlying reason for the alignments of halos and galaxies seen in numerical simulations and in observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Helium reionization simulations. I. Modeling quasars as radiation sources

We introduce a new project to understand helium reionization using fully coupled N-body, hydrodynamics, and radiative transfer simulations. This project aims to capture correctly the thermal history of the intergalactic medium as a result of reionization and make predictions about the Lyα forest and baryon temperature–density relation. The dominant sources of radiation for this transition are quasars, so modeling the source population accurately is very important for making reliable predictions. In this first paper, we present a new method for populating dark matter halos with quasars. Our set of quasar models includes two different light curves, a lightbulb (simple on/off) and symmetric exponential model, and luminosity-dependent quasar lifetimes. Our method self-consistently reproduces an input quasar luminosity function given a halo catalog from an N-body simulation, and propagates quasars through the merger history of halo hosts. After calibrating quasar clustering using measurements from the Baryon Oscillation Spectroscopic Survey, we find that the characteristic mass of quasar hosts is M{sub h}∼2.5×10{sup 12} h{sup −1} M{sub ⊙} for the lightbulb model, and M{sub h}∼2.3×10{sup 12} h{sup −1} M{sub ⊙} for the exponential model. In the latter model, the peak quasar luminosity for a given halo mass is larger than that in the former, typically by a factor of 1.5–2. The effective lifetime for quasars in the lightbulb model is 59 Myr, and in the exponential case, the effective time constant is about 15 Myr. We include semi-analytic calculations of helium reionization, and discuss how to include these quasars as sources of ionizing radiation for full hydrodynamics with radiative transfer simulations in order to study helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

An accretion-jet model for M87: interpreting the spectral energy distribution and faraday rotation measure

M87 is arguably the best supermassive black hole (BH) to explore jet and/or accretion physics, due to its proximity and fruitful high-resolution multi-waveband observations. We model the multi-wavelength spectral energy distribution (SED) of the M87 core that observed at a scale of 0.4 arcsec (∼10{sup 5} R {sub g}, R {sub g} is gravitational radius), as recently presented by Prieto et al. Similar to Sgr A*, we find that the millimeter bump as observed by the Atacama Large Millimeter/submillimeter Array can be modeled by the synchrotron emission of the thermal electrons in an advection-dominated accretion flow (ADAF), while the low-frequency radio emission and X-ray emission may predominantly come from the jet. The millimeter radiation from ADAF predominantly comes from the region within 10R {sub g}, which is roughly consistent with the recent very long baseline interferometry observations at 230 GHz. We further calculate the Faraday rotation measure (RM) from both ADAF and jet models, and find that the RM predicted from the ADAF is roughly consistent with the measured value, while the RM predicted from the jet is much higher if jet velocity close to the BH is low or moderate (e.g., v {sub jet} ≲ 0.6 c). With the constraints from the SED modeling and RM, we find that the accretion rate close to the BH horizon is ∼(0.2--1) ×10{sup −3} M{sub ⊙} yr{sup −1}≪ M-dot {sub B}∼0.2 M{sub ⊙} yr{sup −1} ( M-dot {sub B} is Bondi accretion rate), where the electron density profile, n {sub e} ∝ r {sup ∼−1}, in the accretion flow, is consistent with that determined from X-ray observation inside the Bondi radius and recent numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Magnetic Flux Emergence and Decay Rates for Preceder and Follower Sunspots Observed with HMI

We quantify the emergence and decay rates of preceder (p) and follower (f) sunspots within 10 active regions from 2010 to 2014 using Space-weather Helioseismic Magnetic Imager Active Region Patch data. The sunspots are small to mid-sized regions and contain a signed flux within a single polarity sunspot of (1.1--6.5)×10{sup 21} Mx. The net unsigned flux within the regions, including plage, ranges from (5.1--20)×10{sup 21} Mx. Rates are calculated with and without intensity contours to differentiate between sunspot formation and flux emergence. Signed flux emergence rates, calculated with intensity contours, for the p (f) spots average 6.8(4.9)×10{sup 19} Mx hr{sup −1}, while decay rates are −1.9(−3.4)×10{sup 19} Mx hr{sup −1}. The mean, signed flux emergence rate of the regions, including plage, is 7.1×10{sup 19} Mx hr{sup −1}, for a mean peak flux of 5.9×10{sup 21} Mx. Using a synthesis of these results and others reported previously, there is a clear trend for larger flux regions to emerge faster than smaller ones. Observed emergence rates (dϕ/dt, Mx hr{sup −1}) scale with total signed peak flux, ϕ-tilde {sub max}, as a power law with an exponent of 0.36, i.e., dϕ/dt=A ϕ-tilde {sub max}{sup 0.36}. The observed rates may assist in constraining the boundary and initial conditions in simulations which already demonstrate increased rates for flux tubes with higher buoyancy and twist, or in the presence of a strong upflow. Overall, the observed emergence rates are smaller than those in simulations, which may indicate a slower rise of the flux in the interior than what is captured in simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Is There a Disk of Satellites around the Milky Way?

The “disk of satellites” (DoS) around the Milky Way is a highly debated topic with conflicting interpretations of observations and their theoretical models. We perform a comprehensive analysis of all of the dwarfs detected in the Milky Way and find that the DoS structure depends strongly on the plane identification method and the sample size. In particular, we demonstrate that a small sample size can artificially produce a highly anisotropic spatial distribution and a strong clustering of the angular momentum of the satellites. Moreover, we calculate the evolution of the 11 classical satellites with proper motion measurements and find that the thin DoS in which they currently reside is transient. Furthermore, we analyze two cosmological simulations using the same initial conditions of a Milky-Way-sized galaxy, an N -body run with dark matter only, and a hydrodynamic one with both baryonic and dark matter, and find that the hydrodynamic simulation produces a more anisotropic distribution of satellites than the N -body one. Our results suggest that an anisotropic distribution of satellites in galaxies can originate from baryonic processes in the hierarchical structure formation model, but the claimed highly flattened, coherently rotating DoS of the Milky Way may be biased by the small-number selection effect. These findings may help resolve the contradictory claims of DoS in galaxies and the discrepancy among numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Effective use of DXTRAN in MCNP - 25480

The DXTRAN variance reduction method enables Monte Carlo radiation transport calculations to deterministically place particle tracks on spheres surrounding geo-metric regions of interest. For many problems DXTRAN is the best way to transport particles to small distant regions or to angle bias particles exiting collisions. We provide guidance for optimal selection of parameters when using both single and nested DXTRAN spheres in Monte Carlo N-Particle radiation transport code (MCNP) calculations. Our recommendations are based upon both the underlying mathematics of DXTRAN and two demonstration MCNP simulations. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Weak Convergence of a Mass-Structured Individual-Based Model

We propose a model of chemostat where the bacterial population is individually-based, each bacterium is explicitly represented and has a mass evolving continuously over time. The substrate concentration is represented as a conventional ordinary differential equation. These two components are coupled with the bacterial consumption. Mechanisms acting on the bacteria are explicitly described (growth, division and washout). Bacteria interact via consumption. We set the exact Monte Carlo simulation algorithm of this model and its mathematical representation as a stochastic process. We prove the convergence of this process to the solution of an integro-differential equation when the population size tends to infinity. Finally, we propose several numerical simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Microstructural study and simulation of intrinsic two-way shape memory behavior of functionally graded Ni-rich/NiTiCu thin film

Highlights: • Ni-rich/Ti-rich bi-layer thin film was deposited on Si by using DC magnetron sputtering. • The bi-layer exhibits the gradient structure and mechanical properties through the thickness. • The bi-layer shows intrinsic two-way shape memory effect due to functionally graded structure. • The intrinsic two-way shape memory effect was confirmed by finite element simulation. - Abstract: In this research, the gradient structure and mechanical properties of an austenitic/martensitic (Ni{sub 50.8}Ti/Ni{sub 45}TiCu{sub 5} (at.%)) bi-layer thin film was investigated. The bi-layer was deposited on a Si substrate using DC magnetron sputtering. After crystallization of the film at 773 K for 60 min, the microstructure and mechanical properties of the bi-layer thin film were characterized using X-ray diffraction (XRD), transmission electron microscopy (STEM and HRTEM) and nanoindentation, respectively. The diffraction pattern illustrated that the crystallized bi-layer was combined of martensitic and austenitic layers while secondary ion mass spectroscopy (SIMS) and high resolution transmission electron microscopy analysis demonstrated the existence of a compositional gradient through the thickness and residual strain in the interface of the bi-layer, respectively. Furthermore, the compositional gradient in the bi-layer led to gradual variations in the structure, hardness and Young's modulus through the thickness of the bi-layer. The intrinsic two-way shape memory effect due to its functionally graded structure and the induced stress was confirmed by experimental test and finite element simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Direct imaging of light elements by annular dark-field aberration-corrected scanning transmission electron microscopy

In this report, we show that an annular dark-field detector in an aberration-corrected scanning transmission electron microscope allows the direct observation of light element columns in crystalline lattices. At specific imaging conditions, an enhancement of the intensities of light element columns in the presence of heavy element columns is observed. Experimental results are presented for imaging the nitrogen and carbon atomic columns at the GaN-SiC interface and within the GaN and SiC compounds. The crystal polarity of GaN at the interface is identified. The obtained findings are discussed and are well supported by image simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

A fail-safe magnetorheological energy absorber for shock and vibration isolation

Magnetorheological (MR) energy absorbers (EAs) are an effective adaptive EA technology with which to maximize shock and vibration isolation. However, to realize maximum performance of the semi-active control system, the off-state (i.e., field off) stroking load of the MREA must be minimized at all speeds, and the dynamic range of the MREA must be maximized at high speed. This study presents a fail-safe MREA (MREA-FS) concept that, can produce a greater dynamic range at all piston speeds. A bias damping force is generated in the MREA-FS using permanent magnetic fields, which enables fail-safe behavior in the case of power failure. To investigate the feasibility and capability of the MREA-FS in the context of the semi-active control systems, a single-degree-of-freedom base excited rigid payload is mathematically constructed and simulated with skyhook control.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗