Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational physics”

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 649 records · Page 36

General Improvements to the MCNP Alpha-Eigenvalue Solver

In this document, a new algorithm for computing the α-eigenvalue is implemented and tested in MCNP. The algorithm follows the traditional k-α method in which neutrons are followed from birth until fission in batches. During this simulation, a number of tallies are performed which are derived from integrating the transport equation over all phase space. These new tallies along with the collapsed transport equation are used to compute a physically valid α with or without delayed neutrons. The performance and accuracy of the method is then tested, in which it is found that the new algorithm typically out- performs the current implementation while also representing a more complete physical picture. Through a convergence analysis, it was found that many simple test problems had unexpectedly large biases for reasonable quantities of neutrons per batch. This occurred with both the old and the new algorithm. However, the tally method had lower biases for all problems tested. Special implementation considerations are also discussed.

97 MATHEMATICS AND COMPUTING↗

Multislice does it all—calculating the performance of nanofocusing X-ray optics

Here, we describe an approach to calculating the optical performance of a wide range of nanofocusing X-ray optics using multislice scalar wave propagation with a complex X-ray refractive index. This approach produces results indistinguishable from methods such as coupled wave theory, and it allows one to reproduce other X-ray optical phenomena such as grazing incidence reflectivity where the direction of energy flow is changed significantly. Just as finite element analysis methods allow engineers to compute the thermal and mechanical responses of arbitrary structures too complex to model by analytical approaches, multislice propagation can be used to understand the properties of the real-world optics of finite extent and with local imperfections, allowing one to better understand the limits to nanoscale X-ray imaging.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assigning the absolute configuration of single aliphatic molecules by visual inspection

Deciphering absolute configuration of a single molecule by direct visual inspection is the next step in compound identification, with far-reaching implications for medicinal chemistry, pharmacology, and natural product synthesis. We demonstrate the feasibility of this approach utilizing low temperature atomic force microscopy (AFM) with a CO-functionalized tip to determine the absolute configuration and orientation of a single, adsorbed [123]tetramantane molecule, the smallest chiral diamondoid. We differentiate between single enantiomers on Cu(111) by direct visual inspection, and furthermore identify molecular dimers and molecular clusters. The experimental results are confirmed by a computational study that allowed quantification of the corresponding intermolecular interactions. The unique toolset of absolute configuration determination combined with AFM tip manipulation opens a route for studying molecular nucleation, including chirality-driven assembly or reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement

Alamo is a high-performance scientific code that uses block-structured adaptive mesh refinement to solve such problems as: the ignition and burn of solid rocket propellant, plasticity, damage and fracture in materials undergoing loading, and the interaction of compressible flow with eroding solid materials. Alamo is powered by AMReX, and provides a set of unique methods, models, and algorithms that enable it to solve solid-mechanics problems (coupled to other physical behavior such as fluid flow or thermal diffusion) using the power of block-structured adaptive mesh refinement.

36 MATERIALS SCIENCE↗

Investigation of a model based reconstruction technique for neutron tomography at the PULSTAR reactor

Herein, a model based neutron computed tomography (CT) reconstruction technique is developed and implemented to improve the spatial resolution of the reconstructed images. The spatial resolution of a scintillator and a charge-coupled device (CCD) based digital neutron imaging system is quantified and used as the system response function in the iterative reconstruction process. The response function is included in the forward projection model to provide better estimation of the actual projection, thus improving the spatial resolution of the reconstructed images. Experiments were performed at the neutron imaging facility of the PULSTAR research reactor to corroborate the performance of the developed algorithm with a 50 μm and 250 μa m 6 LiF:ZnS thick scintillator. The reconstruction results of the model based technique have shown an improvement of nearly 40% in achievable spatial resolution as compared to the standard filtered back projection (FBP) technique.

47 OTHER INSTRUMENTATION↗

Multiphysics Modeling of a Single Channel in a Nuclear Thermal Propulsion Grooved Ring Fuel Element

In the past, fuel rods have been used in nuclear propulsion applications. A new fuel element concept that reduces weight and increases efficiency uses a stack of grooved discs. Each fuel element is a flat disc with a hole on the interior and grooves across the top. Many grooved ring fuel elements for use in nuclear thermal propulsion systems have been modeled, and a single flow channel for each design has been analyzed. For increased efficiency, a fuel element with a higher surface-area-to-volume ratio is ideal. When grooves are shallower, i.e., they have a lower surface area, the results show that the exit temperature is higher. By coupling the physics of turbulence with those of heat transfer, the effects on the cooler gas flowing through the grooves of the thermally excited solid can be predicted. Parametric studies were done to show how a pressure drop across the axial length of the channels will affect the exit temperatures of the gas. Geometric optimization was done to show the behaviors that result from the manipulation of various parameters. Temperature profiles of the solid and gas showed that more structural optimization is needed to produce the desired results. Keywords: Nuclear Thermal Propulsion, Fuel Element, Heat Transfer, Computational Fluid Dynamics, Coupled Physics Computations, Finite Element Analysis

Kim, Tony↗

Computation of Transonic Flows Using Potential Methods

The proposed paper will describe the state of the art associated with numerical solution of the full or exact velocity potential equation for solving transonic, external-aerodynamic flows. The presentation will begin with a review of the literature emphasizing research activities of the past decade. Next, the various forms of the full or exact velocity potential equation, the equation's corresponding mathematical characteristics, and the derivation assumptions will be presented and described in detail. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, will be presented and discussed relative to the more complete Euler or Navier-Stokes formulations. The technical presentation will continue with a description of recently developed full potential numerical approach characteristics. This description will include governing equation nondimensionalization, physical-to-computational-domain mapping procedures, a limited description of grid generation requirements, the spatial discretization scheme, numerical implementation of boundary conditions, and the iteration scheme. The next portion of the presentation will present and discuss numerical results for several two- and three-dimensional aerodynamic applications. Included in the results section will be a discussion and demonstration of a typical grid refinement analysis for determining spatial convergence of the numerical solution and level of solution accuracy. Computer timings for a variety of full potential applications will be compared and contrasted with similar results for the Euler equation formulation. Finally. the presentation will end with concluding remarks and recommendations for future work.

Hoist, Terry L.↗

Feasibility study for measurement of insulation compaction in the cryogenic rocket fuel storage tanks at Kennedy Space Center by fast/thermal neutron techniques

The liquid hydrogen and oxygen cryogenic storage tanks at John F. Kennedy Space Center (KSC) use expanded perlite as thermal insulation. Some of the perlite may have compacted over time, compromising the thermal performance and also the structural integrity of the tanks. Neutrons can readily penetrate through the 1.75 cm outer steel shell and through the entire 120 cm thick perlite zone. Neutrons interactions with materials produce characteristic gamma rays which are then detected. In compacted perlite the count rates in the individual peaks in the gamma ray spectrum will increase. Portable neutron generators can produce neutron simultaneous fluxes in two energy ranges: fast (14 MeV) and thermal (25 meV). Fast neutrons produce gamma rays by inelastic scattering which is sensitive to Si, Al, Fe and O. Thermal neutrons produce gamma rays by radiative capture in prompt gamma neutron activation (PGNA), which is sensitive to Si, Al, Na, K and H among others. The results of computer simulations using the software MCNP and measurements on a test article suggest that the most promising approach would be to operate the system in time-of-flight mode by pulsing the neutron generator and observing the subsequent die away curve in the PGNA signal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gate-free state preparation for fast variational quantum eigensolver simulations

Abstract The variational quantum eigensolver is currently the flagship algorithm for solving electronic structure problems on near-term quantum computers. The algorithm involves implementing a sequence of parameterized gates on quantum hardware to generate a target quantum state, and then measuring the molecular energy. Due to finite coherence times and gate errors, the number of gates that can be implemented remains limited. In this work, we propose an alternative algorithm where device-level pulse shapes are variationally optimized for the state preparation rather than using an abstract-level quantum circuit. In doing so, the coherence time required for the state preparation is drastically reduced. We numerically demonstrate this by directly optimizing pulse shapes which accurately model the dissociation of H 2 and HeH + , and we compute the ground state energy for LiH with four transmons where we see reductions in state preparation times of roughly three orders of magnitude compared to gate-based strategies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evolutionary optimization of PAW data-sets for accurate high pressure simulations

We examine the challenge of performing accurate electronic structure calculations at high pressures by comparing the results of all-electron full potential linearized augmented-plane-wave calculations, as implemented in the WIEN2k code, with those of the projector augmented wave (PAW) method, as implemented in Quantum ESPRESSO or Abinit code. In particular, we focus on developing an automated and consistent way of generating transferable PAW data-sets that can closely produce the all electron equation of state defined from zero to arbitrary high pressures. The technique we propose is an evolutionary search procedure that exploits the ATOMPAW code to generate atomic data-sets and the Quantum ESPRESSO software suite for total energy calculations. We demonstrate different aspects of its workability by optimizing PAW basis functions of some elements relatively abundant in planetary interiors. In addition, we introduce a new measure of atomic data-set goodness by considering their performance uniformity over an extended pressure range. - Highlights: • An attempt made to achieve the AE-FLAPW accuracy, at the favorable computational efficiency of the PAW approach for an outspread pressure range. • Provided an automatic and consistent route for developing efficient and transferable PAW data-sets. • Explicated an improved measure for the goodness of data-sets by including uniformity of the data-set performance over an extended pressure range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multi-modal tomographic imaging system for poolside characterization of nuclear test fuels: Design considerations and studies

Testing and qualification of advanced nuclear fuels involves an iterative process of prototyping, in-pile irradiation testing, and in-situ or ex-situ examination. Fuel restructuring and fission product migration during burnup are among the most important aspects of fuel evolution that affect several important performance characteristics such as heat removal, accident tolerance, and fission product retention. Pool-side non-destructive characterization techniques provide fuel developers with tools to understand fuel evolution at different points of burnup. A design for a compact, submersible, and multi-modal gamma-ray tomography system for imaging irradiated nuclear fuel is presented here. Detector selection, collimator geometry and fabrication, mechanical design, imaging protocol and acquisition speed are discussed. Modeling calculations show that sub-millimeter resolution can be achieved in both transmission computed tomography images as well as in emission computed tomography images in a matter of hours. Several design compromises and fabrication challenges are discussed which should be taken into consideration for future submersible gamma-ray tomography instruments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Optimization Studies for ISOL Type High-Powered Targets

The research studied one-step and two-step Isotope Separation on Line (ISOL) targets for future radioactive beam facilities with high driver-beam power through advanced computer simulations. As a target material uranium carbide in the form of foils was used because of increasing demand for actinide targets in rare-isotope beam facilities and because such material was under development in ISAC at TRIUMF when this project started. Simulations of effusion were performed for one-step and two step targets and the effects of target dimensions and foil matrix were studied. Diffusion simulations were limited by availability of diffusion parameters for UC x material at reduced density; however, the viability of the combined diffusion?effusion simulation methodology was demonstrated and could be used to extract physical parameters such as diffusion coefficients and effusion delay times from experimental isotope release curves. Dissipation of the heat from the isotope-producing targets is the limiting factor for high-power beam operation both for the direct and two-step targets. Detailed target models were used to simulate proton beam interactions with the targets to obtain the fission rates and power deposition distributions, which were then applied in the heat transfer calculations to study the performance of the targets. Results indicate that a direct target, with specification matching ISAC TRIUMF target, could operate in 500-MeV proton beam at beam powers up to ~40 kW, producing ~8 10 13 fission/s with maximum temperature in UCx below 2200 C. Targets with larger radius allow higher beam powers and fission rates. For the target radius in the range 9-mm to 30-mm the achievable fission rate increases almost linearly with target radius, however, the effusion delay time also increases linearly with target radius.

43 PARTICLE ACCELERATORS↗

Realizing string-net condensation: Fibonacci anyon braiding for universal gates and sampling chromatic polynomials

Abstract The remarkable complexity of a topologically ordered many-body quantum system is encoded in the characteristics of its anyons. Quintessential predictions emanating from this complexity employ the Fibonacci string net condensate (Fib SNC) and its anyons: sampling Fib-SNC would estimate chromatic polynomials while exchanging its anyons would implement universal quantum computation. However, physical realizations remained elusive. We introduce a scalable dynamical string net preparation (DSNP) that constructs Fib SNC and its anyons on reconfigurable graphs suitable for near-term superconducting processors. Coupling the DSNP approach with composite error-mitigation on deep circuits, we create, measure, and braids Fibonacci anyons; charge measurements show 94% accuracy, and exchanging the anyons yields the expected golden ratioϕwith 98% average accuracy. We then sample the Fib SNC to estimate chromatic polynomial atϕ + 2 for several graphs. Our results establish the proof of principle for using Fib-SNC and its anyons for fault-tolerant universal quantum computation and aim at a classically hard problem.

Science & Technology - Other Topics↗

Electron, Phonon, and Superconducting Properties of Yttrium and Sulfur Hydrides under High Pressures

We have analyzed the electron, phonon, and superconducting properties of binary hydrides YH{sub 6}, H{sub 3}S, YH{sub 10}, UH{sub 10}, and CaH{sub 10}, as compared to the properties of ternary hydride YS{sub 4}H{sub 4} using a unified physical and computational approach to determining the stability range and the superconducting transition temperature. The normal and superconducting characteristics of ternary hydride YS{sub 4}H{sub 4} at high pressure are calculated. The results for ternary hydride YS{sub 4}H{sub 4} are compared with the properties of YH{sub 6}, H{sub 3}S, YH{sub 10}, UH{sub 10}, and CaH{sub 10} compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A computational survey of semiconductors for power electronics

Power electronics (PE) are used to control and convert electrical energy in a wide range of applications from consumer products to large-scale industrial equipment. While Si-based power devices account for the vast majority of the market, wide band gap semiconductors such as SiC, GaN, and Ga 2 O 3 are starting to gain ground. However, these emerging materials face challenges due to either non-negligible defect densities, high synthesis and processing costs, or poor thermal properties. Here, we report on a broad computational survey aimed to identify promising materials for future power electronic devices beyond SiC, GaN, and Ga 2 O 3 . We consider 863 oxides, sulfides, nitrides, carbides, silicides, and borides that exhibit finite calculated band gaps. We utilize ab initio methods in conjunction with models for intrinsic carrier mobility, and critical breakdown field to compute the widely used Baliga figure of merit. We also compute the lattice thermal conductivity and use it as an additional screening parameter. In addition to correctly identifying known PE materials, our survey has revealed a number of promising candidates exhibiting the desirable combination of high figure of merit and high lattice thermal conductivity, which we propose for further experimental investigations.

36 MATERIALS SCIENCE↗

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

Advanced Test Reactor Core Modeling Update Project Annual Report for Fiscal Year 2013

Legacy computational reactor physics software tools and protocols currently used for support of Advanced Test Reactor (ATR) core fuel management and safety assurance, and to some extent, experiment management, are inconsistent with the state of modern nuclear engineering practice, and are difficult, if not impossible, to verify and validate (V&V) according to modern standards. Furthermore, the legacy staff knowledge required for effective application of these tools and protocols from the 1960s and 1970s is rapidly being lost due to staff turnover and retirements. In late 2009, the Idaho National Laboratory (INL) initiated a focused effort, the ATR Core Modeling Update Project, to address this situation through the introduction of modern high-fidelity computational software and protocols. This aggressive computational and experimental campaign will have a broad strategic impact on the operation of the ATR, both in terms of improved computational efficiency and accuracy for support of ongoing DOE programs as well as in terms of national and international recognition of the ATR National Scientific User Facility (NSUF).

97 MATHEMATICS AND COMPUTING↗