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

Improve earthquake hypocenter using adaptive simulated annealing inversion in regional tectonic, volcano tectonic, and geothermal observation

Observation of earthquakes is routinely used widely in tectonic activity observation, and also in local scale such as volcano tectonic and geothermal activity observation. It is necessary for determining the location of precise hypocenter which the process involves finding a hypocenter location that has minimum error between the observed and the calculated travel times. When solving this nonlinear inverse problem, simulated annealing inversion method can be applied to such global optimization problems, which the convergence of its solution is independent of the initial model. In this study, we developed own program codeby applying adaptive simulated annealing inversion in Matlab environment. We applied this method to determine earthquake hypocenter using several data cases which are regional tectonic, volcano tectonic, and geothermal field. The travel times were calculated using ray tracing shooting method. We then compared its results with the results using Geiger’s method to analyze its reliability. Our results show hypocenter location has smaller RMS error compared to the Geiger’s result that can be statistically associated with better solution. The hypocenter of earthquakes also well correlated with geological structure in the study area. Werecommend using adaptive simulated annealing inversion to relocate hypocenter location in purpose to get precise and accurate earthquake location.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A neural approach for the numerical modeling of two-dimensional magnetic hysteresis

This paper deals with a neural network approach to model magnetic hysteresis at macro-magnetic scale. Such approach to the problem seems promising in order to couple the numerical treatment of magnetic hysteresis to FEM numerical solvers of the Maxwell's equations in time domain, as in case of the non-linear dynamic analysis of electrical machines, and other similar devices, making possible a full computer simulation in a reasonable time. The neural system proposed consists of four inputs representing the magnetic field and the magnetic inductions components at each time step and it is trained by 2-d measurements performed on the magnetic material to be modeled. The magnetic induction B is assumed as entry point and the output of the neural system returns the predicted value of the field H at the same time step. A suitable partitioning of the neural system, described in the paper, makes the computing process rather fast. Validations with experimental tests and simulations for non-symmetric and minor loops are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Gyrokinetic studies of core turbulence features in ASDEX Upgrade H-mode plasmas

Gyrokinetic validation studies are crucial for developing confidence in the model incorporated in numerical simulations and thus improving their predictive capabilities. As one step in this direction, we simulate an ASDEX Upgrade discharge with the GENE code, and analyze various fluctuating quantities and compare them to experimental measurements. The approach taken is the following. First, linear simulations are performed in order to determine the turbulence regime. Second, the heat fluxes in nonlinear simulations are matched to experimental fluxes by varying the logarithmic ion temperature gradient within the expected experimental error bars. Finally, the dependence of various quantities with respect to the ion temperature gradient is analyzed in detail. It is found that density and temperature fluctuations can vary significantly with small changes in this parameter, thus making comparisons with experiments very sensitive to uncertainties in the experimental profiles. However, cross-phases are more robust, indicating that they are better observables for comparisons between gyrokinetic simulations and experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Chaotic attractors with separated scrolls

This paper proposes a new behavior of chaotic attractors with separated scrolls while combining Julia's process with Chua's attractor and Lorenz's attractor. The main motivation of this work is the ability to generate a set of separated scrolls with different behaviors, which in turn allows us to choose one or many scrolls combined with modulation (amplitude and frequency) for secure communication or synchronization. This set seems a new class of hyperchaos because each element of this set looks like a simple chaotic attractor with one positive Lyapunov exponent, so the cardinal of this set is greater than one. This new approach could be used to generate more general higher-dimensional hyperchaotic attractor for more potential application. Numerical simulations are given to show the effectiveness of the proposed theoretical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Accurate transport properties for H–CO and H–CO{sub 2}

Transport properties for collisions of hydrogen atoms with CO and CO{sub 2} have been computed by means of quantum scattering calculations. The carbon oxides are important species in hydrocarbon combustion. The following potential energy surfaces (PES’s) for the interaction of the molecule fixed in its equilibrium geometry were employed: for H–CO, the PES was taken from the work of Song et al. [J. Phys. Chem. A 117, 7571 (2013)], while the PES for H–CO{sub 2} was computed in this study by a restricted coupled cluster method that included single, double, and (perturbatively) triple excitations. The computed transport properties were found to be significantly different from those computed by the conventional approach that employs isotropic Lennard-Jones (12-6) potentials. The effect of using the presently computed accurate transport properties in 1-dimensional combustion simulations of methane-air flames was investigated.

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

Compton polarimeter for 10–30 keV x rays

We present a simple and versatile polarimeter for x rays in the energy range of 10–30 keV. It uses Compton scattering in low-Z materials such as beryllium or boron carbide. The azimuthal distribution of the scattered x rays is sampled by an array of 12 silicon PIN diodes operated at room temperature. We evaluated the polarimetry performance using Monte-Carlo simulations and show experimental results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Rubber elasticity for percolation network consisting of Gaussian chains

A theory describing the elastic modulus for percolation networks of Gaussian chains on general lattices such as square and cubic lattices is proposed and its validity is examined with simulation and mechanical experiments on well-defined polymer networks. The theory was developed by generalizing the effective medium approximation (EMA) for Hookian spring network to Gaussian chain networks. From EMA theory, we found that the ratio of the elastic modulus at p, G to that at p = 1, G{sub 0}, must be equal to G/G{sub 0} = (p − 2/f)/(1 − 2/f) if the position of sites can be determined so as to meet the force balance, where p is the degree of cross-linking reaction. However, the EMA prediction cannot be applicable near its percolation threshold because EMA is a mean field theory. Thus, we combine real-space renormalization and EMA and propose a theory called real-space renormalized EMA, i.e., REMA. The elastic modulus predicted by REMA is in excellent agreement with the results of simulations and experiments of near-ideal diamond lattice gels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Design of 3 GeV booster ring lattice

The aim of this study is to design of a 3 GeV booster ring for the 3 GeV storage ring. Electrons are needed to be accelerated to 3.0 GeV from 0.15 GeV energy. In this frame, we studied on two options for booster ring; a compact booster and the booster that shares the same tunnel with the storage ring. The lattice type has been chosen FODO for both options, lattice parameters are calculated, sextupole magnets are used to decrease dynamic aperture problem and dynamic aperture calculations are also made with considering of the necessary conditions. After designing and calculating of the parameters, these designs have been compared with each other. In addition to this comparison, these booster design parameters have been compared with some world centers design parameters and the reliability of the booster design is seen. Beam optics, OPA and Elegant simulation programs have been used in the study calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The importance of the cathode plume and its interactions with the ion beam in numerical simulations of Hall thrusters

Numerical simulations with a 2-D axisymmetric multi-fluid plasma code illustrate the significance of the near-plume interactions in investigations of the anomalous electron transport in Hall thrusters. In our simulations, the transport of electrons is modeled using an anomalous collision frequency, ν{sub anom}, yielding ν{sub anom} ≈ ω{sub ce} (i.e., the electron cyclotron frequency) in the near-plume region. We first show that restricting the anomalous collision frequency in this region to only within the ion beam, where the current density of ions is large, does not alter the plasma discharge in the Hall thruster as long as the interaction between the beam and the cathode plume is captured properly. These simulations suggest that electron transport occurs largely inside the beam. A second finding is on the significance of accounting for the ion acoustic turbulence (IAT), now known to occur in the vicinity of the cathode exit. We have included in our simulations a model of the IAT-driven anomalous collision frequency based on Sagdeev's model for saturation of the ion-acoustic instability. This implementation has allowed us to achieve excellent agreement with experimental measurements in the near plume of the H6 Hall thruster. Low frequency plasma oscillations similar in both magnitude and frequency to those found in the H6 thruster are recovered in our simulations when the model for the anomalous collision frequency in the cathode plume is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Impact of hyperbolicity on chimera states in ensembles of nonlocally coupled chaotic oscillators

In this work we analyse nonlocally coupled networks of identical chaotic oscillators. We study both time-discrete and time-continuous systems (Henon map, Lozi map, Lorenz system). We hypothesize that chimera states, in which spatial domains of coherent (synchronous) and incoherent (desynchronized) dynamics coexist, can be obtained only in networks of chaotic non-hyperbolic systems and cannot be found in networks of hyperbolic systems. This hypothesis is supported by numerical simulations for hyperbolic and non-hyperbolic cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Two dimensional kinetic analysis of electrostatic harmonic plasma waves

Electrostatic harmonic Langmuir waves are virtual modes excited in weakly turbulent plasmas, first observed in early laboratory beam-plasma experiments as well as in rocket-borne active experiments in space. However, their unequivocal presence was confirmed through computer simulated experiments and subsequently theoretically explained. The peculiarity of harmonic Langmuir waves is that while their existence requires nonlinear response, their excitation mechanism and subsequent early time evolution are governed by essentially linear process. One of the unresolved theoretical issues regards the role of nonlinear wave-particle interaction process over longer evolution time period. Another outstanding issue is that existing theories for these modes are limited to one-dimensional space. The present paper carries out two dimensional theoretical analysis of fundamental and (first) harmonic Langmuir waves for the first time. The result shows that harmonic Langmuir wave is essentially governed by (quasi)linear process and that nonlinear wave-particle interaction plays no significant role in the time evolution of the wave spectrum. The numerical solutions of the two-dimensional wave spectra for fundamental and harmonic Langmuir waves are also found to be consistent with those obtained by direct particle-in-cell simulation method reported in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Adaptive smart simulator for characterization and MPPT construction of PV array

Partial shading conditions are among the most important problems in large photovoltaic array. Many works of literature are interested in modeling, control and optimization of photovoltaic conversion of solar energy under partial shading conditions, The aim of this study is to build a software simulator similar to hard simulator and to produce a shading pattern of the proposed photovoltaic array in order to use the delivered information to obtain an optimal configuration of the PV array and construct MPPT algorithm. Graphical user interfaces (Matlab GUI) are built using a developed script, this tool is easy to use, simple, and has a rapid of responsiveness, the simulator supports large array simulations that can be interfaced with MPPT and power electronic converters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Numerical solution of 2D-vector tomography problem using the method of approximate inverse

We propose a numerical solution of reconstruction problem of a two-dimensional vector field in a unit disk from the known values of the longitudinal and transverse ray transforms. The algorithm is based on the method of approximate inverse. Numerical simulations confirm that the proposed method yields good results of reconstruction of vector fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

The preplasma effect on the properties of the shock wave driven by a fast electron beam

Strong shock wave generation by a mono-energetic fast electron beam in a plasma with an increasing density profile is studied theoretically. The proposed analytical model describes the shock wave characteristics for a homogeneous plasma preceded by a low density precursor. The shock pressure and the time of shock formation depend on the ratio of the electron stopping length to the preplasma areal density and on the initial energy of injected electrons. The conclusions of theoretical model are confirmed in numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

On the reliability of microvariability tests in quasars

Microvariations probe the physics and internal structure of quasars. Unpredictability and small flux variations make this phenomenon elusive and difficult to detect. Variance-based probes such as the C and F tests, or a combination of both, are popular methods to compare the light curves of the quasar and a comparison star. Recently, detection claims in some studies have depended on the agreement of the results of the C and F tests, or of two instances of the F-test, for rejecting the non-variation null hypothesis. However, the C-test is a non-reliable statistical procedure, the F-test is not robust, and the combination of tests with concurrent results is anything but a straightforward methodology. A priori power analysis calculations and post hoc analysis of Monte Carlo simulations show excellent agreement for the analysis of variance test to detect microvariations as well as the limitations of the F-test. Additionally, the combined tests yield correlated probabilities that make the assessment of statistical significance unworkable. However, it is possible to include data from several field stars to enhance the power in a single F-test, increasing the reliability of the statistical analysis. This would be the preferred methodology when several comparison stars are available. An example using two stars and the enhanced F-test is presented. These results show the importance of using adequate methodologies and avoiding inappropriate procedures that can jeopardize microvariability detections. Power analysis and Monte Carlo simulations are useful tools for research planning, as they can demonstrate the robustness and reliability of different research approaches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

CARMA OBSERVATIONS OF PROTOSTELLAR OUTFLOWS IN NGC 1333

We present observations of outflows in the star-forming region NGC 1333 using the Combined Array for Research in Millimeter-Wave Astronomy (CARMA). We combined the {sup 12}CO and {sup 13}CO (1-0) CARMA mosaics with data from the 14 m Five College Radio Astronomy Observatory to probe the central, most dense, and active region of this protostellar cluster at scales from 5'' to 7' (or 1000 AU to 0.5 pc at a distance of 235 pc). We map and identify {sup 12}CO outflows, and along with {sup 13}CO data we estimate their mass, momentum, and energy. Within the 7' Multiplication-Sign 7' map, the 5'' resolution allows for a detailed study of morphology and kinematics of outflows and outflow candidates, some of which were previously confused with other outflow emission in the region. In total, we identify 22 outflow lobes, as well as 9 dense circumstellar envelopes marked by continuum emission, of which 6 drive outflows. We calculate a total outflow mass, momentum, and energy within the mapped region of 6 M{sub Sun }, 19 M{sub Sun} km s{sup -1}, and 7 Multiplication-Sign 10{sup 44} erg, respectively. Within this same region, we compare outflow kinematics with turbulence and gravitational energy, and we suggest that outflows are likely important agents for the maintenance of turbulence in this region. In the earliest stages of star formation, outflows do not yet contribute enough energy to totally disrupt the clustered region where most star formation is happening, but have the potential to do so as the protostellar sources evolve. Our results can be used to constrain outflow properties, such as outflow strength, in numerical simulations of outflow-driven turbulence in clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

UPDATE ON THE CETUS POLAR STREAM AND ITS PROGENITOR

We trace the Cetus Polar Stream (CPS) with blue horizontal branch and red giant stars from Data Release 8 of the Sloan Digital Sky Survey. Using a larger data set than was available previously, we are able to refine the measured distance and velocity to this tidal debris star stream in the south Galactic cap. Assuming that the tidal debris traces the progenitor's orbit, we fit an orbit to the CPS and find that the stream is confined between ∼24 and 36 kpc on a rather polar orbit inclined 87° to the Galactic plane. The eccentricity of the orbit is 0.20, and the period is ∼700 Myr. If we instead matched N-body simulations to the observed tidal debris, these orbital parameters would change by 10% or less. The CPS stars travel in the opposite direction to those from the Sagittarius tidal stream in the same region of the sky. Through N-body models of satellites on the best-fitting orbit, and assuming that mass follows light, we show that the stream width, line-of-sight depth, and velocity dispersion imply a progenitor of ∼> 10{sup 8} M{sub ☉}. However, the density of stars along the stream requires either a disruption time on the order of one orbit or a stellar population that is more centrally concentrated than the dark matter. We suggest that an ultrafaint dwarf galaxy progenitor could reproduce a large stream width and velocity dispersion without requiring a very recent deflection of the progenitor into its current orbit. We find that most Cetus stars have metallicities of –2.5 < [Fe/H] <–2.0, similar to the observed metallicities of the ultrafaint dwarfs. Our simulations suggest that the parameters of the dwarf galaxy progenitors, including their dark matter content, could be constrained by observations of their tidal tails through comparison of the debris with N-body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Conservative GRMHD simulations of moderately thin, tilted accretion disks

This paper presents our latest numerical simulations of accretion disks that are misaligned with respect to the rotation axis of a Kerr black hole. In this work, we use a new, fully conservative version of the Cosmos++ general relativistic magnetohydrodynamics (GRMHD) code, coupled with an ad hoc cooling function designed to control the thickness of the disk. Together these allow us to simulate the thinnest tilted accretion disks ever using a GRMHD code. In this way, we are able to probe the regime where the dimensionless stress and scale height of the disk become comparable. We present results for both prograde and retrograde cases. The simulated prograde tilted disk shows no sign of Bardeen-Petterson alignment even in the innermost parts of the disk. The simulated retrograde tilted disk, however, does show modest alignment. The implication of these results is that the parameter space associated with Bardeen-Petterson alignment for prograde disks may be rather small, only including very thin disks. Unlike our previous work, we find no evidence for standing shocks in our simulated tilted disks. We ascribe this to the black hole spin, tilt angle, and disk scale height all being small in these simulations. We also add to the growing body of literature pointing out that the turbulence driven by the magnetorotational instability in global simulations of accretion disks is not isotropic. Finally, we provide a comparison between our moderately thin, untilted reference simulation and other numerical simulations of thin disks in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗