Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference model”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Verification of a Numerical Model of the Offshore Wind Turbine from the Alpha Ventus Wind Farm Within OC5 Phase III: Preprint

The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project, is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project analyzes the Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. This paper shows results of the verification of the OWT models (code-to-code comparison). A subsequent publication will focus on their validation (comparison of simulated results to measured physical system response data). Based on the available data, the participants of Phase III set up numerical models of the OWT in their simulation tools. It was necessary to verify and to tune these models. The verification and tuning were performed against an OWT model available at the University of Stuttgart - Stuttgart Wind Energy (SWE) and documentation provided by Senvion and OWEC Tower. A very good match was achieved between the results from the reference SWE model and models set up by OC5 Phase III participants.

17 WIND ENERGY↗
DOE OSTI2022

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Utah FORGE: 2024 Discrete Fracture Network Model Data

The Utah FORGE 2024 Discrete Fracture Network (DFN) Model dataset provides a set of files representing discrete fracture network modeling for the FORGE site near Milford, Utah. The dataset includes four distinct DFN model file sets, each corresponding to different time frames and modeling approaches in 2024. These models characterize both natural and induced fractures in the geothermal reservoir, which consists of crystalline granitic and metamorphic rock approximately 8,000 feet below the ground surface. The dataset includes a reference DFN model from February 2024 that incorporates planar fractures and well trajectories, as well as upscaled permeability, porosity, compressibility, and storage values on specified grids. Additionally, there are models based on new microseismic (MEQ) data from May and July 2024, including fracture planes fitted to the latest MEQ catalog datasets, tensile fractures from hydraulic stimulation, and an alternative connected DFN for modeling purposes. Coordinate data is provided in both global and local frames, with detailed instructions on the transformations used to align with principal stress orientations. The dataset also includes notes and calculation files for estimating fracture sizes and differences between various fracture sets. There are subfolders for Global Coordinates and Local Coordinates. To move from the global to the local coordinate frame, fractures and wells were a) rotated 20 degrees counterclockwise looking down about the global point (335376.400482041, 4263189.99998761, 250.093546450195) to better align with the principal stresses; and b) translated by (-335408.68, -4263010.9, 1150). Upscaled permeability values using the _XYZ suffix show directions with respect to the global XYZ coordinate frame, while those using the _IJK suffix are aligned with local coordinate frame.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

SU-F-T-144: Analytical Closed Form Approximation for Carbon Ion Bragg Curves in Water

Purpose: Semi-empirical modeling is a powerful computational method in radiation dosimetry. A set of approximations exist for proton ion depth dose distribution (DDD) in water. However, the modeling is more complicated for carbon ions due to fragmentation. This study addresses this by providing and evaluating a new methodology for DDD modeling of carbon ions in water. Methods: The FLUKA, Monte Carlo (MC) general-purpose transport code was used for simulation of carbon DDDs for energies of 100–400 MeV in water as reference data model benchmarking. Based on Thomas Bortfeld’s closed form equation approximating proton Bragg Curves as a basis, we derived the critical constants for a beam of Carbon ions by applying models of radiation transport by Lee et. al. and Geiger to our simulated Carbon curves. We hypothesized that including a new exponential (κ) residual distance parameter to Bortfeld’s fluence reduction relation would improve DDD modeling for carbon ions. We are introducing an additional term to be added to Bortfeld’s equation to describe fragmentation tail. This term accounts for the pre-peak dose from nuclear fragments (NF). In the post peak region, the NF transport will be treated as new beams utilizing the Glauber model for interaction cross sections and the Abrasion- Ablation fragmentation model. Results: The carbon beam specific constants in the developed model were determined to be : p= 1.75, β=0.008 cm-1, γ=0.6, α=0.0007 cm MeV, σmono=0.08, and the new exponential parameter κ=0.55. This produced a close match for the plateau part of the curve (max deviation 6.37%). Conclusion: The derived semi-empirical model provides an accurate approximation of the MC simulated clinical carbon DDDs. This is the first direct semi-empirical simulation for the dosimetry of therapeutic carbon ions. The accurate modeling of the NF tail in the carbon DDD will provide key insight into distal edge dose deposition formation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Confirming a pyrolitic lower mantle using self-consistent pressure scales and new constraints on CaSiO 3 perovskite

In this study, we have examined the lower mantle composition and mineralogy by modeling the density ( ρ ), bulk sound velocity ( V Φ ), and dln ρ /dln V Φ profiles of candidate lower mantle minerals using literature and new experimental equation of state (EoS) results. Additionally, for CaSiO 3 perovskite, complimentary synchrotron X-ray diffraction measurements in a laser-heated diamond anvil cell were conducted up to 156 GPa between 1200 K and 2600 K to provide more reliable constraints on the thermal EoS parameters. These new experimental results as well as literature P-V-T data sets are systematically analyzed using an internally self-consistent pressure scale. We have modeled ρ , V Φ , and dln ρ /dln V Φ profiles of the lower mantle with representative pyrolitic and chondritic compositional models in which the effect of Fe spin transition in ferropericlase is also taken into account. Our modeling results show that a pyrolitic lower mantle with an aggregate mineralogy of 75 vol % bridgmanite, 17 vol % ferropericlase, and 8 vol % CaSiO 3 perovskite produces ρ and V Φ profiles in better agreement with preliminary reference Earth model than a lower mantle with a chondritic composition. The modeled ρ , V Φ , and dln ρ /dln V Φ are mainly affected by the relative ratio of bridgmanite and ferropericlase but are not sensitive to the variation of the CaSiO 3 perovskite content. In addition, the spin crossover of Fe in ferropericlase can greatly raise the value of dln ρ /dln V Φ in the middle lower mantle, which is useful in detecting the presence of ferropericlase in the region. Based on these new mineral physical constraints and radial seismic structure, our study suggests the lower mantle is pyrolitic, which is chemically indistinguishable from the upper mantle.

58 GEOSCIENCES↗
DOE OSTI2025

Sound Velocity of Hematite up to 13 GPa

We performed acoustic measurements on polycrystalline hematite using ultrasonic interferometric techniques in a Kawai-type multi-anvil apparatus. The elastic moduli and their pressure derivatives were determined by fitting the experimental data to third-order finite strain equations, yielding: 𝐾 𝑆,0 = 214.8 (20) GPa, 𝐺 0 = 92.0 (2) GPa, 𝐾$^{′}_{S,0}$ = 5.70 (1), and 𝐺$^{′}_{0}$ = 0.88 (4). Modeling of the velocity profiles of hematite along hot and cold slab geotherms shows that compressional (V P ) and shear (V S ) wave velocities are 22%–24% and 53%–55% lower, respectively, compared to the Preliminary Reference Earth Model at the top of the lower mantle. Furthermore, these findings likely offer a potential explanation for the low-velocity anomalies associated with slab graveyards at the base of the lower mantle.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Comparison of Serpent and HELIOS-2 as applied for the PWR few-group cross section generation

This paper discusses recent modifications to the Serpent Monte Carlo code methodology and related to the calculation of few-group diffusion coefficients and reflector discontinuity factors The new methods were assessed in the following manner. First, few-group homogenized cross sections calculated by Serpent for a reference PWR core were compared with those generated 1 commercial deterministic lattice transport code HELIOS-2. Second, Serpent and HELIOS-2 fe group cross section sets were later employed by nodal diffusion code DYN3D for the modeling the reference PWR core. Finally, the nodal diffusion results obtained using the both cross section sets were compared with the full core Serpent Monte Carlo solution. The test calculations show that Serpent can calculate the parameters required for nodal analyses similar to conventional deterministic lattice codes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

UNIFI's Grid-Forming (GFM) Inverter Reference Design: A Tutorial on Modeling, Control, and Experimental Implementation of GFM Inverters

The UNIFI Consortium's tutorial on grid-forming (GFM) inverters provides a comprehensive guide to the modeling, control, and experimental implementation of GFM inverters. As the integration of renewable energy accelerates, the transition from traditional grid following (GFL) to GFM inverters is crucial to ensure stable and sustainable power systems. This document outlines a reference design for three-phase and single-phase GFM inverters developed at the University of Texas at Austin. The tutorial also provides step-by-step guidance for accessing and using UNIFI’s GitHub repository, enabling users to design, build, and test GFM inverters efficiently. By fostering collaboration and equipping users with accessible resources, this initiative aims to drive widespread adoption of GFM technology across academia, utilities, and industries.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2012

Use of event-level neutrino telescope data in global fits for theories of new physics

We present a fast likelihood method for including event-level neutrino telescope data in parameter explorations of theories for new physics, and announce its public release as part of DarkSUSY 5.0.6. Our construction includes both angular and spectral information about neutrino events, as well as their total number. We also present a corresponding measure for simple model exclusion, which can be used for single models without reference to the rest of a parameter space. We perform a number of supersymmetric parameter scans with IceCube data to illustrate the utility of the method: example global fits and a signal recovery in the constrained minimal supersymmetric standard model (CMSSM), and a model exclusion exercise in a 7-parameter phenomenological version of the MSSM. The final IceCube detector configuration will probe almost the entire focus-point region of the CMSSM, as well as a number of MSSM-7 models that will not otherwise be accessible to e.g. direct detection. Our method accurately recovers the mock signal, and provides tight constraints on model parameters and derived quantities. We show that the inclusion of spectral information significantly improves the accuracy of the recovery, providing motivation for its use in future IceCube analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Development of the Technical Basis of a Unified Viscoplastic Model of 316H Stainless Steel for Incorporation into ASME Division 5

This report describes the development of a unified viscoplastic model describing the deformation of 316 stainless steel. The model was developed for incorporation into Section III, Division 5 of the ASME Boiler and Pressure Vessel Code covering the design and fabrication of Class A high temperature nuclear reactor components. The current version of the Code does not provide reference inelastic models, leaving the responsibility of defining an adequate model up to either the plant owner or designer. This approach represents a substantial barrier to using the ASME design by inelastic analysis methods. As historical practice demonstrates these methods produce much more efficient designs, compared to the design by elastic analysis methods, eliminating this barrier could result in safer, more effective, and more efficient structural designs thereby lowering the cost of future advanced reactors. This report summarizes the process of gathering experimental data, calibrating the candidate inelastic model against that data, and validating the model. The report also summarizes work on a new method of calculating the sensitivities of history dependent nonlinear constitutive models that was used to calibrate the 316H model using gradient-based optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2011

B Mixing in the Standard Model and Beyond: Lattice QCD

We give a brief overview and progress report on our lattice QCD calculation of neutral B mixing hadronic matrix elements needed for Standard Model and Beyond the Standard Model physics. Reference [1] contains more details and results. We have a good start on a large-data-set lattice QCD calculation of the matrix elements that describe neutral-B mixing. Our calculation will cover the operators needed for both Standard Model and Beyond the Standard Model physics. For the Standard Model matrix elements, we expect to halve the error on current, published calculations. In the Beyond the Standard Model case, this will be the first full-QCD (unquenched) calculation and the first update in ten years.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Universal topological quantum computation with strongly correlated Majorana edge modes

Abstract Majorana-based quantum gates are not complete for performing universal topological quantum computation while Fibonacci-based gates are difficult to be realized electronically and hardly coincide with the conventional quantum circuit models. In reference Hu and Kane (2018 Phys. Rev. Lett. 120 066801), it has been shown that a strongly correlated Majorana edge mode in a chiral topological superconductor can be decomposed into a Fibonacci anyon τ and a thermal operator anyon ɛ in the tricritical Ising model. The deconfinement of τ and ɛ via the interaction between the fermion modes yields the anyon collisions and gives the braiding of either τ or ɛ . With these braidings, the complete members of a set of universal gates, the Pauli gates, the Hadamard gate and extra phase gates for one-qubit as well as controlled-NOT (CNOT) gate for two-qubits, are topologically assembled. Encoding quantum information and reading out the computation results can be carried out through electric signals. With the sparse-dense mixed encodings, we set up the quantum circuit where the CNOT gate turns out to be a probabilistic gate and design the corresponding devices with thin films of the chiral topological superconductor. As an example of the universal topological quantum computing, we show the application to Shor’s integer factorization algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Sheath-induced distortions in particle distributions near enhanced polar outflow probe particle sensors

We discuss sheath and kinetic effects on ion and electron distribution functions at the aperture of enhanced Polar Outflow Probe particle sensors. For this purpose, the interaction between the CASSIOPE spacecraft and space environment is simulated fully kinetically using the electrostatic Particle In Cell code PTetra. The simulations account for the geometry of the main features of the spacecraft body, the booms, and the sensors. In addition to the background plasma, the model also accounts for Earth magnetic field. The plasma parameters assumed in the simulations are obtained from the latest version of the International Reference Ionosphere (IRI) model and the value of magnetic field is obtained from the International Geophysical Reference Field model. Our analysis shows significant distortions in the ion distribution function in the plane of the sensor aperture, as well as in the direction along the boom holding the sensor. We argue that significant distortions and asymmetries should also occur at the aperture of the suprathermal electron imager when suprathermal electrons are detected, with energies of 5 eV or more.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Application of PDF mixing models to premixed flames with differential diffusion

Differential diffusion alters the balance of reaction and diffusion in turbulent premixed combustion, affecting the performance and emissions of combustion devices. Modelling combustion devices with Probability or Filtered Density Function (PDF or FDF) methods provides an exact treatment for the change in composition due to chemical reaction, while molecular mixing has to be modelled. Previous PDF molecular mixing models do not account for differential diffusion in a manner which satisfies realizability requirements. A new approach for treating differential diffusion, which ensures realizability, is proposed for pairwise-exchange mixing models in general, and applied in the Interaction by Exchange with the Mean (IEM) model of Dopazo [26], and in the Euclidean Minimum Spanning Tree (EMST) model of Subramaniam and Pope [5]. The new differential diffusion models are referred to as IEM-DD and EMST-DD respectively. Results from two and three-dimensional DNS of turbulent premixed methane–air combustion show that mixing rates and conditional statistics of species mass fractions depend on species diffusivities and the combustion regime. Zero-dimensional PDF model results obtained for the two-dimensional DNS case show that the EMST-DD model best reproduces the features that characterize differential diffusion in the DNS. The essential feature of the EMST-DD model, which accounts for its success in turbulent premixed combustion, is that differential mixing rates are imposed within a model which mixes locally in composition space.

Energy & Fuels↗
DOE OSTI2025

Transport and confinement physics: Chapter 2 of the special issue: on the path to tokamak burning plasma operation

Progress in physics understanding and theoretical model development of plasma transport and confinement (TC) in the ITPA TC Topical Group since the publication of the ITER Physics Basis (IPB) document (Doyle et al 2007 Nucl. Fusion 47 S18) was summarized focusing on the contributions to ITER and burning plasma prediction and control. This paper provides a general and streamlined overview on the advances that were mainly led by the ITPA TC joint experiments and joint activities for the last 15 years (see JEX/JA table in appendix). This paper starts with the scientific strategy and scope of the ITPA TC Topical group and overall picture of the major progress, followed by the progress of each research field: particle transport, impurity transport, ion and electron thermal turbulent transport, momentum transport, impact of 3D magnetic fields on transport, confinement mode transitions, global confinement, and reduced transport modeling. Cross references with other Topical Groups are given in order to highlight overlapped topics, such as the 3D effect on the plasma transport in the edge and L-H transition physics. The increasing overlap between the topical groups is a reflection of the progress on integrating the known physics into comprehensive models that are better and better able to reproduce the plasma transport. In recent years, such integration has become increasingly prevalent when considering transport from the SOL, through the edge pedestal, and into the plasma core. In the near future, increased collaboration also with the magneto-hydrodynamic and energetic particles community will be important as we approach burning plasma conditions in next-step fusion devices. A summary of remaining challenges and next steps for each research field is given in the Summary section.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Collaborative Project. 3D Radiative Transfer Parameterization Over Mountains/Snow for High-Resolution Climate Models. Fast physics and Applications

Under the support of the aforementioned DOE Grant, we have made two fundamental contributions to atmospheric and climate sciences: (1) Develop an efficient 3-D radiative transfer parameterization for application to intense and intricate inhomogeneous mountain/snow regions. (2) Innovate a stochastic parameterization for light absorption by internally mixed black carbon and dust particles in snow grains for understanding and physical insight into snow albedo reduction in climate models. With reference to item (1), we divided solar fluxes reaching mountain surfaces into five components: direct and diffuse fluxes, direct- and diffuse-reflected fluxes, and coupled mountain-mountain flux. “Exact” 3D Monte Carlo photon tracing computations can then be performed for these solar flux components to compare with those calculated from the conventional plane-parallel (PP) radiative transfer program readily available in climate models. Subsequently, Parameterizations of the deviations of 3D from PP results for five flux components are carried out by means of the multiple linear regression analysis associated with topographic information, including elevation, solar incident angle, sky view factor, and terrain configuration factor. We derived five regression equations with high statistical correlations for flux deviations and successfully incorporated this efficient parameterization into WRF model, which was used as the testbed in connection with the Fu-Liou-Gu PP radiation scheme that has been included in the WRF physics package. Incorporating this 3D parameterization program, we conducted simulations of WRF and CCSM4 to understand and evaluate the mountain/snow effect on snow albedo reduction during seasonal transition and the interannual variability for snowmelt, cloud cover, and precipitation over the Western United States presented in the final report. With reference to item (2), we developed in our previous research a geometric-optics surface-wave approach (GOS) for the computation of light absorption and scattering by complex and inhomogeneous particles for application to aggregates and snow grains with external and internal mixing structures. We demonstrated that a small black (BC) particle on the order of 1 μm internally mixed with snow grains could effectively reduce visible snow albedo by as much as 5–10%. Following this work and within the context of DOE support, we have made two key accomplishments presented in the attached final report.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Representation of the Community Earth System Model (CESM1) CAM4-chem within the Chemistry-Climate Model Initiative (CCMI)

The Community Earth System Model (CESM1) CAM4-chem has been used to perform the Chemistry Climate Model Initiative (CCMI) reference and sensitivity simulations. In this model, the Community Atmospheric Model version 4 (CAM4) is fully coupled to tropospheric and stratospheric chemistry. Details and specifics of each configuration, including new developments and improvements are described. CESM1 CAM4-chem is a low-top model that reaches up to approximately 40 km and uses a horizontal resolution of 1.9° latitude and 2.5° longitude. For the specified dynamics experiments, the model is nudged to Modern-Era Retrospective Analysis for Research and Applications (MERRA) reanalysis. We summarize the performance of the three reference simulations suggested by CCMI, with a focus on the last 15 years of the simulation when most observations are available. Comparisons with selected data sets are employed to demonstrate the general performance of the model. We highlight new data sets that are suited for multi-model evaluation studies. Most important improvements of the model are the treatment of stratospheric aerosols and the corresponding adjustments for radiation and optics, the updated chemistry scheme including improved polar chemistry and stratospheric dynamics and improved dry deposition rates. These updates lead to a very good representation of tropospheric ozone within 20 % of values from available observations for most regions. In particular, the trend and magnitude of surface ozone is much improved compared to earlier versions of the model. Furthermore, stratospheric column ozone of the Southern Hemisphere in winter and spring is reasonably well represented. In conclusion, all experiments still underestimate CO most significantly in Northern Hemisphere spring and show a significant underestimation of hydrocarbons based on surface observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Simulations of Seismic Wave Propagation on Mars

In this paper, we present global and regional synthetic seismograms computed for 1D and 3D Mars models based on the spectral-element method. For global simulations, we implemented a radially-symmetric Mars model with a 110 km thick crust. For this 1D model, we successfully benchmarked the 3D seismic wave propagation solver SPECFEM3D_GLOBE against the 2D axisymmetric wave propagation solver AxiSEM at periods down to 10 s. We also present higher-resolution body-wave simulations with AxiSEM down to 1 s in a model with a more complex 1D crust, revealing wave propagation effects that would have been difficult to interpret based on ray theory. For 3D global simulations based on SPECFEM3D_GLOBE, we superimposed 3D crustal thickness variations capturing the distinct crustal dichotomy between Mars’ northern and southern hemispheres, as well as topography, ellipticity, gravity, and rotation. The global simulations clearly indicate that the 3D crust speeds up body waves compared to the reference 1D model, whereas it significantly changes surface waveforms and their dispersive character depending on its thickness. We also perform regional simulations with the solver SES3D based on 3D crustal models derived from surface composition, thereby addressing the effects of various distinct crustal features down to 2 s. The regional simulations confirm the strong effects of crustal variations on waveforms. Finally, we conclude that the numerical tools are ready for examining more scenarios, including various other seismic models and sources.

58 GEOSCIENCES↗