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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Proton decoupling and recoupling under double-nutation irradiation in solid-state NMR

The effect of {sup 1}H decoupling in magic-angle spinning solid-state NMR is studied under radiofrequency irradiation causing simultaneous nutations around a pair of orthogonal axes. Double-nutation with an arbitrary pair of nutation frequencies is implemented through modulation of the amplitude, phase, and frequency of the transmitting pulses. Similarity and difference of double-nutation decoupling and two-pulse phase-modulation decoupling schemes [A. E. Bennett, C. M. Rienstra, M. Auger, K. V. Lakshmi, and R. G. Griffin, J. Chem. Phys. 103, 6951–6958 (1995) and I. Scholz, P. Hodgkinson, B. H. Meier, and M. Ernst, J. Chem. Phys. 130, 114510 (2009)] are discussed. The structure of recoupling bands caused by interference of the {sup 1}H spin nutation with sample spinning is studied by both experiments and numerical simulations.

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

Influence of changes in initial conditions for the simulation of dynamic systems

Chaos theory is a field of study in mathematics, with applications in several disciplines including meteorology, sociology, physics, engineering, economics, biology, and philosophy. Chaos theory studies the behavior of dynamical systems that are highly sensitive to initial conditions—a paradigm popularly referred to as the butterfly effect. Small differences in initial conditions field widely diverging outcomes for such dynamical systems, rendering long-term prediction impossible in general. This happens even though these systems are deterministic, meaning that their future behavior is fully determined by their initial conditions, with no random elements involved. In this paperinfluence of changes in initial conditions will be presented for the simulation of Lorenz system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Theory of runaway electrons in ITER: Equations, important parameters, and implications for mitigation

The plasma current in ITER cannot be allowed to transfer from thermal to relativistic electron carriers. The potential for damage is too great. Before the final design is chosen for the mitigation system to prevent such a transfer, it is important that the parameters that control the physics be understood. Equations that determine these parameters and their characteristic values are derived. The mitigation benefits of the injection of impurities with the highest possible atomic number Z and the slowing plasma cooling during halo current mitigation to ≳40 ms in ITER are discussed. The highest possible Z increases the poloidal flux consumption required for each e-fold in the number of relativistic electrons and reduces the number of high energy seed electrons from which exponentiation builds. Slow cooling of the plasma during halo current mitigation also reduces the electron seed. Existing experiments could test physics elements required for mitigation but cannot carry out an integrated demonstration. ITER itself cannot carry out an integrated demonstration without excessive danger of damage unless the probability of successful mitigation is extremely high. The probability of success depends on the reliability of the theory. Equations required for a reliable Monte Carlo simulation are derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Ternary eutectic dendrites: Pattern formation and scaling properties

Extending previous work [Pusztai et al., Phys. Rev. E 87, 032401 (2013)], we have studied the formation of eutectic dendrites in a model ternary system within the framework of the phase-field theory. We have mapped out the domain in which two-phase dendritic structures grow. With increasing pulling velocity, the following sequence of growth morphologies is observed: flat front lamellae → eutectic colonies → eutectic dendrites → dendrites with target pattern → partitionless dendrites → partitionless flat front. We confirm that the two-phase and one-phase dendrites have similar forms and display a similar scaling of the dendrite tip radius with the interface free energy. It is also found that the possible eutectic patterns include the target pattern, and single- and multiarm spirals, of which the thermal fluctuations choose. The most probable number of spiral arms increases with increasing tip radius and with decreasing kinetic anisotropy. Our numerical simulations confirm that in agreement with the assumptions of a recent analysis of two-phase dendrites [Akamatsu et al., Phys. Rev. Lett. 112, 105502 (2014)], the Jackson-Hunt scaling of the eutectic wavelength with pulling velocity is obeyed in the parameter domain explored, and that the natural eutectic wavelength is proportional to the tip radius of the two-phase dendrites. Finally, we find that it is very difficult/virtually impossible to form spiraling two-phase dendrites without anisotropy, an observation that seems to contradict the expectations of Akamatsu et al. Yet, it cannot be excluded that in isotropic systems, two-phase dendrites are rare events difficult to observe in simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Two-dimensional studies of relativistic electron beam plasma instabilities in an inhomogeneous plasma

Relativistic electron beam propagation in plasma is fraught with several micro instabilities like two stream, filamentation, etc., in plasma. This results in severe limitation of the electron transport through a plasma medium. Recently, however, there has been an experimental demonstration of improved transport of Mega Ampere of electron currents (generated by the interaction of intense laser with solid target) in a carbon nanotube structured solid target [G. Chatterjee et al., Phys. Rev. Lett. 108, 235005 (2012)]. This then suggests that the inhomogeneous plasma (created by the ionization of carbon nanotube structured target) helps in containing the growth of the beam plasma instabilities. This manuscript addresses this issue with the help of a detailed analytical study and 2-D Particle-In-Cell simulations. The study conclusively demonstrates that the growth rate of the dominant instability in the 2-D geometry decreases when the plasma density is chosen to be inhomogeneous, provided the scale length 1/k{sub s} of the inhomogeneous plasma is less than the typical plasma skin depth (c/ω{sub 0}) scale. At such small scale lengths channelization of currents is also observed in simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Higher order mode of a microstripline fed cylindrical dielectric resonator antenna

A microstrip transmission line can be used to excite the broadside radiating mode of a cylindrical dielectric resonator antenna (CDRA). The same is found to excite considerably well a higher order mode (HOM) as well. However unlike the broadside mode, the higher order mode gives distorted radiation pattern which makes this mode less useful for practical applications. The cause of distortion in the HOM radiation and the dependence of HOM coupling on the microstrip feed line are explored using HFSS simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Micro-taper as focusing or scattering optical element

We consider micro-taper (narrow refractive axicon) as optical element which is focusing or scattering in dependence on axicon’s cone angle. The diffraction of laser beam by micro-taper is simulated by two methods: multiply internal ray reflections using geometrical approach and Helmholtz equation solving using finite elements method. Based on ray optics we derive analytic formulas for conical angles values which provide focusing or scattering features of micro-taper. Numerical simulation by finite elements method verifies theoretical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Energy optimization for a wind DFIG with flywheel energy storage

The type of distributed generation unit that is the subject of this paper relates to renewable energy sources, especially wind power. The wind generator used is based on a double fed induction Generator (DFIG). The stator of the DFIG is connected directly to the network and the rotor is connected to the network through the power converter with three levels. The objective of this work is to study the association a Flywheel Energy Storage System (FESS) in wind generator. This system is used to improve the quality of electricity provided by wind generator. It is composed of a flywheel; an induction machine (IM) and a power electronic converter. A maximum power tracking technique « Maximum Power Point Tracking » (MPPT) and a strategy for controlling the pitch angle is presented. The model of the complete system is developed in Matlab/Simulink environment / to analyze the results from simulation the integration of wind chain to networks.

17 WIND ENERGY↗
DOE OSTI2018

Physical mechanism of electron bunch generation by an ultrarelativistic-intensity laser pulse passing through a sharp plasma boundary

The physical mechanism of the interaction process of a laser pulse of ultrarelativistic intensity with semi-bounded plasma having a sharp boundary is studied in one-dimensional geometry. It is shown that electron bunches are generated by the laser pulse due to the multiflow motion of plasma electrons with crossing trajectories. An analytical relation that allows an electron bunch charge to be estimated as a function of plasma characteristics and electron trajectory parameters is derived and confirmed by the results of numerical simulations. (paper)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

CHARACTERISTICS OF SPIRAL ARMS IN LATE-TYPE GALAXIES

We have measured the positions of large numbers of H II regions in four nearly face-on, late-type, spiral galaxies: NGC 628 (M74), NGC 1232, NGC 3184, and NGC 5194 (M51). Fitting log-periodic spiral models to segments of each arm yields local estimates of spiral pitch angle and arm width. While pitch angles vary considerably along individual arms, among arms within a galaxy, and among galaxies, we find no systematic trend with galactocentric distance. We estimate the widths of the arm segments from the scatter in the distances of the H II regions from the spiral model. All major arms in these galaxies show spiral arm width increasing with distance from the galactic center, similar to the trend seen in the Milky Way. However, in the outermost parts of the galaxies, where massive star formation declines, some arms reverse this trend and narrow. We find that spiral arms often appear to be composed of segments of ∼5 kpc length, which join to form kinks and abrupt changes in pitch angle and arm width; these characteristics are consistent with properties seen in the large N-body simulations of D'Onghia et al. and others.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

CONTRIBUTION OF VELOCITY VORTICES AND FAST SHOCK REFLECTION AND REFRACTION TO THE FORMATION OF EUV WAVES IN SOLAR ERUPTIONS

We numerically study the detailed evolutionary features of the wave-like disturbance and its propagation in the eruption. This work is a follow-up to Wang et al., using significantly upgraded new simulations. We focus on the contribution of the velocity vortices and the fast shock reflection and refraction in the solar corona to the formation of the EUV waves. Following the loss of equilibrium in the coronal magnetic structure, the flux rope exhibits rapid motions and invokes the fast-mode shock at the front of the rope, which then produces a type II radio burst. The expansion of the fast shock, which is associated with outward motion, takes place in various directions, and the downward expansion shows the reflection and the refraction as a result of the non-uniform background plasma. The reflected component of the fast shock propagates upward and the refracted component propagates downward. As the refracted component reaches the boundary surface, a weak echo is excited. The Moreton wave is invoked as the fast shock touches the bottom boundary, so the Moreton wave lags the type II burst. A secondary echo occurs in the area where reflection of the fast shock encounters the slow-mode shock, and the nearby magnetic field lines are further distorted because of the interaction between the secondary echo and the velocity vortices. Our results indicate that the EUV wave may arise from various processes that are revealed in the new simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

STRAHL FORMATION IN THE SOLAR WIND ELECTRONS VIA WHISTLER INSTABILITY

This Letter puts forth a possible explanation for the formation of a solar wind strahl electron distribution function by means of local wave–particle interactions. A solar wind electron consists of a core and the hot “halo” electrons, which possess a net drift speed with respect to the core. According to the present model, pitch angle scattering of the initially isotropic drifting halo occurs when the enhanced whistler waves are excited by mildly anisotropic core electrons. The pitch angle scattering primarily affects the halo moving in the anti-sunward direction, resulting in pitch angle diffusion across the 90° gap through a nonlinear scattering process, and consequently leading to a reduction in the net drift speed of halo electrons. The remaining portion of the anti-sunward moving halo, which is not affected by pitch angle scattering, simply appears to form a field-aligned strahl in the electron velocity distribution. The present scenario of local generation of the strahl is demonstrated by the particle-in-cell simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SU-E-J-145: Validation of An Analytical Model for in Vivo Range Verification Using GATE Monte Carlo Simulation in Proton Therapy

Purpose: Predicted PET images on the basis of analytical filtering approach for proton range verification has been successful developed and validated using FLUKA Monte Carlo (MC) codes and phantom measurements. The purpose of the study is to validate the effectiveness of analytical filtering model for proton range verification on GATE/GEANT4 Monte Carlo simulation codes. Methods: In this study, we performed two experiments for validation of predicted β+-isotope by the analytical model with GATE/GEANT4 simulations. The first experiments to evaluate the accuracy of predicting β+-yields as a function of irradiated proton energies. In second experiment, we simulate homogeneous phantoms of different materials irradiated by a mono-energetic pencil-like proton beam. The results of filtered β+-yields distributions by the analytical model is compared with those of MC simulated β+-yields in proximal and distal fall-off ranges. Results: The results investigate the distribution between filtered β+-yields and MC simulated β+-yields distribution in different conditions. First, we found that the analytical filtering can be applied over the whole range of the therapeutic energies. Second, the range difference between filtered β+-yields and MC simulated β+-yields at the distal fall-off region are within 1.5mm for all materials used. The findings validated the usefulness of analytical filtering model on range verification of proton therapy on GATE Monte Carlo simulations. In addition, there is a larger discrepancy between filtered prediction and MC simulated β+-yields using GATE code, especially in proximal region. This discrepancy might Result from the absence of wellestablished theoretical models for predicting the nuclear interactions. Conclusion: Despite the fact that large discrepancies of the distributions between MC-simulated and predicted β+-yields were observed, the study prove the effectiveness of analytical filtering model for proton range verification using GATE Monte Carlo simulation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Application of Backward Differentiation Formulas to Neutron Kinetics Problems

The neutron kinetics equations belong to the class of stiff equations for numerical time integration schemes. In this work, the accuracy and speed of algorithms based on backward differentiation formulas (BDFs) are studied with regard to point and spatial neutron kinetics problems. Using a BDF algorithm with the automatic selection of time step and order, solutions to a number of model problems with both a positive and negative reactivity insertion are analyzed. Plots of numerical cost dependences on local and global errors are presented. The results indicate that algorithms based on BDFs are highly efficient and that their application to nuclear reactor simulation is justified.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Specificities of one-dimensional dissipative magnetohydrodynamics

One-dimensional dynamics of a plane slab of cold (β ≪ 1) isothermal plasma accelerated by a magnetic field is studied in terms of the MHD equations with a finite constant conductivity. The passage to the limit β → 0 is analyzed in detail. It is shown that, at β = 0, the character of the solution depends substantially on the boundary condition for the electric field at the inner plasma boundary. The relationship between the boundary condition for the pressure at β > 0 and the conditions for the electric field at β = 0 is found. The stability of the solution against one-dimensional longitudinal perturbations is analyzed. It is shown that, in the limit β → 0, the stationary solution is unstable if the time during which the acoustic wave propagates across the slab is longer than the time of magnetic field diffusion. The growth rate and threshold of instability are determined, and results of numerical simulation of its nonlinear stage are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Auxiliary-Field Quantum Monte Carlo Simulations of Strongly-Correlated Molecules and Solids

We propose a method of implementing projected wave functions for second-quantized auxiliary-field quantum Monte Carlo (AFQMC) techniques. The method is based on expressing the two-body projector as one-body terms coupled to binary Ising fields. To benchmark the method, we choose to study the two-dimensional (2D) one-band Hubbard model with repulsive interactions using the constrained-path MC (CPMC). The CPMC uses a trial wave function to guide the random walks so that the so-called fermion sign problem can be eliminated. The trial wave function also serves as the importance function in Monte Carlo sampling. As such, the quality of the trial wave function has a direct impact to the efficiency and accuracy of the simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

TOWARD MORE REALISTIC ANALYTIC MODELS OF THE HELIOTAIL: INCORPORATING MAGNETIC FLATTENING VIA DISTORTION FLOWS

Both physical arguments and simulations of the global heliosphere indicate that the tailward heliopause is flattened considerably in the direction perpendicular to both the incoming flow and the large-scale interstellar magnetic field. Despite this fact, all of the existing global analytical models of the outer heliosheath's magnetic field assume a circular cross section of the heliotail. To eliminate this inconsistency, we introduce a mathematical procedure by which any analytically or numerically given magnetic field can be deformed in such a way that the cross sections along the heliotail axis attain freely prescribed, spatially dependent values for their total area and aspect ratio. The distorting transformation of this method honors both the solenoidality condition and the stationary induction equation with respect to an accompanying flow field, provided that both constraints were already satisfied for the original magnetic and flow fields prior to the transformation. In order to obtain realistic values for the above parameters, we present the first quantitative analysis of the heliotail's overall distortion as seen in state-of-the-art three-dimensional hybrid MHD–kinetic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗