Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 631 records · Page 35

Asynchronous domain decomposition methods for nonlinear PDEs

One- and two-level parallel asynchronous methods for the numerical solution of nonlinear systems of equations, especially those arising from (nonlinear) partial differential equations, are studied. The proposed methods are based on domain decomposition techniques. Local convergence theorems are presented in several cases, with appropriate hypotheses. Computational results on a shared memory multiprocessor machine for various problems exhibiting nonlinearities are reported, illustrating the potential of these asynchronous methods, especially for heterogeneous clusters.

97 MATHEMATICS AND COMPUTING↗

Parallel-in-Time Integration for Nonlinear Hyperbolic Problems (Final Report)

The work for the subcontract was situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. Over the past two decades, much progress has been made in this area; however, this progress has largely been limited to diffusion-dominated PDEs, with some recent success in scalar linear hyperbolic PDEs. Given the ubiquity of numerical simulations of hyperbolic PDEs throughout the sciences, in particular, nonlinear hyperbolic systems, there is a strong need to develop efficient parallel-in-time techniques for hyperbolic problems beyond simple scalar and linear cases, which is the main topic of this subcontract. The main focus of the work was to further develop and perfect coarse-grid operators for the Multigrid Reduction-inTime (MGRIT) method applied to hyperbolic PDEs that were recently proposed in PhD thesis, based on a modified semi-Lagrangian approach.

97 MATHEMATICS AND COMPUTING↗

LHC Event Generation in the Exascale Era

MCFM is a dedicated Monte-Carlo simulation program for collider phenomenology at highest energies. Designed during the Tevatron era, it has successfully incorporated the latest developments needed for LHC precision calculations and remained on the forefront of collider phenomenology. The Fortran code includes interfaces to modern PDF and loop reduction libraries but has been unchanged structurally compared to the earlier versions. Parallel computing has been enabled using OpenMP and MPI. MCFM provides numerically highly stable one-loop amplitudes and superior phase-space efficiency, leading to excellent performance in NXLO calculations using jettiness or qT subtraction techniques for IR regularization.

Campbell, John [Fermilab]↗

Modeling the Influenza A NP-vRNA-Polymerase Complex in Atomic Detail

Seasonal flu is an acute respiratory disease that exacts a massive toll on human populations, healthcare systems and economies. The disease is caused by an enveloped Influenza virus containing eight ribonucleoprotein (RNP) complexes. Each RNP incorporates multiple copies of nucleoprotein (NP), a fragment of the viral genome (vRNA), and a viral RNA-dependent RNA polymerase (POL), and is responsible for packaging the viral genome and performing critical functions including replication and transcription. A complete model of an Influenza RNP in atomic detail can elucidate the structural basis for viral genome functions, and identify potential targets for viral therapeutics. In this work we construct a model of a complete Influenza A RNP complex in atomic detail using multiple sources of structural and sequence information and a series of homology-modeling techniques, including a motif-matching fragment assembly method. Our final model provides a rationale for experimentally-observed changes to viral polymerase activity in numerous mutational assays. Further, our model reveals specific interactions between the three primary structural components of the RNP, including potential targets for blocking POL-binding to the NP-vRNA complex. The methods developed in this work open the possibility of elucidating other functionally-relevant atomic-scale interactions in additional RNP structures and other biomolecular complexes.

59 BASIC BIOLOGICAL SCIENCES↗

Design of proton deflectometry with in situ x-ray fiducial for magnetized high-energy-density systems

We report a design and implementation of proton deflectometry with an in situ reference x-ray image of a mesh to precisely measure non-uniform magnetic fields in expanding plasmas at the OMEGA and OMEGA EP laser facilities. The technique has been developed with proton and x-ray sources generated from both directly driven capsule implosions and short pulse laser–solid interactions. The accuracy of the measurement depends on the contrast of both the proton and x-ray images. Here we present numerical and analytic studies to optimize the image contrast using a variety of mesh materials and grid spacings. Our results show clear enhancement of the image contrast by factors of four to six using a high Z mesh with large grid spacing. This leads to further improvement in the accuracy of the magnetic field measurement using this technique in comparison with its first demonstration at the OMEGA laser facility [Rev. Sci. Instrum. 93, 023502 (2022) [CrossRef]].

47 OTHER INSTRUMENTATION↗

The TRANSCEND University Consortium: Theme 2 - Site Decommissioning, Deconstruction and Remediation - 20329

This paper will discuss the aims, objectives and progress to-date under Theme 2 of the TRANSCEND consortium project: Site decommissioning, deconstruction and remediation. Decommissioning nuclear sites involves waste retrieval, decontamination, deconstruction and, where necessary, containment and/or remediation of the remaining structure and surrounding land. Critical to management of these processes is limiting radiation exposure for the workforce, restricting the spread of radionuclides in groundwater, surface water and airborne particulates, and minimising the volume of contaminated waste for disposal. The aim of Theme 2 research is to develop new technologies for monitoring, remediation and containment that serve to minimise the volume of radioactively contaminated waste for disposal, for application prior to, during and after retrieval, deconstruction and decontamination operations. Prior research, conducted under the previous DISTINCTIVE project, demonstrated that colloidal silica grout can penetrate low permeability materials (including cement) for hydraulic barrier formation, and improved sorption capacity. The silica grout can be injected at surface using extremely small (potentially gravity-driven) fluid pressures, without the need for borehole drilling. Current work is investigating the erodibility of silica-grouted soils for inhibition of airborne and water-borne particulates as well as enhancing the grout's sorption capacity by addition of other materials to provide a chemical, as well as hydraulic, barrier to subsurface migration. EK remediation uses low voltage DC current to control migration of contaminants in porous media as well as to remove or degrade them. Researchers in the consortium have already demonstrated that low-energy ex-situ EK techniques can be used to provide remediation and volume minimisation for AWE legacy wastes in the UK. Under the current project researchers are building EK test cells containing simulated site materials at laboratory and intermediate(m)-scales to: remove, focus or degrade contaminants (remediation or waste minimisation); and direct subsurface water, chemical and colloid flow (fencing/containment or forced migration). This lab-based research is being informed by numerical models of EK processes that can subsequently be used to design full-scale on-site applications by nuclear site holders. Research under theme 2 of TRANSCEND will ultimately combine novel non-invasive detection technologies with EK techniques and colloidal silica grout barriers, to optimise the containment of radionuclide contamination in soils. Thus, allowing us to detect in-situ contamination, mobilise it to a selected location and grout it in-situ, prior to the initiation of decommissioning and deconstruction operations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mitigation of distortion of Al/steel part under simulated paint baking condition: Experiment and numerical model studies

Multi-material joining of lightweight structures is essential to reduce vehicle weight for more energy savings and less greenhouse gas emission. However, mismatch of thermal expansion coefficient for dissimilar materials during the paint baking process can induce part distortion and joint failure for adhesive bonding. Here, in the present work, a thermomechanical model based on contact mechanics and large deformation theory was developed for dissimilar high-strength Al alloy and steel components to study the distortion mechanism and influential factors of the residual gap. The established model was used to optimize joint conditions, such as pitch distance and part geometry. When a weld pitch is shorter than 100 mm, the maximum gap between Al and steel part can be greatly reduced to 0.1 mm, and the local stress and plastic strain around the joint during the oven heating and cooling cycle are also substantially reduced compared with the long pitch case (900 mm). The numerical modeling results revealed that a comparable bending stiffness ratio between the steel and Al cross sections is critical to the minimization of gap and distortion under paint baking condition. Digital image correlation technique was used to measure the overall part distortion and local strain distribution that were used to validate the model prediction. Weld bonding (adhesive bonding with friction bit joining) process was successfully employed to join Al to steel component without gap opening in adhesive after the paint baking and cooling.

36 MATERIALS SCIENCE↗

Lifetime renormalization of driven weakly anharmonic superconducting qubits. II. The readout problem

Recent experiments in superconducting qubit systems have shown an unexpectedly strong dependence of the qubit relaxation rate on the readout drive power. This phenomenon limits the maximum measurement strength and thus the achievable readout speed and fidelity. We address this problem here and provide a plausible mechanism for drive-power dependence of relaxation rates. To this end we introduce a two-parameter perturbative expansion in qubit anharmonicity and the drive amplitude through a unitary transformation technique introduced in Part I. This approach naturally reveals number-nonconserving terms in the Josephson potential as a fundamental mechanism through which applied microwave drives can activate additional relaxation mechanisms. Here, we present our results in terms of an effective master equation with renormalized state- and drive-dependent transition frequency and relaxation rates. Comparison of numerical results from this effective master equation to those obtained from a Lindblad master equation which only includes number-conserving terms (i.e., Kerr interactions) shows that number-nonconserving terms can lead to significant drive-power dependence of the relaxation rates. The systematic expansion technique introduced here is of general applicability to obtaining effective master equations for driven-dissipative quantum systems that contain weakly nonlinear degrees of freedom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of flow properties coupled to experimental and numerical analyses of dense, granular flows for solar thermal energy storage

Granular flows of sintered bauxite proppants were researched along an inclined plane for solar thermal energy storage applications. Granular flow properties needed to drive numeric granular models were measured for improved numerical model predictions for Carbobead CP 50/140, 40/100, and 30/60 particles. Particle shape and size distributions were determined by coupling optical microscopy to an in-house image processing algorithm. The impulse excitation technique was used to measure elastic and shear moduli, and compute Poisson’s ratio. The coefficient of static sliding friction was measured using the slip-stick method, and the static rolling friction was determined from measured shear on particles positioned between two hot-pressed plates. The coefficient of restitution was measured by dropping particles on a surface and determining the kinetic energy before and after impact with the surface using high resolution particle tracking velocimetry. Particle size did not significantly impact the coefficients of restitution and static rolling friction, however, particle shape distribution resulted in a large variation in measurements. An inclined flow experiment was performed to characterize granular flows of Carbobead CP 30/60 particles using particle image velocimetry. Numerical models of the experiment using discrete element method were generated with the measured mechanical properties as inputs for comparison with experimental results. A constant directional torque rolling friction model best predicted bulk granular flow behavior. Furthermore, good agreement between the model and experiment was achieved at ambient, steady state conditions, with average velocity differences <10%.

14 SOLAR ENERGY↗

Continuous integration data-driven platform of industrial-scale subsurface storage for real-time analytics

This project helped address the growing need for efficient and scalable models to support geological carbon and energy storage, which are crucial for achieving net-zero emissions. Traditionally accurate high-fidelity numerical models have been used to simulate relevant storage processes under a handful of processes, however such models are computationally demanding, making uncertainty quantification impractical. Consequently, we first developed a machine learning framework, based on Graph Neural Operators (GNOs), to improving the accuracy of model predictions for a fixed computational budget. We then developed an Ensemble of Improved Neural Operators (ENO), which uses bagging and Monte Carlo dropout techniques, to further improve prediction accuracy. Lastly, we developed the way to explain progressive transfer learning methods to reduce the amount of training data and computational cost of training (i.e., reduce trainable parameters) when using our models for multiple storage sites. Our numerical investigation, which used real-world case studies, demonstrated that our framework can significantly improve the safety and efficiency of geological storage operations, with potential applications in other domains such as geothermal reservoirs and climate modeling.

54 ENVIRONMENTAL SCIENCES↗

3D micromechanical simulation of PBX composites

Previous research studied the constitutive response and interface strength of a bonded bicrystal system containing two HMX explosive grains bound with a HTPB polymer binder by performing FE simulations on real geometry imaged with micro Computed Tomography (µCT). The parameters generated from this past study were successfully applied to 2D mesoscale simulations of µCT imaged HMX-HTPB polycrystal samples. Presently, the mechanical response of PBXs containing a Nitro-Plasticized Estane (NPE) binder with similar microstructural geometry to the HMX-HTPB polycrystal samples were studied. Imaging of a HMX-NPE bicrystal sample allowed for simulation of the new material pairing. However, imaging using µCT techniques on the new polycrystal samples was difficult due to the decrease in contrast between the HMX grains and the NPE binder. Mesoscale polycrystal simulations utilizing a single µCT imaged microstructure were performed on a 3D representative volume element (RVE) comparing the delamination behavior and mechanical response of both HMX-HTPB and HMX-NPE PBXs. Furthermore, this technique shows the ability to virtually explore the mechanical response of a range of hypothetical materials that share common microstructural geometric characteristics in addition to studying numerically altered geometries.

36 MATERIALS SCIENCE↗

Audible enclaves crafted by nonlinear self-bending ultrasonic beams

Delivering audible content to a targeted listener without disturbing others is paramount in audio engineering. However, achieving this goal has long been challenging due to the diffraction of low-frequency (long-wavelength) audio waves in linear acoustics. Here, we introduce an approach for creating remote audio spots, dubbed audible enclaves, by harnessing the local nonlinear interaction of two self-bending ultrasonic beams with distinct spectra. The self-bending ultrasonic beams created by acoustic metasurfaces, though inaudible, can bypass obstacles such as human heads. At their intersection behind obstacles, highly localized audible enclaves are formed due to the local nonlinear interactions. Additionally, we demonstrate the ultrabroadband capabilities of our metasurface-based implementation both numerically and experimentally, spanning from 125 Hz to 4 kHz (6 octave bands), covering the majority of the audible frequency range. The practicality of our proposed technique is underscored by its compact implementation size (0.16 m, equivalent to 0.06 wavelengths at 125 Hz), as well as its robust performance under wideband transient audio signal excitation and in a common room with reverberations. Our proposed audible enclaves hold significant potential for various applications in advanced audio engineering, including private speech communications, immersive spatial audio reproduction, and high-resolution sound/quiet zone control.

Science & Technology - Other Topics↗

Adiabatic cooling of planar motion in a Penning-trap ion crystal to sub-millikelvin temperatures

Two-dimensional planar ion crystals in a Penning trap are a platform for quantum information science experiments. However, the low-frequency planar modes of these crystals are not efficiently cooled by laser cooling, which can limit the utility of the drumhead modes for quantum information processing. Recently, it has been shown that nonlinear mode coupling can enhance the cooling of the low-frequency planar modes. Here, in this study, we demonstrate in numerical simulations that this coupling can be dynamically tuned by adiabatically changing the rotation frequency of the ion crystal during experiments. Furthermore, we show that this technique can, in addition, produce lower temperatures for the low-frequency planar modes via an adiabatic cooling process. This result allows cooling of the planar modes to sub-millikelvin temperatures, resulting in improved spectral resolution of the drumhead modes at experimentally relevant rotation frequencies, which is crucial for quantum information processing applications.

Johnson, Wes [Univ. of Colorado, Boulder, CO (Unit↗

Microscopic Theory of Fluctuating Hydrodynamics in Nonlinear Lattices

The theory of fluctuating hydrodynamics has been an important tool for analyzing macroscopic behavior in nonlinear lattices. However, despite its practical success, its microscopic derivation is still incomplete. In this work, we provide the microscopic derivation of fluctuating hydrodynamics, using the coarse-graining and projection technique; the equivalence of ensembles turns out to be critical. The Green-Kubo (GK)-like formula for the bare transport coefficients are presented in a numerically computable form. Our numerical simulations show that the bare transport coefficients exist for a sufficiently large but finite coarse-graining length in the infinite lattice within the framework of the GK-like formula. Finally, this demonstrates that the bare transport coefficients uniquely exist for each physical system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multicarrier Spread Spectrum Communications With Noncontiguous Subcarrier Bands for HF Skywave Links

Existing high-frequency (HF) radio platforms offer robust performance against the volatile HF propagation channel. However, the growing traffic across the band contests the reliability of these systems. While techniques to mitigate the effects of narrowband interference have been thoroughly explored, they are insufficient against wideband interference or when the transmission band is occupied by numerous scattered users. To improve reliability in these congested channel conditions, we propose a filter-bank based multicarrier spread-spectrum waveform with noncontiguous subcarrier bands. Using noncontiguous subcarrier bands enables the system to at once leverage the robustness of a wideband system while retaining the frequency agility of a narrowband system. In this study, we modify a filter-bank transmitter structure to accommodate noncontiguous subcarrier bands and consider several immediate impacts of this change, such as elevated peak-to-average-power ratios (PAPRs). A receiver architecture to process the noncontiguous spread-spectrum signal is also introduced, along with details regarding wideband channel estimation. Finally, we develop efficient transmitter and receiver structures to support practical system implementations. We conclude by comparing the performance of contiguous and noncontiguous systems through both simulation and over-the-air testing. The results show that the noncontiguous system remains robust in typical HF channels while significantly outperforming the contiguous system in congested spectral conditions.

(PAPR↗

Morphogene-assisted transformation of Sorghum bicolor allows more efficient genome editing

Sorghum bicolor (L.) Moench, the fifth most important cereal worldwide, is a multi-use crop for feed, food, forage and fuel. To enhance the sorghum and other important crop plants, establishing gene function is essential for their improvement. For sorghum, identifying genes associated with its notable abiotic stress tolerances requires a detailed molecular understanding of the genes associated with those traits. The limits of this knowledge became evident from our earlier in-depth sorghum transcriptome study showing that over 40% of its transcriptome had not been annotated. Here, we describe a full spectrum of tools to engineer, edit, annotate and characterize sorghum’s genes. Efforts to develop those tools began with a morphogene-assisted transformation (MAT) method that led to accelerated transformation times, nearly half the time required with classical callus-based, non-MAT approaches. These efforts also led to expanded numbers of amenable genotypes, including several not previously transformed or historically recalcitrant. Another transformation advance, termed altruistic, involved introducing a gene of interest in a separate Agrobacterium strain from the one with morphogenes, leading to plants with the gene of interest but without morphogenes. The MAT approach was also successfully used to edit a target exemplary gene, phytoene desaturase. To identify single-copy transformed plants, we adapted a high-throughput technique and also developed a novel method to determine transgene independent integration. These efforts led to an efficient method to determine gene function, expediting research in numerous genotypes of this widely grown, multi-use crop.

54 ENVIRONMENTAL SCIENCES↗