Results for “SIMULATION METHODS” — Page 34

DOE OSTI2016

Efficient implicit LES method for the simulation of turbulent cavitating flows

We present a numerical method for efficient large-eddy simulation of compressible liquid flows with cavitation based on an implicit subgrid-scale model. Phase change and subgrid-scale interface structures are modeled by a homogeneous mixture model that assumes local thermodynamic equilibrium. Unlike previous approaches, emphasis is placed on operating on a small stencil (at most four cells). The truncation error of the discretization is designed to function as a physically consistent subgrid-scale model for turbulence. We formulate a sensor functional that detects shock waves or pseudo-phase boundaries within the homogeneous mixture model for localizing numerical dissipation. In smooth regions of the flow field, a formally non-dissipative central discretization scheme is used in combination with a regularization term to model the effect of unresolved subgrid scales. The new method is validated by computing standard single- and two-phase test-cases. Comparison of results for a turbulent cavitating mixing layer obtained with the new method demonstrates its suitability for the target applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Development of a Smart Signal Detection Method for Cyclic Voltammetry via Artificial Neural Intelligence

The electrorefiner (ER) is the heart of pyroprocessing technology which is a high-temperature method for separating uranium from Used Nuclear Fuel (UNF). It is important to improve this technology with respect to nuclear materials detection and accountability. Artificial Neural Intelligence (ANI) is a novel data analysis and simulation method that can be applied to electrochemical data sets. This computational code, which has been performed using the commercial software MATLAB, can be trained to generalize adequate electrical current and potential simulated data sets for the unseen data with a high accuracy of prediction. For this purpose, a massive collection of cyclic voltammetry (CV) data sets by Hoover (2014), for 0.5, 1, 2.5, and 5 wt% of zirconium chloride in LiCl-KCl molten salt with different scan rates at 773 K has been considered. The computer is trained via ANI to predict the unseen data after providing suitable hidden layers and validation numbers. In addition, this work can trace the CV plot for a blind condition by interpolating between two simulated data set. The different hidden layers with various neurons (from 5 to 30) at several validation numbers (from 5 to 30) has been studied and the average percent error between experimental and theoretical data for 0.5 wt% with 200 and 450 mV/s has been calculated. Preliminary results demonstrate that if the number of hidden layers increases from one to three, the average error falls down from 44% to around 8%. The best condition which gives a minimum average percent error has been discovered and the simulated CV graph has been compared with the experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2015

SU-E-J-163: A Biomechanical Lung Model for Respiratory Motion Study

Purpose: This work presents a biomechanical model to investigate the complex respiratory motion for the lung tumor tracking in radiosurgery by computer simulation. Methods: The models include networked massspring-dampers to describe the tumor motion, different types of surrogate signals, and the force generated by the diaphragm. Each mass-springdamper has the same mechanical structure and each model can have different numbers of mass-spring-dampers. Both linear and nonlinear stiffness parameters were considered, and the damping ratio was tuned in a range so that the tumor motion was over-damped (no natural tumor oscillation occurs without force from the diaphragm). The simulation was run by using ODE45 (ordinary differential equations by Runge-Kutta method) in MATLAB, and all time courses of motions and inputs (force) were generated and compared. Results: The curvature of the motion time courses around their peaks was sensitive to the damping ratio. Therefore, the damping ratio can be determined based on the clinical data of a high sampling rate. The peak values of different signals and the time the peaks occurred were compared, and it was found that the diaphragm force had a time lead over the tumor motion, and the lead time (0.1–0.4 seconds) depended on the distance between the tumor and the diaphragm. Conclusion: We reported a model based analysis approach for the spatial and temporal relation between the motion of the lung tumor and the surrogate signals. Due to the phase lead of the diaphragm in comparing with the lung tumor motion, the measurement of diaphragm motion (or its electromyography signal) can be used as a beam gating signal in radiosurgery, and it can also be an additional surrogate signal for better tumor motion tracking. The research is funded by the American Cancer Society (ACS) grant. The grant name is: Frameless SRS Based on Robotic Head Motion Cancellation. The grant number is: RSG-13-313-01-CCE.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Scalable, Physical Effects Measureable Microgridfor Cyber Resilience Analysis

SPEMMCRA provides a low-cost simulated method to provide for testing of control theory and cybersecurity approaches to microgrid implementations. Specifically, this physics-based approach provides realistic responses to compare control strategies in simulation, but also to demonstrate actions that can or cannot be taken in taking cybersecurity actions . As communications interfaces or other actions are taken that affect the flow of data, SPEMMCRA can show whether the microgrid operation is adversely affected by what is intended to mitigate a cybersecurity event.

Rieger, CraigG↗

NASA NTRS2001

Transcatheter Microwave Antenna

A method, simulation, and apparatus are provided that are highly suitable for treatment of benign prostatic hyperplasia (BPH). A catheter is disclosed that includes a small diameter disk loaded monopole antenna surrounded by fusion material having a high heat of fusion and a melting point preferably at or near body temperature. Microwaves from the antenna heat prostatic tissue to promote necrosing of the prostatic tissue that relieves the pressure of the prostatic tissue against the urethra as the body reabsorbs the necrosed or dead tissue. The fusion material keeps the urethra cool by means of the heat of fusion of the fusion material. This prevents damage to the urethra while the prostatic tissue is necrosed. A computer simulation is provided that can be used to predict the resulting temperature profile produced in the prostatic tissue. By changing the various control features of the catheter and method of applying microwave energy a temperature profile can be predicted and produced that is similar to the temperature profile desired for the particular patient.

Dickey G Arndt↗

NASA NTRS2003

Method for Selective Thermal Ablation

A method, simulation, and apparatus are provided that are highly suitable for treatment of benign prostatic hyperplasia (BPH). A catheter is disclosed that includes a small diameter disk loaded monopole antenna surrounded by fusion material having a high heat of fusion and a melting point preferably at or near body temperature. Microwaves from the antenna heat prostatic tissue to promote necrosing of the prostatic tissue that relieves the pressure of the prostatic tissue against the urethra as the body reabsorbs the necrosed or dead tissue. The fusion material keeps the urethra cool by means of the heat of fusion of the fusion material. This prevents damage to the urethra while the prostatic tissue is necrosed. A computer simulation is provided that can be used to predict the resulting temperature profile produced in the prostatic tissue. By changing the various control features of the catheter and method of applying microwave energy a temperature profile can be predicted and produced that is similar to the temperature profile desired for the particular patient.

Arndt, G. Dickey↗

NASA NTRS2003

Method of Constructing a Microwave Antenna

A method, simulation, and apparatus are provided that are highly suitable for treatment of benign prostatic hyperplasia (BPH). A catheter is disclosed that includes a small diameter disk loaded monopole antenna surrounded by fusion material having a high heat of fusion and a melting point preferably at or near body temperature. Microwaves from the antenna heat prostatic tissue to promote necrosing of the prostatic tissue that relieves the pressure of the prostatic tissue against the urethra as the body reabsorbs the necrosed or dead tissue. The fusion material keeps the urethra cool by means of the heat of fusion of the fusion material. This prevents damage to the urethra while the prostatic tissue is necrosed. A computer simulation is provided that can be used to predict the resulting temperature profile produced in the prostatic tissue. By changing the various control features of the catheter and method of applying microwave energy a temperature profile can be predicted and produced that is similar to the temperature profile desired for the particular patient.

Arndt, G. Dickey↗

NASA NTRS2003

Method for selective thermal ablation

A method, simulation, and apparatus are provided that are highly suitable for treatment of benign prostatic hyperplasia (BPH). A catheter is disclosed that includes a small diameter disk loaded monopole antenna surrounded by fusion material having a high heat of fusion and a melting point preferably at or near body temperature. Microwaves from the antenna heat prostatic tissue to promote necrosing of the prostatic tissue that relieves the pressure of the prostatic tissue against the urethra as the body reabsorbs the necrosed or dead tissue. The fusion material keeps the urethra cool by means of the heat of fusion of the fusion material. This prevents damage to the urethra while the prostatic tissue is necrosed. A computer simulation is provided that can be used to predict the resulting temperature profile produced in the prostatic tissue. By changing the various control features of the catheter and method of applying microwave energy a temperature profile can be predicted and produced that is similar to the temperature profile desired for the particular patient.

Arndt, G. Dickey↗

NASA NTRS2003

Method of constructing a microwave antenna

A method, simulation, and apparatus are provided that are highly suitable for treatment of benign prostatic hyperplasia (BPH). A catheter is disclosed that includes a small diameter disk loaded monopole antenna surrounded by fusion material having a high heat of fusion and a melting point preferably at or near body temperature. Microwaves from the antenna heat prostatic tissue to promote necrosing of the prostatic tissue that relieves the pressure of the prostatic tissue against the urethra as the body reabsorbs the necrosed or dead tissue. The fusion material keeps the urethra cool by means of the heat of fusion of the fusion material. This prevents damage to the urethra while the prostatic tissue is necrosed. A computer simulation is provided that can be used to predict the resulting temperature profile produced in the prostatic tissue. By changing the various control features of the catheter and method of applying microwave energy a temperature profile can be predicted and produced that is similar to the temperature profile desired for the particular patient.

Arndt, G. Dickey↗

NASA NTRS2012

Characterization and Analysis of Triaxially Braided Polymer Composites under Static and Impact Loads

In order to design impact resistant aerospace components made of triaxially-braided polymer matrix composite materials, a need exists to have reliable impact simulation methods and a detailed understanding of the material behavior. Traditional test methods and specimen designs have yielded unrealistic material property data due to material mechanisms such as edge damage. To overcome these deficiencies, various alternative testing geometries such as notched flat coupons have been examined to alleviate difficulties observed with standard test methods. The results from the coupon level tests have been used to characterize and validate a macro level finite element-based model which can be used to simulate the mechanical and impact response of the braided composites. In the analytical model, the triaxial braid unit cell is approximated by using four parallel laminated composites, each with a different fiber layup, which roughly simulates the braid architecture. In the analysis, each of these laminated composites is modeled as a shell element. Currently, each shell element is considered to be a smeared homogeneous material. Simplified micromechanics techniques and lamination theory are used to determine the equivalent stiffness properties of each shell element, and results from the coupon level tests on the braided composite are used to back out the strength properties of each shell element. Recent improvements to the model include the incorporation of strain rate effects into the model. Simulations of ballistic impact tests have been carried out to investigate and verify the analysis approach.

Goldberg, Robert K.↗

DOE OSTI2014

SU-E-T-455: Characterization of 3D Printed Materials for Proton Beam Therapy

Purpose: The widespread availability of low cost 3D printing technologies provides an alternative fabrication method for customized proton range modifying accessories such as compensators and boluses. However the material properties of the printed object are dependent on the printing technology used. In order to facilitate the application of 3D printing in proton therapy, this study investigated the stopping power of several printed materials using both proton pencil beam measurements and Monte Carlo simulations. Methods: Five 3–4 cm cubes fabricated using three 3D printing technologies (selective laser sintering, fused-deposition modeling and stereolithography) from five printers were investigated. The cubes were scanned on a CT scanner and the depth dose curves for a mono-energetic pencil beam passing through the material were measured using a large parallel plate ion chamber in a water tank. Each cube was measured from two directions (perpendicular and parallel to printing plane) to evaluate the effects of the anisotropic material layout. The results were compared with GEANT4 Monte Carlo simulation using the manufacturer specified material density and chemical composition data. Results: Compared with water, the differences from the range pull back by the printed blocks varied and corresponded well with the material CT Hounsfield unit. The measurement results were in agreement with Monte Carlo simulation. However, depending on the technology, inhomogeneity existed in the printed cubes evidenced from CT images. The effect of such inhomogeneity on the proton beam is to be investigated. Conclusion: Printed blocks by three different 3D printing technologies were characterized for proton beam with measurements and Monte Carlo simulation. The effects of the printing technologies in proton range and stopping power were studied. The derived results can be applied when specific devices are used in proton radiotherapy.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2016

WE-FG-207B-06: Plaque Composition Measurement with Dual Energy Computed Tomography

Purpose: To investigate the feasibility of characterizing arterial plaque composition in terms of water, lipid and protein or calcium using dual energy computed tomography. Characterization of plaque composition can potentially help distinguish vulnerable from stable plaques. Methods: Simulations studies were performed by the CT simulator based on ASTRA tomography toolbox. The beam energy for dual energy images was selected to be 80 kVp and 135 kVp. The radiation dose and energy spectrum for the CT simulator were carefully calibrated with respect to a 320-slice CT scanner. A digital chest phantom was constructed using Matlab for calibration and plaque measurement. Pure water, lipid, protein or calcium was used for calibration and a mixture of different volume percentages of these materials were used for validation purposes. Non-calcified plaque was simulated using water, lipid and protein with volumetric percentage range of 35%∼65%, 5%∼60% and 5%∼40%, respectively. Calcified plaque was simulated using water, lipid and calcium with volumetric percentage range of 50%∼80%, 8%∼45% and 3%∼13%, respectively. We employed iterative sinogram processing (ISP) to reduce the beam hardening effect in the simulation to improve the decomposition results. Results: The simulated known composition and dual energy decomposition results were in good agreement. Water, lipid and protein (calcium) mixtures were decomposed into water, lipid and protein (calcium) contents. The RMS errors of volumetric percentage for the water, lipid and protein (non-calcified plaque) decomposition, as compared to known values, were estimated to be approximately 5.74%, 2.54%, and 0.95% respectively. The RMS errors of volumetric percentage for the water, lipid and Calcium (calcified plaque) decomposition, as compared to known values, were estimated to be approximately 7.4%, 8.64%, and 0.08% respectively. Conclusion: The results of this study suggest that the dual energy decomposition can potentially be used to quantify the water, lipid, and protein or calcium composition of a plaque with relatively good accuracy. Grant funding from Toshiba Medical Systems and Philips Medical Systems.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Generalized Metropolis acceptance criterion for hybrid non-equilibrium molecular dynamics—Monte Carlo simulations

A family of hybrid simulation methods that combines the advantages of Monte Carlo (MC) with the strengths of classical molecular dynamics (MD) consists in carrying out short non-equilibrium MD (neMD) trajectories to generate new configurations that are subsequently accepted or rejected via an MC process. In the simplest case where a deterministic dynamic propagator is used to generate the neMD trajectories, the familiar Metropolis acceptance criterion based on the change in the total energy ΔE, min[1, exp( − βΔE)], guarantees that the hybrid algorithm will yield the equilibrium Boltzmann distribution. However, the functional form of the acceptance probability is more complex when the non-equilibrium switching process is generated via a non-deterministic stochastic dissipative propagator coupled to a heat bath. Here, we clarify the conditions under which the Metropolis criterion remains valid to rigorously yield a proper equilibrium Boltzmann distribution within hybrid neMD-MC algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Experimental and Monte Carlo evaluation of Eclipse treatment planning system for effects on dose distribution of the hip prostheses

The present study aimed to investigate the effects of titanium, titanium alloy, and stainless steel hip prostheses on dose distribution based on the Monte Carlo simulation method, as well as the accuracy of the Eclipse treatment planning system (TPS) at 6 and 18 MV photon energies. In the present study the pencil beam convolution (PBC) method implemented in the Eclipse TPS was compared to the Monte Carlo method and ionization chamber measurements. The present findings show that if high-Z material is used in prosthesis, large dose changes can occur due to scattering. The variance in dose observed in the present study was dependent on material type, density, and atomic number, as well as photon energy; as photon energy increased back scattering decreased. The dose perturbation effect of hip prostheses was significant and could not be predicted accurately by the PBC method for hip prostheses. The findings show that for accurate dose calculation the Monte Carlo-based TPS should be used in patients with hip prostheses.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2015

SU-E-T-667: Radiosensitization Due to Gold Nanoparticles: A Monte Carlo Cellular Dosimetry Investigation of An Expansive Parameter Space

Purpose: To investigate dose enhancement to cellular compartments following gold nanoparticle (GNP) uptake in tissue, varying cell and tissue morphology, intra and extracellular GNP distribution, and source energy using Monte Carlo (MC) simulations. Methods: Models of single and multiple cells are developed for normal and cancerous tissues; cells (outer radii 5–10 µm) are modeled as concentric spheres comprising the nucleus (radii 2.5–7.5 µm) and cytoplasm. GNP distributions modeled include homogeneous distributions throughout the cytoplasm, variable numbers of GNP-containing endosomes within the cytoplasm, or distributed in a spherical shell about the nucleus. Gold concentrations range from 1 to 30 mg/g. Dose to nucleus and to cytoplasm for simulations including GNPs are compared to simulations without GNPs to compute Nuclear and Cytoplasm Dose Enhancement Factors (NDEF, CDEF). Photon source energies are between 20 keV and 1.25 MeV. Results: DEFs are highly sensitive to GNP intracellular distribution; for a 2.5 µm radius nucleus irradiated by a 30 keV source, NDEF varies from 1.2 for a single endosome containing all GNPs to 8.2 for GNPs distributed about the nucleus (7 mg/g). DEFs vary with cell dimensions and source energy: NDEFs vary from 2.5 (90 keV) to 8.2 (30 keV) for a 2.5 µm radius nucleus and from 1.1 (90 keV) to 1.3 (30 keV) for a 7.5 µm radius nucleus, both with GNPs in a spherical shell about the nucleus (7 mg/g). NDEF and CDEF are generally different within a single cell. For multicell models, the presence of gold within intervening tissues between source and target perturbs the fluence reaching cellular targets, resulting in DEF inhomogeneities within a population of irradiated cells. Conclusion: DEFs vary by an order of magnitude for different cell models, GNP distributions, and source energies, demonstrating the importance of detailed modelling for advancing GNP development for radiotherapy. Funding provided by the Natural Sciences and Engineering Council of Canada (NSERC), and the Canada Research Chairs Program (CRC)

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2026

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗

DOE OSTI2018

Dynamic Potential of Interaction between Nitrogen Atoms and Aluminum Crystal Surface

Specific features of the angular distributions of accelerated neutral nitrogen atoms at the grazing angles of incidence on the Al(001) crystal surface have been investigated by the computer simulation method. The N–Al pair interaction potential is approximated by the three-parameter Morse potential with the energy dependent coefficients. The angular distributions of scattered atoms have been simulated taking into account the interaction between atoms and several atomic layers in the lattice and the atomic displacement during thermal oscillations. The parameters of the pair potential of accelerated neutral nitrogen atoms in the energy range from 10 to 70 keV have been determined according to the best agreement between the calculated dependence of the rainbow scattering angle on the energy of particles incident on the crystal surface and the available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2013

Uncertainty Quantification of GEOS-5 L-band Radiative Transfer Model Parameters Using Bayesian Inference and SMOS Observations

Uncertainties in L-band (1.4 GHz) radiative transfer modeling (RTM) affect the simulation of brightness temperatures (Tb) over land and the inversion of satellite-observed Tb into soil moisture retrievals. In particular, accurate estimates of the microwave soil roughness, vegetation opacity and scattering albedo for large-scale applications are difficult to obtain from field studies and often lack an uncertainty estimate. Here, a Markov Chain Monte Carlo (MCMC) simulation method is used to determine satellite-scale estimates of RTM parameters and their posterior uncertainty by minimizing the misfit between long-term averages and standard deviations of simulated and observed Tb at a range of incidence angles, at horizontal and vertical polarization, and for morning and evening overpasses. Tb simulations are generated with the Goddard Earth Observing System (GEOS-5) and confronted with Tb observations from the Soil Moisture Ocean Salinity (SMOS) mission. The MCMC algorithm suggests that the relative uncertainty of the RTM parameter estimates is typically less than 25 of the maximum a posteriori density (MAP) parameter value. Furthermore, the actual root-mean-square-differences in long-term Tb averages and standard deviations are found consistent with the respective estimated total simulation and observation error standard deviations of m3.1K and s2.4K. It is also shown that the MAP parameter values estimated through MCMC simulation are in close agreement with those obtained with Particle Swarm Optimization (PSO).

MCMC↗