SEARCH · Search NASA
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.
This is my PhD dissertation that describes two computational methods, the finite-difference frequency-domain method for general bianisotropic media and the spatially-variant lattice algorithm with the finite element method.
We present an algorithm for the numerical simulation of seismic wave propagation in models with a complex near surface part and free surface topography. The approach is based on the combination of finite differences with the discontinuous Galerkin method. The discontinuous Galerkin method can be used on polyhedral meshes; thus, it is easy to handle the complex surfaces in the models. However, this approach is computationally intense in comparison with finite differences. Finite differences are computationally efficient, but in general, they require rectangular grids, leading to the stair-step approximation of the interfaces, which causes strong diffraction of the wavefield. In this research we present a hybrid algorithm where the discontinuous Galerkin method is used in a relatively small upper part of the model and finite differences are applied to the main part of the model.
The NESTLE code is a few-group, diffusion-based core simulation code that employs the Nodal Expansion Method (NEM). NESTLE has the capability to solve both steady-state and transient problems, using either Cartesian or hexagonal geometry, and to apply thermal-hydraulic feedback including single and two-phase flow. A modern version of the NESTLE nodal-diffusion simulator is currently maintained and developed at the University of Tennessee. The original version of the code was limited to only two or four neutron energy groups, therefore, this article describes recent work carried out to expand NESTLE into a general multi-group solver. Generalization of a fixed group finite-difference diffusion theory into arbitrary N groups is a well-known problem [3], however, this paper describes the intricacies associated with the generalization of the Nodal Expansion Method, and the methodical and computational challenges encountered and managed along the way.
We developed a low order Mimetic Finite Difference (MFD) discretization for the equations of magnetohydrodynamics (MHD) in two space dimensions. These equations describe the evolution of the electric and magnetic fields in the presence of prescribed velocity field.The method is designed to work on general polygonal meshes and preserves the divergence-free condition on the magnetic field. The electric field is discretized at the vertices/nodes and the magnetic field uses edge-based discretization. The method reconstructs the magnetic field to extract nodal values necessary to approximate some terms present in Ohm’s law. We test the robustness of our numerical scheme on three different types of meshes: with triangular elements, quadrilateral and unstructured polyhedrons obtain from a Voronoi tesselation. Analysis of the convergence for each of the aforementioned mesh types is presented. We finish with a test problem that shows the method is capable of modelling magnetic reconnection.
The Coarse Mesh Finite Difference (CMFD) method and the Generalized Coarse Mesh Rebalance (GCMR) method have been applied to many reactor physics codes and these methods contribute for improvement of calculation speed. On the other hand, the numerical instability of these methods have been a crucial issue. In order to overcome this issue, some stabilization techniques such as damping factor and correction techniques of diffusion coefficient have been developed. The systematic comparison of the calculation stability and efficiency between GCMR and CMFD with various stabilization techniques has not been reported and these results are important for the improvement of nuclear reactor analysis. The present summary focuses on this point. Namely, calculation stability and efficiency of the CMFD and the GCMR are compared with investigation of the effect of stabilization techniques on these methods. We applied these methods and techniques to 2D geometries with and without voided regions. (authors)
The oxidation of spent uranium carbide fuel, a candidate fuel for Generation IV nuclear reactors, is an important process in its potential reprocessing cycle. However, the oxidation of uranium carbide in air is highly exothermic. A model has therefore been developed to predict the temperature rise, as well as other useful information such as reaction completion times, under different reaction conditions in order to help in deriving safe oxidation conditions. Finite difference-methods are used to model the heat and mass transfer processes occurring during the reaction in two dimensions and are coupled to kinetics found in the literature.
Here, the linear and nonlinear benchmarks between the CLT code and the M3D-C1 code for the 2/1 resistive tearing mode and the 1/1 resistive kink mode are presented. CLT is an explicit finite difference code, while M3D-C1 is an implicit finite element code. Although the implementations of CLT and M3D-C1 are totally different, we find that the simulation results of the resistive-kink mode and the m/n=2/1 tearing mode from M3D-C1 and CLT are almost the same, including the linear and nonlinear growth rates, the mode structures, the nonlinear saturation levels, the Poincare plots, and the scaling laws. This confirms that the nonlinear results for the 1/1 resistive-kink mode and 2/1 tearing mode are accurate and reliable.
The Landau–Lifshitz equation describes the dynamics of the magnetization inside ferromagnetic materials. This equation is highly nonlinear and has a non-convex constraint (the magnitude of the magnetization is constant) which poses interesting challenges in developing numerical methods. We develop and analyze explicit and implicit mimetic finite difference schemes for this equation. These schemes work on general polytopal meshes which provide enormous flexibility to model magnetic devices with various shapes. A projection on the unit sphere is used to preserve the magnitude of the magnetization. We also provide a proof that shows the exchange energy is decreasing in certain conditions. The developed schemes are tested on general meshes that include distorted and randomized meshes. As a result, the numerical experiments include a test proposed by the National Institute of Standard and Technology and a test showing formation of domain wall structures in a thin film.
Synchronizations of processing elements (PEs) in massively parallel simulations, which arise due to communication or load imbalances between PEs, significantly affect the scalability of scientific applications. We have recently proposed a method based on finite-difference schemes to solve partial differential equations in an asynchronous fashion – synchronization between PEs is relaxed at a mathematical level. While standard schemes can maintain their stability in the presence of asynchrony, their accuracy is drastically affected. In this work, we present a general methodology to derive asynchrony-tolerant (AT) finite difference schemes of arbitrary order of accuracy, which can maintain their accuracy when synchronizations are relaxed. We show that there are several choices available in selecting a stencil to derive these schemes and discuss their effect on numerical and computational performance. We provide a simple classification of schemes based on the stencil and derive schemes that are representative of different classes. Their numerical error is rigorously analyzed within a statistical framework to obtain the overall accuracy of the solution. Results from numerical experiments are used to validate the performance of the schemes.
This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.
The quest to understand three-body dynamics from first-principle QCD includes the study of nonresonant and resonant systems. The isospin I = 2 system is of particular interest having no three-body resonance but featuring a resonance in a subchannel, while also being a coupled-channel problem. In this study, we calculate the finite-volume spectrum from lattice QCD at two different pion masses, map the amplitude to the infinite volume through a generalized Finite-Volume Unitarity three-body quantization condition, investigate the limit of a narrow ρ , and compare with an effective Lagrangian prediction at leading order. Chiral extrapolations between different pion masses are performed.
The growth dynamics of gold-silver core-shell (Au@Ag) nanoparticles are studied using in situ time-dependent second harmonic generation (SHG) and extinction spectroscopy to investigate the nanoparticle shell formation. The silver shell is grown by reduction of silver cations onto a 14 nm gold core using ascorbic acid in colloidal aqueous solution under varying reaction concentrations producing Au@Ag nanoparticles of final sizes ranging from 51 to 78 nm in diameter. The in situ extinction spectra show a rapid increase in intensity on the timescale of 5–6 s with blue shifting and narrowing of the plasmonic peak during the silver shell formation. The in situ SHG signals show an abrupt rise at early times of the reaction, followed by a time-dependent biexponential decrease, where the faster SHG lifetime corresponds to the timescale of the shell growth, and where the slower SHG lifetime is attributed to changes in the nanoparticle surface charge density. A large enhancement in the SHG signal at early stages of the reaction is caused by plasmonic hot spots due to the nanoparticle surface morphology, which becomes smoother as the reaction proceeds. The final extinction spectra are compared to finite-difference time-domain (FDTD) calculations, showing general agreement with experiment, where the plasmon peak red shifts and increases in spectral width as the silver shell thickness increases. These in situ SHG and extinction spectroscopy results, combined with FDTD calculations, help characterize the complicated processes involved in colloidal nanoparticle shell formation in real time for developing potential plasmon-enhanced nanomaterial applications.
Here, we present the arbitrary order mimetic finite difference (MFD) discretization for the diffusion equation with non-symmetric tensorial diffusion coefficient in a mixed formulation on general polygonal meshes. The diffusion tensor is assumed to be positive definite. The asymmetry of the diffusion tensor requires changes to the standard MFD construction. We also present new approach for the construction that guarantees positive definiteness of the non-symmetric mass matrix in the space of discrete velocities. The numerically observed convergence rate for the scalar quantity matches the predicted one in the case of the lowest order mimetic scheme. For higher orders schemes, we observed super-convergence by one order for the scalar variable which is consistent with the previously published result for a symmetric diffusion tensor. The new scheme was also tested on a time-dependent problem modeling the Hall effect in the resistive magnetohydrodynamics.
An added-mass partitioned (AMP) algorithm is described for solving fluid–structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully second-order accurate and stable, without sub-time-step iterations, even for very light structures when added-mass effects are strong. The fluid, governed by the incompressible Navier–Stokes equations, is solved in velocity-pressure form using a fractional-step method; large deformations are treated with a mixed Eulerian-Lagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler–Bernoulli beam model, and these equations are solved in a Lagrangian frame using two approaches, one based on finite differences and the other on finite elements. The key AMP interface condition is a generalized Robin (mixed) condition on the fluid pressure. This condition, which is derived at a continuous level, has no adjustable parameters and is applied at the discrete level to couple the partitioned domain solvers. Special treatment of the AMP condition is required to couple the finite-element beam solver with the finite-difference-based fluid solver, and two coupling approaches are described. A normal-mode stability analysis is performed for a linearized model problem involving a beam separating two fluid domains, and it is shown that the AMP scheme is stable independent of the ratio of the mass of the fluid to that of the structure. A traditional partitioned (TP) scheme using a Dirichlet–Neumann coupling for the same model problem is shown to be unconditionally unstable if the added mass of the fluid is too large. A series of benchmark problems of increasing complexity are considered to illustrate the behavior of the AMP algorithm, and to compare the behavior with that of the TP scheme. The results of all these benchmark problems verify the stability and accuracy of the AMP scheme. Results for one benchmark problem modeling blood flow in a deforming artery are also compared with corresponding results available in the literature.
Many geologic materials and minerals are seismically anisotropic, with the most general anisotropic material having up to 21 independent elastic coefficients. This report outlines the development of a 3-D, generally anisotropic, linear elastic full waveform finite-difference solver. First, a mathematical description of the solution equations will be described. The finite-difference implementation of these equations will then be shown. Finally, a comparison of results from this new solver to other solutions will be provided as verification that the new algorithm can accurately replicate these solutions.
Underground explosions can produce infrasound in the atmosphere, and the wavefield characteristics are often governed by the ground surface motions. Finite-difference methods are popular for infrasound simulation as their generality and robustness allow for complex atmospheric structures and surface topography. A simple point-source approximation is often used because infrasound wavelengths tend to be large relative to the source dimensions. However, this assumption may not be able to capture the complexity of explosion-induced ground motions if the surface area is not compact, and appropriate source models must be incorporated into the finite-difference simulations for accurate infrasound prediction. In this study, we develop a point source representation of the complex ground motions for infrasound sources. Instead of a single point source, we use a series of point sources distributed over the source area. These distributed point sources can be equivalent to air volume changes produced by the ground motions in the atmosphere. Further, we apply the distributed point-source method to a series of buried chemical explosions conducted during the Source Physics Experiment Phase I. Epicentral ground-motion measurements during the experiments provide a way to calculate accurate distributed point sources. We validate and evaluate the accuracy of distributed point source approach for infrasound simulations by direct comparison with acoustic observations in the field experiment.
Charging of two conducting spheres in a weakly ionized collisional plasma flow is considered. The spheres are arranged along the flow, and the plasma is assumed to consist of two ion species with the charges equal in magnitude but opposite in sign. The problem is analyzed with allowance for the external electric field, charging of the spheres due to the sedimentation of plasma ions on them, the fields of the sphere charges, the space charge field, and the processes of recombination and molecular diffusion. The interaction between the spheres and plasma is studied by numerically solving a time-dependent problem in a bispherical coordinate system by the finite difference method. The steady-state values of the sphere charges and the distributions of the space charge and ion densities in the ambient plasma are found as functions of the plasma parameters and the distance between the spheres. The electrostatic forces acting on the spheres are determined, and the effects of the external field, the space charge fields, and the fields of the sphere charges are comparatively analyzed. It is shown that, for the considered plasma parameters, the main electrostatic effect in the interaction between two spheres is their mutual approach in the external field due to the difference in their charges (one sphere catches up with the other). Due to the friction force with the neutral gas, this mutual approach is much slower than all other processes in the system. For widely spaced spheres, the results coincide with the solution obtained previously for a solitary sphere.
In this work we consider a modification to an arbitrary order mixed mimetic finite difference method (MFD) for a diffusion equation on general polygonal meshes [1]. The modification is based on moving some degrees of freedom (DoF) for a flux variable from edges to vertices. We showed that for a non-degenerate element this transformation is locally equivalent, i.e. there is a one-to-one map between the new and the old DoF. Globally, on the other hand, this transformation leads to a reduction of the total number of degrees of freedom (by up to 40%) and additional continuity of the discrete flux.