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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

OBSERVATIONAL SIGNATURES OF A MASSIVE DISTANT PLANET ON THE SCATTERING DISK

The orbital element distribution of trans-Neptunian objects (TNOs) with large pericenters has been suggested to be influenced by the presence of an undetected, large planet at >200 au from the Sun. To find additional observables caused by this scenario, we present here the first detailed emplacement simulation in the presence of a massive ninth planet on the distant Kuiper Belt. We perform 4 Gyr N -body simulations with the currently known solar system planetary architecture, plus a 10 M {sub ⊕} planet with similar orbital parameters to those suggested by Trujillo and Sheppard or Batygin and Brown, and 10{sup 5} test particles in an initial planetesimal disk. We find that including a distant super-Earth-mass planet produces a substantially different orbital distribution for the scattering and detached TNOs, raising the pericenters and inclinations of moderate semimajor axis (50 < a < 500 au) objects. We test whether this signature is detectable via a simulator with the observational characteristics of four precisely characterized TNO surveys. We find that the qualitatively very distinct solar system models that include a ninth planet are essentially observationally indistinguishable from an outer solar system produced solely by the four giant planets. We also find that the mass of the Kuiper Belt’s current scattering and detached populations is required to be 3–10 times larger in the presence of an additional planet. We do not find any evidence for clustering of orbital angles in our simulated TNO population. Wide-field, deep surveys targeting inclined high-pericenter objects will be required to distinguish between these different scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Time-resolved Rotational Velocities in the Upper Atmosphere of WASP-33 b

While steady empirical progress has been made in understanding the structure and composition of hot-planet atmospheres, direct measurements of velocity signatures, including winds, rotation, and jets, have lagged behind. Quantifying atmospheric dynamics of hot planets is critical for a complete understanding of their atmospheres, and such measurements may even illuminate other planetary properties, such as magnetic field strengths. In this manuscript we present the first detection of the Balmer lines Hα and Hβ in the atmosphere of the ultra-hot Jupiter WASP-33 b. Using atmospheric models that include the effects of atmospheric dynamics, we show that the shape of the average Balmer line transmission spectrum is consistent with rotational velocities in the planet’s thermosphere of v{sub rot}=10.1{sub −1.0}{sup +0.8} km s{sup −1}. We also measure a low-significance day- to nightside velocity shift of −4.6{sub −3.4}{sup +3.4} km s{sup −1} in the transmission spectrum, which is naturally explained by a global wind across the planet’s terminator. In a separate analysis the time-resolved velocity centroids of individual transmission spectra show unambiguous evidence of rotation, with a best-fit velocity of 10.0{sub −2.0}{sup +2.4} km s{sup −1}, consistent with the value of v {sub rot} derived from the shape of the average Balmer line transmission spectrum. Our observations and analysis confirm the power of time-resolved transmission spectra with a high signal-to-noise ratio to measure the velocity structures in exoplanet atmospheres. The high rotational and wind velocities we measure highlight the need for more detailed 3D global climate simulations of the rarefied upper atmospheres of ultra-hot gas giants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

X-Ray Spectro-polarimetry with Photoelectric Polarimeters

We derive a generalization of forward fitting for X-ray spectroscopy to include linear polarization of X-ray sources, appropriate for the anticipated next generation of space-based photoelectric polarimeters. We show that the inclusion of polarization sensitivity requires joint fitting to three observed spectra, one for each of the Stokes parameters, I(E), U(E), and Q(E). The equations for Stokes’ I(E) (the total intensity spectrum) are identical to the familiar case with no polarization sensitivity, and for which the model-predicted spectrum is obtained by a convolution of the source spectrum, F(E′), with the familiar energy response function, ϵ(E′) R(E′, E), where ϵ(E′) and R(E′, E) are the effective area and energy redistribution matrix, respectively. In addition to the energy spectrum, the two new relations for U(E) and Q(E) include the source polarization fraction and position angle versus energy, a(E′), and ψ{sub 0}{sup ′}(E{sup ′}), respectively, and the model-predicted spectra for these relations are obtained by a convolution with the “modulated” energy response function, μ(E′) ϵ (E′)R(E, E′), where μ(E′) is the energy-dependent modulation fraction that quantifies a polarimeter’s angular response to 100% polarized radiation. We present results of simulations with response parameters appropriate for the proposed PRAXyS Small Explorer observatory to illustrate the procedures and methods, and we discuss some aspects of photoelectric polarimeters with relevance to understanding their calibration and operation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Triggering Process of the X1.0 Three-ribbon Flare in the Great Active Region NOAA 12192

The solar magnetic field in a flare-producing active region (AR) is much more complicated than theoretical models, which assume a very simple magnetic field structure. The X1.0 flare, which occurred in AR 12192 on 2014 October 25, showed a complicated three-ribbon structure. To clarify the trigger process of the flare and to evaluate the applicability of a simple theoretical model, we analyzed the data from Hinode /Solar Optical Telescope and the Solar Dynamics Observatory /Helioseismic and Magnetic Imager, Atmospheric Imaging Assembly. We investigated the spatio-temporal correlation between the magnetic field structures, especially the non-potentiality of the horizontal field, and the bright structures in the solar atmosphere. As a result, we determined that the western side of the positive polarity, which is intruding on a negative polarity region, is the location where the flare was triggered. This is due to the fact that the sign of the magnetic shear in that region was opposite that of the major shear of the AR, and the significant brightenings were observed over the polarity inversion line (PIL) in that region before flare onset. These features are consistent with the recently proposed flare-trigger model that suggests that small reversed shear (RS) magnetic disturbances can trigger solar flares. Moreover, we found that the RS field was located slightly off the flaring PIL, contrary to the theoretical prediction. We discuss the possibility of an extension of the RS model based on an extra numerical simulation. Our result suggests that the RS field has a certain flexibility for displacement from a highly sheared PIL, and that the RS field triggers more flares than we expected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Photospheric Velocity Structures during the Emergence of Small Active Regions on the Sun

We study the plasma flows in the solar photosphere during the emergence of two small active regions, NOAA 9021 and 10768. Using Solar and Heliospheric Observatory /Michelson Doppler Imager data, we find that the strong plasma upflows appear at the initial stage of active region formation, with maximum upflow velocities of −1650 and −1320 m s{sup −1}. The structures with enhanced upflows have size ∼8 Mm in diameter, and they exist for 1–2 hr. The parameters of the enhanced upflows are consistent with those of the large active region NOAA 10488, which may suggest the possibility that the elementary emerging magnetic loops that appear at the earliest phase of active region formation have similar properties, irrespective of scales of active regions. Comparison between the observations and a numerical simulation of magnetic flux emergence shows a striking consistency. We find that the driving force of the plasma upflow is at first the gas pressure gradient and later the magnetic pressure gradient.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Experimental study on the guillotine break accident of the main steam line with Atlas

A main steam line break (MSLB) accident is defined as a pipe break in the main steam system of a nuclear power plant. The MSLB accident is characterized as cool-down event due to the increased steam flow rate, which causes excessive energy removal from the steam generators and eventually reactor coolant system. This results in a decrease in reactor coolant temperature. The cool-down causes an increase in core reactivity due to the negative moderator and Doppler reactivity coefficients. The thermalhydraulic integral effect test was performed to simulate a double-ended guillotine break accident of the main steam line of a pressurized water reactor. With an aim of simulating the MSLB accident as realistically as possible, a pertinent scaling approach was taken, especially focusing on a break flow. The main objectives of this test were not only to provide physical insight into the system response of a nuclear power plant during the MSLB accident but also to produce an integral effect test data to validate the safety analysis codes. In the present test, asymmetric thermal-hydraulic behavior was clearly observed between the loops. In the cold leg-1A and -1B where the affected steam generator is connected, the flow rates were higher and the fluid temperatures were lower than those in the cold leg-2A and -2B. This could be attributed to the excessive heat removal in the cold leg-1A and -1B through the affected steam generator. This integral effect test data will be used to evaluate the prediction capability of existing safety analysis codes. Furthermore, this data can be utilized to identify any code deficiency for a MSLB simulation, especially for DVI-adapted plants. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Measurement Concept for a Possible Clearance of Mercury Waste from Nuclear Facilities - 19138

This paper describes the development of a measurement concept for a clearance of mercury waste from the decommissioning of nuclear facilities. A decontamination of mercury is performed after which large parts of the waste could possibly be released into reuse or conventional mercury disposal. The focus of this work is on the development and validation of a gamma-spectrometric measurement setup that is used for decision-making measurements for a possible clearance procedure. The measurement system is validated by a comparison of experimental measurements with nuclear simulations using Monte Carlo N-Particle{sup R} Transport Code (MCNP{sup R}). The validated MCNP{sup R} simulation model is used to determine the energy-dependent photopeak efficiencies and subsequently the sensitivity of the measurements. Preliminary results of gamma spectrometry in combination with MCNP{sup R} simulations indicate that a clearance of elemental mercury should be possible after decontamination. (authors)

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

Benefits of digitalizing and employing simulation to increase plant system performance and ensure compliance with technical specifications - 109

The essential services water system in a BWR nuclear power plant was affected by corrosion, soiling and blockages due to ageing. Those episodes affected safety margins hindering the Technical Specifications compliance during periodic system testing. Tecnatom designed, implemented and tested a design modification for the on-line monitoring of this system by integrating it in the existing Digital Control System (DCS). Key system parameters were acquired in real time to be displayed in the Human System Interface, used in performing calculations and stored their historical evolution. The objective was to optimize monitoring and surveillance of the essential service water system. Furthermore, Tecnatom developed an engineering simulator ('What if' simulator) consisting of a hydraulic model of the system. This simulator takes as inputs each heat exchanger performance parameters from the plant DCS, allows the user to change the essential services water system configuration (valves position, UHS level...) and calculates theoretical process values predicting the system real behavior. In a second stage, more system instruments were wired to the DCS for Technical Specifications compliance automatic surveillance regarding opening and closing times of system valves. Operating experience shows important benefits after the digitalization. It enhances on-line information of system operation and automates calculation of load loss factor of each heat exchanger. It eases the evaluation of blockages and soiling in the system allowing optimization of cleaning proceedings. System Surveillance Requirements and Test Requirements reports have been automated optimizing the operation workload. The simulator use enables to identify optimal system configuration increasing safety margins of the system and ensuring Technical Specifications fulfillment. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

The Virtual Environment for Reactor Applications (VERA): Design and architecture

VERA, the Virtual Environment for Reactor Applications, is the system of physics capabilities being developed and deployed by the Consortium for Advanced Simulation of Light Water Reactors (CASL). CASL was established for the modeling and simulation of commercial nuclear reactors. VERA consists of integrating and interfacing software together with a suite of physics components adapted and/or refactored to simulate relevant physical phenomena in a coupled manner. VERA also includes the software development environment and computational infrastructure needed for these components to be effectively used. We describe the architecture of VERA from both software and numerical perspectives, along with the goals and constraints that drove major design decisions, and their implications. We explain why VERA is an environment rather than a framework or toolkit, why these distinctions are relevant (particularly for coupled physics applications), and provide an overview of results that demonstrate the use of VERA tools for a variety of challenging applications within the nuclear industry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Modulational instability and soliton trains in a model for two-mode fiber ring lasers

A model of orthogonally polarized two-field fiber ring laser with a linear gain is considered, with emphasis on the continuous-wave stability and the existence of soliton trains. The continuous-wave stability analysis is carried out within the framework of the modulational-instability approach, the variations of the gain spectrum with the modulation frequency and characteristic parameters of the model give rise to a rich variety of stability features including single-band and multiband stability regions. Seeking for pulse structures of the model, the two coupled cubic complex Ginzburg–Landau equations describing individual mode propagations are transformed into a set of coupled, first-order nonlinear ordinary-differential equations for the amplitudes and phases of the two modes. Numerical simulations of the last set of coupled equations indicate that in the anomalous dispersion regime, envelopes of the two fields are periodic trains of pulses the amplitudes of which are affected by the linear gain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Periodically driven ergodic and many-body localized quantum systems

We study dynamics of isolated quantum many-body systems whose Hamiltonian is switched between two different operators periodically in time. The eigenvalue problem of the associated Floquet operator maps onto an effective hopping problem. Using the effective model, we establish conditions on the spectral properties of the two Hamiltonians for the system to localize in energy space. We find that ergodic systems always delocalize in energy space and heat up to infinite temperature, for both local and global driving. In contrast, many-body localized systems with quenched disorder remain localized at finite energy. We support our conclusions by numerical simulations of disordered spin chains. We argue that our results hold for general driving protocols, and discuss their experimental implications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Isentropic expansion of copper plasma in Mbar pressure range at “Luch” laser facility

We present experimental results on thermodynamic properties of dense copper plasma in Mbar pressure range. The laser facility “Luch” with laser intensity 10{sup 14} W/cm{sup 2} is used to compress copper up to ∼8 Mbar by a strong shock wave; subsequent expansion of copper plasma into Al, Ti, Sn allows us to obtain release isentropes of copper by the impedance–matching method. A theoretical analysis and quantum simulations show that in our experiments strongly coupled quantum plasma is generated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Stationary solid particle attractors in standing waves

The present analysis extends earlier theories on patterns formed by the spontaneous accumulation and ordering of solid particles in certain types of flow by considering the case in which the particle carrier flow has the typical features of a “standing wave.” For the first time an explanation for this phenomenon is elaborated through arguments based on the interplay between vorticity and wave-interference dynamics (following a deductive approach after the so-called phase-locking or “resonance” model originally introduced by Pushkin et al. [Phys. Rev. Lett. 106, 234501 (2011)] and later variants developed by Lappa [Phys. Fluids 25(1), 012101 (2013) and Lappa, Chaos 23(1), 013105 (2013)]). The results of dedicated numerical simulations are used in synergy with available experimental work. An interesting analogy is proposed with the famous Chladni's series of experiments on patterns formed by sand on vibrating plates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Simulation of dilute polymeric fluids in a three-dimensional contraction using a multiscale FENE model

The multiscale FENE model is applied to a 3D square-square contraction flow problem. For this purpose, the stochastic Brownian configuration field method (BCF) has been coupled with our fully parallelized three-dimensional Navier-Stokes solver NaSt3DGPF. The robustness of the BCF method enables the numerical simulation of high Deborah number flows for which most macroscopic methods suffer from stability issues. The results of our simulations are compared with that of experimental measurements from literature and show a very good agreement. In particular, flow phenomena such as a strong vortex enhancement, streamline divergence and a flow inversion for highly elastic flows are reproduced. Due to their computational complexity, our simulations require massively parallel computations. Using a domain decomposition approach with MPI, the implementation achieves excellent scale-up results for up to 128 processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Adsorption of PTCDA on Si(001) − 2 × 1 surface

Adsorption structures of the 3,4,9,10-perylenetetracarboxylic dianhydride (PTCDA) molecule on the clean Si(001) − 2 × 1 surface were investigated using scanning tunneling microscopy (STM) experiments in conjunction with first principles theoretical calculations. Four dominant adsorption structures were observed in the STM experiments and their atomic coordinates on the Si(001) surface were determined by comparison between the experimental STM images and the theoretical simulations. Maximizing the number of the Si—O bonds is more crucial than that of the Si—C bonds in the PTCDA adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Time delay and noise explaining the behaviour of the cell growth in fermentation process

This paper proposes to investigate the interplay between time delay and external noise in explaining the behaviour of the microbial growth in batch fermentation process. Time delay and noise are modelled jointly via stochastic delay differential equations (SDDEs). The typical behaviour of cell concentration in batch fermentation process under this model is investigated. Milstein scheme is applied for solving this model numerically. Simulation results illustrate the effects of time delay and external noise in explaining the lag and stationary phases, respectively for the cell growth of fermentation process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Heteroclinic tangle phenomena in nanomagnets subject to time-harmonic excitations

Magnetization dynamics in uniformly magnetized nanomagnets excited by time-harmonic (AC) external fields or spin-polarized injected currents is considered. The analysis is focused on the behaviour of the AC-excited dynamics near saddle equilibria. It turns out that this dynamics has a chaotic character at moderately low power level. This chaotic and fractal nature is due to the phenomenon of heteroclinic tangle which is produced by the combined effect of AC-excitations and saddle type dynamics. By using the perturbation technique based on Melnikov function, analytical formulas for the threshold AC excitation amplitudes necessary to create the heteroclinic tangle are derived. Both the cases of AC applied fields and AC spin-polarized injected currents are treated. Then, by means of numerical simulations, we show how heteroclinic tangle is accompanied by the erosion of the safe basin around the stable regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Coalescence-controlled and coalescence-free growth regimes during deposition of pulsed metal vapor fluxes on insulating surfaces

The morphology and physical properties of thin films deposited by vapor condensation on solid surfaces are predominantly set by the processes of island nucleation, growth, and coalescence. When deposition is performed using pulsed vapor fluxes, three distinct nucleation regimes are known to exist depending on the temporal profile of the flux. These regimes can be accessed by tuning deposition conditions; however, their effect on film microstructure becomes marginal when coalescence sets in and erases morphological features obtained during nucleation. By preventing coalescence from being completed, these nucleation regimes can be used to control microstructure evolution and thus access a larger palette of film morphological features. Recently, we derived the quantitative criterion to stop coalescence during continuous metal vapor flux deposition on insulating surfaces—which typically yields 3-dimensional growth—by describing analytically the competition between island growth by atomic incorporation and the coalescence rate of islands [Lü et al., Appl. Phys. Lett. 105, 163107 (2014)]. Here, we develop the analytical framework for entering a coalescence-free growth regime for metal vapor deposition on insulating substrates using pulsed vapor fluxes, showing that there exist three distinct criteria for suppressing coalescence that correspond to the three nucleation regimes of pulsed vapor flux deposition. The theoretical framework developed herein is substantiated by kinetic Monte Carlo growth simulations. Our findings highlight the possibility of using atomistic nucleation theory for pulsed vapor deposition to control morphology of thin films beyond the point of island density saturation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗