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

A hybrid finite element-spectral boundary integral approach: Applications to dynamic rupture modeling in unbounded domains

The finite element method (FEM) and the spectral boundary integral method (SBI) have both been widely used in the study of dynamic rupture simulations along a weak interface. Here, we present a hybrid method that combines FEM and SBI through the consistent exchange of displacement and traction boundary conditions, thereby benefiting from the flexibility of FEM in handling problems with nonlinearities or small-scale heterogeneities and from the superior performance and accuracy of SBI. We validate the hybrid method using a benchmark problem from the Southern California Earthquake Center's dynamic rupture simulation validation exercises.We further demonstrate the capability and computational efficiency of the hybrid scheme for resolving off-fault heterogeneities by studying a 2D in-plane shear crack in two different settings: one where the crack is embedded in a high-velocity zone and another where it is embedded in a low-velocity zone. Finally, we discuss the potential of the hybrid method for addressing a wide range of problems in geophysics and engineering.

58 GEOSCIENCES↗
DOE OSTI2018

Phase diagrams and Hofstadter butterflies in the strongly correlated bosonic systems on the lattices with Dirac points

Gauge potentials with different configurations have been recently realized in the optical lattice experiments. It is remarkable that one of the simplest gauge potential can generate particle energy spectrum with the self-similar structure known as a Hofstadter butterfly. We investigate theoretically the impact of strong on-site interaction on such a spectrum in the bosonic Mott insulator within Bose–Hubbard model. In particular, it is shown that the fractal structure is encoded in the quasi-particle and hole bosonic branches for different lattice backgrounds. For example a square lattice and other structures (brick-wall and staggered magnetic flux lattice) which contain Dirac points in energy dispersions are considered. This shows that single-particle physics is still present even in the strong interaction limit for whole Hofstadter spectrum. Additionally we observe, that although in brick-wall and staggered flux lattices the quasi-particle densities of states look qualitatively similar, the corresponding Hofstadter butterfly assumes different forms. In particular, we use a superposition of two different synthetic gauge fields which appears to be a generator of non-trivial phenomena in the optical lattice systems. We also discuss the consequences of these phenomena on the phase diagrams between bosonic Mott insulator and superfluid phase. The analysis is carried out within the strong coupling expansion method on the finite size lattices and also at finite temperatures which are relevant for the currently made experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Fourier Analysis-Based Evolutionary Multi-Objective Multiphysics Optimization of Liquid-Cooled Heat Sinks

Optimal heat dissipation in power modules can significantly increase their power density. Removing the generated heat is critical for capturing the benefits of advanced semiconductor materials and improving the reliability of the device operation. This article proposes a design optimization method for liquid-cooled heat sinks that use a Fourier analysis-based tool and an evolutionary optimization algorithm to optimize the heat sink geometry for specified objectives. The optimized heat sink geometry is then compared with state-of-the-art solutions in literature based on finite element analysis of different designs. The proposed methodology can develop complex geometries that outperform the conventional heat sink geometries.

Sahu, Raj↗
DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗
DOE OSTI2022

Quantum spin liquids bootstrapped from Ising criticality in Rydberg arrays

Arrays of Rydberg atoms constitute a highly tunable, strongly interacting venue for the pursuit of exotic states of matter. We develop a strategy for accessing a family of fractionalized phases known as quantum spin liquids in two-dimensional Rydberg arrays. We specifically use effective field theory methods to study arrays assembled from Rydberg chains tuned to an Ising phase transition that famously hosts emergent fermions propagating within each chain. This highly entangled starting point allows us to naturally access spin liquids familiar from Kitaev's honeycomb model—albeit from an entirely different framework. In particular, we argue that finite-range repulsive Rydberg interactions, which frustrate nearby symmetry-breaking orders, can enable coherent propagation of emergent fermions between the chains in which they were born. Delocalization of emergent fermions across the full two-dimensional Rydberg array yields a gapless Z 2 spin liquid with a single massless Dirac cone. Here, the Rydberg occupation numbers exhibit universal power-law correlations that provide a straightforward experimental diagnostic of this phase. We further show that explicitly breaking symmetries perturbs the gapless spin liquid into gapped, topologically ordered descendants: Breaking lattice symmetries generates toric-code topological order, whereas introducing Floquet-mediated chirality generates non-Abelian Ising topological order. In the toric-code phase, we analytically construct microscopic incarnations of non-Abelian defects, which can be created and transported by dynamically controlling the atom positions in the array. Our work suggests that appropriately tuned Rydberg arrays provide a cold-atoms counterpart of solid-state “Kitaev materials” and, more generally, it spotlights a different angle for pursuing experimental platforms for Abelian and non-Abelian fractionalization.

36 MATERIALS SCIENCE↗
DOE OSTI2018

SELF-LIMITING TRANSIENT PULSE SIMULATION METHOD EXHIBITING TIME LAG PHENOMENON USING MAMMOTH

Coupling of multi-physics and multi-scale allows the exploration of phenomenon not previously studied with high fidelity. The heat diffusion time lag in TREAT-like fuel requires both multi-physics and multi-scale coupling to model the full transient behavior. The reactor physics code MAMMOTH, which has the ability to couple the heat conduction equation and the neutron diffusion equation, simulated a self-limiting transient pulse meeting these requirements. The MOOSE MultiApp system allowed the coupling of many micro-scale simulations with a single macro-scale simulation. The micro-scale geometry was a 0.0044 cm diameter UO2 fuel grain surrounded by a graphite moderator. The macro-scale geometry was a homogeneous fuel with the same isotope ratio of the micro-scale material surrounded by a graphite reflector. Two temperature-dependent two energy group cross-section libraries were computed using Serpent 2 for use in MAMMOTH. Both the homogeneous unit cell and heterogeneous unit cell cross-section library had infinite reactor boundary conditions. The macro-scale simulation had finite boundary conditions. By varying the feedback temperature, macro or micro, given to the cross-sections, the effect of the time lag was explored using a Picard coupling method. The time lag increases the peak power density by 9 % to 10 % and the energy deposited by 8 % to 9 %. The effect of differing grain and moderator temperature at the micro-scale was also explored showing differences in time lag behavior. Future research will improve the fidelity of the method through representative macro-scale geometries, temperature-dependent cross-sections, and greater energy group resolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Excited-state uncertainties in lattice-QCD calculations of multi-hadron systems

Excited-state effects lead to hard-to-quantify systematic uncertainties in lattice quantum chromodynamics (LQCD) spectroscopy calculations when computationally accessible imaginary times are smaller than inverse excitation gaps, as often arises for multi-hadron systems with signal-to-noise problems. Lanczos residual bounds address this by providing two-sided constraints on energies that do not require assumptions beyond Hermiticity, but often give very conservative systematic uncertainty estimates. Here, a more-constraining set of gap bounds is introduced for hadron spectroscopy. These bounds provide tighter constraints whose validity requires an explicit assumption about an energy gap. Exactly solvable lattice field theory correlators are used to test the utility of residual and gap bounds at finite and infinite statistics. Two-sided bounds and other analysis methods are then applied to a high-statistics LQCD calculation of nucleon-nucleon scattering at $m_π\sim 800$ MeV. Generalized eigenvalue problem (GEVP) and Lanczos energy estimators are compatible when applied to the same correlator data, but analyses including different interpolating operators show statistically significant inconsistencies. However, two-sided bounds from all operators are consistent. Under the assumption that the number of energy levels below $NΔ$ and $ΔΔ$ thresholds is the same as for non-interacting nucleons, gap bounds are sufficient to constrain nucleon-nucleon scattering amplitudes at phenomenologically relevant precision. Lanczos methods further reveal that energy-eigenstate estimates from previously studied asymmetric correlators have not converged over accessible imaginary times. Nevertheless, data-driven examples demonstrate why assumptions are required to draw conclusions about the natures of two-nucleon ground states at these masses.

Detmold, William [MIT, Cambridge, CTP]↗
DOE OSTI2019

Uncertainty quantification methodology for hyperbolic systems with application to blood flow in arteries

We present a Stochastic Finite Volume - ADER (SFV-ADER) methodology for Uncertainty Quantification (UQ) in the general framework of systems of hyperbolic balance laws with uncertainty in parameters. The resulting method is weakly intrusive, meaning that deterministic solvers need minor modifications to include random contributions, and has no theoretical accuracy barrier. An illustration of the second-order version for the viscous Burgers equation with uncertain viscosity coefficient is first given in detail, under different deterministic initial conditions; the attainment of the theoretically expected convergence rate is demonstrated empirically; results and features of our scheme are discussed. We then extend the SFV-ADER method to a non-linear hyperbolic system with source terms that models one-dimensional blood flow in arteries, assuming uncertainties in a set of parameters of the problem. Results are then compared with measurements in a well-defined 1:1 experimental replica of the network of largest arteries in the human systemic circulation. For a 99%-confidence level, severe influence of parameter fluctuations is seen in pressure profiles, while secondary effects on flow rate profiles become visible, particularly in terminal branches. In conclusion, results suggest that the proposed methodology could successfully be applied to other physical problems involving hyperbolic balance laws.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Elastic geothermobarometry: Corrections for the geometry of the host-inclusion system

Elastic geothermobarometry on inclusions is a method to determine pressure-temperature conditions of mineral growth independent of chemical equilibrium. Because of the difference in their elastic properties, an inclusion completely entrapped inside a host mineral will develop a residual stress upon exhumation, from which one can back-calculate the entrapment pressure. Current elastic geobarometric models assume that both host and inclusion are elastically isotropic and have an ideal geometry (the inclusion is spherical and isolated at the center of an infinite host). These conditions do not commonly occur in natural rocks, and the consequences for inclusion pressures can only be quantified with numerical approaches. In this paper, we report the results of numerical simulations of inclusions with the finite element method on elastically isotropic systems. We define and determine a geometrical factor (Γ) that allows measured residual pressures to be corrected for the effects of non-ideal geometry. We provide simple guidelines as to which geometries can safely be used for elastic geobarometry without correcting for the geometry. We also show that the discrepancies between elastic and conventional geobarometry reported in literature are not due to geometrical effects, and therefore result from other factors not yet included in current models.

58 GEOSCIENCES↗
DOE OSTI2017

Merging symmetry projection methods with coupled cluster theory: Lessons from the Lipkin model Hamiltonian

Coupled cluster and symmetry projected Hartree-Fock are two central paradigms in electronic structure theory. However, they are very different. Single reference coupled cluster is highly successful for treating weakly correlated systems but fails under strong correlation unless one sacrifices good quantum numbers and works with broken-symmetry wave functions, which is unphysical for finite systems. Symmetry projection is effective for the treatment of strong correlation at the mean-field level through multireference non-orthogonal configuration interaction wavefunctions, but unlike coupled cluster, it is neither size extensive nor ideal for treating dynamic correlation. We here examine different scenarios for merging these two dissimilar theories. We carry out this exercise over the integrable Lipkin model Hamiltonian, which despite its simplicity, encompasses non-trivial physics for degenerate systems and can be solved via diagonalization for a very large number of particles. We show how symmetry projection and coupled cluster doubles individually fail in different correlation limits, whereas models that merge these two theories are highly successful over the entire phase diagram. In conclusion, despite the simplicity of the Lipkin Hamiltonian, the lessons learned in this work will be useful for building an ab initio symmetry projected coupled cluster theory that we expect to be accurate in the weakly and strongly correlated limits, as well as the recoupling regime.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Merging symmetry projection methods with coupled cluster theory: Lessons from the Lipkin model Hamiltonian

Coupled cluster and symmetry projected Hartree-Fock are two central paradigms in electronic structure theory. However, they are very different. Single reference coupled cluster is highly successful for treating weakly correlated systems, but fails under strong correlation unless one sacrifices good quantum numbers and works with broken-symmetry wave functions, which is unphysical for finite systems. Symmetry projection is effective for the treatment of strong correlation at the mean-field level through multireference non-orthogonal configuration interaction wavefunctions, but unlike coupled cluster, it is neither size extensive nor ideal for treating dynamic correlation. We here examine different scenarios for merging these two dissimilar theories. We carry out this exercise over the integrable Lipkin model Hamiltonian, which despite its simplicity, encompasses non-trivial physics for degenerate systems and can be solved via diagonalization for a very large number of particles. We show how symmetry projection and coupled cluster doubles individually fail in different correlation limits, whereas models that merge these two theories are highly successful over the entire phase diagram. Despite the simplicity of the Lipkin Hamiltonian, the lessons learned in this work will be useful for building an ab initio symmetry projected coupled cluster theory that we expect to be accurate in the weakly and strongly correlated limits, as well as the recoupling regime.

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

Effect of Weld Consumable Conditioning on the Diffusible Hydrogen and Subsequent Residual Stress and Flexural Strength of Multipass Welded P91 Steels

P91 steel weld joint was prepared using the shielded metal arc welding process and four different conditions of weld consumable that provide the different levels of diffusible hydrogen in deposited metal. In the present research, the effects of diffusible hydrogen content on the flexural strength, lower critical stress, and tensile strength of P91 steel welds were also determined with respect to different electrode conditions. To investigate the effect of diffusible hydrogen on multipass welding, top and root side flexural tests were performed. The residual stresses (axial stress and transverse stress) were also measured using the blind hole drilling method for different conditions of welding consumable. The peak value of residual stresses was measured at the center of the weld fusion zone. The maximum value of transverse stress was measured to be 355 MPa for case II (6.21 mL/100 g of diffusible hydrogen), while the maximum axial stress was about 218 MPa for case IV (12.43 mL/100 g of diffusible hydrogen). A three-dimensional finite element simulation was also performed to predict the residual stress distribution and thermal profile along the welded joint. The experimentally determined residual stresses correlated well with the numerically estimated residual stresses. The diffusible hydrogen content was not observed to have any significant effect on the residual stresses. The corrected residual stress values were also predicted by considering the plasticity-induced error. However, the flexural performance of the welded joint was affected by the diffusible hydrogen content. The top and root flexural strength was measured to be optimum for the low level of diffusible hydrogen content, and the values decreased with an increase in diffusible hydrogen content.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Influence of Forced Flow on the Dendritic Growth of Fe-C Alloy: 3D vs 2D Simulation

A 3D parallel cellular automaton-finite volume method (CA-FVM) model was used to simulate the equiaxed dendritic growth of an Fe-0.82 wt pct C alloy with xy-in-out and xyz-in-out type forced flows and the columnar dendritic growth with y-in-out type forced flow. In addition, the similarities and differences between the results of the 3D and 2D models are discussed and summarized in detail. The capabilities of the 3D and 2D CA-FVM models to predict the dendritic growth of the alloy with forced flow are validated through comparison with the boundary layer correction and Oseen–Ivanstov models, respectively. Because the forced flow can pass around perpendicular arms of the dendrites, the secondary arms at the sides upstream from the perpendicular arms are more developed than those on the upstream side of the upstream arms, especially at higher inlet velocities. In addition, compared to the xy-in-out case, the growth of the downstream arms is less inhibited and the secondary arms are more developed in the xyz-in-out case because of the greater lateral flow around their tips. Compared to the 3D case, the 2D equiaxed dendrites are more asymmetrical and lack secondary arms because of the thicker solute envelope. In the 3D case, the columnar dendrites on the upstream side (left one) are promoted, while the middle and downstream dendrites are inhibited in sequence. However, the sequential inhibition starts on the upstream side in the 2D case. This is mainly because the melt can pass around the upstream branch in 3D space. However, it can only climb over the upstream tip in 2D space. Additionally, the secondary arms show upstream development, which is more significant with increasing inlet velocity. The level of development of the secondary arms is also affected by the decay of the forced flow in the flow direction.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Experimental Validation of Numerical Simulations of a New-Generation NiTi Endodontic File Under Bending

The superelasticity of NiTi shape memory alloy has been used in endodontics since the 1990s. To study the mechanical behavior of endodontic instruments, a traditional approach consists in experimental investigations. However, finite element analysis constitutes another way to assess their mechanical behavior and to facilitate their design. The main aim of this study is to compare experimental and numerical bending results on different structures (NiTi wire, spreader, and instruments) to estimate the reliability of the finite element simulations. These investigations were carried out as follows. Firstly, experimental material parameters identification was performed using NiTi wires. These parameters were implemented in an appropriate NiTi model. Bending was numerically applied to the meshed structures generated by the finite element method. Experimental tests were performed on real structures under bending, with exactly the same loading, in order to compare experimental and numerical results. These results were in good agreement for each of the considered structures. This enabled the validation of the simulation results and the use of simulations to design new endodontic instruments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An $\bar{F}$ Meshfree Treatment for Nearly Incompressible Materials

The computational modeling of nearly incompressible materials is a difficult task for many numerical methods, and even after several decades of investigation, it is still an active research area. This report seeks to address the treatment of incompressible materials in meshfree methods using a synergistic combination of two treatments. The first treatment is an $\bar{F}$ method, where the decomposed dilatational and deviatoric parts are calculated over different smoothing domains. The second treatment “activates” additional nodes throughout the domain to increase the flexibility of the model. We implement this synergistic combination in the context of the reproducing kernel particle method (RKPM) and present results for the Cook’s membrane benchmark problem. The results are compared with those using the composite tet10 finite element with a volume-averaged J formulation. We show that the combined treatment is an effective way to deal with nearly incompressible materials in a meshfree framework and compares well with other highly-effective treatments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗