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,603 records · Page 89

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗

The evolution history of superalloy powders during hot consolidation and plastic deformation

Highlights: • Element segregations in a powder metallurgy superalloy at different states are identified. • The surface energy drives precipitation of oxides, carbides, and γ' particles on powder surface during hot consolidation. • Prior particle boundary structure is distorted and more uniformly refined grains are obtained after hot extrusion. - Abstract: Powder metallurgy superalloys are widely used to produce critical rotating parts in turbine engines. However, the issues like what is the driving force of precipitation on prior particle boundary (PPB) during consolidation and how the PPB structure evolves during plastic deformation still remain to be resolved. In this study, with the characterization of the microstructure evolution from powder, consolidation, to plastic deformation, the segregation of alloying elements on the particle boundaries at elevated temperature was analyzed by equilibrium thermodynamic calculations, and the distortion of PPB structure under large plastic deformation was quantitively simulated by finite element (FE) methods. The results indicated that surface energy serves as the main driving force of PPB precipitation and hot extrusion could distort the PPB structure and yield more uniformly refined grains.

36 MATERIALS SCIENCE↗

Differential analysis of incompressibility in neutron-rich nuclei

Both the incompressibility K A of a finite nucleus of mass A and that (K ∞ ) of infinite nuclear matter are fundamentally important for many critical issues in nuclear physics and astrophysics. While some consensus has been reached about K ∞ , accurate theoretical predictions and experimental extractions of K τ characterizing the isospin dependence of K A have been very difficult. We propose a differential approach to extract K τ and K ∞ independently from the K A data of any two nuclei in a given isotope chain. Applying this method to the K A data from isoscalar giant monopole resonances (ISGMR) in even-even Pb, Sn, Cd, and Ca isotopes taken by Garg et al. at the Research Center for Nuclear Physics (RCNP), Osaka University, Japan, we find that the 106 Cd– 116 Cd and 112 Sn– 124 Sn pairs having the largest differences in isospin asymmetries in their respective isotope chains measured so far provide consistently the most accurate up-to-date K τ value of K τ = –616 ± 59 MeV and K τ =–623 ± 86 MeV, respectively, largely independent of the remaining uncertainties of the surface and Coulomb terms in expanding K A , while the K ∞ values extracted from different isotopes chains are all well within the current uncertainty range of the community consensus for K ∞ . Moreover, the size and origin of the “soft Sn puzzle” is studied with respect to the “stiff Pb phenomenon.” Furthermore, it is found that the latter is favored due to a much larger (by ≈ 380 MeV) K τ for Pb isotopes than for Sn isotopes, while K ∞ from analyzing the K A data of Sn isotopes is only about 5 MeV less than that from analyzing the Pb data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Surface-Enhanced Raman Scattering (SERS) Studies of Disc-on-Pillar (DOP) Arrays: Contrasting Enhancement Factor with Analytical Performance

The use of nanomachining methods capable of reproducible construction of nano-arrayed devices have revolutionized the field of plasmonic sensing by the introduction of a diversity of rationally engineered designs. Significant strides have been made to fabricate plasmonic platforms with tailored interparticle gaps to improve their performance for surface-enhanced Raman scattering (SERS) applications. Over time, a dichotomy has emerged in the implementation of SERS for analytical applications, the construction of substrates, optimization of interparticle spacing as a means to optimize electromagnetic field enhancement at the localized surface plasmon level, and the substrate sensitivity over extended areas to achieve quantitative performance. Here, we assessed the enhancement factor of plasmonic Ag/SiO 2 /Si disc-on-pillar (DOP) arrays of variable pitch with its analytical performance for quantitative applications. Experimental data were compared with those from finite-difference time-domain (FDTD) simulations used in the optimization of the array dimensions. A self-assembled monolayer (SAM) of benzenethiol rendered highly reproducible signals (RSD ~4–10%) and SERS substrate enhancement factor (SSEF) values in the orders of 10 6 –10 8 for all pitches. Furthermore, spectra corresponding to rhodamine 6G (R6G) and 4-aminobenzoic acid demonstrated the advantages of using the more densely packed DOP arrays with a 160 nm pitch (gap= 40 nm) for quantitation in spite of the strongest SSEF was attained for a pitch of 520 nm corresponding to a 400 nm gap.

36 MATERIALS SCIENCE↗

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗

Probing Surface Plasmon Dynamics in Periodic Nanostructures through Ultrafast Electron Microscopy

Surface plasmon polaritons (SPPs) can be manipulated to localize and guide light in subwavelength distances, enabling them to find applications in a wide range of areas, from sensing to quantum computing. Among several methods of SPP excitation, periodic arrays of nano- and microstructures are of particular interest, as they enable engineering SPP properties through structural parameters. Here, in this study, using the photon-induced near-field electron microscopy (PINEM) technique, we investigated the mode formation, coupling, interference, and decay of SPPs in square and hexagonal arrays of circular nanoholes under both visible and near-infrared excitation. Polarization-resolved analysis revealed the key factors governing SPP localization and interference patterns, showing that the periodicity and symmetry of the array primarily determine the SPP interference patterns and their orientation, while pump polarization mainly modulates their intensity. Time-resolved PINEM measurements demonstrated the spatial dependence of the SPP temporal characteristics. In addition, cathodoluminescence (CL) spectroscopy was employed to examine the intrinsic plasmonic characteristics of the structure. Finite difference time domain (FDTD) simulations showed strong agreement with both PINEM and CL measurements on the spatial and spectral behavior of SPPs. Understanding the spatiotemporal dynamics of SPPs on nanostructures beyond the diffraction limit is crucial for optimizing plasmonic structures for advanced photonic and quantum technologies.

Plasmonics↗

Remarks on mechanical approach to observable Universe

We consider the Universe deep inside the cell of uniformity. At these scales, the Universe is filled with inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies), which perturb the background Friedmann model. Here, the mechanical approach (Eingorn and Zhuk, 2012) is the most appropriate to describe the dynamics of the inhomogeneities which is defined, on the one hand, by gravitational potentials of inhomogeneities and, on the other hand, by the cosmological expansion of the Universe. In this paper, we present additional arguments in favor of this approach. First, we estimate the size of the cell of uniformity. With the help of the standard methods of statistical physics and for the galaxies of the type of the Milky Way and Andromeda, we get that it is of the order of 190 Mpc which is rather close to observations. Then, we show that the nonrelativistic approximation (with respect to the peculiar velocities) is valid for z∼<10, i.e. approximately for 13 billion years from the present moment. We consider scalar perturbations and, within the ΛCDM model, justify the main equations. Moreover, we demonstrate that radiation can be naturally incorporated into our scheme. This emphasizes the viability of our approach. This approach gives a possibility to analyze different cosmological models and compare them with the observable Universe. For example, we indicate some problematic aspects of the spatially flat models. Such models require a rather specific distribution of the inhomogeneities to get a finite potential at any points outside gravitating masses. We also criticize the application of the Schwarzschild-de Sitter solution to the description of the motion of test bodies on the cosmological background.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-dimensional gyrokinetic simulation of the relaxation of a magnetized temperature filament

An electromagnetic, 3D gyrokinetic particle code is used to study the relaxation of a magnetized electron temperature filament embedded in a large, uniform plasma of lower temperature. The study provides insight into the role played by unstable drift-Alfvén waves observed in a basic electron heat transport experiment [D. C. Pace et al., Phys. Plasmas 15, 122304 (2008)] in which anomalous cross-field transport has been documented. The simulation exhibits the early growth of temperature-gradient-driven, drift-Alfvén fluctuations that closely match the eigenmodes predicted by linear theory. At the onset of saturation, the unstable fluctuations display a spiral spatial pattern, similar to that observed in the laboratory, which causes the rearrangement of the temperature profile. After saturation of the linear instability, the system exhibits a markedly different behavior depending on the inclusion in the computation of modes without variation along the magnetic field, i.e., k{sub z} = 0. In their absence, the initial filament evolves into a broadened temperature profile, self-consistent with undamped, finite amplitude drift-Alfvén waves. But the inclusion of k{sub z} = 0 modes causes the destruction of the filament and damping of the drift-Alfvén modes leading to a final state consisting of undamped convective cells and multiple, smaller-scale filaments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of wave damping and finite perpendicular scale on three-dimensional Alfvén wave parametric decay in low-beta plasmas

Shear Alfvén wave parametric decay instability (PDI) provides a potential path toward significant wave dissipation and plasma heating. However, fundamental questions regarding how PDI is excited in a realistic three-dimensional (3D) open system and how the finite perpendicular wave scale—as found in both laboratory and space plasmas—affects the excitation remain poorly understood. Here, we present the first 3D, open-boundary, hybrid kinetic-fluid simulations of kinetic Alfvén wave PDI in low-beta plasmas. Key findings are that the PDI excitation is strongly limited by the wave damping present, including electron–ion collisional damping (represented by a constant resistivity) and geometrical attenuation associated with the finite-scale Alfvén wave, and ion Landau damping of the child acoustic wave. The perpendicular wave scale alone, however, plays no discernible role: waves of different perpendicular scales exhibit similar instability excitation as long as the magnitude of the parallel ponderomotive force remains unchanged. These findings are corroborated by theoretical analysis and estimates. This new understanding of 3D kinetic Alfvén wave PDI physics is essential for laboratory study of the basic plasma process and may also aid future evaluation of the relevance/role of PDI in low-beta space plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flow-based likelihoods for non-Gaussian inference

We investigate the use of data-driven likelihoods to bypass a key assumption made in many scientific analyses, which is that the true likelihood of the data is Gaussian. In particular, we suggest using the optimization targets of flow-based generative models, a class of models that can capture complex distributions by transforming a simple base distribution through layers of nonlinearities. We call these flow-based likelihoods (FBL). We analyze the accuracy and precision of the reconstructed likelihoods on mock Gaussian data, and show that simply gauging the quality of samples drawn from the trained model is not a sufficient indicator that the true likelihood has been learned. We nevertheless demonstrate that the likelihood can be reconstructed to a precision equal to that of sampling error due to a finite sample size. We then apply FBLs to mock weak lensing convergence power spectra, a cosmological observable that is significantly non-Gaussian (NG). We find that the FBL captures the NG signatures in the data extremely well, while other commonly used data-driven likelihoods, such as Gaussian mixture models and independent component analysis, fail to do so. This suggests that works that have found small posterior shifts in NG data with data-driven likelihoods such as these could be underestimating the impact of non-Gaussianity in parameter constraints. By introducing a suite of tests that can capture different levels of NG in the data, we show that the success or failure of traditional data-driven likelihoods can be tied back to the structure of the NG in the data. Here, unlike other methods, the flexibility of the FBL makes it successful at tackling different types of NG simultaneously. Because of this, and consequently their likely applicability across datasets and domains, we encourage their use for inference when sufficient mock data are available for training.

79 ASTRONOMY AND ASTROPHYSICS↗

Nucleon resonance structure in the finite volume of lattice QCD

An approach for relating the nucleon resonances extracted from πN reaction data to lattice QCD calculations has been developed by using the finite-volume Hamiltonian method. Within models of πN reactions, bare states are introduced to parametrize the intrinsic excitations of the nucleon. We show that the resonance can be related to the probability P N* (E) of finding the bare state, N*, in the πN scattering states in infinite volume. We further demonstrate that the probability P V N* (E) of finding the same bare states in the eigenfunctions of the underlying Hamiltonian in finite volume approaches P N* (E) as the volume increases. Our findings suggest that the comparison of P N* (E) and P V N* (E) can be used to examine whether the nucleon resonances extracted from the πN reaction data within the dynamical models are consistent with lattice QCD calculation. We also discuss the measurement of P V N* (E) directly from lattice QCD. Furthermore, the practical differences between our approach and the approach using the Lüscher formalism to relate LQCD calculations to the nucleon resonance poles embedded in the data are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗

Hybrid Optimization Software Suite

The hybrid optimization software suite (HOSS) is a general purpose fully 2D/3D parallel combined finite discrete element (FDEM) code that can be used to simulate problems involving fracture and fragmentation processes, large deformation and large rotations, discrete particle systems, solid-fluid interaction in hydro fracture problems, etc. HOSS uses a hybrid approach that combines finite-element and discrete-element methods with a novel computational fluid dynamics solver. The finite-element method is often used to analyze a material or object and how it responds to stress. The discrete-element method analyzes stresses and displacements in a volume containing a large number of particles, such as grains of sand. Fluid dynamics analyzes the fluid flow inside, around, or through solid domains. With these processes combined, HOSS represents a paradigm shift when it comes to generating accurate simulations of material deformations and failure. HOSS can perform the following: Resolves problems that consist of millions of deforming fracturing, and interacting solid particles; Requires no coupling—it naturally integrates with all regimes of fluid flow; Resolves all regimes of fluid flow, such as Stokes, low/high Reynolds, subsonic/transonic/supersonic/hypersonic, compressible/incompressible, viscid/inviscid, Newtonian/Non-Newtonian and turbulent/laminar flow; Combines all flow regimes in the same solver; Numerically stable for all flow regimes; Naturally combines different flow regimes in the same problem; Possesses material library for fluid behavior; Naturally matches fluid solver to solid solver; Includes non-inertial Eulerian formulation for fluids; Uses highly efficient parallel computing based on the Virtual Parallel Machine, thus enabling easy porting between different parallel architectures, including possible future architectures; Designed specifically for multi-physics problems in research and industry in virtual experimentation format, thereby complementing both theoretical-computational and experimental tools; Can be extended to address millions of deforming solid particles immersed in fluid flowing at all regimes; Can be extended to address biological applications, such as modeling bone fracture or red blood cells in blood plasma

Knight, Earl E.↗

A parametric study on leading-edge slots used on wind turbine airfoils at various angles of attack

Here in the present work, a narrow span-wise rectangular channel (referred to as a slot) is introduced and drilled near the leading edge of a finite-span cambered airfoil to study its impact on the overall aerodynamic performance. This slot starts from the vicinity of the leading edge and exits from the pressure-side of the airfoil after making a slight turn. NACA 4412 is used as the baseline airfoil profile, and the influence of several geometrical parameters of the slot at the chord-based Reynolds number of 1.6 M, and different angles of attack (AoAs) ranging from 0 to 16° on the lift and drag coefficients are investigated. The effect of slot's width (thickness), inlet angle, and vertical position are studied, and the appropriate range of values for each variable is determined. For one of the best cases examined, a lift coefficient improvement as large as 30% is observed, while the drag penalty is insignificant. Computational fluid dynamics (CFD) method is used for modeling and analysis, and the results are validated against both published data and the data from a set of wind-tunnel experiments.

17 WIND ENERGY↗

Evaluation of Nanoscale Accessible Pore Structures for Improved Prediction of Gas Production Potential in Chinese Marine Shales

The Lower Cambrian Niutitang and Lower Silurian Longmaxi shales in the Upper Yangtze Platform (UYP) are the most promising strata for shale gas exploration in China. Knowledge of the nanoscale pore structure may improve the prediction of the gas production potential in Chinese marine shales. A systematic investigation of the pore accessibility and its impact on methane adsorption capacity has been conducted on shale samples using various techniques including geochemical and mineralogical analyses, field-emission scanning electron microscopy (FE-SEM), small-angle neutron scattering (SANS), helium porosimetry, and methane adsorption. The results show that organic matter (OM) pores with various shapes dominate the pore systems of these shales. OM tended to mix with clay minerals and converted to organoclay complexes, developing plentiful micro- and mesopores. A unified fit model with two pore structures, fractal pores and finite pores, was used to model the SANS data to characterize the pore structure of the shales. Both mass and surface fractals are identified for each pore structure. The total porosity estimated by the Porod invariant method ranges between 2.35 and 16.40%, of which the porosity for finite pores ranges between 0.35 and 6.36%, and the porosity for the fractal pores ranges between 2.07 and 8.51%. The fraction of open pores was evaluated by comparing the porosities estimated by He porosimetry and SANS. We find that the fraction of open pores is higher than 64% for most of these shales. Correlation analyses suggest that clay and total organic carbon (TOC) have opposite effects on pore structure and methane adsorption capacity. Samples with higher clay contents have higher pore accessibility and lower total porosity, surface area, and maximum methane adsorption, whereas samples with higher TOC content show the inverse relationships. The high percentage of open pores may reduce methane adsorption capacity in these shales, whereas low pore accessibility may reduce methane production at specific pressure differences. Furthermore, both TOC and pore accessibility may be essential controlling factors in methane production from shale gas reservoirs.

03 NATURAL GAS↗

Emergence of many-body quantum chaos via spontaneous breaking of unitarity

It is suggested that many-body quantum chaos appears as the spontaneous symmetry breaking of unitarity in interacting quantum many-body systems. It has been shown that many-body level statistics, probed by the spectral form factor (SFF) defined as K(η,t)=$\langle$|Tr exp(-ηH+itH)| 2 $\rangle$, is dominated by a diffuson-type mode in a field theory analysis. The key finding of this Letter is that the “unitary” η = 0 case is different from the η → 0 ± limit, with the latter leading to a finite mass of these modes due to interactions. This mass suppresses a rapid exponential ramp in the SFF, which is responsible for the fast emergence of Poisson statistics in the noninteracting case, and gives rise to a nontrivial random matrix structure of many-body levels. The interaction-induced mass in the SFF shares similarities with the dephasing rate in the theory of weak localization and the Lyapunov exponent of the out-of-time-ordered correlators.

36 MATERIALS SCIENCE↗

Study of 3D Pores and Its Relationship with Crack Initiation Factors of Aluminum Alloy Die Castings

Observing different types of pores in aluminum alloy from A356 high-pressure die casting by X-ray computerized tomography and scanning electron microscopy, we found that equivalent spherical diameters of pores obeyed the 3-parameter lognormal distribution rather than the 2-parameter one. Statistically, the sphericity of most round gas pores with smooth inner surfaces consisting of oxide films was > 0.58, while the sphericity of shrinkage pores with unfractured dendrites was < 0.40. Sphericity range for the gas-shrinkage pores, composed of gas pores and shrinkage pores with long tails, was 0.40-0.58. Finite element analysis with simulation of the Von Mises stress distribution around the actual 3D pores among crack generation was performed. Tensile test simulation showed that cracks often originated from the shrinkage pores and the gas-shrinkage pores, rather than from the gas pores. This method also can be used to predict fracture strength of samples with pores. For our samples, results obtained using this method showed that simulated fracture strength was slightly higher than experimentally determined strength.

36 MATERIALS SCIENCE↗

Predicting effective thermal and elastic properties of cementitious composites containing polydispersed hollow and core-shell micro-particles

Micro- and nano-sized hollow and core-shell particles (CSPs) have attracted tremendous interests in developing functional cementitious composites and concretes. In this paper, a general method is developed to predict the thermal and elastic properties of cementitious composites containing core-shell and hollow micro-particles. The model follows a two-stage homogenization process – CSP inclusions together with their surrounding interfacial transition zone (ITZ) are first treated as equivalent solid particles, and then the homogenized properties of the composite system are obtained using numerical (i.e., finite element) or analytical (Mori-Tanaka) approaches. Compared with other homogenization methods, the proposed model is easy-to-use and versatile. The numerical model, validated by experimental results, is then used to guide the design of functional cementitious composites. The results show that CSP particle size (or relative shell thickness), volume fraction, shell property, and the size of the ITZ have significant impacts on the effective thermal and elastic properties of cementitious composites; whereas the particle size distribution pattern has relatively minor influence. In addition, the effective property of the composite is more sensitive to the particle property change when CSP inclusions are much stiffer or more thermally conductive than the matrix. Lastly, a scoping study is conducted to elucidate the relationship between the effective thermal conductivity and effective elastic moduli for cementitious composites containing CSP additives with different shell materials, particle sizes, and volume concentrations etc.

36 MATERIALS SCIENCE↗