Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗
DOE OSTI2014

Vertically aligned crystalline silicon nanowires with controlled diameters for energy conversion applications: Experimental and theoretical insights

Vertically orientated single crystalline silicon nanowire (SiNW) arrays with controlled diameters are fabricated via a metal-assisted chemical etching method. The diameter of the fabricated nanowires is controlled by simply varying the etching time in HF/H{sub 2}O{sub 2} electrolytes. The fabricated SiNWs have diameters ranging from 117 to 650 nm and lengths from 8 to 18 μm. The optical measurements showed a significant difference in the reflectance/absorption of the SiNWs with different diameters, where the reflectance increases with increasing the diameter of the SiNWs. The SiNWs showed significant photoluminescence (PL) emission spectra with peaks lying between 380 and 670 nm. The PL intensity increases as the diameter increases and shows red shift for peaks at ∼670 nm. The increase or decrease of reflectivity is coincident with PL intensity at wavelength ∼660 nm. The x-ray diffraction patterns confirm the high crystallinity of the fabricated SiNWs. In addition, the Raman spectra showed a shift in the first order transverse band toward lower frequencies compared to that usually seen for c-Si. Finite difference time domain simulations have been performed to confirm the effect of change of diameter on the optical properties of the nanowires. The simulation results showed good agreement with the experimental results for the SiNWs of different diameters.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Large eddy simulation/dynamic thickened flame modeling of a high Karlovitz number turbulent premixed jet flame

Due to the complex multiscale interaction between intense turbulence and relatively weak flames, turbulent premixed flames in the thin and broken reaction zones regimes exhibit strong finite-rate chemistry and strain effects and are hence challenging to model. Here, in this paper, a laboratory premixed jet flame in the broken reaction zone, which has recently been studied using direct numerical simulation (DNS), is modeled using a large eddy simulation (LES)/dynamic thickened flame (DTF) approach with detailed chemistry. The presence of substantial flame thickening due to strong turbulence-chemistry interactions, which can be characterized by a high Karlovitz number ( Ka ), requires the DTF model to thicken the flame in an adaptive way based on the local resolution of flame scales. Here, an appropriate flame sensor and strain-sensitive flame thickness are used to automatically determine the thickening location and thickening factor, respectively. To account for finite-rate chemistry and strain effects, the chemistry is described in two different ways: (1) detailed chemistry denoted as full transport and chemistry (FTC), and (2) tabulated chemistry based on a strained premixed flamelet (SPF) model. The performance of the augmented LES/DTF approach for modeling the high Ka premixed flame is assessed through detailed a posteriori comparisons with DNS of the same flame. It is found that the LES/DTF/FTC model is capable of reproducing most features of the high Ka turbulent premixed flame including accurate CO and NO prediction. The LES/DTF/SPF model has the potential to capture the impact of strong turbulence on the flame structure and provides reasonable prediction of pollutant emissions at a reasonable computational cost. In order to identify the impact of aerodynamic strain, the turbulent flame structure is analyzed and compared with unstrained and strained premixed flamelet solutions. Lastly, the results indicate that detailed strain effects should be considered when using tabulated methods to model high Ka premixed flames.

42 ENGINEERING↗
DOE OSTI2015

On the scalability of the Albany/FELIX first-order Stokes approximation ice sheet solver for large-scale simulations of the Greenland and Antarctic ice sheets

We examine the scalability of the recently developed Albany/FELIX finite-element based code for the first-order Stokes momentum balance equations for ice flow. We focus our analysis on the performance of two possible preconditioners for the iterative solution of the sparse linear systems that arise from the discretization of the governing equations: (1) a preconditioner based on the incomplete LU (ILU) factorization, and (2) a recently-developed algebraic multigrid (AMG) preconditioner, constructed using the idea of semi-coarsening. A strong scalability study on a realistic, high resolution Greenland ice sheet problem reveals that, for a given number of processor cores, the AMG preconditioner results in faster linear solve times but the ILU preconditioner exhibits better scalability. A weak scalability study is performed on a realistic, moderate resolution Antarctic ice sheet problem, a substantial fraction of which contains floating ice shelves, making it fundamentally different from the Greenland ice sheet problem. Here, we show that as the problem size increases, the performance of the ILU preconditioner deteriorates whereas the AMG preconditioner maintains scalability. This is because the linear systems are extremely ill-conditioned in the presence of floating ice shelves, and the ill-conditioning has a greater negative effect on the ILU preconditioner than on the AMG preconditioner.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Alternative modeling methods for plasma-based Rf ion sources

Rf-driven ion sources for accelerators and many industrial applications benefit from detailed numerical modeling and simulation of plasma characteristics. For instance, modeling of the Spallation Neutron Source (SNS) internal antenna H{sup −} source has indicated that a large plasma velocity is induced near bends in the antenna where structural failures are often observed. This could lead to improved designs and ion source performance based on simulation and modeling. However, there are significant separations of time and spatial scales inherent to Rf-driven plasma ion sources, which makes it difficult to model ion sources with explicit, kinetic Particle-In-Cell (PIC) simulation codes. In particular, if both electron and ion motions are to be explicitly modeled, then the simulation time step must be very small, and total simulation times must be large enough to capture the evolution of the plasma ions, as well as extending over many Rf periods. Additional physics processes such as plasma chemistry and surface effects such as secondary electron emission increase the computational requirements in such a way that even fully parallel explicit PIC models cannot be used. One alternative method is to develop fluid-based codes coupled with electromagnetics in order to model ion sources. Time-domain fluid models can simulate plasma evolution, plasma chemistry, and surface physics models with reasonable computational resources by not explicitly resolving electron motions, which thereby leads to an increase in the time step. This is achieved by solving fluid motions coupled with electromagnetics using reduced-physics models, such as single-temperature magnetohydrodynamics (MHD), extended, gas dynamic, and Hall MHD, and two-fluid MHD models. We show recent results on modeling the internal antenna H{sup −} ion source for the SNS at Oak Ridge National Laboratory using the fluid plasma modeling code USim. We compare demonstrate plasma temperature equilibration in two-temperature MHD models for the SNS source and present simulation results demonstrating plasma evolution over many Rf periods for different plasma temperatures. We perform the calculations in parallel, on unstructured meshes, using finite-volume solvers in order to obtain results in reasonable time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Exotic phases in finite-density ℤ 3 theories

Lattice ℤ 3 theories with complex actions share many key features with finite- density QCD including a sign problem and $\mathcal{CK}$ symmetry. Complex ℤ 3 spin and gauge models exhibit a generalized Kramers-Wannier duality mapping them onto chiral ℤ 3 spin and gauge models, which are simulatable with standard lattice methods in large regions of parameter space. The Migdal-Kadanoff real-space renormalization group (RG) preserves this duality, and we use it to compute the approximate phase diagram of both spin and gauge ℤ 3 models in dimensions one through four. Chiral ℤ 3 spin models are known to exhibit a Devil’s Flower phase structure, with inhomogeneous phases that can be thought of as ℤ 3 analogues of chiral spirals. Out of the large class of models we study, we find that only chiral spin models and their duals have a Devil’s Flower structure with an infinite set of inhomogeneous phases, a result we attribute to Elitzur’s theorem. We also find that different forms of the Migdal-Kadanoff RG produce different numbers of phases, a violation of the expectation for universal behavior from a real-space RG. We discuss extensions of our work to ℤ N models, SU(N) models and nonzero temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Investigation of a passively Q-switched Raman laser at 1176 nm with Nd{sup 3+}:YAG/Cr{sup 4+}:YAG/YAG composite crystal and a coupled cavity

To reduce the thermal lens effect and increase the efficiency, a composite crystal of Nd{sup 3+}:YAG/Cr{sup 4+}:YAG/YAG and a coupled cavity structure were employed in our passively Q-switched YVO{sub 4}-based Raman laser system. Both the simulation and the experiment results showed that the undoped YAG crystal bonded in the end-face of Cr{sup 4+}:YAG crystal can speed up the heat transfer of Cr{sup 4+}:YAG crystal and reduce the temperature difference between its side surfaces and the center significantly, and the coupled cavity structure can improve the laser performance, especially at high pump power. The maximum output power of Raman laser of 0.92 W at 1176 nm was obtained with an optical-to-optical efficiency of 11.51%. The corresponding repetition rate and the pulse width is 53.5 kHz and 3.22 ns, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Revisiting Multi-Material Composite Structures with Homogenized Composite Properties

Composite structures inherently develop residual stresses during their curing process. Driven predominately by mismatched thermal strains between differing materials or ply orientations, but also affected by curing process phenomena like polymer shrinkage, these residual stresses can lead to failure within composite structures. There are several methods varying in complexity that can be used to model the development of residual stresses, all of which are capable of capturing sufficient detail to understand the residual stress state at the ply level. However, explicitly modeling all plies of a layup in a composite structure can be prohibitively expensive based on the number of plies, structure size, and required element size. The computational cost can be reduced through the homogenization of the composite layup without losing much fidelity of the overall response of the structure. The homogenization process reduces the many plies of a laminate to a single lamina that reduces complexity and increases the mesh size where a single element can span multiple plies. This report focuses on verification and validation efforts for a homogenization process using a suite of finite element simulations rather than an analytic solution derived from classical laminate theory. Initial verification using representative element volumes indicated there was minimal error in the homogenization process; however, this compounded to a small, but acceptable error in strip and split ring experimental composite structures. The error does under predict the residual stress state in the strip and split ring and should be accounted for when simulating composite structures with homogenized properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗
DOE OSTI2014

Inflection points of microcanonical entropy: Monte Carlo simulation of q state Potts model on a finite square lattice

Traditional definition of phase transition involves an infinitely large system in thermodynamic limit. Finite systems such as biological proteins exhibit cooperative behavior similar to phase transitions. We employ recently discovered analysis of inflection points of microcanonical entropy to estimate the transition temperature of the phase transition in q state Potts model on a finite two dimensional square lattice for q=3 (second order) and q=8 (first order). The difference of energy density of states (DOS) Δ ln g(E) = ln g(E+ ΔE) −ln g(E) exhibits a point of inflexion at a value corresponding to inverse transition temperature. This feature is common to systems exhibiting both first as well as second order transitions. While the difference of DOS registers a monotonic variation around the point of inflexion for systems exhibiting second order transition, it has an S-shape with a minimum and maximum around the point of inflexion for the case of first order transition.

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

Physical-mass calculation of ρ ( 770 ) and K * ( 892 ) resonance parameters via π π and K π scattering amplitudes from lattice QCD

We present our study of the ρ ( 770 ) and K * ( 892 ) resonances from lattice quantum chromodynamics (QCD) employing domain-wall fermions at physical quark masses. We determine the finite-volume energy spectrum in various momentum frames and obtain phase-shift parametrizations via the Lüscher formalism and as a final step the complex resonance poles of the π π and K π elastic scattering amplitudes via an analytical continuation of the models. By sampling a large number of representative sets of underlying energy-level fits, we also assign a systematic uncertainty to our final results. This is a significant extension to data-driven analysis methods that have been used in lattice QCD to date, due to the two-step nature of the formalism. Our final pole positions, M + i Γ / 2 , with all statistical and systematic errors exposed, are M K * = 893 ( 2 ) ( 8 ) ( 54 ) ( 2 ) MeV and Γ K * = 51 ( 2 ) ( 11 ) ( 3 ) ( 0 ) MeV for the K * ( 892 ) resonance and M ρ = 796 ( 5 ) ( 15 ) ( 48 ) ( 2 ) MeV and Γ ρ = 192 ( 10 ) ( 28 ) ( 12 ) ( 0 ) MeV for the ρ ( 770 ) resonance. The four differently grouped sources of uncertainties are, in the order of occurrence: statistical, data-driven systematic, an estimation of systematic effects beyond our computation (dominated by the fact that we employ a single lattice spacing), and the error from the scale-setting uncertainty on our ensemble. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Grating-based x-ray differential phase contrast imaging with twin peaks in phase-stepping curves—phase retrieval and dewrapping

Purpose: X-ray differential phase contrast CT implemented with Talbot interferometry employs phase-stepping to extract information of x-ray attenuation, phase shift, and small-angle scattering. Since inaccuracy may exist in the absorption grating G{sub 2} due to an imperfect fabrication, the effective period of G{sub 2} can be as large as twice the nominal period, leading to a phenomenon of twin peaks that differ remarkably in their heights. In this work, the authors investigate how to retrieve and dewrap the phase signal from the phase-stepping curve (PSC) with the feature of twin peaks for x-ray phase contrast imaging. Methods: Based on the paraxial Fresnel–Kirchhoff theory, the analytical formulae to characterize the phenomenon of twin peaks in the PSC are derived. Then an approach to dewrap the retrieved phase signal by jointly using the phases of the first- and second-order Fourier components is proposed. Through an experimental investigation using a prototype x-ray phase contrast imaging system implemented with Talbot interferometry, the authors evaluate and verify the derived analytic formulae and the proposed approach for phase retrieval and dewrapping. Results: According to theoretical analysis, the twin-peak phenomenon in PSC is a consequence of combined effects, including the inaccuracy in absorption grating G{sub 2}, mismatch between phase grating and x-ray source spectrum, and finite size of x-ray tube’s focal spot. The proposed approach is experimentally evaluated by scanning a phantom consisting of organic materials and a lab mouse. The preliminary data show that compared to scanning G{sub 2} over only one single nominal period and correcting the measured phase signal with an intuitive phase dewrapping method that is being used in the field, stepping G{sub 2} over twice its nominal period and dewrapping the measured phase signal with the proposed approach can significantly improve the quality of x-ray differential phase contrast imaging in both radiograph and CT. Conclusions: Using the phase retrieval and dewrapping methods proposed to deal with the phenomenon of twin peaks in PSCs and phase wrapping, the performance of grating-based x-ray differential phase contrast radiography and CT can be significantly improved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Two-phase Flow Regime Transition Study using Level-Set Method

In Boiling Water Reactors (BWRs) the coolant exiting the core often has void fractions as high as 50%, resulting in high gas volume fraction two-phase flow regimes, such as annular or churn-turbulent flow. As the coolant boils along the subchannels it will experience flow regime transitions from bubbly to slug, slug to churn, and, finally, churn to annular. While these transition points are known to occur as a higher void develops along the flow path, the driving physics of the transition phenomenon is not well understood. Currently, experiment based correlations are the most reliable method for predicting where in the core the flow regime transitions will occur. However, experimentally, it's ever difficult to observe the details of the two-phase flow at prototypical pressures/temperatures, or, as with high fidelity experiments, the conditions differ from those in a reactor. Being able to predict these flow regime transitions is critical to safe and efficient reactor operations as well as maintaining fuel rod integrity by avoiding dryout. The characterization of flow regimes has been extensively studied, producing a wide set of flow regime maps based on characteristic parameters such as vapor and liquid superficial velocity, time dependent pressure drop, interfacial area density, etc. While these studies are useful in predicting the flow regime for certain geometry and flow conditions, they do not provide a universal predictive capability for any geometry under any conditions. Additionally, it is difficult to achieve reactor realistic pressure and temperature conditions in experiments. Direct Numerical Simulation (DNS) coupled with Interface Tracking Methods (ITMs) becomes advantageous in this case because the simulations are not limited by the material properties or to the real world conditions. DNS allows for analysis of flows that would be impossible to achieve, or observe in detail, experimentally. The finite element based DNS code, PHASTA (parallel, hierarchic, higher-order accurate, adaptive, stabilized transient analysis flow solver), used in this work utilizes the level set ITM and has been previously shown to properly resolve bubbly, slug, and annular flows, but not transitions between the two-phase flow regimes. The present work uses PHASTA to analyze flow in a pipe for various void fractions and superficial velocities. Flow regime transitions were induced to demonstrate the transient behavior of this phenomenon. With the unique advanced analysis capabilities developed in PHASTA that allow for individual bubble tracking, it becomes possible to extract new localized information about the flow behavior that is currently not attainable from experiments. The potential of this capability is to uncover the underlying physics of the flow regime transition phenomenon. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Temporal breakdown and Borel resummation in the complex Langevin method

We reexamine the Parisi-Klauder conjecture for complex e{sup i{theta}/2}{phi}{sup 4} measures with a Wick rotation angle 0{<=}{theta}/2{<=}{pi}/2 interpolating between Euclidean signature and Lorentzian signature. Our main result is that the asymptotics for short stochastic times t encapsulates information also about the equilibrium aspects. The moments evaluated with the complex measure and with the real measure defined by the stochastic Langevin equation have the same t{yields}0 asymptotic expansion which is shown to be Borel summable. The Borel transform correctly reproduces the time dependent moments of the complex measure for all t, including their t{yields}{infinity} equilibrium values. On the other hand the results of a direct numerical simulation of the Langevin moments are found to disagree from the 'correct' result for t larger than a finite t{sub c}. The breakdown time t{sub c} increases powerlike for decreasing strength of the noise's imaginary part but cannot be excluded to be finite for purely real noise. To ascertain the discrepancy we also compute the real equilibrium distribution for complex noise explicitly and verify that its moments differ from those obtained with the complex measure. - Highlights: Black-Right-Pointing-Pointer The Parisi-Klauder conjecture is reexamined for complex e{sup i{theta}/2}{phi}{sup 4} measures. Black-Right-Pointing-Pointer The time dependent moments are evaluated by temporal Borel resummation. Black-Right-Pointing-Pointer The results disagree with the Langevin simulations beyond a critical time t{sub c}. Black-Right-Pointing-Pointer t{sub c} increases with decreasing strength of the noise's imaginary part. Black-Right-Pointing-Pointer The technical reason for the breakdown is identified.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Investigations of α-helix↔β-sheet transition pathways in a miniprotein using the finite-temperature string method

A parallel implementation of the finite-temperature string method is described, which takes into account the invariance of coordinates with respect to rigid-body motions. The method is applied to the complex α-helix↔β-sheet transition in a β-hairpin miniprotein in implicit solvent, which exhibits much of the complexity of conformational changes in proteins. Two transition paths are considered, one derived from a linear interpolant between the endpoint structures and the other derived from a targeted dynamics simulation. Two methods for computing the conformational free energy (FE) along the string are compared, a restrained method, and a tessellation method introduced by E. Vanden-Eijnden and M. Venturoli [J. Chem. Phys. 130, 194103 (2009)]. It is found that obtaining meaningful free energy profiles using the present atom-based coordinates requires restricting sampling to a vicinity of the converged path, where the hyperplanar approximation to the isocommittor surface is sufficiently accurate. This sampling restriction can be easily achieved using restraints or constraints. The endpoint FE differences computed from the FE profiles are validated by comparison with previous calculations using a path-independent confinement method. The FE profiles are decomposed into the enthalpic and entropic contributions, and it is shown that the entropy difference contribution can be as large as 10 kcal/mol for intermediate regions along the path, compared to 15–20 kcal/mol for the enthalpy contribution. This result demonstrates that enthalpic barriers for transitions are offset by entropic contributions arising from the existence of different paths across a barrier. The possibility of using systematically coarse-grained representations of amino acids, in the spirit of multiple interaction site residue models, is proposed as a means to avoid ad hoc sampling restrictions to narrow transition tubes.

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

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A combined experimental and crystal plasticity study of grain size effects in magnesium alloys

This work presents a method to incorporate the micro Hall-Petch equation into the crystal plasticity finite element (CPFE) framework accounting for the microstructural features to understand the coupling between grain size, texture, and loading direction in magnesium alloys. The effect of grain size and texture is accounted for by modifying the slip resistances of individual basal and prismatic systems based on the micro Hall-Petch equation. The modification based on the micro Hall-Petch equation endows every slip system at each microstructural point with a slip system-level grain size and maximum compatibility factor, which are in turn used to modify the slip resistance. While the slip-system level grain size is a measure of the grain size, the maximum compatibility factor encodes the effect of the grain boundary on the slip system resistance modification and is computed based on the Luster-Morris factor. The model is calibrated using experimental stress-strain curves of Mg-4Al samples with three different grain sizes from which the Hall-Petch coefficients are extracted and compared with Hall-Petch coefficients predicted using original parameters from previous work. The predictability of the model is then evaluated for a Mg-4Al sample with different texture and three grain sizes subjected to loading in different directions. The calibrated parameters are then used for some parametric studies to investigate the variation of Hall-Petch slope for different degrees of simulated spread in basal texture, variation of Hall-Petch slope with loading direction relative to basal poles for a microstructure with strong basal texture, and variation of yield strength with change in grain morphology. The proposed approach to incorporate the micro Hall-Petch equation into the CPFE framework provides a foundation to quantitatively model more complicated scenarios of coupling between grain size, texture and loading direction in the plasticity of Mg alloys.

Crystal plasticity↗
DOE OSTI2019

Process Modeling of U-10wt% Mo Alloys Using Integrated Computational Materials Engineering

Low-enriched uranium alloyed with 10 wt% molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium as fuel for the nation’s high-performance research reactors. Manufacturing the U-10Mo alloy consists of multiple complex thermomechanical processes, which is highly challenging for computational modeling. The integrated computational materials engineering (ICME) concept supports building a microstructure-based framework to investigate the effect of manufacturing processes on the material microstructure. In this report, an ICME model that combines homogenization, hot rolling, annealing, and cold rolling is presented as an example of implementing the ICME concept for modeling the manufacturing processes of the U-10Mo alloy. The integrated model enables information to be passed between the different model components, leading to an improved prediction and a better understanding of microstructure across multiple processes. The integrated model was also implemented to investigate the variation of chemical composition, thickness of the zirconium interlayer, grain growth, and carbide redistribution and fracture of U-10Mo during the thermomechanical processes. With the homogenization model, a molybdenum concentration map can be reconstructed from any given microstructure, and the homogenization time needed to reach a desired homogenization level can also be approximated. The microstructure-based finite-element rolling model, shows that the homogenization process reduces the nonuniformity in the thickness of the zirconium interlayers on the U-10Mo. With the output from the rolling model, a statistical analysis of second-phase hard particles was conducted to quantify the particle redistribution and interparticle spacing between carbides during the deformation. Also, new stringer identification criteria were proposed to predict the correlation between stringer volume fractions and hot-rolling reduction. A macro-finite-element method hot-rolling model showed that higher reduction passes can be achieved within the mill separation force range, and a new “aggressive” rolling schedule was proposed. The model also showed the importance of can material on the quality of rolled foils, with stronger can materials resulting in hot-rolled U-10Mo free of defects (waviness, dog-boning). The deformation-induced recrystallization model demonstrated capability to use inputs from deformation simulations of strain accumulation in U-10Mo samples to predict grain growth behavior of the fuel during thermal treatment stages. Finally, a plane-strain compression mode simulated particle fracture and calculates the void volume fraction generated during cold rolling. The macro-finite-element method cold-rolling model showed that smaller diameter rolls are necessary to achieve the desired foil thickness.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗