Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite size effects”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Real-time dynamics of the Schwinger model as an open quantum system with Neural Density Operators

Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ultrafast valley polarization of graphene nanorings

Here, we study theoretically electron dynamics of a graphene nanoring placed in the field of an ultrashort optical pulse. We describe the graphene nanoring within an effective model with infinite mass boundary conditions. For an optical pulse with a duration of just a few femtoseconds, the electron dynamics is coherent and is described by a time-dependent Schrödinger equation. If the optical pulse is circularly polarized, then two valleys of graphene are populated differently, resulting in a finite valley polarization of the system after the pulse. Such a valley polarization is a unique property of graphene nanoscale systems, while for a graphene monolayer, a circularly polarized pulse does not produce any valley polarization. The valley polarization of the graphene nanoring depends on parameters of the system, such as inner and outer radii. With the system's size increasing, the valley polarization monotonically decreases, converging to its zero value for the infinite graphene monolayer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A micromorphic filter for determining stress and deformation measures from direct numerical simulations of lower length scale behavior

A micromorphic filter is presented for the extraction of stress and deformation measures from underlying direct numerical simulation (DNS) of microstructural mechanical response. The filter is consistent with the micromorphic continuum theory of Eringen and Suhubi (1964), and allows the interrogation of a DNS to generate higher order constitutive models. The filter adopts aspects of the overlap coupling techniques for atomistic-continuum scale-bridging in order to update the nodal degrees of freedom of the finite-element-based filter. Here, two methods for the extraction of quantities of interest from the underlying DNS are presented. The first is a least squares fitting technique which is appropriate for finite-element-style macro-scale domains, and the second is a variationally based approach which generalizes the method to any macro-scale domain which has a variational form. No restrictions on the form of the micro-scale simulation are inherent to the approach, provided that the required quantities (stress, position, etc.) can be provided. The filter is applied to a simple homogeneous medium undergoing homogeneous deformation to demonstrate the recovery of a classical continuum response, what we call the null case. We then apply the framework to the same homogeneous medium with a central and edge penny-shaped crack to show micromorphic effects. We demonstrate that varying the size of the filter causes the micromorphic behaviors to arise and then dissipate as the filter becomes smaller than the crack's length scale. Additionally we show the effect of filter size on a chiral element network along with a plate perforated with evenly spaced and perturbed holes through the thickness. The micromorphic effects are shown to be present for the chiral elements under any loading but become most pronounced when a single element is located within each averaging domain. Micromorphic effects become pronounced for the plate upon the onset of localization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Efficient backward x-ray emission in a finite-length plasma irradiated by a laser pulse of picosecond duration

Motivated by experiments employing picosecond-long, kilojoule laser pulses, we examined x-ray emission in a finite-length underdense plasma irradiated by such a pulse using two-dimensional particle-in-cell simulations. We found that, in addition to the expected forward emission, the plasma also efficiently emits in the backward direction. Our simulations reveal that the backward emission occurs when the laser exits the plasma. The longitudinal plasma electric field generated by the laser at the density down-ramp turns around some of the laser-accelerated electrons and re-accelerates them in the backward direction. As the electrons collide with the laser, they emit hard x rays. The energy conversion efficiency is comparable to that for the forward emission, but the effective source size is smaller. We show that the picosecond laser duration is required for achieving a spatial overlap between the laser and the backward energetic electrons. At peak laser intensity of 1.4×1020 W/cm2, backward-emitted photons (energies above 100 keV and 10° divergence angle) account for 2×10−5 of the incident laser energy. This conversion efficiency is three times higher than that for similarly selected forward-emitted photons. The source size of the backward photons (5 μm) is three times smaller than the source size of the forward photons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Precision study of the massive Schwinger model near quantum criticality

We perform a numerical analysis of the massive Schwinger model in the presence of a background electric field. Using the Density Matrix Renormalization Group approach, we efficiently compute the spectrum of the Schwinger model on a staggered lattice with up to 3000 qubits. As a result, we achieve a precise computation of the critical mass of the massive Schwinger model to five digits using four different “criticality criteria,” observing perfect agreement among them Additionally, we discuss the effect of a four-fermion operator deformation of the Schwinger model and compute the critical mass for various values of the deformation parameter.

Critical phenomena↗
DOE OSTI2020

An embedded 3D fracture modeling approach for simulating fracture-dominated fluid flow and heat transfer in geothermal reservoirs

In this paper, we describe an efficient modeling approach, named embedded discrete fracture method (EDFM), for incorporating arbitrary 3D, discrete fractures, such as hydraulic fractures or faults, into modeling fracture-dominated fluid flow and heat transfer in fractured geothermal reservoirs. This technique allows 3D discrete fractures to be discretized independently from surrounding rock volume and inserted explicitly into a primary fracture/matrix grid, generated without including 3D discrete fractures in prior. An effective computational algorithm is developed to discretize these 3D discrete fractures and construct local connections between 3D fractures and fracture/matrix grid blocks representing the surrounding rock volume. The constructed gridding information on 3D fractures is then added to the primary grid. This embedded fracture modeling approach can be directly implemented into a developed geothermal reservoir simulator via the integral finite difference (IFD) method or with TOUGH2 technology. This embedded fracture modeling approach is very promising and computationally efficient to handle realistic 3D discrete fractures with complicated geometries, connections, and spatial distributions. Compared with other fracture modeling approaches, it avoids cumbersome 3D unstructured, local refining procedures, and increases computational efficiency by simplifying Jacobian matrix size and sparsity, while maintaining enough accuracy. Several numeral simulations are presented to demonstrate the utility and robustness of the proposed technique. Our numerical experiments show that this approach captures all the key patterns about fluid flow and heat transfer dominated by fractures in these cases. Thus, this approach is readily available to the simulation of fractured geothermal reservoirs with both artificial and natural fractures.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Specimen Size Artifacts Associated with a Glovebox Deployable Laser Flash Diffusivity Measurement System

To accelerate nuclear fuel qualification and deployment efforts, Oak Ridge National Laboratory has led the world in developing vehicles for accelerating burnup accumulation in nuclear fuel by testing miniature samples in the High Flux Isotope Reactor. However, a challenge remains to establish parity between post irradiation characterization methods for miniature fuel specimens in comparison to conventional fuel elements. This work aims to identify artifacts associated with the thermal analysis of miniature disk specimens in comparison to samples of standardized dimensions through a combination of experimental and computational tools. Using a newly procured Netzsch 427 laser flash diffusivity system, preliminary data on mild steel specimens were collected, and finite element techniques were used to model heat transfer using 1D conditions to assess issues associated with pulse width and data collection frequency effects. These preliminary investigations found that significant variability in measured thermal diffusivity exist within datasets. These variations are not necessarily attributed to limitations with respect to data collection rate; nor are they associated with limitations on pulse width. Instead, these artifacts may be a result of limitations in detector intensity. Future work aims to design a focusing lens to better measure the smaller surface area of miniature specimens analyzed using the laser flash analysis system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Optimization of overhang printing in ceramic binder jetting additive manufacturing

Cylindrical-shaped horizontally-hanging overhangs with different heights on a solid cubic body were printed via binder jetting additive manufacturing with sililcon carbide powder feedstock. Cracking on the overhang neck and sagging on the overhang body were observed both pre- and post-curing, especially on printed design with a smaller particle size of 10 μm. Here, it was found that the width of cracking and depth of sagging were dependent on the height relative to the solid body in the powder bed, and elevated position of overhang with large volume of underneath loose powder led to increased extend of cracking and sagging. Overhang diameter and weight were varied, which were found with no apparent effects on its deformation extend. To mitigate cracking and sagging, contactless support was designed and tested, which could significantly improve the overhang stability and avoid any overhang deformations. Support gap and support distance were studied for their effects on support capability. The printing processes with and without contactless support were simulated by a finite element modeling approach to understand the stress distribution and evolution of the overhang structures, which revealed stress concentrations on the upper neck point of the designed overhang.

Binder jetting↗
DOE OSTI2025

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗
DOE OSTI2020

Hall viscosity of the composite-fermion Fermi seas for fermions and bosons

The Hall viscosity has been proposed as a topological property of incompressible fractional quantum Hall states and can be evaluated as Berry curvature. This paper reports on the Hall viscosities of composite-fermion Fermi seas at ν = 1/m, where m is even for fermions and odd for bosons. A well-defined value for the Hall viscosity is not obtained by viewing the 1/m composite-fermion Fermi seas as the n → ∞ limit of the Jain ν = n/(nm ± 1) states, whose Hall viscosities (±n + m)ℏρ/4 (ρ is the two-dimensional density) approach ±∞ in the limit n → ∞. A direct calculation shows that the Hall viscosities of the composite-fermion Fermi sea states are finite and also relatively stable with system size variation, although they are not topologically quantized in the entire τ space. Here, I find that the ν = 1/2 composite-fermion Fermi sea wave function for a square torus yields a Hall viscosity that is expected from particle-hole symmetry and is also consistent with the orbital spin of 1/2 for Dirac composite fermions. I compare my numerical results with some theoretical conjectures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Structure and dynamics of hydrodynamically interacting finite-size Brownian particles in a spherical cavity: Spheres and cylinders

The structure and dynamics of confined suspensions of particles of arbitrary shape are of interest in multiple disciplines from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an immersed-boundary general geometry Ewald-like method to capture lubrication and long-range hydrodynamics and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction, and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than 10%, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion toward the walls, and a sub-diffusive regime-caused by crowding-in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗
DOE OSTI2026

On the roles of welding residual stresses in determination of fracture toughness in austenitic stainless steel SUS 304 pipeline girth welds

Welding residual stresses especially the high tensile stresses are proved to have negative impacts on the fatigue and fracture behaviors of welded structures. In this study, a virtual fabrication of test specimens from welding process to specimen preparation was carried out by numerical simulation. An austenitic stainless steel multi-pass pipe welding was simulated by transient thermal–mechanical finite element analysis, the residual stresses were then mapped into the test specimen to evaluate fracture toughness. The findings in this study confirmed that, residual stress can be high in a sub-sized compact tensile specimen, which may accelerate or hinder the crack propagation during actual fatigue and fracture tests as reported in recent years. The influence of the cutting location and orientation of the specimen on fracture performance was investigated systematically to provide a fundamental understanding of welding residual stress and necessary insights into the specimen preparation procedure. Considering the limitation of measuring techniques and the complexity of the stress distribution, the developed numerical model can be a very useful tool to elucidate the stress evolution and quantify the effect of remaining welding stress on fracture toughness.

Fracture behavior↗
DOE OSTI2021

Topological protection of coherence in disordered open quantum systems

Here, we consider topological protection mechanisms in dissipative quantum systems in the presence of quenched disorder, with the intent to prolong the coherence time of a fiducial qubit. The qubit is part of a network of other qubits and dissipative cavities whose coupling parameters are tunable, such that topological edge states can be stabilized. The evolution of the fiducial qubit is entirely determined by a non-Hermitian Hamiltonian which thus emerges from a bona fide physical process. Even in the presence of disorder, a winding number W can be defined and evaluated in real space, as long as certain symmetries are preserved. Hence we can construct the topological phase diagrams of noisy open quantum models, such as the non-Hermitian disordered Su-Schrieffer-Heeger dimer model and a trimer model that includes longer-range couplings. For finite-size systems we find that there are precisely W modes localized at one end of the chain. In such topological phases the qubit's coherence lifetime is exponentially large in the system size. In the presence of competing disorder parameters, interesting reentrance phenomena of topologically nontrivial sectors are observed. This means that in certain parameter regions, increasing disorder drastically increases the coherence time of the fiducial qubit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mechanical analysis of heterogeneous materials with higher-order parameters

Abstract Even though heterogeneous porous materials are widely used in a variety of engineering and scientific fields, such as aerospace, energy-storage technology, and bio-engineering, the relationship between effective material properties of porous materials and their underlying morphology is still not fully understood. To contribute to this knowledge gap, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we use the second-order scheme that is an extension of the first-order computational homogenization framework, where a generalized continuum enables us to introduce length scale into the material constitutive law and capture both pore size and pore distribution. Through several numerical case studies with different combinations of porosity, pore shapes, and distributions, we systematically studied the relationship between the underlying morphology and effective mechanical properties. The results highlight the necessity of higher-order homogenization in understanding the mechanical properties and reveal that higher-order parameters are required to capture the role of realistic pore morphologies on effective mechanical properties. Furthermore, for specific pore shapes, higher-order parameters exhibit dominant influence over the first-order continuum.

42 ENGINEERING↗
DOE OSTI2022

The Study on the Shock Wave Propagation Rule of a Gas Explosion in a Gas Compartment

Combined with the k-ε turbulence model of general application, a refined finite element model of a utility tunnel’s gas compartment filled with the methane/air mixture is developed. A series of analyses are made by using the powerful industry-leading computational fluid dynamics (CFD) software flame acceleration simulator (FLACS) to study the shock wave propagation rule in the gas compartment. The longitudinal and transversal distribution laws of the explosion shock wave are gained taking into consideration the spatial characteristics of the gas compartment. The influences of a few parameters, such as initial conditions and section size of the gas compartment, on the shock wave propagation rule are further discussed. The basic procedure for predicting the peak pressure of the blast wave is provided by considering the initial conditions and the gas compartment, and the corresponding injury effect of the explosion wave on the living beings is assessed. The investigation demonstrates that the peak pressure by the coupled effect between the initial conditions is significantly influenced, especially at the upper and lower gas explosion limits. The peak pressure increases gradually as the width or height increases, and both basically meet the linear relation. The proposed method can forecast the peak pressure of the explosion shock wave in the gas compartment accurately. According to the peak pressure longitudinal and transversal distributions of the blast wave, the peak pressure is far greater than the killing pressure threshold in the underground and closed space; consequently, it is not safe for the living beings in the gas compartment.

Lin, Shu-Chao↗
DOE OSTI2021

Driven-dissipative phase transition in a Kerr oscillator: From semiclassical PT symmetry to quantum fluctuations

We study a minimal model that has a driven-dissipative quantum phase transition, namely a Kerr non-linear oscillator subject to driving and dissipation. Using mean-field theory, exact diagonalization, and the Keldysh formalism, we analyze the critical phenomena in this system, showing which aspects can be captured by each approach and how the approaches complement each other. Then critical scaling and finite-size scaling are calculated analytically using the quantum Langevin equation. Furthermore, the physics contained in this simple model is surprisingly rich: it includes a continuous phase transition, Z 2 symmetry breaking, PT symmetry, state squeezing, and critical fluctuations. Due to its simplicity and solvability, this model can serve as a paradigm for exploration of open quantum many-body physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗