A new theory of macular organization based on computer-assisted 3-D reconstruction, Monte Carlo simulation and symbolic modeling of vestibular maculas
No abstract available
Acoustic Maculae/anatomy & histology/physiology↗
SEARCH · Engineering Papers
Results for “Monte Carlo Modelling”
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.
Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification
CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.
36 MATERIALS SCIENCE↗
Dust Morphology of Comet Hale-Bopp (C/1995 01). II. Introduction of a Working Model
A Monte Carlo image simulation code for dust features in comets is applied to comet Hale-Bopp in order to model the object's persistent porcupine-like appearance on high-resolution images taken between May 11 and November 2, 1996. In this scenario, the dust-emission profile is dominated by several brief flare-ups, or puffs in the production of dust from one of the sources.
comet↗
An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices
In this paper, we report a random sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for under- standing the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbor hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbor hopping and 2D lattices. In addition, constraints on real-space con gurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.
74 ATOMIC AND MOLECULAR PHYSICS↗
The modulation of galactic protons by the solar wind - A Monte Carlo approach
Solar wind modulation effects on galactic protons and cosmic rays, using radial magnetic field mathematical model and Monte Carlo method
Giuli, R. T.↗
Final Technical Report for project entitled Highly Active, Durable, and Ultra-Low PGM NSTF Thin Film ORR Catalysts and Support
In this project, the objective was to develop new oxygen reduction reaction (ORR) electrocatalysts for proton exchange membrane fuel cells (PEMFCs) which could exceed all of the Department of Energy (DOE) 2020 targets listed in DE-FOA-0001224, Subtopic 1b, Table I. The expected outcome was development of one or more electrocatalysts which are substantially improved over the current state-of-the-art in terms of overall activity, durability, and cost, and are suitable for automotive traction and stationary fuel cell applications. This report summarizes this project’s progress towards meeting the stated objectives and includes a summary of key findings and conclusions. The development activities towards new highly active and durable thin film electrocatalysts has led to the discovery of several electrocatalysts which approached or exceeded several DOE 2020 targets for activity, durability, and performance in proton exchange membrane fuel cell (membrane electrode assemblies (MEAs). In this project, electrocatalyst development focused on systematic physical and electrochemical characterization of electrocatalyst activity, durability, and performance in MEAs as a function of electrocatalyst fabrication, compositional and structural variables. The project approach was to develop relationships between the catalyst functional responses (activity, durability and performance) with catalyst physical properties and catalyst fabrication methods. Additionally, the thin film catalysts were integrated onto the unique and durable 3M Nanostructured Thin Film (NSTF) support, consisting of arrays of self-assembled organic crystalline whiskers. The catalysts were evaluated against the project targets via extensive electrochemical characterization in MEA format at 3M, advanced structural and compositional microscopy at Oak Ridge National Laboratory (ORNL), and atomic structure analysis via XAFS at Argonne National Laboratory (ANL). The extensive electrochemical and physical characterization resulted in development of several trends which correlate the catalysts’ electrochemical properties to their physical properties, including composition, structure, and method of fabrication. The development was guided by density functional theory (DFT) modeling at Purdue University and kinetic Monte Carlo (kMC) modeling at Johns Hopkins University. The catalyst simulations provided key insights into the observed experimental activity and durability trends, and additionally were utilized to assess new electrocatalyst concepts prior to or coincident with physical catalyst development. This development has led to several electrocatalyst candidates with activity and durability which approach or exceed DOE targets. One class of catalysts, based on nm-scale thin layers of Pt on Ir, met or exceeded 6 of the 6 DOE targets the project addressed.
08 HYDROGEN↗
Node-avoiding Levy flight - A numerical test of the epsilon expansion
A study is conducted of an extension of Levy flight to include self-repulsion in the path of the walk. The extension is called node-avoiding Levy flight and its equivalence to the n approaches 0 limit of a statistical mechanical model for a magnetic system with long-range interactions between the spins is shown. By use of this equivalence it is possible to make a detailed comparison beween the results of the epsilon expansion for the magnetic model, a Monte Carlo simulation of the Levy flight model, and the results of a Flory-type argument. This is the first comparison of the epsilon expansion for epsilon much less than 1 with a numerical simulation for any model. Some speculations are made on applications of the model of node-avoiding Levy flight.
Halley, J. W.↗
Composition/Structure/Dynamics of comet and planetary satellite atmospheres
This research program addresses two cases of tenuous planetary atmospheres: comets and Io. The comet atmospheric research seeks to analyze a set of spatial profiles of CN in comet Halley taken in a 7.4-day period in April 1986; to apply a new dust coma model to various observations; and to analyze observations of the inner hydrogen coma, which can be optically thick to the resonance scattering of Lyman-alpha radiation, with the newly developed approach that combines a spherical radiative transfer model with our Monte Carlo H coma model. The Io research seeks to understand the atmospheric escape from Io with a hybrid-kinetic model for neutral gases and plasma given methods and algorithms developed for the study of neutral gas cometary atmospheres and the earth's polar wind and plasmasphere. Progress is reported on cometary Hydrogen Lyman-alpha studies; time-series analysis of cometary spatial profiles; model analysis of the dust comae of comets; and a global kinetic atmospheric model of Io.
Combi, Michael R.↗
Setup and verification of a SCALE/KENO platform for generic FHR benchmark calculations
The work presented in this article is preliminary to downstream analysis of a generic fluoride salt-cooled high-temperature reactor (gFHR) core performed by the University of Tennessee in collaboration with Kairos Power (KP). A Monte Carlo transport model of the publicly available gFHR equilibrium core is developed in SCALE/KENO with multigroup energy treatment. Several output quantities of interest are used to verify the simulations against a benchmark model developed by KP using the continuous energy Monte Carlo code Serpent 2. Good agreement is seen in flux and fission rate profiles with a maximum relative difference of 1.4% and 2.8% respectively. Furthermore, an effective multiplication factor bias of 44 pcm was observed between the two simulations. The fuel temperature reactivity coefficient calculated with SCALE is within uncertainty to the reference model. This verification acts as a publicly reproducible benchmark for the gFHR in SCALE/KENO. A simplified depletion model is also presented where a single fuel pebble is depleted to discharge burnup through the equilibrium core while the equilibrium core is assumed to be invariant. This method produces results that intercept the equilibrium core concentrations in every case, however, an interesting artifact of this particular depletion model is uncovered. The phenomenon is shown to be a fundamental feature of the differential rate equations and inspires questions about how this system behaves when the time evolution of the equilibrium core is considered. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Atomic Oxygen Energy in Low Frequency Hyperthermal Plasma Ashers
Experimental and analytical analysis of the atomic oxygen erosion of pyrolytic graphite as well as Monte Carlo computational modeling of the erosion of Kapton H (DuPont, Wilmington, DE) polyimide was performed to determine the hyperthermal energy of low frequency (30 to 35 kHz) plasma ashers operating on air. It was concluded that hyperthermal energies in the range of 0.3 to 0.9 eV are produced in the low frequency air plasmas which results in texturing similar to that in low Earth orbit (LEO). Monte Carlo computational modeling also indicated that such low energy directed ions are fully capable of producing the experimentally observed textured surfaces in low frequency plasmas.
Oxygen atoms↗
Magnetic multilayer edges in Bernal-stacked hexagonal boron nitride
Single-layer h-BN is known to have edges with unique magnetism; however, in the commonly fabricated multilayer AA'-h-BN, edge relaxations occur that create interlayer bonds and eliminate the unpaired electrons at the edge. Recently, a robust method of growing the unconventional Bernal-stacked h-BN (AB-h-BN) has been reported. Here, we use theoretical approaches to investigate the nitrogen-terminated zigzag edges in AB-h-BN that can be formed in a controlled fashion using a high-energy electron beam. We find that these "open"edges remain intact in bilayer and multilayer AB-h-BN, enabling researchers potentially to investigate these edge states experimentally. We also investigate the thermodynamics of the spin configurations at the edge by constructing a lattice model that is based on parameters extracted from a set of first-principles calculations. We find that the edge spins in neighboring layers interact very weakly, resulting in a sequence of independent spin chains in multilayer samples. By solving this model using Monte Carlo simulations, we can determine nm-scale correlation lengths at liquid-N 2 temperatures and lower. At low temperatures, these edges may be utilized in magnetoresistance and spintronics applications.
1-dimensional spin chains↗
BATSE software for the analysis of the gamma ray burst spatial distribution
The Burst and Transient Source Experiment (BATSE) on the Gamma Ray Observatory (GRO) is designed to study astronomical gamma ray sources and to provide better positional, spectral, and time resolution about these objects than has previously been possible from one experiment. The procedure to be used in the analysis of the gamma ray burst spatial distribution is presented. Data is input from BATSE via the Gamma Ray Burst Catalog (listing individual burst positions, flux values, and associated errors) and the Sky Sensitivity Map (which summarizes observational selection effects in table format). A FORTRAN program generates Monte Carlo burst catalogs, which are models to be compared to the actual distribution. The Monte Carlo models are then filtered through the Sky Sensitivity Map so that they suffer from the same selection effects as the actual catalog data. Additionally, each burst position is converted into a probability distribution to mimic BATSE positional sensitivity. The Burst Catalog, Monte Carlo burst catalog, and Sky Sensitivity Map are then passed onto an IDL program that compares the catalogs for statistical significance. The Sky Sensitivity Map is used to estimate how often each sky area is observed above the minimum flux level in question. Each burst found in this sky area is then weighted according to the frequency with which this sky area is observed. The catalogs are then compared via tests of homogeneity (based on their radial distributions) and isotropy (based upon their angular distributions). The results of the statistical comparisons along with graphs and charts of the summaries, are output from the IDL program for study.
Hakkila, Jon↗
Monte Carlo analysis of lobular gas-surface scattering in tubes applied to thermal transpiration
A model of rarefied gas flow in tubes was developed which combines a lobular distribution with diffuse reflection at the wall. The model with Monte Carlo techniques was used to explain previously observed deviations in the free molecular thermal transpiration ratio which suggest molecules can have a greater tube transmission probability in a hot-to-cold direction than in a cold-to-hot direction. The model yields correct magnitudes of transmission probability ratios for helium in Pyrex tubing (1.09 to 1.14), and some effects of wall-temperature distribution, tube surface roughness, tube dimensions, gas temperature, and gas molecular mass.
Smith, J. D.↗
Fragile Topology and Flat-Band Superconductivity in the Strong-Coupling Regime
In flat bands, superconductivity can lead to surprising transport effects. The superfluid “mobility”, in the form of the superfluid weight D s , does not draw from the curvature of the band but has a purely band-geometric origin. In a mean-field description, a nonzero Chern number or fragile topology sets a lower bound for D s , which, via the Berezinskii-Kosterlitz-Thouless mechanism, might explain the relatively high superconducting transition temperature measured in magic-angle twisted bilayer graphene (MATBG). For fragile topology, relevant for the bilayer system, the fate of this bound for finite temperature and beyond the mean-field approximation remained, however, unclear. Here, we numerically use exact Monte Carlo simulations to study an attractive Hubbard model in flat bands with topological properties akin to those of MATBG. We find a superconducting phase transition with a critical temperature that scales linearly with the interaction strength. Then, we investigate the robustness of the superconducting state to the addition of trivial bands that may or may not trivialize the fragile topology. Our results substantiate the validity of the topological bound beyond the mean-field regime and further stress the importance of fragile topology for flat-band superconductivity.
2-dimensional systems↗
Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal
Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Integrating chromosome conformation and DNA repair in a computational framework to assess cell radiosensitivity
Objective. The arrangement of chromosomes in the cell nucleus has implications for cell radiosensitivity. The development of new tools to utilize Hi-C chromosome conformation data in nanoscale radiation track structure simulations allows for in silico investigation of this phenomenon. We have developed a framework employing Hi-C-based cell nucleus models in Monte Carlo radiation simulations, in conjunction with mechanistic models of DNA repair, to predict not only the initial radiation-induced DNA damage, but also the repair outcomes resulting from this damage, allowing us to investigate the role chromosome conformation plays in the biological outcome of radiation exposure. Approach. In this study, we used this framework to generate cell nucleus models based on Hi-C data from fibroblast and lymphoblastoid cells and explore the effects of cell type-specific chromosome structure on radiation response. The models were used to simulate external beam irradiation including DNA damage and subsequent DNA repair. The kinetics of the simulated DNA repair were compared with previous results. Main results. We found that the fibroblast models resulted in a higher rate of inter-chromosome misrepair than the lymphoblastoid model, despite having similar amounts of initial DNA damage and total misrepairs for each irradiation scenario. Significance. This framework represents a step forward in radiobiological modeling and simulation allowing for more realistic investigation of radiosensitivity in different types of cells.
59 BASIC BIOLOGICAL SCIENCES↗
Coupled Monte Carlo and thermal-fluid modeling of high temperature gas reactors using Cardinal
Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computa-tional fluid dynamics to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. Here, we provide an introduction to Cardinal's software design, data mapping, and multi -physics coupling strategy to highlight our approach to overcoming common challenges in multiphysics simulation. We then describe an application of Cardinal to prismatic High Temperature Gas Reactors (HTGRs) with various combinations of NekRS, OpenMC, BISON, and THM. A high-resolution coupling of NekRS, OpenMC, and BISON provides a reference solution at the unit cell level and shows excellent agree-ment with a lower-resolution coupling of THM, OpenMC, and BISON. A full core coupling of THM, OpenMC, and BISON resolving the three-dimensional conjugate heat transfer and sub-pin power distri-bution then provides detailed predictions of HTGR temperatures and the fission distribution.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
E-Area Low-Level Waste Facility GoldSim System Model
This report documents the development of the E-Area Low-Level Waste Facility (ELLWF) trench system model. The GoldSim® Monte Carlo simulation software is utilized to model the release and transport of radiological inventory disposed (both currently and in the future) within Engineered and Slit Trenches. This model is in support of the sensitivity and uncertainty analysis for the ELLWF Performance Assessment. The ELLWF system model utilizes a hybrid-approach to accurately describe the disposal system. The Hydrologic Evaluation of Landfill Performance model provides the infiltration data to both PORFLOW and GoldSim. PORFLOW is used to calibrate the GoldSim model to ensure confidence in the stochastic results. Finally, the concentrations from GoldSim transport simulations are fed into the SRNL Dose Toolkit to calculate dose impacts and assess plume interaction
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗