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

Numerical analysis of experiments on the generation of shock waves in aluminium under indirect (X-ray) action on the Iskra-5 facility

The dynamics of laser and X-ray radiation fields in experiments with cylindrical converter boxes (illuminators), which had earlier been carried out on the Iskra-5 laser facility (the second harmonic of iodine laser radiation, {lambda} = 0.66 {mu}m) was investigated in a sector approximation using the SND-LIRA numerical technique. In these experiments, the X-ray radiation temperature in the box was determined by measuring the velocity of the shock wave generated in the sample under investigation, which was located at the end of the cylindrical illuminator. Through simulations were made using the SND-LIRA code, which took into account the absorption of laser driver radiation at the box walls, the production of quasithermal radiation, as well as the formation and propagation of the shock wave in the sample under investigation. An analysis of the experiments permits determining the electron thermal flux limiter f: for f = 0.03 it is possible to match the experimental scaling data for X-ray in-box radiation temperature to the data of our simulations. The shock velocities obtained from the simulations are also consistent with experimental data. In particular, in the experiment with six laser beams (and a laser energy E{sub L} = 1380 J introduced into the box) the velocity of the shock front (determined from the position of a laser mark) after passage through a 50-{mu}m thick base aluminium layer was equal to 35{+-}1.6 km s{sup -1}, and in simulations to 36 km s{sup -1}. In the experiment with four laser beams (for E{sub L} = 850 J) the shock velocity (measured from the difference of transit times through the base aluminium layer and an additional thin aluminium platelet) was equal to 30{+-}3.6 km s{sup -1}, and in simulations to 30 km s{sup -1}. (interaction of laser radiation with matter)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Modulation instability and short-pulse generation in media with relaxing Kerr nonlinearity and high self-steepening

The modulation instability in waveguides with high Kerr nonlinearity, characterised by a delayed nonlinear response, has been investigated with allowance for the self-steepening parameter and third-order dispersion. General expressions for the modulation gain are obtained. The influence of the waveguide parameters on the gain is analysed. It is shown that the joint effect of the delayed nonlinear response and negative nonlinearity dispersion leads to an increase in the modulation gain. The relations obtained are confirmed by numerical simulation. The results of this study can be used to design compact generators of high-frequency pulse trains. (nonlinear optical phenomena)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Simulation of picosecond pulse propagation in fibre-based radiation shaping units

We have performed a numerical simulation of picosecond pulse propagation in a combined stretcher consisting of a segment of a telecommunication fibre and diffraction holographic gratings. The process of supercontinuum generation in a nonlinear photonic-crystal fibre pumped by picosecond pulses is simulated by solving numerically the generalised nonlinear Schrödinger equation; spectral and temporal pulse parameters are determined. Experimental data are in good agreement with simulation results. The obtained results are used to design a high-power femtosecond laser system with a pulse repetition rate of 1 kHz.

42 ENGINEERING↗
DOE OSTI2011

ASSESSING QUANTITATIVE RESULTS IN ACCRETION SIMULATIONS: FROM LOCAL TO GLOBAL

Discretized numerical simulations are a powerful tool for the investigation of nonlinear MHD turbulence in accretion disks. However, confidence in their quantitative predictions requires a demonstration that further refinement of the spatial grid scale would not result in any significant change. This has yet to be accomplished, particularly for global disk simulations. In this paper, we combine data from previously published stratified shearing box simulations and new global disk simulations to calibrate several quantitative diagnostics by which one can estimate progress toward numerical convergence of the magnetic field. Using these diagnostics, we find that the established criterion for an adequate numerical description of linear growth of the magneto-rotational instability (the number of cells across a wavelength of the fastest-growing vertical wavenumber mode) can be extended to a criterion for the adequate description of nonlinear MHD disk turbulence, but the standard required is more stringent. We also find that azimuthal resolution, which has received little attention in previous studies, can significantly affect the evolution of the poloidal magnetic field. We further analyze the comparative resolution requirements of a small sample of initial magnetic field geometries; not surprisingly, more complicated initial field geometries require higher spatial resolution. Otherwise, they tend to evolve to qualitatively similar states if evolved for sufficient time. Applying our quantitative resolution criteria to a sample of previously published global simulations, we find that, with perhaps a single exception, they are significantly underresolved, and therefore underestimate the magnetic turbulence and resulting stress levels throughout the accretion flow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

DEEP ABSORPTION LINE STUDIES OF QUIESCENT GALAXIES AT z {approx} 2: THE DYNAMICAL-MASS-SIZE RELATION AND FIRST CONSTRAINTS ON THE FUNDAMENTAL PLANE

We present dynamical and structural scaling relations of quiescent galaxies at z = 2, including the dynamical-mass-size relation and the first constraints on the fundamental plane (FP). The backbone of the analysis is a new, very deep Very Large Telescope/X-shooter spectrum of a massive, compact, quiescent galaxy at z = 2.0389. We detect the continuum between 3700 and 22,000 A and several strong absorption features (Balmer series, Ca H+K, G band) from which we derive a stellar velocity dispersion of 318 {+-} 53 km s{sup -1}. We perform detailed modeling of the continuum emission and line indices and derive strong simultaneous constraints on the age, metallicity, and stellar mass. The galaxy is a dusty (A{sub V} = 0.77{sup +0.36}{sub -0.32}) solar metallicity (log(Z/Z{sub Sun }) = 0.02{sup +0.20}{sub -0.41}) post-starburst galaxy, with a mean-luminosity-weighted log(age/yr) of 8.9 {+-} 0.1. The galaxy formed the majority of its stars at z > 3 and currently has little or no ongoing star formation. We compile a sample of three other z {approx} 2 quiescent galaxies with measured velocity dispersions, two of which are also post-starburst like. Their dynamical-mass-size relation is offset significantly less than the stellar-mass-size relation from the local early-type relations, which we attribute to a lower central dark matter fraction. Recent cosmological merger simulations agree qualitatively with the data, but cannot fully account for the evolution in the dark matter fraction. The z {approx} 2 FP requires additional evolution beyond passive stellar aging to be in agreement with the local FP. The structural evolution predicted by the cosmological simulations is insufficient, suggesting that additional, possibly non-homologous, structural evolution is needed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A DERIVATION OF (HALF) THE DARK MATTER DISTRIBUTION FUNCTION

All dark matter structures appear to follow a set of universalities, such as phase-space density or velocity anisotropy profiles; however, the origin of these universalities remains a mystery. Any equilibrated dark matter structure can be fully described by two functions, namely the radial and tangential velocity distribution functions (VDFs), and once these two are understood we will understand all the observed universalities. Here, we demonstrate that if we know the radial VDF then we can derive and understand the tangential VDF. This is based on simple dynamical arguments about properties of collisionless systems. We use a range of controlled numerical simulations to demonstrate the accuracy of this result. We therefore boil the question of the dark matter structural properties down to understanding the radial VDF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

INTERPLAY OF NEUTRINO OPACITIES IN CORE-COLLAPSE SUPERNOVA SIMULATIONS

We have conducted a series of numerical experiments using spherically symmetric, general relativistic, neutrino radiation hydrodynamics with the code Agile-BOLTZTRAN to examine the effects of modern neutrino opacities on the development of supernova simulations. We test the effects of opacities by removing opacities or by undoing opacity improvements for individual opacities and groups of opacities. We find that improvements to electron capture (EC) on nuclei, namely EC on an ensemble of nuclei using modern nuclear structure models rather than the simpler independent-particle approximation (IPA) for EC on a mean nucleus, plays the most important role during core collapse of all tested neutrino opacities. Low-energy neutrinos emitted by modern nuclear EC preferentially escape during collapse without the energy downscattering on electrons required to enhance neutrino escape and deleptonization for the models with IPA nuclear EC. During shock breakout the primary influence on the emergent neutrinos arises from non-isoenergetic scattering (NIS) on electrons. For the accretion phase, NIS on free nucleons and pair emission by e {sup +} e {sup -} annihilation have the largest impact on the neutrino emission and shock evolution. Other opacities evaluated, including nucleon-nucleon bremsstrahlung and especially neutrino-positron scattering, have little measurable impact on neutrino emission or shock dynamics. Modern treatments of nuclear EC, e {sup +} e {sup -}-annihilation pair emission, and NIS on electrons and free nucleons are critical elements of core-collapse simulations of all dimensionality.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE GALEX ARECIBO SDSS SURVEY. VII. THE BIVARIATE NEUTRAL HYDROGEN-STELLAR MASS FUNCTION FOR MASSIVE GALAXIES

We present the bivariate neutral atomic hydrogen (H I)-stellar mass function (HISMF) φ(M{sub H{sub I}}, M{sub *}) for massive (log M{sub *}/M{sub ☉} \gt 10) galaxies derived from a sample of 480 local (0.025 < z < 0.050) galaxies observed in H I at Arecibo as part of the GALEX Arecibo SDSS Survey. We fit six different models to the HISMF and find that a Schechter function that extends down to a 1% H I gas fraction, with an additional fractional contribution below that limit, is the best parameterization of the HISMF. We calculate Ω{sub H{sub I,{sub M{sub *>10{sup 1}{sup 0}}}}} and find that massive galaxies contribute 41% of the H I density in the local universe. In addition to the binned HISMF, we derive a continuous bivariate fit, which reveals that the Schechter parameters only vary weakly with stellar mass: M{sub H{sub I}{sup *}}, the characteristic H I mass, scales as M{sub *}{sup 0.39}; α, the slope of the HISMF at moderate H I masses, scales as M{sub *}{sup 0.07}; and f, the fraction of galaxies with H I gas fraction greater than 1%, scales as M{sub *}{sup -0.24}. The variation of f with stellar mass should be a strong constraint for numerical simulations. To understand the physical mechanisms that produce the shape of the HISMF, we redefine the parameters of the Schechter function as explicit functions of stellar mass and star formation rate (SFR) to produce a trivariate fit. This analysis reveals strong trends with SFR. While M{sub H{sub I}{sup *}} varies weakly with stellar mass and SFR (M{sub H{sub I}{sup *}} ∝ M{sub *}{sup 0.22}, M{sub H{sub I}{sup *}} ∝ SFR{sup –0.03}), α is a stronger function of both stellar mass and especially SFR (α ∝ M{sub *}{sup 0.47}, α ∝ SFR{sup 0.95}). The HISMF is a crucial tool that can be used to constrain cosmological galaxy simulations, test observational predictions of the H I content of populations of galaxies, and identify galaxies whose properties deviate from average trends.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

The pre-merger impact velocity of the binary cluster A1750 from X-ray, lensing, and hydrodynamical simulations

Since the discovery of the 'Bullet Cluster', several similar cases have been uncovered that suggest relative velocities well beyond the tail of high speed collisions predicted by the concordance ΛCDM model. However, quantifying such post-merger events with hydrodynamical models requires a wide coverage of possible initial conditions. Here, we show that it is simpler to interpret pre-merger cases, such as A1750, where the gas between the colliding clusters is modestly affected, so that the initial conditions are clear. We analyze publicly available Chandra data confirming a significant increase in the projected X-ray temperature between the two cluster centers in A1750 consistent with our expectations for a merging cluster. We model this system with a self-consistent hydrodynamical simulation of dark matter and gas using the FLASH code. Our simulations reproduce well the X-ray data and the measured redshift difference between the two clusters in the phase before the first core passage viewed at an intermediate projection angle. The deprojected initial relative velocity derived using our model is 1460 km s{sup –1}, which is considerably higher than the predicted mean impact velocity for simulated massive haloes derived by recent ΛCDM cosmological simulations, but is within the allowed range. Our simulations demonstrate that such systems can be identified using a multi-wavelength approach and numerical simulations, for which the statistical distribution of relative impact velocities may provide a definitive examination of a broad range of dark matter scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

On the outer edges of protoplanetary dust disks

The expectation that aerodynamic drag will force the solids in a gas-rich protoplanetary disk to spiral in toward the host star on short timescales is one of the fundamental problems in planet formation theory. The nominal efficiency of this radial drift process is in conflict with observations, suggesting that an empirical calibration of solid transport mechanisms in a disk is highly desirable. However, the fact that both radial drift and grain growth produce a similar particle size segregation in a disk (such that larger particles are preferentially concentrated closer to the star) makes it difficult to disentangle a clear signature of drift alone. We highlight a new approach, by showing that radial drift leaves a distinctive 'fingerprint' in the dust surface density profile that is directly accessible to current observational facilities. Using an analytical framework for dust evolution, we demonstrate that the combined effects of drift and (viscous) gas drag naturally produce a sharp outer edge in the dust distribution (or, equivalently, a sharp decrease in the dust-to-gas mass ratio). This edge feature forms during the earliest phase in the evolution of disk solids, before grain growth in the outer disk has made much progress, and is preserved over longer timescales when both growth and transport effects are more substantial. The key features of these analytical models are reproduced in detailed numerical simulations, and are qualitatively consistent with recent millimeter-wave observations that find gas/dust size discrepancies and steep declines in dust continuum emission in the outer regions of protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

AN EMBEDDED ACTIVE NUCLEUS IN THE OH MEGAMASER GALAXY IRAS16399–0937

We present a multiwavelength study of the OH megamaser galaxy IRAS16399–0937, based on new Hubble Space Telescope (HST)/Advanced Camera for Surveys F814W and Hα+[N II] images and archive data from HST, Two Micron All Sky Survey, Spitzer, Herschel and the Very Large Array. This system has a double nucleus, whose northern (IRAS16399N) and southern (IRAS16399S) components have a projected separation of ∼6'' (3.4 kpc) and have previously been identified based on optical spectra as a low ionization nuclear emission line region (LINER) and starburst nucleus, respectively. The nuclei are embedded in a tidally distorted common envelope, in which star formation is mostly heavily obscured. The infrared spectrum is dominated by strong polycyclic aromatic hydrocarbon, but deep silicate and molecular absorption features are also present, and are strongest in the IRAS16399N nucleus. The 0.435-500 μm spectral energy distribution was fitted with a model including stellar, interstellar medium and active galactic nucleus (AGN) torus components using our new Markov Chain Monte Carlo code, CLUMPYDREAM. The results indicate that the IRAS16399N contains an AGN (L {sub bol} ∼ 10{sup 44} erg s{sup –1}) deeply embedded in a quasi-spherical distribution of optically thick clumps with a covering fraction ≈1. We suggest that these clumps are the source of the OHM emission in IRAS16399–0937. The high torus covering fraction precludes AGN photoionization as the origin of the LINER spectrum, however, the spectrum is consistent with shocks (v ∼ 100-200 km s{sup –1}). We infer that the ∼10{sup 8} M {sub ☉} black hole in IRAS16399N is accreting at a small fraction (∼1%) of its Eddington rate. The low accretion rate and modest nuclear star formation rates suggest that while the gas-rich major merger forming the IRAS16399–0937 system has triggered widespread star formation, the massive gas inflows expected from merger simulations have not yet fully developed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SPIRAL DENSITY WAVES IN M81. II. HYDRODYNAMIC SIMULATIONS OF THE GAS RESPONSE TO STELLAR SPIRAL DENSITY WAVES

The gas response to the underlying stellar spirals is explored for M81 using unmagnetized hydrodynamic simulations. Constrained within the uncertainty of observations, 18 simulations are carried out to study the effects of self-gravity and to cover the parameter space comprising three different sound speeds and three different arm strengths. The results are confronted with the data observed at wavelengths of 8 μm and 21 cm. In the outer disk, the ring-like structure observed in the 8 μm image is consistent with the response of cold neutral medium with an effective sound speed 7 km s{sup –1}. For the inner disk, the presence of spiral shocks can be understood as a result of 4:1 resonances associated with the warm neutral medium with an effective sound speed 19 km s{sup –1}. Simulations with a single effective sound speed alone cannot simultaneously explain the structures in the outer and inner disks. Instead this justifies the coexistence of cold and warm neutral media in M81. The anomalously high streaming motions observed in the northeast arm and the outward shifted turning points in the iso-velocity contours seen along the southwest arm are interpreted as signatures of interactions with companion galaxies. The level of simulated streaming motions narrows down the uncertainty of the observed arm strength toward larger amplitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

NUMERICAL SIMULATION OF FAST-MODE MAGNETOSONIC WAVES EXCITED BY PLASMOID EJECTIONS IN THE SOLAR CORONA

The Atmospheric Imaging Assembly instrument on board the Solar Dynamics Observatory has directly imaged the fast-propagating magnetosonic waves (FMWs) successively propagating outward along coronal magnetic funnels. In this study we perform a numerical investigation of the excitation of FMWs in the interchange reconnection scenario, with footpoint shearing flow being used to energize the system and drive the reconnection. The modeling results show that as a result of magnetic reconnection, the plasma in the current sheet is heated up by Joule dissipation to ∼10 MK and is ejected rapidly, developing the hot outflows. Meanwhile, the current sheet is torn into plasmoids, which are shot quickly both upward and downward. When the plasmoids reach the outflow regions, they impact and collide with the ambient magnetic field there, which consecutively launches FMWs. The FMWs propagate outward divergently away from the impact regions, with a phase speed of the Alfvén speed of ∼1000 km s{sup –1}. In the k – ω diagram of the Fourier wave power, the FMWs display a broad frequency distribution with a straight ridge that represents the dispersion relation. With the WKB approximation, at the distance of 15 Mm from the wave source region, we estimate the energy flux of FMWs to be E ∼ 7.0 × 10{sup 6} erg cm{sup –2} s{sup –1}, which is ∼50 times smaller than the energy flux related to the tube-channeled reconnection outflow. These simulation results indicate that energetically and dynamically the outflow is far more important than the waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Estimating CDM particle trajectories in the mildly non-linear regime of structure formation. Implications for the density field in real and redshift space

We obtain approximations for the CDM particle trajectories starting from Lagrangian Perturbation Theory. These estimates for the CDM trajectories result in approximations for the density in real and redshift space, as well as for the momentum density that are better than what standard Eulerian and Lagrangian perturbation theory give. For the real space density, we find that our proposed approximation gives a good cross-correlation ( > 95%) with the non-linear density down to scales almost twice smaller than the non-linear scale, and six times smaller than the corresponding scale obtained using linear theory. This allows for a speed-up of an order of magnitude or more in the scanning of the cosmological parameter space with N-body simulations for the scales relevant for the baryon acoustic oscillations. Possible future applications of our method include baryon acoustic peak reconstruction, building mock galaxy catalogs, momentum field reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Measuring global monopole velocities, one by one

We present an estimation of the average velocity of a network of global monopoles in a cosmological setting using large numerical simulations. In order to obtain the value of the velocity, we improve some already known methods, and present a new one. This new method estimates individual global monopole velocities in a network, by means of detecting each monopole position in the lattice and following the path described by each one of them. Using our new estimate we can settle an open question previously posed in the literature: velocity-dependent one-scale (VOS) models for global monopoles predict two branches of scaling solutions, one with monopoles moving at subluminal speeds and one with monopoles moving at luminal speeds. Previous attempts to estimate monopole velocities had large uncertainties and were not able to settle that question. Our simulations find no evidence of a luminal branch. We also estimate the values of the parameters of the VOS model. With our new method we can also study the microphysics of the complicated dynamics of individual monopoles. Finally we use our large simulation volume to compare the results from the different estimator methods, as well as to asses the validity of the numerical approximations made.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Symmetry-breaking instability of quadratic soliton bound states

We study both numerically and experimentally two-dimensional soliton bound states in quadratic media and demonstrate their symmetry-breaking instability. The experiment is performed in a potassium titanyl phosphate crystal in a type-II configuration. The bound state is generated by the copropagation of the antisymmetric fundamental beam locked in phase with the symmetrical second harmonic one. Experimental results are in good agreement with numerical simulations of the nonlinear wave equations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Determination of neutron spectra within the energy of 1 keV to 1 MeV by means of reactor dosimetry

The standard procedure for neutron reactor dosimetry is based on neutron irradiation of a target and its post-irradiation analysis by Gamma and/or X-ray spectrometry. Nowadays, the neutron spectra can be easily characterized for thermal and fast energies (respectively 0.025 eV and >1 MeV). In this work we propose a new target and an innovating post-irradiation technique of analysis in order to detect the neutron spectra within the energy of 1 keV to 1 MeV. This article will present the calculations performed for the selection of a suitable nuclear reaction and isotope, the results predicted by simulations, the irradiation campaign that is proposed and the post-irradiation technique of analysis. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Monte Carlo simulation of breast imaging using synchrotron radiation

Purpose: Synchrotron radiation (SR), being the brightest artificial source of x-rays with a very promising geometry, has raised the scientific expectations that it could be used for breast imaging with optimized results. The ''in situ'' evaluation of this technique is difficult to perform, mostly due to the limited available SR facilities worldwide. In this study, a simulation model for SR breast imaging was developed, based on Monte Carlo simulation techniques, and validated using data acquired in the SYRMEP beamline of the Elettra facility in Trieste, Italy. Furthermore, primary results concerning the performance of SR were derived. Methods: The developed model includes the exact setup of the SR beamline, considering that the x-ray source is located at almost 23 m from the slit, while the photon energy was considered to originate from a very narrow Gaussian spectrum. Breast phantoms, made of Perspex and filled with air cavities, were irradiated with energies in the range of 16-28 keV. The model included a Gd{sub 2}O{sub 2}S detector with the same characteristics as the one available in the SYRMEP beamline. Following the development and validation of the model, experiments were performed in order to evaluate the contrast resolution of SR. A phantom made of adipose tissue and filled with inhomogeneities of several compositions and sizes was designed and utilized to simulate the irradiation under conventional mammography and SR conditions. Results: The validation results of the model showed an excellent agreement with the experimental data, with the correlation for contrast being 0.996. Significant differences only appeared at the edges of the phantom, where phase effects occur. The initial evaluation experiments revealed that SR shows very good performance in terms of the image quality indices utilized, namely subject contrast and contrast to noise ratio. The response of subject contrast to energy is monotonic; however, this does not stand for contrast to noise ratio, since there is a range of optimal performance for SR (18-21 keV). In comparison to conventional mammography, SR shows improved subject contrast for energies lower than the mean energy of each spectrum. Conclusions: The comparison of the results of the two models, conventional and SR, proved that SR exhibits better performance in the majority of cases. The proposed simulation model offers the possibility to perform exhaustive search to evaluate the performance of SR in clinical applications such as breast imaging.

62 RADIOLOGY AND NUCLEAR MEDICINE↗