Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer experiment”

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 973 records · Page 54

Polymer Compositions in Critical Experiments: Possibly Not What You Think

The Chlorine Worth Study (CWS) Experiment was performed at the National Criticality Experiments Research Center (NCERC) in December 2021. Its goal was to provide a validation benchmark experiment with chlorine and plutonium in the thermal neutron spectrum. This purpose is necessary to reduce the margin of subcriticality in aqueous chloride operations at the Los Alamos National Laboratory Plutonium Facility (PF-4). Reducing the margin of subcriticality will enable higher throughput, required to meet NNSA mission needs. The experiment used layers of plutonium plates, polyethylene (HDPE), aluminum, and polyvinyl chloride (PVC) or chlorinated polyvinyl chloride (CPVC). They were optimized to match plutonium-chloride solutions of 30 g/L, 300 g/L, and 600 g/L. Upon completion of the experiment, a International Criticality Safety Benchmark Evaluation Program (ICSBEP) report was immediately started. The compositions of all materials were assumed pure unless additional information was known (such as for the plutonium plates). During the benchmark analysis, the assumed CPVC composition was questioned. The follow-on work led to lessons learned on compositions in benchmarks. Materials, and specifically polymers, are often much more complex than a basic chemical formula. The CWS experiment is used as an example in this paper to document the lessons learned.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Sandia Electromagnetic Environment Simulator Facilities

Sandia National Laboratories’ (Sandia’s) Electrical Sciences Group (1350) provides a spectrum of solutions for electromagnetic (EM) environment effects on electrical systems to advance physical understanding including experiments, analytical and numerical model development leveraging first-principles physics, and high-performance computational modeling and simulation (M&S). The six departments in the Electrical Sciences Group provide analysis, design guidance, and experiments in support of nuclear weapons qualification in normal, abnormal, and hostile environments, as well as research and development for advanced electrical systems that can operate through these environments. This document is intended to provide a quick summary of Sandia’s EM experimental facilities, most of which provide unique national capabilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular-level design of alternative media for energy-saving pilot-scale fibrillation of nanocellulose

The outstanding mechanical properties, light weight, and biodegradability of cellulose nanofibrils (CNFs) make them promising components of renewable and sustainable next-generation reinforced composite biomaterials and bioplastics. Manufacturing CNFs at a pilot scale requires disc-refining fibrillation of dilute cellulose fibers in aqueous pulp suspensions to shear the fibers apart into their nanodimensional forms, which is, however, an energy-intensive process. Here, in this study, we used atomistic molecular dynamics (MD) simulation to examine media that might facilitate the reduction of interactions between cellulose fibers, thereby reducing energy consumption in fibrillation. The most suitable medium found by the simulations was an aqueous solution with 0.007:0.012 wt.% NaOH:urea, and indeed this was found in pilot-scale experiments to reduce the fibrillation energy by ~21% on average relative to water alone. The NaOH:urea-mediated CNFs have similar crystallinity, morphology, and mechanical strength to those formed in water. The NaOH and urea act synergistically on CNFs to aid fibrillation but at different length scales. NaOH deprotonates hydroxyl groups leading to mesoscale electrostatic repulsion between fibrils, whereas urea forms hydrogen bonds with protonated hydroxyl groups thus disrupting interfibril hydrogen bonds. This suggests a general mechanism in which an aqueous medium that contains a strong base and a small organic molecule acting as a hydrogen-bond acceptor and/or donor may be effectively employed in materials processes where dispersion of deprotonable polymers is required. The study demonstrates how atomic-detail computer simulation can be integrated with pilot-scale experiments in the rational design of materials processes for the circular bioeconomy.

36 MATERIALS SCIENCE↗

Charge Conservation beyond Uniformity: Spatially Inhomogeneous Electromagnetic Response in Periodic Solids

Nonlinear electromagnetic response functions have reemerged as a crucial tool for studying quantum materials, due to recently appreciated connections between optical response functions, quantum geometry, and band topology. Most attention has been paid to responses to spatially uniform electric fields, relevant to low-energy optical experiments in conventional solid state materials. However, magnetic and magnetoelectric phenomena are naturally connected by responses to spatially varying electric fields due to Maxwell’s equations. Furthermore, in the emerging field of moiré materials, characteristic lattice scales are much longer, allowing spatial variation of optical electric fields to potentially have a measurable effect in experiments. In order to address these issues, we develop a formalism for computing linear and nonlinear responses to spatially inhomogeneous electromagnetic fields. Starting with the continuity equation, we derive an expression for the second-quantized current operator that is manifestly conserved and model independent. Crucially, our formalism makes no assumptions on the form of the microscopic Hamiltonian and so is applicable to model Hamiltonians derived from tight-binding or calculations. We then develop a diagrammatic Kubo formalism for computing the wave vector dependence of linear and nonlinear conductivities, using Ward identities to fix the value of the diamagnetic current order by order in the vector potential. We apply our formula to compute the magnitude of the Kerr effect at oblique incidence for a model of a moiré-Chern insulator and demonstrate the experimental relevance of spatially inhomogeneous fields in these systems. We further show how our formalism allows us to compute the (orbital) magnetic multipole moments and magnetic susceptibilities in insulators. Turning to nonlinear response, we use our formalism to compute the second-order transverse response to spatially varying transverse electric fields in our moiré-Chern insulator model, with an eye toward the next generation of experiments in these systems. Published by the American Physical Society 2024

Physics↗

In situ compression artifact removal in scientific data using deep transfer learning and experience replay

The massive amount of data produced during simulation on high-performance computers has grown exponentially over the past decade, exacerbating the need for streaming compression and decompression methods for efficient storage and transfer of this data---key to realizing the full potential of large-scale computational science. Lossy compression approaches such as JPEG when applied to scientific simulation data realized as a stream of images can achieve good compression rates but at the cost of introducing compression artifacts and loss of information. This paper develops a unified framework for in situ compression artifact removal in which the fully convolutional neural network architectures are combined with scalable training, transfer learning, and experience replay to achieve superior accuracy and efficiency while significantly decreasing the storage footprint as compared with the traditional optimization-based approaches. We demonstrate the proposed approach and compare it with compressed sensing postprocessing and other baseline deep learning models using climate simulations and nuclear reactor simulations, both of which are driven by hyperbolic partial differential equations. Our approach when applied to remove the compression artifacts on the JPEG-compressed nuclear reactor simulation data (using a transfer-trained model that was pretrained on the climate simulation data and updated incrementally as the nuclear reactor simulation progressed), achieved a significant improvement---mean peak signal-to-noise ratio of 42.438 as compared with 27.725 obtained with the compressed sensing approach.

97 MATHEMATICS AND COMPUTING↗

High-pressure phase transition of olivine-type Mg 2 GeO 4 to a metastable forsterite-III type structure and their equations of state

Germanates are often used as structural analogs of planetary silicates. We have explored the high-pressure phase relations in Mg 2 GeO 4 using diamond-anvil cell experiments combined with synchrotron X-ray diffraction and computations based on density functional theory. Upon room temperature compression, forsterite-type Mg 2 GeO 4 remains stable up to 30 GPa. At higher pressures, a phase transition to a forsterite-III type (Cmc2 1 ) structure was observed, which remained stable to the peak pressure of 105 GPa. Using a third-order Birch Murnaghan fit to the experimental data, we obtained V 0 = 305.1(3) Å3, K 0 = 124.6(14) GPa, and $K'_0$ = 3.86 (fixed) for forsterite-type Mg 2 GeO 4 and V 0 = 263.5(15) Å 3 , K 0 = 175(7) GPa, and $K'_0$ = 4.2 (fixed) for the forsterite-III type phase. The forsterite-III type structure was found to be metastable when compared to the stable assemblage of perovskite/post-perovskite + MgO, as observed during laser-heating experiments. Understanding the phase relations and physical properties of metastable phases is crucial for studying the mineralogy of impact sites, understanding metastable wedges in subducting slabs, and interpreting the results of shock compression experiments.

58 GEOSCIENCES↗

A multifidelity Bayesian optimization method for inertial confinement fusion design

Due to their cost, experiments for inertial confinement fusion (ICF) heavily rely on numerical simulations to guide design. As simulation technology progresses, so too can the fidelity of models used to plan for new experiments. However, these high-fidelity models are by themselves insufficient for optimal experimental design, because their computational cost remains too high to efficiently and effectively explore the numerous parameters required to describe a typical experiment. Therefore, traditionally, ICF design has relied on low-fidelity modeling to initially identify potentially interesting design regions, which are then subsequently explored via selected high-fidelity modeling. In this paper, we demonstrate that this two-step approach can be insufficient: even for simple design problems, a two-step optimization strategy can lead high-fidelity searching toward incorrect regions and consequently waste computational resources on parameter regimes far away from the true optimal solution. We reveal that a primary cause of this behavior in ICF design problems is the presence of low-fidelity optima in different regions of the parameter space far away from high-fidelity optima. To address this issue, we propose an iterative multifidelity Bayesian optimization method based on Gaussian Process Regression that leverages both low- and high-fidelity models simultaneously. We demonstrate, using both two- and eight-dimensional ICF test problems, that our algorithm can effectively utilize both low-fidelity and high-fidelity models to refine the designs. This approach proves to be more efficient than relying solely on high-fidelity modeling for optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimizing accuracy and efficacy in data-driven materials discovery for the solar production of hydrogen

The production of hydrogen fuels, via water splitting, is of practical relevance for meeting global energy needs and mitigating the environmental consequences of fossil-fuel-based transportation. Water photoelectrolysis has been proposed as a viable approach for generating hydrogen, provided that stable and inexpensive photocatalysts with conversion efficiencies over 10% can be discovered, synthesized at scale, and successfully deployed. While a number of first-principles studies have focused on the data-driven discovery of photocatalysts, in the absence of systematic experimental validation, the success rate of these predictions may be limited. We address this problem by developing a screening procedure with co-validation between experiment and theory to expedite the synthesis, characterization, and testing of the computationally predicted, most desirable materials. Here, starting with 70 150 compounds in the Materials Project database, the proposed protocol yielded 71 candidate photocatalysts, 11 of which were synthesized as single-phase materials. Experiments confirmed hydrogen generation and favorable band alignment for 6 of the 11 compounds, with the most promising ones belonging to the families of alkali and alkaline-earth indates and orthoplumbates. This study shows the accuracy of a nonempirical, Hubbard-corrected density-functional theory method to predict band gaps and band offsets at a fraction of the computational cost of hybrid functionals, and outlines an effective strategy to identify photocatalysts for solar hydrogen generation.

08 HYDROGEN↗

Propagation of positive discharges in an air bubble having an embedded water droplet

Low temperature plasmas generated inside gas bubbles immersed in water is an effective method of rapidly transferring plasma generated reactive species to the water for applications in biomedicine, agriculture and environment. Reactive species are generally produced in the gas phase plasma and then solvate into the liquid. The large surface-to-volume ratio (SVR) of the bubble accelerates this process. In generating bubbles in water, aerosols and droplets are also contained within the bubble. These droplets also have a large SVR and so can be rapidly plasma activated. However, the presence of the droplets can also impact the propagation of the plasma in the bubble. In this paper, results are discussed from computational and experimental investigations of the formation and evolution of discharges in an air bubble immersed in water with an embedded water droplet. The computations were performed with a two-dimensional plasma hydrodynamics model. Experiments were performed with a quasi-2D bubble apparatus. In bubbles having a droplet, a plasma filament typically bridges from the powered electrode to the droplet, and then from the droplet to the bubble surface. A surface-hugging streamer also occurs on the inner bubble surface and on the surface of the droplet. Both surface streamers result in part from surface charge accumulation and can dominate the formation of reactive species that transport into the droplet. Increasing droplet conductivity suppresses propagation of the surface discharge and leads to a lower density of aqueous reactive species. Increasing conductivity of the surrounding water does not change the overall structure of the discharge but does slightly elevate the discharge intensity. Finally, the size and shape of the embedded droplet can significantly affect the formation and propagation of the streamer.

2D fluid modelling↗

Reactor physics benchmark experiments at the JSI TRIGA MARK II reactor - Current status and future outlook

Full text of publication follows. With the development of new high-fidelity computational methods, improvement of nuclear data, and multiphysics modelling, there is an increased need for benchmark experiments to experimentally validate the models, methods and input data. Many of the nuclear facilities designed to perform reactor physics benchmark experiments have been shut down. Therefore, research reactors offer a great opportunity for benchmark experiments, if they are well designed and performed with great care and accuracy. In this presentation we provide an overview of the past and ongoing activities related to benchmark experiments at the Jozef Stefan Institute TRIGA Mark II research reactor. The following experiments have been performed: criticality with fresh fuel, {sup 197}Au(n,γ) and {sup 27}Al(n,α) reaction rates in irradiation channels, absolute and relative {sup 197}Au(n,γ), {sup 235}U(n,f) and {sup 238}U(n,f) reaction rates in the core, burnup, kinetic parameters, control rod worth, isothermal reactivity coefficient, self-shielding, slow and fast (pulse) transients, nuclear heating, delayed and prompt gamma ray production, temperature profiles for multi-physics. Since the existing fleet of research reactors is ageing very rapidly and new experiments are needed, new research reactors should be designed and built to meet the needs of future advanced reactors, education and training, and other technologies in the coming years. We will review planned activities at the JSI TRIGA reactors and plans for the new research reactor in Slovenia. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental measurements of gas-particle flows in large-scale strippers

Measurements of gas–solid particle flow in a large-scale stripper unit are reported. The experiments target industrial (pilot) scale measurements of gas-particle flow that can be used to validate numerical methods for modeling multiphase flows, such as computational fluid dynamics coupled to the discrete element method (CFD-DEM). Specifically, experiments were performed with Geldart Group B (533 µm) glass beads in a stripper unit. In general, previous CFD-DEM validation datasets were performed in benchtop experimental systems and limited to O(10 5 ) particles. Here, in this work, the number of particles in the isolated section of the stripper is O(10 7 ) for all tests. Measurements are reported for a ~2 m tall section of the ~1 m diameter stripper. Radial mass flux profiles measured at five axial positions and three axial pressure drop measurements are reported for six operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Requirements, Status, and Plans for Track Reconstruction at the sPHENIX Experiment

sPHENIX is a new experiment that is being constructed at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory. The primary physics goals of sPHENIX are to measure jets, their substructure, and the upsilon resonances in p+p, p+Au and Au+Au collisions. To realize these goals, a tracking system composed of a time projection chamber and several silicon detectors will be used to identify tracks that correspond to jets and upsilon decays. However, the sPHENIX experiment will collect approximately 200 PB of data utilizing a finite-sized computing center; thus, performing track reconstruction in a timely manner is a challenge due to the large occupancy of heavy-ion collisions. The sPHENIX experiment, its track reconstruction, and the need for implementing faster track-fitting algorithms, such as that provided by the A Common Tracking Software package, into the sPHENIX software stack are discussed.

Osborn, Joe↗

Integration of the NCRC Database and Other INL Databases

The Nuclear Computational Resource Center provides a portal by which industry professionals, educational staff, students, national laboratory employees, and others may request access to certain engineering software tools. As the tools provided through the Nuclear Computational Resource Center portal are not open-source and freely available, a set of approvals are necessary before access is granted. All code recipients must be associated with an institution that has a license with Idaho National Laboratory for the code requested. Information about these licenses is controlled by Idaho National Laboratory’s Technology Deployment organization and housed in a Technology Deployment database. Those requesting code access who are not citizens of the United States must also have a security plan, mandated by Idaho National Laboratory policy. Security plans are managed by the International Access Program and are stored in an International Access Program database known as IFacts. Granting access to software thus depends on information stored in the Technology Deployment database and IFacts. In the past, no connection between the Nuclear Computational Resource Center portal and these databases existed, making checking the status of license agreements and security plans time consuming and error prone. This report demonstrates that the Nuclear Computational Resource Center portal now connects to both the Technology Deployment database and IFacts, greatly improving the ease of use of the Nuclear Computational Resource Center system for administrators, which leads to a better overall experience for those requesting code access.

99 GENERAL AND MISCELLANEOUS↗

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗

Localization of fluorescent targets in deep tissue with expanded beam illumination for studies of cancer and the brain

Imaging fluorescence through millimeters or centimeters of tissue has important in vivo applications, such as guiding surgery and studying the brain. Often, the important information is the location of one of more optical reporters, rather than the specifics of the local geometry, motivating the need for a localization method that provides this information. We present an optimization approach based on a diffusion model for the fast localization of fluorescent inhomogeneities in deep tissue with expanded beam illumination that simplifies the experiment and the reconstruction. We show that the position of a fluorescent inhomogeneity can be estimated while assuming homogeneous tissue parameters and without having to model the excitation profile, reducing the computational burden and improving the utility of the method. We perform two experiments as a demonstration. First, a tumor in a mouse is localized using a near infrared folate-targeted fluorescent agent (OTL38). This result shows that localization can quickly provide tumor depth information, which could reduce damage to healthy tissue during fluorescence-guided surgery. Second, another near infrared fluorescent agent (ATTO647N) is injected into the brain of a rat, and localized through the intact skull and surface tissue. This result will enable studies of protein aggregation and neuron signaling.

60 APPLIED LIFE SCIENCES↗

Evaluating Portable Parallelization Strategies for Heterogeneous Architectures in High Energy Physics

High-energy physics (HEP) experiments have developed millions of lines of code over decades that are optimized to run on traditional x86 CPU systems. However, we are seeing a rapidly increasing fraction of floating point computing power in leadership-class computing facilities and traditional data centers coming from new accelerator architectures, such as GPUs. HEP experiments are now faced with the untenable prospect of rewriting millions of lines of x86 CPU code, for the increasingly dominant architectures found in these computational accelerators. This task is made more challenging by the architecture-specific languages and APIs promoted by manufacturers such as NVIDIA, Intel and AMD. Producing multiple, architecture-specific implementations is not a viable scenario, given the available person power and code maintenance issues. The Portable Parallelization Strategies team of the HEP Center for Computational Excellence is investigating the use of Kokkos, SYCL, OpenMP, std::execution::parallel and alpaka as potential portability solutions that promise to execute on multiple architectures from the same source code, using representative use cases from major HEP experiments, including the DUNE experiment of the Long Baseline Neutrino Facility, and the ATLAS and CMS experiments of the Large Hadron Collider. This cross-cutting evaluation of portability solutions using real applications will help inform and guide the HEP community when choosing their software and hardware suites for the next generation of experimental frameworks. We present the outcomes of our studies, including performance metrics, porting challenges, API evaluations, and build system integration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗