Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Regularization”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Strongly hyperbolic quasilinear systems revisited, with applications to relativistic fluid dynamics

We revisit the theory of first-order quasilinear systems with diagonalizable principal part and only real eigenvalues, what is commonly referred to as strongly hyperbolic systems. We provide a self-contained and simple proof of local well-posedness, in the Hadamard sense, of the Cauchy problem. Our regularity assumptions are very minimal. As an application, we apply our results to systems of ideal and viscous relativistic fluids, where the theory of strongly hyperbolic equations has been systematically used to study several systems of physical interest.

Mathematics↗
DOE OSTI2024

Structural Effects on Solubility and Crystallinity in Polyamide Ionomers

Although polyamides have been an industrial staple for decades, their solubility and processability remain limited due to the strong hydrogen bonding between amide groups. Here, we report a family of polyamides in which aryl rings, alkyl spacers, sulfonate groups, and amide linkages are regularly spaced along the polymer backbone, allowing us to probe the structural elements that affect processability. We accessed these polymers through a synthetic route based on an expanded diamine monomer and characterized them through a range of physical and spectroscopic techniques. Our results show that the combination of sulfonate groups with increased amide content results in highly soluble polymers, which can dissolve in polar solvents such as water, methanol, and N,N-dimethylformamide. Infrared spectroscopy on the solid polymers shows that the aliphatic amides engage in hydrogen bonding with the sulfonate groups, thus inhibiting ion aggregation, crystallization, and microphase separation. These findings expand the avenue to sulfonate-containing polyamide chemistry and provide design rules for the synthesis of more processable ionic polyamides.

amides↗
DOE OSTI2021

Modelling a Heaving Point-Absorber with a Closed-Loop Control System Using the DualSPHysics Code

The present work addresses the need for an efficient, versatile, accurate and open-source numerical tool to be used during the design stage of wave energy converters (WECs). The device considered here is the heaving point-absorber developed and tested by Sandia National Laboratories. The smoothed particle hydrodynamics (SPH) method, as implemented in DualSPHysics, is proposed since its meshless approach presents some important advantages when simulating floating devices. The dynamics of the power take-off system are also modelled by coupling DualSPHysics with the multi-physics library Project Chrono. A satisfactory matching between experimental and numerical results is obtained for: (i) the heave response of the device when forced via its actuator; (ii) the vertical forces acting on the fixed device under regular waves and; (iii) the heave response of the WEC under the action of both regular waves and the actuator force. This proves the ability of the numerical approach proposed to simulate accurately the fluid–structure interaction along with the WEC’s closed-loop control system. In addition, radiation models built from the experimental and WAMIT results are compared with DualSPHysics by plotting the intrinsic impedance in the frequency domain, showing that the SPH method can be also employed for system identification.

42 ENGINEERING↗
DOE OSTI2023

(U) Correlated Sampling Using Batch Statistics to Reduce the Uncertainty of Combinations of KSEN Outputs with MCNP6

The relative sensitivity of k eff to the densities of nuclides in a material are combined to compute relative sensitivities to other inputs. When computed in a single Monte Carlo run, the nuclide density sensitivities are correlated, and the statistical uncertainties propagated to other inputs will be incorrect unless those correlations are accounted for. Equations are presented to apply correlated sampling using batch statistics for the sum of an arbitrary number of random tallies and the difference of two random tallies when each is multiplied by a different constant. When correlated sampling is used on a recent benchmark evaluation, the correct statistical uncertainties for certain combinations of sensitivities are dramatically smaller than the incorrect uncertainties. A user-controlled, regular output of MCNP6’s KSEN sensitivities that allows batch statistics to be applied to combinations would be extremely valuable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Structural response reconstruction using a system-equivalent singular vector basis

Here, this paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.

42 ENGINEERING↗
DOE OSTI2024

Asymptotic Freedom at the Berezinskii-Kosterlitz-Thouless Transition without Fine-Tuning Using a Qubit Regularization

We propose a two-dimensional hard-core loop-gas model as a way to regularize the asymptotically free massive continuum quantum field theory that emerges at the Berezinskii-Kosterlitz-Thouless transition. Without fine-tuning, our model can reproduce the universal step-scaling function of the classical lattice X Y model in the massive phase as we approach the phase transition. This is achieved by lowering the fugacity of Fock-vacuum sites in the loop-gas configuration space to zero in the thermodynamic limit. Some of the universal quantities at the Berezinskii-Kosterlitz-Thouless transition show smaller finite size effects in our model as compared to the traditional X Y model. Our model is a prime example of qubit regularization of an asymptotically free massive quantum field theory in Euclidean space-time and helps understand how asymptotic freedom can arise as a relevant perturbation at a decoupled fixed point without fine-tuning. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

Plug-and-Play Methods for Integrating Physical and Learned Models in Computational Imaging: Theory, algorithms, and applications

Plug-and-play (PnP) priors constitute one of the most widely used frameworks for solving computational imaging problems through the integration of physical models and learned models. PnP leverages high-fidelity physical sensor models and powerful machine learning methods for prior modeling of data to provide state-of-the-art reconstruction algorithms. PnP algorithms alternate between minimizing a data fidelity term to promote data consistency and imposing a learned regularizer in the form of an image denoiser. Recent highly successful applications of PnP algorithms include biomicroscopy, computerized tomography (CT), magnetic resonance imaging (MRI), and joint ptychotomography. This article presents a unified and principled review of PnP by tracing its roots, describing its major variations, summarizing main results, and discussing applications in computational imaging. Additionally, we also point the way toward further developments by discussing recent results on equilibrium equations that formulate the problem associated with PnP algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

MOOSE Optimization Module Overview: Application to Residual Stress Inversion in Nuclear Fuel Plates

Manufacturing methods, material selections and unique plate geometries used in low enriched uranium plate fuels for the High Performance Research Reactor (HPRR) have led to concerns about the residual stresses effect on bond strength of the fuel/cladding interfaces and overall fuel performance under irradiation. The current HPRR plate fuel test program performed residual stress calculations using the contour method with experimental displacement data obtained by incrementally cutting the plate fuel. In this work, we use a gradient based inverse optimization module developed within the MOOSE framework to estimate normal and shear residual stresses within the plate fuels using displacement data taken from a cut placed in the HPRR fuel plate. These results are verified against residual stress calculations obtained from the contour method. The inversion protocol used in this work utilizes two meshes to separately discretize the physics and parameter space. This allows a fine mesh to be used to fully resolve the physics and geometry of the forward problem. A coarser mesh is used to resolve the parameter space with spatially varying mesh density determined from a sensitivity analysis. The adaptive resolution of the parameter space essentially has the effect of regularization as it is designed to achieve the bias-variance tradeoff.

42 ENGINEERING↗
DOE OSTI2020

Influence of crystal structure on constitutive anisotropy of silica sand at particle-scale

In spite of recent advanced geotechnical research regarding the influence of particle-scale morphology and fabric on the constitutive behavior of silica sands, there is still no published literature that investigated the influence of crystal structure of sand particles on their constitutive behavior. The main constituent of silica sand particles is quartz mineral, which has a trigonal crystal system and hexagonal lattice. This paper investigates the influence of quartz crystal structure on the constitutive response of silica sand particles using 3D X-ray diffraction (3DXRD), synchrotron micro-computed tomography (SMT), and 3D finite element (FE) analysis. Here, unconfined uniaxial compression experiments on single particles of synthetic silica cubes exposed the notion of constitutive anisotropy caused by the crystal structure of quartz. Since the silica cubes have simple regular geometry (1 mm) 3 , they were loaded in different directions with respect to their crystal local planes. The crystal structure was also identified for a natural silica sand specimen composed of 2705 particles and loaded in confined uniaxial compression. Physical limitations of the experimental measurements made a strong case for the use of 3D FE analysis to study how the change in crystal local orientation of individual sand particles can fundamentally produce anisotropy in the constitutive response of natural sand specimens. The 3D FE analysis was validated to accurately model the crystal-based constitutive anisotropy in silica particles using the results of the silica cube experiments. The analysis was then executed on 3D FE meshes that closely matched the physical shape and fabric of the 2705 particles in the natural sand specimen. Both the experiments and 3D FE analysis showed that the crystal structure of quartz essentially causes a directional anisotropy in the constitutive behavior of silica particles.

58 GEOSCIENCES↗
DOE OSTI2021

Dynamics of a single polyampholyte chain

Polymers that feature both positive and negative charges along chains, known as polyampholytes, represent a class of materials that hold promise for a new generation of energy storage devices, the design of which will require knowledge of the underlying structure and dynamics. Here, we develop a theory based on the Rouse model for the dynamic structure factor of a single polyampholyte chain in the weak coupling regime (negligible intramolecular electrostatics) or subjected to weak external electric fields (governed by linear response). Neglecting effects of small ions, we find deviations in scaling from the classic Rouse theory and make predictions for scattering experiments performed on polyampholytes. We find that, under weak coupling with arbitrarily strong fields, the dynamics are highly dependent on the charge distribution and consequently look at two representative examples—random charge densities and periodic charge densities—with different scaling properties. Additionally, under weak fields, the dynamics are largely independent of charge distribution. Finally, we investigate the influence of hydrodynamic effects and the implications of including inertial effects in the model.

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

Fuzzy gauge theory for quantum computers

Continuous gauge theories, because of their bosonic degrees of freedom, have an infinite-dimensional local Hilbert space. Encoding these degrees of freedom on qubit-based hardware demands some sort of “qubitization” scheme, where one approximates the behavior of a theory while using only finitely many degrees of freedom. We propose a novel qubitization strategy for gauge theories, called “fuzzy gauge theory,” building on the success of the fuzzy σ -model in earlier work. We provide arguments that the fuzzy gauge theory lies in the same universality class as regular gauge theory, in which case its use would obviate the need of any further limit besides the usual spatial continuum limit. Furthermore, we demonstrate that these models are relatively resource-efficient for quantum simulations. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2023

Quantum simulations of SO(5) many-fermion systems using qudits

The structure and dynamics of many-body systems are the result of a delicate interplay between underlying interactions. Fermionic pairing, for example, plays a central role in various physical systems, ranging from condensed matter to nuclear systems, where it can lead to collective phenomena such as superconductivity and superfluidity. In atomic nuclei, the interplay between pairing and particle-hole interactions leads to a high degree of complexity and intricate entanglement structures. Despite this apparent complexity, symmetries emerge and manifest themselves in observable regular patterns. These symmetries and their breakings have long been used to determine relevant degrees of freedom and simplify classical descriptions of many-body systems. Here, this work explores the potential utility of quantum computers with arrays of qudits in simulating interacting fermionic systems, when the qudits can naturally map the relevant degrees of freedom determined by an underlying symmetry group. The Agassi model of fermions interacting via particle-hole and pairing interactions is based on an underlying so(5) algebra. Such systems can intuitively be partitioned into pairs of modes with five basis states, which thus naturally map to arrays of d = 5 qudits (qu5its). Classical noiseless simulations of the time evolution of systems with up to twelve qu5its are performed, by implementing quantum circuits that are developed herein, using PYTHON codes invoking Google's CIRQ software. The resource requirements of the qu5it circuits are analyzed and compared with two different mappings to qubit systems: a physics-aware Jordan-Wigner mapping requiring four qubits per mode pair and a state-to-state mapping requiring three qubits per mode pair. While the dimensionality of Hilbert spaces in mappings to qu5it systems are less than those for the corresponding qubit systems, the number of entangling operations, depending on the available hardware, can either be greater or smaller than for the physics-aware Jordan-Wigner mapping. The state-to-state mapping, while having a smaller Hilbert space than Jordan-Wigner mappings, appears to be the least efficient in gate counts. Further, a previously unknown sign problem has been identified from Trotterization errors in time evolving high-energy excitations. There appear to be advantages in employing quantum computers with arrays of qudits to perform simulations of many-body dynamics that exploit the role of underlying symmetries, specifically in lowering the required quantum resources and in reducing anticipated errors that take the simulation out of the physical space. If the necessary entangling gates are not directly supported by the hardware, physics-aware mappings to qubits may, however, be advantageous for other aspects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Repositioning Quantum Cellular Automata for Dependable Quantum-Classical Systems

Quantum Cellular Automata (QCA) provides a structured model of distributed quantum computation with inherent locality and regularity properties that are suited to dependable execution. However, QCA remain largely absent from discussions on reproducibility, fault management, and orchestration in heterogeneous quantum-classical systems. We propose a dual-axis framework that situates QCA within both computation and physical realizability, revealing regions where robust, scalable, and hardware-constrained quantum dynamics may reside. By revisiting prior results through the lens of reproducibility and architecture resilience, we suggest that QCA offers a potential substrate for benchmarking and system-level co-design.

Stapleton, Nicholas [ORNL] (ORCID:0000000335305325↗
DOE OSTI2024

Fusion Neutron Generator

The proposed code, named FROG (Fusion neutron Generator) is built upon the open-source particle transport Monte Carlo toolkit Geant4. Geant4 provides C++ classes that can be leveraged to build application-specific codes dealing with the transport of particles through matter. Geant4-based codes are applied in high-energy particle physics experiments, medical applications, shielding, and space applications for example. The FROG code allows the user to define the geometry of a neutron converter device shaped as a hollow cylinder, where a neutron breeding material such as lithium deuteride (LiD) is cladded by two concentric cylinders. Such neutron converter is then placed inside a regular nuclear fission reactor, where thermal neutrons will react with the neutron breeder material (typically, Lithium 6), and through a series of reactions, will generate high-energy neutrons – neutrons whose kinetic energy are around 14 MeV. The hollowed central portion can hold a specimen that will be bombarded by high-energy neutrons created inside the neutron breeding material. Figuratively speaking, this type of device transforms neutrons from thermal (~0.625 eV) to fusion (~14 MeV) energies and is sometimes termed “fusion-to-thermal neutron converters” in the literature. The code consists of C++ source file compiled and linked to generate an executable. The user can select the dimensions of the converter (radius, length, and thickness of the breeder material), the breeder material type, the cladding material, and the specimen material that will be activated or irradiated. As input, the neutron flux for a specific location inside a reactor, for instance, positions in ATR, is required. As output, the code predicts the number of high-energy neutrons produced, the total neutron flux and fluence as well as its detailed spectrum. The physics involved in such device is very complex, as it requires modeling neutron transport, light-ion (tritons) transport, as well as fusion reactions. The Geant4 toolkit provides the required physical models.

Martin, NicholasP. [Idaho National Laboratory (INL↗
DOE OSTI2024

One-shot omnidirectional pressure integration through matrix inversion

In this work, we present a method to perform 2D and 3D omnidirectional pressure integration from velocity measurements with a single-iteration matrix inversion approach. This work builds upon our previous work, where the rotating parallel ray approach was extended to the limit of infinite rays by taking continuous projection integrals of the ray paths and recasting the problem as an iterative matrix inversion problem. This iterative matrix equation is now 'fast-forwarded' to the 'infinity' iteration, leading to a different matrix equation that can be solved in a single step, thereby presenting the same computational complexity as the Poisson equation. We observe computational speedups of ~10 6 when compared to brute-force omnidirectional integration methods, enabling the treatment of grids of ~10 9 points and potentially even larger in a desktop setup at the time of publication. Further examination of the boundary conditions of our one-shot method shows that omnidirectional pressure integration implements a boundary condition where the boundary points are treated as interior points to the extent that information is available. Finally, we show how the method can be extended from the regular grids typical of particle image velocimetry to the unstructured meshes characteristic of particle tracking velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Characterizing juvenile salmon predation risk during early marine residence

Predation mortality can influence the distribution and abundance of fish populations. While predation is often assessed using direct observations of prey consumption, potential predation can be predicted from co-occurring predator and prey densities under varying environmental conditions. Juvenile Pacific salmon Oncorhynchus spp. (i.e., smolts) from the Columbia River Basin experience elevated mortality during the transition from estuarine to ocean habitat, but a thorough understanding of the role of predation remains incomplete. We used a Holling type II functional response to estimate smolt predation risk based on observations of piscivorous seabirds (sooty shearwater [ Ardenna griseus ] and common murre [ Uria aalge ]) and local densities of alternative prey fish including northern anchovy ( Engraulis mordax ) in Oregon and Washington coastal waters during May and June 2010–2012. We evaluated predation risk relative to the availability of alternative prey and physical factors including turbidity and Columbia River plume area, and compared risk to returns of adult salmon. Seabirds and smolts consistently co-occurred at sampling stations throughout most of the study area (mean = 0.79 ± 0.41, SD), indicating that juvenile salmon are regularly exposed to avian predators during early marine residence. Predation risk for juvenile coho ( Oncorhynchus kisutch ), yearling Chinook salmon ( O . tshawytscha ), and subyearling Chinook salmon was on average 70% lower when alternative prey were present. Predation risk was greater in turbid waters, and decreased as water clarity increased. Juvenile coho and yearling Chinook salmon predation risk was lower when river plume surface areas were greater than 15,000 km 2 , while the opposite was estimated for subyearling Chinook salmon. These results suggest that plume area, turbidity, and forage fish abundance near the mouth of the Columbia River, all of which are influenced by river discharge, are useful indicators of potential juvenile salmon mortality that could inform salmonid management.

Phillips, Elizabeth M. (ORCID:0000000327752563)↗
DOE OSTI2021

Synthesis and Characterization of Macrocyclic Ionic Liquids for CO 2 Separation

In this work, a series of ionic liquids with symmetric, macrocyclic cations composed of N-heterocycles connected by alkyl or ether linkages coupled with bis(trifluoromethylsulfonyl)imide or bis(pentafluoroethanesulfonyl)amide anions were synthesized to explore their application in CO 2 separation. These novel ionic liquids showed selective physical absorption of CO 2 compared to N 2 making them potential candidates for carbon capture. We examined the effect of cation structure on the physicochemical and separations properties of the ionic liquid and showed that these ionic liquids exhibit an inverse relationship between molar volume and CO 2 sorption capacity following the trend predicted by models based on regular solution theory. The fractional free volumes of the ionic liquids were lower than typical alkyl-imidazolium ionic liquids despite the cyclic structure of the cations which could be related to the flexibility of the alkyl and ether linkages. These results provide added insight into our understanding of structure–property relationships in CO 2 solubility in ionic liquids and afford a path toward enhancing the separations properties based on structure.

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

Metal-assisted etch combined with regularizing etch

In an aspect of the disclosure, a process for forming nanostructuring on a silicon-containing substrate is provided. The process comprises (a) performing metal-assisted chemical etching on the substrate, (b) performing a clean, including partial or total removal of the metal used to assist the chemical etch, and (c) performing an isotropic or substantially isotropic chemical etch subsequently to the metal-assisted chemical etch of step (a). In an alternative aspect of the disclosure, the process comprises (a) performing metal-assisted chemical etching on the substrate, (b) cleaning the substrate, including removal of some or all of the assisting metal, and (c) performing a chemical etch which results in regularized openings in the silicon substrate.

36 MATERIALS SCIENCE↗