Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATOR”

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 325 records · Page 18

SU-D-BRF-06: A Brachytherapy Simulator with Realistic Haptic Force Feedback and Real-Time Ultrasounds Image Simulation for Training and Teaching

Purpose: Surgical procedures require dexterity, expertise and repetition to reach optimal patient outcomes. However, efficient training opportunities are usually limited. This work presents a simulator system with realistic haptic force-feedback and full, real-time ultrasounds image simulation. Methods: The simulator is composed of a custom-made Linear-DELTA force-feedback robotic platform. The needle tip is mounted on a force gauge at the end effector of the robot, which responds to needle insertion by providing reaction forces. 3D geometry of the tissue is using a tetrahedral finite element mesh (FEM) mimicking tissue properties. As the needle is inserted/retracted, tissue deformation is computed using a mass-tensor nonlinear visco-elastic FEM. The real-time deformation is fed to the L-DELTA to take into account the force imparted to the needle, providing feedback to the end-user when crossing tissue boundaries or needle bending. Real-time 2D US image is also generated synchronously showing anatomy, needle insertion and tissue deformation. The simulator is running on an Intel I7 6- core CPU at 3.26 MHz. 3D tissue rendering and ultrasound display are performed on a Windows 7 computer; the FEM computation and L-DELTA control are executed on a similar PC using the Neutrino real-time OS. Both machines communicate through an Ethernet link. Results: The system runs at 500 Hz for a 8333-tetrahedron tissue mesh and a 100-node angular spring needle model. This frame rate ensures a relatively smooth displacement of the needle when pushed or retracted (±20 N in all directions at speeds of up to 2 m/s). Unlike commercially-available haptic platforms, the oblong workspace of the L-DELTA robot complies with that required for brachytherapy needle displacements of 0.1m by 0.1m by 0.25m. Conclusion: We have demonstrated a real-life, realistic brachytherapy simulator developed for prostate implants (LDR/HDR). The platform could be adapted to other sites or training for other types of needle-based procedures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Simulated Polarimetric Fields of Ice Vapor Growth Using the Adaptive Habit Model. Part I: Large-Eddy Simulations

Abstract The bulk adaptive habit model (AHM) explicitly predicts ice particle aspect ratio, improving the representation of microphysical processes and properties, including ice–liquid-phase partitioning. With the unique ability to predict ice particle shape and density, the AHM is combined with an offline forward operator to produce fields of simulated polarimetric variables. An evaluation of AHM-forward-simulated dual-polarization radar signatures in an idealized Arctic mixed-phase cloud is presented. Interpretations of those signatures are provided through microphysical model output using the large-eddy simulation mode of the Weather Research and Forecasting Model. Vapor-grown ice properties are associated with distinct observable signatures in polarimetric radar variables, with clear sensitivities to the simulated ice particle properties, including ice number, size, and distribution shape. In contrast, the liquid droplet number has little influence on both polarimetric and microphysical variables in the case presented herein. Polarimetric quantities are sensitive to the dominating crystal habit type in a volume, with enhancements for aspect ratios much lower or higher than unity. This synthesis of a microphysical model and a polarimetric forward simulator is a first step in the evaluation of detailed AHM microphysics.

54 ENVIRONMENTAL SCIENCES↗

Simulations of the pairwise kinematic Sunyaev-Zel'dovich signal

The pairwise kinematic Sunyaev–Zel’dovich (kSZ) signal from galaxy clusters is a probe of their line of sight momenta, and thus a potentially valuable source of cosmological information. In addition to the momenta, the amplitude of the measured signal depends on the properties of the intracluster gas and observational limitations such as errors in determining cluster centers and redshifts. In this work, we simulate the pairwise kSZ signal of clusters at z<1, using the output from a cosmological N-body simulation and including the properties of the intracluster gas via a model that can be varied in post-processing. We find that modifications to the gas profile due to star formation and feedback reduce the pairwise kSZ amplitude of clusters by ∼50%, relative to the naive “gas traces mass” assumption. We demonstrate that miscentering can reduce the overall amplitude of the pairwise kSZ signal by up to 10%, while redshift errors can lead to an almost complete suppression of the signal at small separations. We confirm that a high-significance detection is expected from the combination of data from current generation, high-resolution cosmic microwave background experiments, such as the South Pole Telescope, and cluster samples from optical photometric surveys, such as the Dark Energy Survey. Furthermore, we forecast that future experiments such as Advanced ACTPol in conjunction with data from the Dark Energy Spectroscopic Instrument will yield detection significances of at least 20σ, and up to 57σ in an optimistic scenario. Our simulated maps are publicly available at http://www.hep.anl.gov/cosmology/ksz.html.

79 ASTRONOMY AND ASTROPHYSICS↗

New Cryogenic Method for Combining Lunar Regolith Simulant and Frozen Volatiles to Generate Icy Lunar Simulated Regolith

There is a requirement within the lunar science and exploration communities to develop lunar simulants containing volatiles that are solids at the cryogenic temperatures found in the permanently shadowed regions (PSRs), such as those found near the Lunar South Pole. Icy regolith simulants would be used for curation training, as well as for various research activities. One of the most critical aspects of developing a regolith simulant that is more physically and chemically like icy lunar regolith is minimizing any form of modification driven by elevated temperatures. Here we document our ongoing efforts to combine materials at ultralow temperatures, designed to minimize any chemical reactions or other physical changes during the production of the icy regolith. Our goal is to document how one can create an essentially “unreacted” icy regolith that can serve as an effective “time zero” start point from which advanced curation research will proceed. This is done using commercial off-the-shelf equipment as much as possible, along with creating a custom spray plate that can be adapted to a wide variety of methods, all for the planetary simulant community. This method creates crystals of various volatiles such as water, methanol, along with CO 2 ice and these components are combined with lunar regolith that is at -196°C temperature to avoid chemical reactions, and/or phase changes thus creating a more chemically relevant icy lunar regolith. The rules of chemistry stay the same regardless of location whether it’s in the lab, or directly on the lunar surface, and therefore we aim for creating a more authentic icy lunar simulant in the laboratory by operating at ultralow temperatures. It is also envisioned that this method will lead to advanced materials testing in the future. In addition, this method is directly applicable to prior missions such as LCROSS.

Cryogenic Lunar Simulant and Icy Volatile Producti↗

New Cryogenic Method for Combining Lunar Regolith Simulant and Frozen Volatiles to Generate Icy Lunar Simulated Regolith

There is a requirement within the lunar science and exploration communities to develop lunar simulants containing volatiles that are solids at the cryogenic temperatures found in the permanently shadowed regions (PSRs), such as those found near the Lunar South Pole. Icy regolith simulants would be used for curation training, as well as for various research activities. One of the most critical aspects of developing a regolith simulant that is more physically and chemically like icy lunar regolith is minimizing any form of modification driven by elevated temperatures. Here we document our ongoing efforts to combine materials at ultralow temperatures, designed to minimize any chemical reactions or other physical changes during the production of the icy regolith. Our goal is to document how one can create an essentially “unreacted” icy regolith that can serve as an effective “time zero” start point from which advanced curation research will proceed. This is done using commercial off-the-shelf equipment as much as possible, along with creating a custom spray plate that can be adapted to a wide variety of methods, all for the planetary simulant community. This method creates crystals of various volatiles such as water, methanol, along with CO 2 ice and these components are combined with lunar regolith that is at -196°C temperature to avoid chemical reactions, and/or phase changes thus creating a more chemically relevant icy lunar regolith. The rules of chemistry stay the same regardless of location whether it’s in the lab, or directly on the lunar surface, and therefore we aim for creating a more authentic icy lunar simulant in the laboratory by operating at ultralow temperatures. It is also envisioned that this method will lead to advanced materials testing in the future. In addition, this method is directly applicable to prior missions such as LCROSS.

Cryogenic Lunar Simulant and Icy Volatile Producti↗

A numerical simulation of a dry-out process for CO 2 sequestration in heterogeneous deep saline aquifers: Original Research Article: Dry-out problem: Multifactorial numerical simulation

Salt precipitation can be generated near the injection well when dry supercritical CO 2 (scCO 2 ) is injected into deep saline aquifers. Traditional laboratory experiments and numerical simulations focus on the generation of salt precipitation and its impact on core permeability. Current laboratory experimental studies have shown that the heterogeneity of the core and water content on the surface of the particles plays an important role in the formation of salt precipitation. The aim of this study is to explore the effects of brine back-flow, heterogeneity, and low water content on salt precipitation. Numerical simulations were performed using core-scale and site-scale models. Simulation results showed the following: 1. Back-flow plays an important role in the salt accumulation near the well – an homogeneous laboratory experiment used with a small core may not provide a full picture of the back-flow phenomenon. 2. The continuous dry-out process reduces the water content near the injection well to zero. As a result, the full range of saturation flow equations should be used to simulate the dry-out process accurately. 3. The heterogeneity of the matrix not only affects the distribution of salt precipitation but also increases the amount of salt precipitation near the well. Homogeneous simulation underestimates the amount of salt precipitation. 4. The reduction in porosity due to salt precipitation could extend the dry-out process in the low permeability zone, further expanding the salt precipitation area.

58 GEOSCIENCES↗

Assessing the mechanisms governing the daytime evolution of marine stratocumulus using large‐eddy simulation

Large‐eddy simulation (LES) output for a case of thin stratocumulus off the coast of California is examined in a mixed‐layer analysis framework to identify the specific mechanisms responsible for governing the evolution of the cloud system. An equation for cloud‐base height tendency isolates the individual cloud‐modulating mechanisms that control the evolution of boundary‐layer liquid‐water static energy ( S l ) and total water mixing ratio ( q T ). With a suitable spin‐up procedure, the control simulation performs admirably compared with observed estimates of liquid water content, vertical velocity variance, and radiative fluxes sampled during an aircraft field campaign. Investigation of the cloud response to various environmental forcing scenarios was addressed through a suite of sensitivity simulations, including variations in subsidence velocity, surface fluxes, wind shear near the inversion, and radiative forcing. In the control simulation, rising cloud‐base tendencies are associated with entrainment warming/drying and short‐wave absorption, whereas lowering cloud‐base tendencies are driven by long‐wave cooling. Even in the presence of substantial afternoon solar heating, entrainment fluxes remained active. The thin cloud demonstrated unexpected resiliency, with mixed‐layer analysis indicating that, as the short‐wave flux decreases later in the afternoon, the relative contribution of long‐wave cooling often becomes large enough to offset entrainment warming/drying and result in a reversal of cloud‐base tendency. The evolution of cloud‐base tendency is found to be insensitive to the net radiative flux divergence for most of the simulations (liquid water path ranging from ~10–50 g/m 2 ). Error analysis in comparison with LES S l and q T budgets suggests that our method of entrainment flux calculation could be improved by a more complete understanding of entrainment‐layer physics.

McMichael, Lucas A.↗

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗

Simulation of the BEAVRS benchmark using VERA

The Virtual Environment for Reactor Applications (VERA) is being developed by the Consortium for the Advanced Simulation of Light Water Reactors (CASL). Here, VERA is used to perform the Benchmark for Evaluation and Validation of Reactor Simulations (BEAVRS), which provides two cycles’ worth of operating power history, along with a full, detailed description of the geometry, and measured data. Cycle 1 and 2 are simulated with VERA and the results are compared to the measured zero power physics tests, critical boron concentration, and flux maps.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying the unknown impact of segmentation uncertainty on image-based simulations

Abstract Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.

42 ENGINEERING↗

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗

Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity

In this report we compare time sampling and ensemble averaging as two different methods available for phase space sampling. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium molecular dynamics. We introduce two different schemes for the ensemble averaging approach, and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical molecular dynamics, the ensemble generation approaches may find their greatest utility in computationally expensive simulations such as first principles molecular dynamics. For such simulations, where each time step is costly, time sampling can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. On the other hand, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each ensemble is independent. For this reason, particularly when using massively parallel architectures, ensemble sampling can result in much shorter simulation times and exhibits similar overall computational effort.

74 ATOMIC AND MOLECULAR PHYSICS↗