Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

Learning nonlinear level sets for dimensionality reduction in function approximation

We developed a Nonlinear Level-set Learning (NLL) method for dimensionality reduction in high-dimensional function approximation with small data. This work is motivated by a variety of design tasks in real-world engineering applications, where practitioners would replace their computationally intensive physical models (e.g., high-resolution fluid simulators) with fast-to-evaluate predictive machine learning models, so as to accelerate the engineering design processes. There are two major challenges in constructing such predictive models: (a) high-dimensional inputs (e.g., many independent design parameters) and (b) small training data, generated by running extremely time-consuming simulations. Thus, reducing the input dimension is critical to alleviate the over-fitting issue caused by data insufficiency. Existing methods, including sliced inverse regression and active subspace approaches, reduce the input dimension by learning a linear coordinate transformation; our main contribution is to extend the transformation approach to a nonlinear regime. Specifically, we exploit reversible networks (RevNets) to learn nonlinear level sets of a high-dimensional function and parameterize its level sets in low-dimensional spaces. A new loss function was designed to utilize samples of the target functions' gradient to encourage the transformed function to be sensitive to only a few transformed coordinates. The NLL approach is demonstrated by applying it to three 2D functions and two 20D functions for showing the improved approximation accuracy with the use of nonlinear transformation, as well as to an 8D composite material design problem for optimizing the buckling-resistance performance of composite shells of rocket inter-stages.

Zhang, Guannan↗

Component-wise reduced-order model design optimization such as for lattice design optimization

Systems and methods for optimizing a lattice structure design are disclosed herein. In some embodiments, a method for optimizing a lattice structure design can include (i) modeling the lattice structure with a component-wise reduced-order model (CWROM) and (ii) optimizing the CWROM based on a selected criterion using a topology optimization algorithm for lattice design. The selected criterion can include a boundary condition and a load applied to the lattice structure. By modeling the lattice structure as a CWROM, the optimization process can be very fast while still permitting the accurate computation of physical quantities of the lattice structure.

Choi, Youngsoo↗

Topological terms with qubit regularization and relativistic quantum circuits

Qubit regularization provides a rich framework to explore quantum field theories. The freedom to choose how the important symmetries of the theory are embedded in the qubit regularization scheme allows us to construct new lattice models with rich phase diagrams. Some of the phases can contain topological terms which lead to critical phases. In this work we introduce and study the SU(3)-F qubit regularization scheme to embed the SO(3) spin-symmetry. We argue that qubit models in this regularization scheme contain several phases including a critical phase which describes the k = 1 Wess-Zumino-Witten (WZW) conformal field theory (CFT) at long distances, and two massive phases one of which is trvially gapped and the other which breaks the lattice translation symmetry. We construct a simple space-time Euclidean lattice model with a single coupling U and study it using the Monte Carlo method. We show the model has a critical phase at small U and a trivially massive phase at large U with a first order transition separating the two. Another feature of our model is that it is symmetric under space-time rotations, which means the temporal and spatial lattice spacing are connected to each other. The unitary time evolution operator obtained by a Wick rotation of the transfer matrix of our model can help us compute the physics of the k = 1 WZW CFT in real time without the need for tuning the temporal lattice spacing to zero. We use this idea to introduce the concept of a relativistic quantum circuit on a discrete space-time lattice.

Bhattacharya, Tanmoy↗

Recent advances in the Mercury Monte Carlo particle transport code

We review recent physics and computational science advances in the Mercury Monte Carlo particle transport code under development at Lawrence Livermore National Laboratory. We describe recent efforts to enable a nuclear resonance fluorescence capability in the Mercury photon transport. We also describe recent work to implement a probability of extinction capability into Mercury. We review the results of current parallel scaling and threading efforts that enable the code to run on millions of MPI processes. (authors)

97 MATHEMATICS AND COMPUTING↗

Development of Fusion Neutronics Methodology, Software, and Testing Facility in FDS Team

Fusion neutronics including methodology, software, experiments, etc. is studied concerning the generation, utilizing and protection of neutrons in fusion facilities. Achievements have been made in the past decades to solve the new problems fusion neutronics. However there is still some way to go such as following. Automatic modeling method based on CAD has been developed to avoid the heavy burden of manual geometry description. The pre-processing costs much time to simplify the CAD model for converting to Constructive Solid Geometry (CSG), 2D or 3D mesh based geometry. This problem has become a new focus. Deep penetration of radiation shielding is notable problem for fusion. Adaptive variance reduction methods such as Consistent Adjoint Driven Importance Sampling (CADIS) were developed to accelerate convergence. However the gap still exists between engineering requirements and current methodology. The geometry simplification, homogenizing based tally, error in nuclear data and error propagation in coupled calculation, etc., causes considerable uncertainty, which leads to 'safety factor' applied in simulation results to achieve required safety margin. It can significantly raise the cost and even makes it impractical to initiate the fusion facilities. Uncertainty qualification, error analysis and error propagation problem should be studied. New methodology and software should be validated. For example the (n,p), (n,a), (n,t) cross section of Cu, Fe, W, Ni which is important for evaluation of material radiation damage is absent. The validation of neutronics properties of special components such as first wall, blanket and coil are needed. The D-T fusion neutron generator can directly represent the neutron environment in fusion reactors. However the neutron yield is several orders of magnitudes lower than requirements. To resolve obstacles in fusion neutronics, broad improvements of methodology, software and testing facility have been achieved by FDS Team in recent years and the main progress was introduced in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics

A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.

Process Modeling↗

Earthquake Ground Motion Simulations on the Sierra High-Performance Computing System

Physics-based three-dimensional (3D) numerical simulations of earthquake ground motions provide an alternative to reliance on limited and highly scattered empirical data and models for seismic hazard assessment. We performed fully 3D earthquake ground motion simulations on the Sierra and Lassen high-performance computing (HPC) platforms at Livermore Computing. These new HPC platforms feature graphics processing unit (GPU) acceleration for improved performance. Ground motions were computed with SW4, a summation-by-parts finite difference code for anelastic wave propagation. SW4 was recently enhanced to run efficiently on GPU-accelerated platforms. We found that SW4 runs almost fifty (50) times faster in terms of node-hours on new GPU-accelerated platforms compared to conventional central processing unit (CPU) platforms with agreement of computed ground motion time-series to the full single precision which they are written. We focused on moment magnitude, M W , 7.0 earthquakes on the Hayward Fault. We were able to run three simulations resolving ground motions from static displacements (0 Hz) to frequencies of 10 Hz (wavelengths as short as 50 meters) making these Sierra calculations the highest resolution regional-scale earthquake simulations ever made for California. Results show good agreement between simulated ground motion intensities (peak ground velocity and acceleration, spectral response) and empirical models, including frequencies above 10 Hz. On Lassen we were able to run a suite of 52 M W 7.0 Hayward Fault earthquake simulations resolved 0-5 Hz for both plane-layered, one-dimensional (1D) and 3D Earth models. These simulations were used to develop a methodology for characterizing and removing systematic path and site effects which bias ground motion intensities. These simulations using the newly enhanced SW4 code allow us to run higher resolution seismic simulations with shorter run times, providing a new capability for seismic hazard and risk studies. Results were presented at scientific conferences and workshops and two publications are being prepared based on work performed in this project.

58 GEOSCIENCES↗

In Search of the Physics: The Interplay of Experiment and Computation in Airframe Noise Research: Flap-Edge Noise

The LaRC investigative process for airframe noise has proven to be a useful guide for elucidation of the physics of flow-induced noise generation over the last five years. This process, relying on a close interplay between experiment and computation, is described and demonstrated here on the archetypal problem of flap-edge noise. Some detailed results from both experiment and computation are shown to illustrate the process, and a description of the multi-source physics seen in this problem is conjectured.

Streett, C. L.↗

Data Storage for HEP Experiments in the Era of High-Performance Computing

As particle physics experiments push their limits on both the energy and the intensity frontiers, the amount and complexity of the produced data are also expected to increase accordingly. With such large data volumes, next-generation efforts like the HL-LHC and DUNE will rely even more on both high-throughput (HTC) and high-performance (HPC) computing clusters. Full utilization of HPC resources requires scalable and efficient data-handling and I/O. For the last few decades, ROOT has been used by most HEP experiments to store data. However, other storage technologies like HDF5 may perform better in HPC environments. Initial explorations with HDF5 have begun using ATLAS, CMS and DUNE data; the DUNE experiment has also adopted HDF5 for its data-acquisition system. This paper presents the future outlook of the HEP computing and the role of HPC, and a summary of ongoing and future works to use HDF5 as a possible data storage technology for the HEP experiments to use in HPC environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Conference on Real-Time Computer Applications in Nuclear, Particle and Plasma Physics, 6th, Williamsburg, VA, May 15-19, 1989, Proceedings

Papers on real-time computer applications in nuclear, particle, and plasma physics are presented, covering topics such as expert systems tactics in testing FASTBUS segment interconnect modules, trigger control in a high energy physcis experiment, the FASTBUS read-out system for the Aleph time projection chamber, a multiprocessor data acquisition systems, DAQ software architecture for Aleph, a VME multiprocessor system for plasma control at the JT-60 upgrade, and a multiasking, multisinked, multiprocessor data acquisition front end. Other topics include real-time data reduction using a microVAX processor, a transputer based coprocessor for VEDAS, simulation of a macropipelined multi-CPU event processor for use in FASTBUS, a distributed VME control system for the LISA superconducting Linac, a distributed system for laboratory process automation, and a distributed system for laboratory process automation. Additional topics include a structure macro assembler for the event handler, a data acquisition and control system for Thomson scattering on ATF, remote procedure execution software for distributed systems, and a PC-based graphic display real-time particle beam uniformity.

Pordes, Ruth↗

Optimization and supervised machine learning methods for fitting numerical physics models without derivatives

Here, we address the calibration of a computationally expensive nuclear physics model for which derivative information with respect to the fit parameters is not readily available. Of particular interest is the performance of optimization-based training algorithms when dozens, rather than millions or more, of training data are available and when the expense of the model places limitations on the number of concurrent model evaluations that can be performed. As a case study, we consider the Fayans energy density functional model, which has characteristics similar to many model fitting and calibration problems in nuclear physics. We analyze hyperparameter tuning considerations and variability associated with stochastic optimization algorithms and illustrate considerations for tuning in different computational settings.

97 MATHEMATICS AND COMPUTING↗