Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Modeling”

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 829 records · Page 46

SOLAR SPECTRAL IRRADIANCE, SOLAR ACTIVITY, AND THE NEAR-ULTRA-VIOLET

The previous calculations of the Solar Spectral Irradiance (SSI) by the Solar Radiation Physical Modeling, version 2 system, are updated in this work by including new molecular photodissociation cross-sections of important species, and many more levels and lines in its treatment of non-LTE radiative transfer. The current calculations including the new molecular photodissociation opacities produce a reduced over-ionizaton of heavy elements in the lower chromosphere and solve the problems with prior studies of the UV SSI in the wavelength range 160–400 nm and now reproduce the available observations with much greater accuracy. Calculations and observations of the near-UV at 0.1 nm resolution and higher are compared. The current set of physical models includes four quiet-Sun and five active-region components, from which radiance is computed for ten observing angles. These radiances are combined with images of the solar disk to obtain the SSI and Total Solar Irradiance and their variations. The computed SSI is compared with measurements from space at several nm resolution and agreement is found within the accuracy level of these measurements. An important result is that the near-UV SSI increase with solar activity is significant for the photodissociation of ozone in the terrestrial atmosphere because a number of highly variable upper chromospheric lines overlap the ozone Hartley band.

79 ASTRONOMY AND ASTROPHYSICS↗

A halo-independent lower bound on the dark matter capture rate in the Sun from a direct detection signal

We show that a positive signal in a dark matter (DM) direct detection experiment can be used to place a lower bound on the DM capture rate in the Sun, independent of the DM halo. For a given particle physics model and DM mass we obtain a lower bound on the capture rate independent of the local DM density, velocity distribution, galactic escape velocity, as well as the scattering cross section. We illustrate this lower bound on the capture rate by assuming that upcoming direct detection experiments will soon obtain a significant signal. When comparing the lower bound on the capture rate with limits on the high-energy neutrino flux from the Sun from neutrino telescopes, we can place upper limits on the branching fraction of DM annihilation channels leading to neutrinos. With current data from IceCube and Super-Kamiokande non-trivial limits can be obtained for spin-dependent interactions and direct annihilations into neutrinos. In some cases also annihilations into ττ or bb start getting constrained. For spin-independent interactions current constraints are weak, but they may become interesting for data from future neutrino telescopes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A modern Monte Carlo investigation of the TG-43 dosimetry parameters for an {sup 125}I seed already having AAPM consensus data

Purpose: To investigate potential causes for differences in TG-43 brachytherapy dosimetry parameters in the existent literature for the model IAI-125A{sup 125}I seed and to propose new standard dosimetry parameters. Methods: The MCNP5 code was used for Monte Carlo (MC) simulations. Sensitivity of dose distributions, and subsequently TG-43 dosimetry parameters, was explored to reproduce historical methods upon which American Association of Physicists in Medicine (AAPM) consensus data are based. Twelve simulation conditions varying{sup 125}I coating thickness, coating mass density, photon interaction cross-section library, and photon emission spectrum were examined. Results: Varying{sup 125}I coating thickness, coating mass density, photon cross-section library, and photon emission spectrum for the model IAI-125A seed changed the dose-rate constant by up to 0.9%, about 1%, about 3%, and 3%, respectively, in comparison to the proposed standard value of 0.922 cGy h{sup −1} U{sup −1}. The dose-rate constant values by Solberg et al. [“Dosimetric parameters of three new solid core {sup 125}I brachytherapy sources,” J. Appl. Clin. Med. Phys. 3, 119–134 (2002)], Meigooni et al. [“Experimental and theoretical determination of dosimetric characteristics of IsoAid ADVANTAGE™ {sup 125}I brachytherapy source,” Med. Phys. 29, 2152–2158 (2002)], and Taylor and Rogers [“An EGSnrc Monte Carlo-calculated database of TG-43 parameters,” Med. Phys. 35, 4228–4241 (2008)] for the model IAI-125A seed and Kennedy et al. [“Experimental and Monte Carlo determination of the TG-43 dosimetric parameters for the model 9011 THINSeed™ brachytherapy source,” Med. Phys. 37, 1681–1688 (2010)] for the model 6711 seed were +4.3% (0.962 cGy h{sup −1} U{sup −1}), +6.2% (0.98 cGy h{sup −1} U{sup −1}), +0.3% (0.925 cGy h{sup −1} U{sup −1}), and −0.2% (0.921 cGy h{sup −1} U{sup −1}), respectively, in comparison to the proposed standard value. Differences in the radial dose functions between the current study and both Solberg et al. and Meigooni et al. were <10% for r ≤ 5 cm, and increased for r > 5 cm with a maximum difference of 29% at r = 9 cm. In comparison to Taylor and Rogers, these differences were lower (maximum of 2% at r = 9 cm). For the similarly designed model 6711 {sup 125}I seed, differences did not exceed 0.5% for 0.5 ≤ r ≤ 10 cm. Radial dose function values varied by 1% as coating thickness and coating density were changed. Varying the cross-section library and source spectrum altered the radial dose function by 25% and 12%, respectively, but these differences occurred at r = 10 cm where the dose rates were very low. The 2D anisotropy function results were most similar to those of Solberg et al. and most different to those of Meigooni et al. The observed order of simulation condition variables from most to least important for influencing the 2D anisotropy function was spectrum, coating thickness, coating density, and cross-section library. Conclusions: Several MC radiation transport codes are available for calculation of the TG-43 dosimetry parameters for brachytherapy seeds. The physics models in these codes and their related cross-section libraries have been updated and improved since publication of the 2007 AAPM TG-43U1S1 report. Results using modern data indicated statistically significant differences in these dosimetry parameters in comparison to data recommended in the TG-43U1S1 report. Therefore, it seems that professional societies such as the AAPM should consider reevaluating the consensus data for this and others seeds and establishing a process of regular evaluations in which consensus data are based upon methods that remain state-of-the-art.

61 RADIATION PROTECTION AND DOSIMETRY↗

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An X-ray and Neutron Scattering Study of the Formation of Core-Shell Type Polyoxometalates

A typical type of core-shell polyoxometalates can be obtained through the Keggin-type polyoxometalatetemplated growth of a layer of spherical shell structure of {Mo72Fe30}. Small angle X-ray scattering is used to study the structural features and stability of the core-shell structures in aqueous solutions. Timeresolved small angle X-ray scattering is applied to monitor the synthetic reactions and a three-stage formation mechanism is proposed to describe the synthesis of the core-shell polyoxometalates based on the monitoring results. New protocols have been developed by fitting the X-ray data with customed physical models, which provide more convincing, objective, and completed data interpretation. Quasielastic and inelastic neutron scattering are used to probe the dynamics of water molecules in the core-shell structures and two different types of water molecules, the confined and structured water, are observed. These water molecules play an important role in bridging core and shell structures and stabilizing the cluster structures.

Yin, Panchao↗

A covariance methodology for theoretically generated thermal neutron inelastic scattering data - 14502

The The S(α,β) double-differential thermal neutron scattering law tabulated in ENDF File 7 is, by convention, produced theoretically through fundamental scattering physics models for material-specific interatomic structures. Differential and integral inelastic scattering cross sections are calculated by numerically integrating over S(α,β). This information is supplied to nuclear system simulation codes. For systems with a significant thermal neutron flux, calculated system response parameters can be sensitive to this thermal cross section data. Formats and procedures are currently established for representing uncertainty information, or covariances, for various types of ENDF reaction data. However, no published ENDF evaluations contain covariance data for S(α,β) or associated differential and integral inelastic scattering cross sections. Furthermore, no accepted methodology exists for quantifying or representing covariances for this thermal scattering data. For most solids, the influence of a material's interatomic structure and dynamics on inelastic scattering cross sections can be adequately described through the phonon density of states (DOS). The phonon DOS is a probability density function of the population of available vibrational modes in energy, and it is typically the fundamental input for calculating S(α,β). Probable variation in the spectral shape of the phonon DOS may be established which characterizes uncertainties in the physics models and methodology employed in its production. Through Monte Carlo sampling of perturbations from the reference phonon spectrum, an S(α,β) covariance matrix may be generated. With appropriate sensitivity information, this can be propagated to generate covariance data for differential and integral inelastic scattering cross sections. In this work, hexagonal graphite is used as an example material for the purpose of demonstrating procedures for analyzing, calculating and representing uncertainty information for thermal neutron inelastic scattering data. (authors)

42 ENGINEERING↗

POLARIZATION SWINGS REVEAL MAGNETIC ENERGY DISSIPATION IN BLAZARS

The polarization signatures of blazar emissions are known to be highly variable. In addition to small fluctuations of the polarization angle around a mean value, large (≳180°) polarization angle swings are sometimes observed. We suggest that such phenomena can be interpreted as arising from light travel time effects within an underlying axisymmetric emission region. We present the first simultaneous fitting of the multi-wavelength spectrum, variability, and time-dependent polarization features of a correlated optical and gamma-ray flaring event of the prominent blazar 3C279, which was accompanied by a drastic change in its polarization signatures. This unprecedented combination of spectral, variability, and polarization information in a coherent physical model allows us to place stringent constraints on the particle acceleration and magnetic field topology in the relativistic jet of a blazar, strongly favoring a scenario in which magnetic energy dissipation is the primary driver of the flare event.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-dimensional atomic packing in amorphous solids with liquid-like structure

Liquids and solids are two fundamental states of matter. However, our understanding of their three-dimensional atomic structure is mostly based on physical models. Here we use atomic electron tomography to experimentally determine the three-dimensional atomic positions of monatomic amorphous solids, namely a Ta thin film and two Pd nanoparticles. We observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Instead of forming icosahedra, the majority of pentagonal bipyramids arrange into pentagonal bipyramid networks with medium-range order. Molecular dynamics simulations further reveal that pentagonal bipyramid networks are prevalent in monatomic metallic liquids, which rapidly grow in size and form more icosahedra during the quench from the liquid to the glass state. Furthermore, these results expand our understanding of the atomic structures of amorphous solids and will encourage future studies on amorphous–crystalline phase and glass transitions in non-crystalline materials with three-dimensional atomic resolution.

36 MATERIALS SCIENCE↗

Comparisons between tokamak fueling of gas puffing and supersonic molecular beam injection in 2D simulations

Plasma fueling with high efficiency and deep injection is very important to enable fusion power performance requirements. It is a powerful and efficient way to study neutral transport dynamics and find methods of improving the fueling performance by doing large scale simulations. Two basic fueling methods, gas puffing (GP) and supersonic molecular beam injection (SMBI), are simulated and compared in realistic divertor geometry of the HL-2A tokamak with a newly developed module, named trans-neut, within the framework of BOUT++ boundary plasma turbulence code [Z. H. Wang et al., Nucl. Fusion 54, 043019 (2014)]. The physical model includes plasma density, heat and momentum transport equations along with neutral density, and momentum transport equations. Transport dynamics and profile evolutions of both plasma and neutrals are simulated and compared between GP and SMBI in both poloidal and radial directions, which are quite different from one and the other. It finds that the neutrals can penetrate about four centimeters inside the last closed (magnetic) flux surface during SMBI, while they are all deposited outside of the LCF during GP. It is the radial convection and larger inflowing flux which lead to the deeper penetration depth of SMBI and higher fueling efficiency compared to GP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nuclear Reaction Data Evaluations for Nuclear Applications Physics and modeling [Slides]

This presentation provides an evaluation workflow scheme for the project. It also discusses nuclear data libraries, including an overview of worldwide nuclear data libraries and current funding agencies. The presentation also discusses the theoretical framework commonly used in the evaluated work, including the R-matrix theory (≈ up to a few keVs for heavy nuclei) as well as pre-equilibrium (PE) and Hauser-Feshbach (HF) theory (⪆1 MeV for medium/heavy nuclei). Inclusion of experimental effects are also discussed. Examples of evaluated work, including theory, experiments, and applications are discussed, specifically n+ 28,29,30 Si evaluation (direct capture component), n+ 53 Cr evaluation (strong multiple scattering effects), n+ 239 Pu(n,γf) reaction (fluctuations in neutron multiplicities), and α+ 17,18 O evaluations (emerging neutron spectra).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗

Physics-based model of irradiation creep for ferritic materials under fusion energy operation conditions

Irradiation creep is known to be an important process for structural materials in nuclear environments, potentially leading to creep failure at temperatures where thermal creep is generally negligible. While there is a great deal of data for irradiation creep in steels and zirconium alloys in light water reactor conditions, much less is known for first wall materials under fusion energy conditions. Lacking suitable fusion neutron sources for detailed experimentation, modeling, and simulation can help bridge the dose-rate and spectral-effects gap and produce quantifiable expectations for creep deformation of first wall materials under standard fusion conditions. In this paper, we develop a comprehensive model for irradiation creep created from merging a crystal plasticity representation of the dislocation microstructure and a defect evolution simulator that accounts for the entire cluster dimensionality space. Both approaches are linked by way of a climb velocity that captures dislocation-biased defect absorption and a dislocation strengthening term that reflects the accumulation of defect clusters in the system. Here, we carry out our study in Fe under first wall fusion reactor conditions, characterized by a fusion neutron spectrum with average recoil energies of 20 keV and a damage dose rate of ≈3×10 –7 dpa/s at temperatures between 300 and 800 K.

36 MATERIALS SCIENCE↗

A Fast and Accurate Transient Stability Assessment Method Based on Deep Learning: WECC Case Study

Transient stability is one of the critical aspects of power system stability assessment. The increasing integration of inverter-based resources and the retirement of conventional synchronous generators result in the decreasing system inertia and growing complexity of system operating conditions. Using a few selected typical operating conditions cannot guarantee system transient stability in all operating conditions, and the time-domain simulation of all operating conditions requires tremendous time and is often infeasible. This paper proposes a more efficient transient stability assessment method based on deep learning. The binary search method is used to determine the critical clearing time (CCT) in creating training databased by time-domain simulation. This method is fast and accurate with 1 ms resolution. The buses whose CCTs are lower than 200 ms are considered critical buses. Buses close to each other are grouped based on their mutual admittance matrix to reduce the search space of the critical buses. This paper also proposes the generator feature normalization based on the physical model. Case study on the reduced 240-bus WECC system model demonstrates that the proposed method can predict CCT accurately and efficiently.

critical clearing time↗

Bottom-up estimation and top-down prediction: Solar energy prediction combining information from multiple sources

Accurately forecasting solar power using the data from multiple sources is an important but challenging problem. Our goal is to combine two different physics model forecasting outputs with real measurements from an automated monitoring network so as to better predict solar power in a timely manner. To this end, we propose a new approach of analyzing large-scale multilevel models with great computational efficiency requiring minimum monitoring and intervention. This approach features a division of the large scale data set into smaller ones with manageable sizes, based on their physical locations, and fit a local model in each area. The local model estimates are then combined sequentially from the specified multilevel models using our novel bottom-up approach for parameter estimation. The prediction, on the other hand, is implemented in a top-down matter. The proposed method is applied to the solar energy prediction problem for the U.S. Department of Energy's SunShot Initiative.

14 SOLAR ENERGY↗

Noise Optimization for MKIDs with Different Design Geometries and Material Selections

The separation and optimization of noise components is critical to microwave-kinetic inductance detector (MKID) development. We analyze the effect of several changes to the lumped-element inductor and interdigitated capacitor geometry on the noise performance of a series of MKIDs intended for millimeter-wavelength experiments. We extract the contributions from two-level system noise in the dielectric layer, the generation-recombination noise intrinsic to the superconducting thin-film, and system white noise from each detector noise power spectrum and characterize how these noise components depend on detector geometry, material, and measurement conditions such as driving power and temperature. We observe a reduction in the amplitude of two-level system noise with both an elevated sample temperature and an increased gap between the fingers within the interdigitated capacitors for both aluminum and niobium detectors. We also verify the expected reduction of the generation-recombination noise and associated quasiparticle lifetime with reduced inductor volume. This study also iterates over different materials, including aluminum, niobium, and aluminum manganese, and compares the results with an underlying physical model.

generation-recombination noise↗

Retrieval of Snow Freeboard of Antarctic Sea Ice Using Waveform Fitting of CryoSat-2 Returns

In this paper we develop a CryoSat-2 algorithm to retrieve the surface elevation of the air–snow interface over Antarctic sea ice. This algorithm utilizes a two-layer physical model that accounts for scattering from a snow layer atop sea ice as well as scattering from below the snow surface. The model produces waveforms that are fit to CryoSat-2 level 1B data through a bounded trust region least-squares fitting process. These fit waveforms are then used to track the air–snow interface and retrieve the surface elevation at each point along the CryoSat-2 ground track, from which the snow freeboard is computed. To validate this algorithm, we compare retrieved surface elevation measurements and snow surface radar return power levels with those from Operation IceBridge, which flew along a contemporaneous CryoSat-2 orbit in October 2011 and November 2012. Average elevation differences (standard deviations) along the flight lines (IceBridge Airborne Topographic Mapper, ATM – CryoSat-2) are found to be 0.016 cm (29.24 cm) in 2011 and 2.58 cm (26.65 cm) in 2012. The spatial distribution of monthly average pan-Antarctic snow freeboard found using this method is similar to what was observed from NASA's Ice, Cloud, and land Elevation Satellite (ICESat), where the difference (standard deviation) between October 2011–2017 CryoSat-2 mean snow freeboard and spring 2003–2007 mean freeboard from ICESat is 1.92 cm (9.23 cm). While our results suggest that this physical model and waveform fitting method can be used to retrieve snow freeboard from CryoSat-2, allowing for the potential to join laser and radar altimetry data records in the Antarctic, larger (∼30 cm) regional differences from ICESat and along-track differences from ATM do exist, suggesting the need for future improvements to the method. Snow–ice interface elevation retrieval is also explored as a potential to obtain snow depth measurements. However, it is found that this retrieval method often tracks a strong scattering layer within the snow layer instead of the actual snow–ice interface, leading to an overestimation of ice freeboard and an underestimation of snow depth in much of the Southern Ocean but with promising results in areas such as the East Antarctic sector.

Fons, Steven W.↗

BRISC v.1.0

BRISC is a developmental prototype for a nextgeneration “systems-level” integrated performance and safety code (IPSC) for nuclear reactors. Its development served to demonstrate how a lightweight multi-physics coupling approach can be used to tightly couple the physics models in several different physics codes (written in a variety of languages) into one integrated package for simulating accident scenarios in a liquid sodium cooled “burner” nuclear reactor. For example, the RIO Fluid Flow and Heat transfer code developed at Sandia (SNL: Chris Moen, Dept. 08005) is used in BRISC to model fluid flow and heat transfer, as well as conduction heat transfer in solids. Because BRISC is a prototype, its most practical application is as a foundation or starting point for developing a true production code. The sub-codes and the associated models and correlations currently employed within BRISC were chosen to cover the required application space and demonstrate feasibility, but were not optimized or validated against experimental data within the context of their use in BRISC.

Source record↗

Measurements of the parallel wavenumber of lower hybrid waves in the scrape-off layer of a high-density tokamak

In lower hybrid current drive (LHCD) experiments on tokamaks, the parallel wavenumber of lower hybrid waves is an important physics parameter that governs the wave propagation and absorption physics. However, this parameter has not been experimentally well-characterized in the present-day high density tokamaks, despite the advances in the wave physics modeling. In this paper, we present the first measurement of the dominant parallel wavenumber of lower hybrid waves in the scrape-off layer (SOL) of the Alcator C-Mod tokamak with an array of magnetic loop probes. The electric field strength measured with the probe in typical C-Mod plasmas is about one-fifth of that of the electric field at the mouth of the grill antenna. The amplitude and phase responses of the measured signals on the applied power spectrum are consistent with the expected wave energy propagation. At higher density, the observed k{sub ||} increases for the fixed launched k{sub ||}, and the wave amplitude decreases rapidly. This decrease is correlated with the loss of LHCD efficiency at high density, suggesting the presence of loss mechanisms. Evidence of the spectral broadening mechanisms is observed in the frequency spectra. However, no clear modifications in the dominant k{sub ||} are observed in the spectrally broadened wave components, as compared to the measured k{sub ||} at the applied frequency. It could be due to (1) the probe being in the SOL and (2) the limited k{sub ||} resolution of the diagnostic. Future experiments are planned to investigate the roles of the observed spectral broadening mechanisms on the LH density limit problem in the strong single pass damping regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗