Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed 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 793 records · Page 44

Simulation of Inrush Dynamics for Unbalanced Distribution Systems using Phasor Dynamics

Microgrids are a possible solution to mitigate the extensive interruptions in service to end-use customers caused by abnormal events, such as major natural disasters. Microgrids can be formed by energizing portions of the distribution system to interconnect generator and critical loads. In this scenario, inrush dynamic currents can present a barrier for forming microgrids by preventing the necessary switching operations and/or damaging equipment. Modeling inrush dynamics, through electromagnetic transient tools, presents a challenge when the microgrid represents a large portion of a distribution feeder. This paper uses a phasor dynamics approach to approximate the envelope of inrush dynamics, thus allowing simulations of large systems. The two main contributions of this paper are: a) a new formulation to solve a set of phasor-dynamic differential equations representing an unbalanced distribution system; and b) modeling of transformer saturation in the phasor-dynamics frame. These two contributions allow the in-rush phenomena to be simulated on full size microgrid and/or distribution feeder models. The new formulation and model are implemented in GridLAB-D, and tested on IEEE test systems.

inrush, GridLAB-D, distribution analysis, microgri↗

EGS Collab Testbed 1: Temperature profile

This submission includes an input file, plot file, Fortran conversion file, and 3D data file from the simulation of the temperature profile within the Test Bed #1 of the EGS Collab project. The simulation was executed with PNNL's STOMP-GT simulator, which reads the input file, and produces the STOMP-GT Simulated 2D Temperature Distribution Data (plot.001302 file). The Fortran conversion file, converts the STOMP-GT Simulated 2D Temperature Distribution Data (plot.001302) 2D results in local coordinates to 3D results in Homestake coordinates.

15 GEOTHERMAL ENERGY↗

GenASiS Basics: Object-oriented utilitarian functionality for large-scale physics simulations (Version 4)

GenASiS Basics provides modern Fortran classes furnishing extensible object-oriented utilitarian functionality for large-scale physics simulations on distributed memory supercomputers. This functionality includes physical units and constants; display to the screen or standard output device; message passing; I/O to disk; and runtime parameter management and usage statistics. Herein, this revision—Version 4 of Basics—includes a name change and additions to functionality, including the facilitation of direct communication between GPUs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THE CENTRAL SLOPE OF DARK MATTER CORES IN DWARF GALAXIES: SIMULATIONS VERSUS THINGS

We make a direct comparison of the derived dark matter (DM) distributions between hydrodynamical simulations of dwarf galaxies assuming a {Lambda}CDM cosmology and the observed dwarf galaxies sample from the THINGS survey in terms of (1) the rotation curve shape and (2) the logarithmic inner density slope {alpha} of mass density profiles. The simulations, which include the effect of baryonic feedback processes, such as gas cooling, star formation, cosmic UV background heating, and most importantly, physically motivated gas outflows driven by supernovae, form bulgeless galaxies with DM cores. We show that the stellar and baryonic mass is similar to that inferred from photometric and kinematic methods for galaxies of similar circular velocity. Analyzing the simulations in exactly the same way as the observational sample allows us to address directly the so-called cusp/core problem in the {Lambda}CDM model. We show that the rotation curves of the simulated dwarf galaxies rise less steeply than cold dark matter rotation curves and are consistent with those of the THINGS dwarf galaxies. The mean value of the logarithmic inner density slopes {alpha} of the simulated galaxies' DM density profiles is {approx}-0.4 {+-} 0.1, which shows good agreement with {alpha} = -0.29 {+-} 0.07 of the THINGS dwarf galaxies. The effect of non-circular motions is not significant enough to affect the results. This confirms that the baryonic feedback processes included in the simulations are efficiently able to make the initial cusps with {alpha} {approx}-1.0 to -1.5 predicted by DM-only simulations shallower and induce DM halos with a central mass distribution similar to that observed in nearby dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved Collective Thomson Scattering measurements of fast ions at ASDEX upgrade

Understanding the behaviour of the confined fast ions is important in both current and future fusion experiments. These ions play a key role in heating the plasma and will be crucial for achieving conditions for burning plasma in next-step fusion devices. Microwave-based Collective Thomson Scattering (CTS) is well suited for reactor conditions and offers such an opportunity by providing measurements of the confined fast-ion distribution function resolved in space, time and 1D velocity space. We currently operate a CTS system at ASDEX Upgrade using a gyrotron which generates probing radiation at 105 GHz. A new setup using two independent receiver systems has enabled improved subtraction of the background signal, and hence the first accurate characterization of fast-ion properties. Here we review this new dual-receiver CTS setup and present results on fast-ion measurements based on the improved background characterization. These results have been obtained both with and without NBI heating, and with the measurement volume located close to the centre of the plasma. The measurements agree quantitatively with predictions of numerical simulations. Hence, CTS studies of fast-ion dynamics at ASDEX Upgrade are now feasible. The new background subtraction technique could be important for the design of CTS systems in other fusion experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermal pressure in the cold neutral medium of nearby galaxies

Dynamic and thermal processes regulate the structure of the multiphase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure (P{sub th}) of the diffuse, neutral gas. By combining [C ii] 158 μm, H i, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of ∼1 kiloparsec. We find a distribution of thermal pressures in the P{sub th}/k∼10{sup 3}--10{sup 5} K cm{sup −3} range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of P{sub th}/k ≈ 3600 K cm{sup −3}. These results are consistent with thermal pressure measurements using other observational methods. We find that P{sub th} increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral media could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about ∼30% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

79 ASTRONOMY AND ASTROPHYSICS↗

Mechanical Properties and Strain-Induced Martensite Transformation in Cold Rolling of 304L Stainless Steel Plate

In this work, cold rolling behavior of stainless steel 304L is investigated and the effects of different process parameters on the occurrence of strain-induced martensite and mechanical properties are studied. The rolling experiments are conducted under different rolling speeds and reductions in which a set of samples is deformed at room temperature and the other set is first cooled to the temperature of − 10 °C and then rolled. Afterward, the developed microstructures and mechanical properties of the rolled steel are evaluated employing different testing techniques. In order to justify the results, mathematical modeling of cold rolling operation is also performed using Abaqus/Explicit to estimate the distribution of equivalent plastic strain and temperature rise during rolling. The results show that strain-induced martensite may occur during cold rolling and it efficiently alters microstructure and mechanical performance of the rolled plate. The amount of reduction, initial temperature, rolling layout, and, to a lesser extent, rolling speed can affect the kinetics of martensite transformation as well as microstructure and mechanical properties of the rolled steel.

36 MATERIALS SCIENCE↗

SN 2013ej: A type IIL supernova with weak signs of interaction

We present optical photometric and spectroscopic observations of SN 2013ej. It is one of the brightest Type II supernovae (SNe II) exploded in a nearby (∼10 Mpc) galaxy, NGC 628. The light-curve characteristics are similar to SNe II, but with a relatively shorter (∼85 days) and steeper (∼1.7 mag (100 days){sup −1} in V) plateau phase. The SN shows a large drop of 2.4 mag in V-band brightness during the plateau-to-nebular transition. The absolute ultraviolet (UV) light curves are identical to SN 2012aw, showing a similar UV-plateau trend extending up to 85 days. The radioactive {sup 56}Ni mass estimated from the tail luminosity is 0.02 M{sub ⊙}, which is significantly lower than typical SNe IIP. The characteristics of spectral features and evolution of line velocities indicate that SN 2013ej is a Type II event. However, light-curve characteristics and some spectroscopic features provide strong support in classifying it as a Type IIL event. A detailed synow modeling of spectra indicates the presence of some high-velocity components in Hα and Hβ profiles, implying a possible ejecta–circumstellar medium interaction. The nebular phase spectrum shows an unusual notch in the Hα emission, which may indicate bipolar distribution of {sup 56}Ni. Modeling of the bolometric light curve yields a progenitor mass of ∼14 M{sub ⊙} and a radius of ∼450 R{sub ⊙}, with a total explosion energy of ∼2.3×10{sup 51} erg.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating the Feasibility of Isolated Pancreatic Perfusion for Chemotherapy Using Computed Tomography: An Experimental Study in Pig Models

PurposePercutaneous isolated pancreatic perfusion (PIPP) is performed along with interventional radiology techniques to obtain high drug concentration by occluding the arterial inlet and venous outlet of the pancreas. The experimental study aimed to evaluate the contrast distribution in PIPP under different flow rates with or without anterior mesenteric artery (AMA) occlusion.Materials and MethodsThis study was approved by a local animal experiment ethics committee. Nine pigs were divided into Groups 1, 2, and 3, by infusion rates of 12, 24, and 36 mL/min. Groups 4 and 5 (3 pigs each) and Group 6 (2 pigs) underwent PIPP at the same respective infusion rates with and without AMA occlusion. Computed tomography (CT) arteriography was performed during PIPP with nonionic contrast media. The enhanced volume was calculated by adding the enhanced area in each slice using 1.25-mm axial images. The percent enhanced volume to the whole pancreas (%eV) was used to simulate drug distribution; the result was compared among groups.ResultsWithout AMA occlusion, a larger %eV was obtained with high infusion rates (P = 0.039). The median %eV in Groups 1, 2, and 3 were 57.7, 74.2, and 90.5%, respectively. With AMA occlusion, CT demonstrated duodenal enhancement at an infusion rate of 36 mL/min, and the median %eV in Groups 4, 5, and 6 were 92.8, 95.4, and 98.5%, respectively. A significantly larger %eV was obtained after AMA occlusion (P = 0.031).ConclusionA higher infusion rate or AMA occlusion increases the enhanced volume in PIPP in pig models. Level of evidence: No level of evidence.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Distribution of annihilation luminosities in dark matter substructure

We calculate the probability distribution function (PDF) of the expected annihilation luminosities of dark matter subhalos as a function of subhalo mass and distance from the Galactic center using a semianalytical model of halo evolution. We find that the PDF of luminosities is relatively broad, exhibiting a spread of as much as an order of magnitude at fixed subhalo mass and halo-centric distance. The luminosity PDF allows for simple construction of mock samples of {gamma}-ray luminous subhalos and assessment of the variance in the predicted {gamma}-ray signals from dark matter annihilation. Other applications include quantifying the variance among the expected luminosities of dwarf spheroidal galaxies, assessing the level at which dark matter annihilation can be a contaminant in the expected {gamma}-ray signal from other astrophysical sources, as well as estimating the level at which nearby subhalos can contribute to the antimatter flux.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The physical properties of Li-doped g-C{sub 3}N{sub 4} monolayer sheet investigated by the first-principles

Highlights: • Systematically research on Li-doped g-C{sub 3}N{sub 4} monolayer sheets by first-principles calculation. • Optimal dopant concentration for optical absorption is 7.12%. • Thermodynamics stability of the doped substrate g-C{sub 3}N{sub 4} decreased with Li dopant concentration increasing. • The values of work function Φ decreased monotonously with the increasing of Li dopant concentration. - Abstract: The geometric, electronic, optical properties, thermodynamic stability, and work function of Li-doped g-C{sub 3}N{sub 4} monolayer were investigated by the first-principles calculation. It was found that the Li atoms were preferentially substituted the open-hollow sites of g-C{sub 3}N{sub 4}. Interestingly, the “odd” number of Li doped g-C{sub 3}N{sub 4} showed metallic properties, while the “even” number of Li atoms widened the band gap of g-C{sub 3}N{sub 4}. The HOMO and LUMO distributions reveal that the active sites located at edge N and C atoms for both pristine and the Li-doped g-C{sub 3}N{sub 4}. In addition, thermodynamic analysis showed that the doped Li atoms reduced the thermodynamic stability of g-C{sub 3}N{sub 4} monolayer sheets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

STELLAR POPULATION SYNTHESIS BASED MODELING OF THE MILKY WAY USING ASTEROSEISMOLOGY OF 13,000 KEPLER RED GIANTS

With current space-based missions it is now possible to obtain age-sensitive asteroseismic information for tens of thousands of red giants. This provides a promising opportunity to study the Galactic structure and evolution. We use asteroseismic data of red giants, observed by Kepler, to test the current theoretical framework of modeling the Galaxy based on population synthesis modeling and the use of asteroseismic scaling relations for giants. We use the open source code Galaxia to model the Milky Way and find the distribution of the masses predicted by Galaxia to be systematically offset with respect to the seismically inferred observed masses. The Galactic model overestimates the number of low-mass stars, and these stars are predominantly old and of low metallicity. Using corrections to the Δν scaling relation suggested by stellar models significantly reduces the disagreement between predicted and observed masses. For a few cases where non-seismic mass estimates are available, the corrections to Δν also improve the agreement between seismic and non-seismic mass estimates. The disagreement between predictions of the Galactic model and the observations is most pronounced for stars with [Fe/H]<−0.5 and [Fe/H]>0 or for T{sub eff}>4700 K. Altering the star formation rate in order to suppress stars older than 10 Gyr improves the agreement for mass but leads to inconsistent color distributions. We also tested the predictions of the TRILEGAL Galactic model. However, unlike Galaxia, it had difficulties in reproducing the photometric properties of the Kepler Input Catalog because it overestimates the number of blue stars. We conclude that either the scaling relations and/or the Galactic models need to be revised to reconcile predictions of theory with asteroseismic observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Calculation of Transfer Functions of Multilayer Biotissues in the Problems of Correction of Their Fluorescence Spectra

A method for rapid calculation of a flux of stimulated fluorescence of a multilayer optically dense medium with inhomogeneous distribution of the fluorophore has been developed. The light field in the medium at the excitation wavelength of fluorescence is represented by a superposition of incident collimated, incident diffuse, and reflected diffuse fluxes. A two-stream approximation is used to describe the light field in the medium at the wavelength of emission of the fluorescence. Fluxes in adjacent elementary layers of the medium and on its surface are connected by simple matrix operators that are obtained using a combination of engineering approaches of radiation-transfer theory and single-scattering approximation. The calculations of fluorescence fluxes of a four-layer biotissue that are excited and recorded at 400–800 nm are compared with their Monte Carlo simulation with a discrepancy of 1%. The effect of the propagation medium on the fluorescence spectra of 5-ALA-induced protoporphyrin IX that are recorded from human skin was studied, and a technique for their correction that is based on measurements and quantitative analysis of the diffuse reflectance spectrum of skin was proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Surface donor states distribution post SiN passivation of AlGaN/GaN heterostructures

In this paper, we present a physics based analytical model to describe the effect of SiN passivation on two-dimensional electron gas density and surface barrier height in AlGaN/GaN heterostructures. The model is based on an extraction technique to calculate surface donor density and surface donor level at the SiN/AlGaN interface. The model is in good agreement with the experimental results and promises to become a useful tool in advanced design and characterization of GaN based heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Emission of photons by positrons channeled in single crystals near an energy of 100 GeV

Emission by 120-GeV positrons in the channeling regime in the (011) plane of a silicon single crystal has been considered. Trajectories of positrons under different initial conditions have been calculated within the theory of nonlinear oscillations. The amplitude distribution function of channeled particles has been determined taking into account the nonlinearity of their motion. The intensity of radiation under various initial conditions has been calculated by two different methods. These results can be useful for comparison with experimental data at energies of positrons beginning with 100 GeV and higher.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Historical (1850-2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols:methodology and application

We present and discuss a new dataset of gridded emissions covering the historical period (1850-2000) in decadal increments at a horizontal resolution of 0.5° in latitude and longitude. The primary purpose of this inventory is to provide consistent gridded emissions of reactive gases and aerosols for use in chemistry model simulations needed by climate models for the Climate Model Intercomparison Program #5 (CMIP5) in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment report. Our best estimate for the year 2000 inventory represents a combination of existing regional and global inventories to capture the best information available at this point; 40 regions and 12 sectors were used to combine the various sources. The historical reconstruction of each emitted compound, for each region and sector, was then forced to agree with our 2000 estimate, ensuring continuity between past and 2000 emissions. Application of these emissions into two chemistry-climate models is used to test their ability to capture long-term changes in atmospheric ozone, carbon monoxide and aerosols distributions. The simulated long-term change in the Northern mid-latitudes surface and mid-troposphere ozone is not quite as rapid as observed. However, stations outside this latitude band show much better agreement in both present-day and long-term trend. The model simulations consistently underestimate the carbon monoxide trend, while capturing the long-term trend at the Mace Head station. The simulated sulfate and black carbon deposition over Greenland is in very good agreement with the ice-core observations spanning the simulation period. Finally, aerosol optical depth and additional aerosol diagnostics are shown to be in good agreement with previously published estimates.

09 BIOMASS FUELS↗

The coupling of fluids, dynamics, and controls on advanced architecture computers

This grant provided for the demonstration of coupled controls, body dynamics, and fluids computations in a workstation cluster environment; and an investigation of the impact of peer-peer communication on flow solver performance and robustness. The findings of these investigations were documented in the conference articles.The attached publication, 'Towards Distributed Fluids/Controls Simulations', documents the solution and scaling of the coupled Navier-Stokes, Euler rigid-body dynamics, and state feedback control equations for a two-dimensional canard-wing. The poor scaling shown was due to serialized grid connectivity computation and Ethernet bandwidth limits. The scaling of a peer-to-peer communication flow code on an IBM SP-2 was also shown. The scaling of the code on the switched fabric-linked nodes was good, with a 2.4 percent loss due to communication of intergrid boundary point information. The code performance on 30 worker nodes was 1.7 (mu)s/point/iteration, or a factor of three over a Cray C-90 head. The attached paper, 'Nonlinear Fluid Computations in a Distributed Environment', documents the effect of several computational rate enhancing methods on convergence. For the cases shown, the highest throughput was achieved using boundary updates at each step, with the manager process performing communication tasks only. Constrained domain decomposition of the implicit fluid equations did not degrade the convergence rate or final solution. The scaling of a coupled body/fluid dynamics problem on an Ethernet-linked cluster was also shown.

Atwood, Christopher↗

Simulating PV System Performance with Component Reliability Distributions

A feature developed by Sandia National Laboratories (SNL) for simulating PV system component reliability has been integrated into the National Renewable Energy Laboratory's System Advisor Model (SAM). The PV Reliability Performance Model (PV-RPM) uses fault/failure and repair distributions to estimate impacts to power, energy production and operational costs. Realization results include the number of failures per component, failure rate, O&M costs incurred based on maintenance strategies, and energy lost based on the failure estimates. As the inputs are inherently uncertain and based on probability distributions, the model realizations include statistical confidence intervals to help bracket the probability that the sample interval contains the true mean value of the variable in question. This paper presents the current validation efforts and provides detail on how this feature is utilized in a PV performance model.

14 SOLAR ENERGY↗