Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical 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 991 records · Page 55

HTR-PROTEUS Pebble Bed Experimental Program Cores 1, 1A, 2, and 3: Hexagonal Close Packing with a 1:2 Moderator-to-Fuel Pebble Ratio

In its deployment as a pebble bed reactor (PBR) critical facility from 1992 to 1996, the PROTEUS facility was designated as HTR-PROTEUS. This experimental program was performed as part of an International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on the Validation of Safety Related Physics Calculations for Low Enriched HTGRs. Within this project, critical experiments were conducted for graphite moderated LEU systems to determine core reactivity, flux and power profiles, reaction-rate ratios, the worth of control rods, both in-core and reflector based, the worth of burnable poisons, kinetic parameters, and the effects of moisture ingress on these parameters. Four benchmark experiments were evaluated in this report: Cores 1, 1A, 2, and 3. These core configurations represent the hexagonal close packing (HCP) configurations of the HTR-PROTEUS experiment with a moderator-to-fuel pebble ratio of 1:2. Core 1 represents the only configuration utilizing ZEBRA control rods. Cores 1A, 2, and 3 use withdrawable, hollow, stainless steel control rods. Cores 1 and 1A are similar except for the use of different control rods; Core 1A also has one less layer of pebbles (21 layers instead of 22). Core 2 retains the first 16 layers of pebbles from Cores 1 and 1A and has 16 layers of moderator pebbles stacked above the fueled layers. Core 3 retains the first 17 layers of pebbles but has polyethylene rods inserted between pebbles to simulate water ingress. The additional partial pebble layer (layer 18) for Core 3 was not included as it was used for core operations and not the reported critical configuration. Cores 1, 1A, 2, and 3 were determined to be acceptable benchmark experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Bidirectional reflectance spectroscopy. III - Correction for macroscopic roughness

A mathematically rigorous formalism is derived by which an arbitrary photometric function for the bidirectional reflectance of a smooth surface may be corrected to include effects of general macroscopic roughness. The correction involves only one arbitrary parameter, the mean slope angle, and is applicable to surfaces of any albedo. Using physically reasonable assumptions and mathematical approximations, the correction expressions are evaluated analytically to second order in the mean slope angle. The correction is applied to the bidirectional-reflectance function of Hapke (1981). Expressions for both the differential and integral brightnesses are obtained. Photometric profiles on hypothetical smooth and rough planets of low and high albedo are shown to illustrate the effects of macroscopic roughness. The theory is applied to observations of Mercury and predicts the integral phase function, the apparent polar darkening, and the lack of limb brightness surge on the planet. The roughness-corrected bidirectional-reflectance function is sufficiently simple that it can be conveniently evaluated on a programmable hand-held calculator.

Hapke, B.↗

Fiducial-cosmology-dependent systematics for the DESI 2024 full-shape analysis

We assess the impact of the fiducial cosmology choice on cosmological inference from full-shape (FS) fits of the galaxy power spectrum in the DESI 2024 Data Release 1 (DR1). Using a suite of AbacusSummit DR1 mock catalogues based on the Planck 2018 best-fit cosmology, we quantify potential systematic shifts introduced by analysing the data under five secondary cosmologies — featuring variations in matter density, thawing dark energy, higher effective number of neutrino species, reduced clustering amplitude, and the DESI DR1 BAO best-fit w 0 w a CDM cosmology — relative to DESI's baseline Planck 2018 cosmology. We investigate two complementary FS analysis approaches: full-modelling (FM) and ShapeFit (SF), each with distinct sensitivities to the assumed fiducial model. Across all tracers, we find for FM that systematic shifts induced by fiducial cosmology mismatches remain well below the DESI DR1 statistical uncertainties, with maximum deviations of 0.22σ DR1 in ΛCDM scenarios and 0.12σ DR1+SN when including SN Ia mock data in extended w 0 w a CDM fits. For SF, the shifts in the compressed parameters remain below 0.45σ DR1 for all tracers and cosmologies.

dark energy experiments↗

Real space iterative reconstruction for vector tomography (RESIRE-V)

Tomography has had an important impact on the physical, biological, and medical sciences. To date, most tomographic applications have been focused on 3D scalar reconstructions. However, in some crucial applications, vector tomography is required to reconstruct 3D vector fields such as the electric and magnetic fields. Over the years, several vector tomography methods have been developed. Here, we present the mathematical foundation and algorithmic implementation of REal Space Iterative REconstruction for Vector tomography, termed RESIRE-V. RESIRE-V uses multiple tilt series of projections and iterates between the projections and a 3D reconstruction. Each iteration consists of a forward step using the Radon transform and a backward step using its transpose, then updates the object via gradient descent. Incorporating with a 3D support constraint, the algorithm iteratively minimizes an error metric, defined as the difference between the measured and calculated projections. The algorithm can also be used to refine the tilt angles and further improve the 3D reconstruction. To validate RESIRE-V, we first apply it to a simulated data set of the 3D magnetization vector field, consisting of two orthogonal tilt series, each with a missing wedge. Our quantitative analysis shows that the three components of the reconstructed magnetization vector field agree well with the ground-truth counterparts. We then use RESIRE-V to reconstruct the 3D magnetization vector field of a ferromagnetic meta-lattice consisting of three tilt series. Our 3D vector reconstruction reveals the existence of topological magnetic defects with positive and negative charges. We expect that RESIRE-V can be incorporated into different imaging modalities as a general vector tomography method. To make the algorithm accessible to a broad user community, we have made our RESIRE-V MATLAB source codes and the data freely available at https://github.com/minhpham0309/RESIRE-V.

47 OTHER INSTRUMENTATION↗

Entanglement-enhanced ac magnetometry in the presence of Markovian noise

Entanglement is a resource to improve the sensitivity of quantum sensors. In an ideal case, using an entangled state as a probe to detect target fields, we can beat the standard quantum limit by which all classical sensors are bounded. However, since entanglement is fragile against decoherence, it is unclear whether entanglement-enhanced metrology is useful in a noisy environment. Its benefit is indeed limited when estimating the amplitude of dc magnetic fields under the effect of parallel Markovian decoherence, where the noise operator is parallel to the target field. In this paper, on the contrary, we show an advantage to using an entanglement over the classical strategy under the effect of parallel Markovian decoherence when we try to detect ac magnetic fields. We consider a scenario to induce a Rabi oscillation of the qubits with the target ac magnetic fields. Although we can, in principle, estimate the amplitude of the ac magnetic fields from the Rabi oscillation, the signal becomes weak if the qubit frequency is significantly detuned from the frequency of the ac magnetic field. We show that, by using the Greenberger-Horne-Zeilinger (GHZ) states, we can significantly enhance the signal of the detuned Rabi oscillation even under the effect of parallel Markovian decoherence. Further, our method is based on the fact that the interaction time between the GHZ states and ac magnetic fields scales as 1/L to mitigate the decoherence effect, where L is the number of qubits, which contributes to improving the bandwidth of the detectable frequencies of the ac magnetic fields. Our results pave the way for new applications of entanglement-enhanced ac magnetometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Kaluza-Klein models: Can we construct a viable example?

In Kaluza-Klein models with toroidal compactification of the extra dimensions, we investigate soliton solutions of Einstein equation. The nonrelativistic gravitational potential of these solitons exactly coincides with the Newtonian one. We obtain the formulas for perihelion shift, deflection of light, time delay of radar echoes and post-Newtonian (PPN) parameters. Using the constraint on PPN parameter {gamma}, we find that the solitonic parameter k should be very big: |k|{>=}2.3x10{sup 4}. We define a soliton solution which corresponds to a pointlike mass source. In this case the soliton parameter k=2, which is clearly contrary to this restriction. A similar problem with the observations takes place for static spherically symmetric perfect fluid with the dustlike equation of state in all dimensions. The common for both of these models is the same (dustlike) equations of state in our three dimensions and in the extra dimensions. All dimensions are treated at equal footing. This is the crucial point. To be in agreement with observations, it is necessary to break the symmetry (in terms of equations of state) between the external/our and internal spaces. It takes place for black strings which are particular examples of solitons with k{yields}{infinity}. For such k, black strings are in concordance with the observations. Moreover, we show that they are the only solitons which are at the same level of agreement with the observations as in general relativity. Black strings can be treated as perfect fluid with dustlike equation of state p{sub 0}=0 in the external/our space and very specific equation of state p{sub 1}=-(1/2){epsilon} in the internal space. The latter equation is due to negative tension in the extra dimension. We also demonstrate that dimension 3 for the external space is a special one. Only in this case we get the latter equation of state. We show that the black string equations of state satisfy the necessary condition of the internal space stabilization. Therefore, black strings are good candidates for a viable model of astrophysical objects (e.g., Sun) if we can provide a satisfactory explanation of negative tension for particles constituting these objects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron stars for holographic metallic criticality

We refer to the ground state of a gravitating, charged ideal fluid of fermions held at a finite chemical potential as an ''electron star.'' In a holographic setting, electron stars are candidate gravity duals for strongly interacting finite fermion density systems. We show how electron stars develop an emergent Lifshitz scaling at low energies. This IR scaling region is a consequence of the two-way interaction between emergent quantum critical bosonic modes and the finite density of fermions. By integrating from the IR region to an asymptotically AdS{sub 4} spacetime, we compute basic properties of the electron stars, including their electrical conductivity. We emphasize the challenge of connecting UV and IR physics in strongly interacting finite density systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Capacitated arc routing problem and its extensions in waste collection

Capacitated arc routing problem (CARP) is the youngest generation of graph theory that focuses on solving the edge/arc routing for optimality. Since many years, operational research devoted to CARP counterpart, known as vehicle routing problem (VRP), which does not fit to several real cases such like waste collection problem and road maintenance. In this paper, we highlighted several extensions of capacitated arc routing problem (CARP) that represents the real-life problem of vehicle operation in waste collection. By purpose, CARP is designed to find a set of routes for vehicles that satisfies all pre-setting constraints in such that all vehicles must start and end at a depot, service a set of demands on edges (or arcs) exactly once without exceeding the capacity, thus the total fleet cost is minimized. We also addressed the differentiation between CARP and VRP in waste collection. Several issues have been discussed including stochastic demands and time window problems in order to show the complexity and importance of CARP in the related industry. A mathematical model of CARP and its new version is presented by considering several factors such like delivery cost, lateness penalty and delivery time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Filament capturing with the Multimaterial Moment-of-Fluid method

A novel method for capturing two-dimensional, thin, under-resolved material configurations, known as “filaments,” is presented in the context of interface reconstruction. This technique uses a partitioning procedure to detect disconnected regions of material in the advective preimage of a cell (indicative of a filament) and makes use of the existing functionality of the Multimaterial Moment-of-Fluid interface reconstruction method to accurately capture the under-resolved feature, while exactly conserving volume. An algorithm for Adaptive Mesh Refinement in the presence of filaments is developed so that refinement is introduced only near the tips of filaments and where the Moment-of-Fluid reconstruction error is still large. Comparison to the standard Moment-of-Fluid method is made. It is demonstrated that using filament capturing at a given resolution yields gains in accuracy comparable to introducing an additional level of mesh refinement at significantly lower cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Focusing of active particles in a converging flow

We consider active particles swimming in a convergent fluid flow in a trapezoid nozzle with no-slip walls. We use mathematical modeling to analyze trajectories of these particles inside the nozzle. By extensive Monte Carlo simulations, we show that trajectories are strongly affected by the background fluid flow and geometry of the nozzle leading to wall accumulation and upstream motion (rheotaxis). In particular, we describe the non-trivial focusing of active rods depending on physical and geometrical parameters. It is also established that the convergent component of the background flow leads to stability of both downstream and upstream swimming at the centerline. The stability of downstream swimming enhances focusing, and the stability of upstream swimming enables rheotaxis in the bulk.

42 ENGINEERING↗

Problem of two-level hierarchical minimax program control the final state of regional social and economic system in the presence of risks

This article discusses a discrete-time dynamical system consisting of a set a controllable objects (region and forming it municipalities). The dynamics each of these is described by the corresponding vector nonlinear discrete-time recurrent vector equations and its control system consist from two levels: basic (control level I) that is dominating and subordinate level (control level II). Both levels have different criterions of functioning and united a priori by determined informational and control connections defined in advance. In this paper we study the problem of optimization of guaranteed result for program control by the final state of regional social and economic system in the presence of risks. For this problem we proposed in this work an economical and mathematical model of two-level hierarchical minimax program control the final state of regional social and economic system in the presence of risks and the general scheme for its solving.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generalized Nuclear Data: A New Structure (with Supporting Infrastructure) for Handling Nuclear Data

The Evaluated Nuclear Data File (ENDF) format was designed in the 1960s to accommodate neutron reaction data to support nuclear engineering applications in power, national security and criticality safety. Over the years, the scope of the format has been extended to handle many other kinds of data including charged particle, decay, atomic, photo-nuclear and thermal neutron scattering. Although ENDF has wide acceptance and support for many data types, its limited support for correlated particle emission, limited numeric precision, and general lack of extensibility mean that the nuclear data community cannot take advantage of many emerging opportunities. More generally, the ENDF format provides an unfriendly environment that makes it difficult for new data evaluators and users to create and access nuclear data. The Cross Section Evaluation Working Group (CSEWG) has begun the design of a new Generalized Nuclear Data (or 'GND') structure, meant to replace older formats with a hierarchy that mirrors the underlying physics, and is aligned with modern coding and database practices. In support of this new structure, Lawrence Livermore National Laboratory (LLNL) has updated its nuclear data/reactions management package Fudge to handle GND structured nuclear data. Fudge provides tools for converting both the latest ENDF format (ENDF-6) and the LLNL Evaluated Nuclear Data Library (ENDL) format to and from GND, as well as for visualizing, modifying and processing (i.e., converting evaluated nuclear data into a form more suitable to transport codes) GND structured nuclear data. GND defines the structure needed for storing nuclear data evaluations and the type of data that needs to be stored. But unlike ENDF and ENDL, GND does not define how the data are to be stored in a file. Currently, Fudge writes the structured GND data to a file using the eXtensible Markup Language (XML), as it is ASCII based and can be viewed with any text editor. XML is a meta-language, meaning that it has a primitive set of definitions for representing hierarchical data/text in a file. Other meta-languages, like HDF5 which stores the data in binary form, can also be used to store GND in a file. In this paper, we will present an overview of the new GND data structures along with associated tools in Fudge.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multimatrix models and tri-Sasaki Einstein spaces

Localization methods reduce the path integrals in N{>=}2 supersymmetric Chern-Simons gauge theories on S{sup 3} to multimatrix integrals. A recent evaluation of such a two-matrix integral for the N=6 superconformal U(N)xU(N) Aharony-Bergman-Jafferis-Maldacena theory produced detailed agreement with the AdS/CFT correspondence, explaining, in particular, the N{sup 3/2} scaling of the free energy. We study a class of p-matrix integrals describing N=3 superconformal U(N){sup p} Chern-Simons gauge theories. We present a simple method that allows us to evaluate the eigenvalue densities and the free energies in the large N limit keeping the Chern-Simons levels k{sub i} fixed. The dual M-theory backgrounds are AdS{sub 4}xY, where Y are seven-dimensional tri-Sasaki Einstein spaces specified by the k{sub i}. The gravitational free energy scales inversely with the square root of the volume of Y. We find a general formula for the p-matrix free energies that agrees with the available results for volumes of the tri-Sasaki Einstein spaces Y, thus providing a thorough test of the corresponding AdS{sub 4}/CFT{sub 3} dualities. This formula is consistent with the Seiberg duality conjectured for Chern-Simons gauge theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reactor Physics Measurements and Benchmark Specifications for Oak Ridge Highly Enriched Uranium Sphere (ORSphere)

In the early 1970s Dr. John T. Mihalczo (team leader), J.J. Lynn, and J.R. Taylor performed experiments at the Oak Ridge Critical Experiments Facility (ORCEF) with highly enriched uranium (HEU) metal (called Oak Ridge Alloy or ORALLOY) in an effort to recreate GODIVA I results with greater accuracy than those performed at Los Alamos National Laboratory in the 1950s. The purpose of the Oak Ridge ORALLOY Sphere (ORSphere) experiments was to estimate the unreflected and unmoderated critical mass of an idealized sphere of uranium metal corrected to a density, purity, and enrichment such that it could be compared with the GODIVA I experiments. Additionally, various material reactivity worths, the surface material worth coefficient, the delayed neutron fraction, the prompt neutron decay constant, relative fission density, and relative neutron importance were all measured. The critical assembly, material reactivity worths, the surface material worth coefficient, and the delayed neutron fraction were all evaluated as benchmark experiment measurements. The reactor physics measurements are the focus of this paper; although for clarity the critical assembly benchmark specifications are briefly discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cylindrical Kernel Density Estimators For Monte Carlo Neutron Transport Reactor Physics Problems

Kernel Density Estimators (KDEs) have been investigated as an alternative to histogram tallies in Monte Carlo radiation transport simulations for capturing spatially-resolved quantities such as scalar flux and reaction rate densities. KDEs are non-parametric estimators that show potential for producing smooth estimates of quantities of interest with reduced variance compared to traditional histogram tallies. Furthermore, KDEs obtain estimates of the underlying distribution at user-defined points without requiring an underlying mesh, thus making KDEs an attractive potential alternative to histogram tallies. Recent work has introduced the Mean Free Path based KDE (MFP KDE) in order to capture distributions with large curvature at material interfaces, e.g. absorption reaction rates in UO{sub 2} fuel in light water reactors. However, the multivariate MFP KDE produces inaccuracies in the form of localized spikes in problems with material cross sections that contain large resonances. Since the MFP KDE works well in one dimension, it is hypothesized that a different MFP KDE formulation where the normalization coefficient only contains one factor of the cross section would be able to capture these distributions with large curvature without producing spikes in the estimate. This can be accomplished by using a cylindrical KDE whose 2-D support area is wedge-shaped and applying the MFP KDE to the radial component of the cylindrical KDE. These cylindrical KDEs use 2-D (r,θ) cylindrical coordinates with no dependence in axial dimension z. This paper derives the cylindrical KDE and the cylindrical MFP KDE and tests them on a reactor physics pin-cell benchmark problem, producing accurate estimates of flux and reaction rate densities with differences less than 0.5% when compared to a reference histogram. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analog neutron transport for nuclear instrumentation applications with the Monte-Carlo code TRIPOLI-4 - Paper 3

Monte-Carlo codes simulating particle transport usually make use of well-known variance reduction techniques, leading to more efficient simulations. Even in configurations which do not deal with deep-penetration, some of these techniques, such as the implicit capture, are actually often used by default. The statistical weight of the simulated particles is then modified in order to take into account the modification of the sampling without introducing any bias in the estimation process. For instance, the implicit capture decreases the statistical weight of the particles at each collision, whereas the sampling of a single outgoing particle increases this statistical weight in case of scattering reactions producing multiple neutrons. In addition, Russian roulette and splitting operations, which are usually also played at each collision, lead to other modifications of the statistical weights. As a result, Monte-Carlo codes provide particle histories that are valid on average, but collect statistical weights which differ from 1, when tallying the estimation of the quantities of interest. Yet, some nuclear instrumentation calculations inherently need to work with statistical weights equal to 1, since either individual events have to be collected (such as the arrival of a particle at a detector, the measurement of a time of flight, etc), or the coincidence of events or also correlations between individual events have to be taken into account. To address these specific applications with Monte-Carlo simulation codes, users need a full analog simulation mode. The continuous-energy Monte-Carlo code TRIPOLI4{sup R} is mainly dedicated to shielding, reactor physics and criticality safety calculations. Its capabilities have recently been extended to widen its use to nuclear instrumentation with the possibility of γ-spectrometry calculations and the simulation of photonuclear reactions for instance. Other features needed by nuclear instrumentation applications have been available in TRIPOLI-4 for a long time, such as the sampling of the time variable and the possibility of tallying in a time grid for most of the available tallies (fluxes, reaction rates, dose equivalent rates, energy depositions, etc). In TRIPOLI-4, the transport between two consecutive collisions is analog if no variance reduction technique is requested by the user to handle deep-penetration configurations. However, up to the version 9 of the code, the collision itself is sampled in a non-analog way. This paper presents a first implementation of a full analog neutron transport mode in TRIPOLI-4, preliminary verification results and an application of this functionality to the simulation of the Nucifer detector located in Saclay, France. It should be noted that only neutron transport is dealt with in this work and that all nuclear data come from JEFF-3.1.1 or ENDF/B-VII.0 evaluations, including fission data when needed. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Production of Entanglement Entropy by Decoherence

We examine here the dynamics of entanglement entropy of all parts in an open system consisting of a two-level dimer interacting with an environment of oscillators. The dimer-environment interaction is almost energy conserving. We find the precise link between decoherence and production of entanglement entropy. We show that not all environment oscillators carry significant entanglement entropy and we identify the oscillator frequency regions which contribute to the production of entanglement entropy. For energy conserving dimer-environment interactions the models are explicitly solvable and our results hold for all dimer-environment coupling strengths. We carry out a mathematically rigorous perturbation theory around the energy conserving situation in the presence of small non-energy conserving interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Big Data in Reciprocal Space: Sliding Fast Fourier Transforms for Determining Periodicity

Significant advances in atomically resolved imaging of crystals and surfaces have occurred in the last decade allowing unprecedented insight into local crystal structures and periodicity. Yet, the analysis of the long-range periodicity from the local imaging data, critical to correlation of functional properties and chemistry to the local crystallography, remains a challenge. Here, we introduce a Sliding Fast Fourier Transform (FFT) filter to analyze atomically resolved images of in-situ grown La5/8Ca3/8MnO3 films. We demonstrate the ability of sliding FFT algorithm to differentiate two sub-lattices, resulting from a mixed-terminated surface. Principal Component Analysis (PCA) and Independent Component Analysis (ICA) of the Sliding FFT dataset reveal the distinct changes in crystallography, step edges and boundaries between the multiple sub-lattices. The method is universal for images with any periodicity, and is especially amenable to atomically resolved probe and electron-microscopy data for rapid identification of the sub-lattices present.

74 ATOMIC AND MOLECULAR PHYSICS↗