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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Accelerated finite-volume calculation of transonic potential flows

A fully conservative finite-volume algorithm is used to calculate transonic potential flows past isolated airfoils and through two-dimensional channels. The difference equations are solved by a multigrid technique which uses an alternating direction implicit method as a smoothing algorithm. The finite-volume formulation provides a framework within which to treat flows past complicated geometries, while the multigrid/alternating direction scheme provides rapid convergence of the solution to very small residuals.

Jameson, A.↗
DOE OSTI2022

Toward a resolution of the NN controversy

Lattice QCD calculations of two-nucleon interactions have been underway for about a decade, but still haven't reached the pion mass regime necessary for matching onto effective field theories and extrapolating to the physical point. Furthermore, results from different methods, including the use of the Lüscher formalism with different types of operators, as well as the HALQCD potential method, do not agree even qualitatively at very heavy pion mass. We investigate the role that different operators employed in the literature may play on the extraction of spectra for use within the Lüscher method. We first explore expectations from Effective Field Theory solved within a finite volume, for which the exact spectrum may be computed given different physical scenarios. We then present preliminary lattice QCD results for two-nucleon spectra calculated using different operators on a common lattice ensemble.

Nicholson, A↗
DOE OSTI2016

Predicting the flow stress of Zircaloy-4 under in-reactor accident conditions - 2016-0075

During loss-of-coolant accident (LOCA) and reactivity-initiated accident (RIA), nuclear fuel rods experience high heating rates that change the microstructure and properties of zirconium cladding materials. We sought to determine how different fast heating rates affect the kinetics of transformation and whether their consequences on the material yield stresses in the dual-phase region could be predicted. Phase fraction and texture of cold-rolled Zircaloy-4 were thus measured at fast heating rates similar to those experienced during LOCA and RIA by in situ high-energy synchrotron X-ray diffraction; the yield stress was also measured at different heating rates to the same dual-phase temperature by the electrical resistance method. A crystal plasticity finite-element model was used to simulate the flow stresses with the inputs of the measured phase fraction and texture. The model was calibrated using the flow stress at the heating rate of 10 deg. C s{sup -1} and used to predict the flow stress at the heating rate of 50 deg. C s{sup -1} at the same temperature. A significant shift of the β transus temperature was identified for fast heating rates in the range of 10 to 100 deg. C s{sup -1}. There was great similarity in texture at different heating rates. In the dual-phase temperature region, the α texture was almost identical to the starting texture, and the b texture was nearly random. Yield stresses of cold-rolled Zircaloy-4 decreased at a faster heating rate at the same temperature of 920 deg. C. A good agreement was found between the experimental and predicted flow stress at different heating rates, indicating this methodology can potentially be used to predict flow behavior in a fast transient regime. (authors)

36 MATERIALS SCIENCE↗
NASA NTRS1983

Finite-element methods for steady solidification problems

Four Galerkin finite-element methods are tested for solving the free-boundary problem that describes steady solidification. The formulations differ in the solution method used to account for the unknown shape of the melt/solid interface, in the interphase condition (either balance of heat flux or equilibrium of temperature) distinguished for locating the interface, and in the technique used for solving the systems of algebraic equations that result from the finite-element approximations. Methods that use the melting point isotherm to locate the melt/solid interface are found more accurate and efficient than formulations based on the interfacial energy balance. Solution by a Galerkin-Newton algorithm of the free-boundary problem transformed to a fixed domain is most efficient when the field problem in each phase is made nonlinear by including radiation from the melt and solid to the surroundings.

Ettouney, H. M.↗
NASA NTRS1981

A numerical simulation of two-dimensional inlet flow fields

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two-dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time-marching finite-difference procedure to solve the thin-layer Navier-Stokes equations formulated in body-fitted coordinates. Because the method can be used for a flow field with both subsonic and supersonic regions, it is applicable to subcritical as well as supercritical inlet operation.

Biringen, S.↗
DOE OSTI2019

Diffraction of waves by a perfect electromagnetic conductor half-plane between isorefractive media

In this study, diffraction of waves by a perfectly electromagnetic conductor half-plane which is lying at the boundary of isorefractive media is investigated. The diffracted field expressions are constructed by using the relation between scattered geometrical optics fields and the transition boundaries. The high frequency asymptotic expressions of diffracted fields reach infinite magnitude values at the transition boundaries. The uniform expressions of fields are obtained by using the uniform theory of diffraction method in order to obtain the finite magnitude values at the transition boundaries and the results are analyzed numerically for different set of parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2016

The Space-Time Conservative Schemes for Large-Scale, Time-Accurate Flow Simulations with Tetrahedral Meshes

Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.

Venkatachari, Balaji Shankar↗
DOE OSTI2016

Uncertainty quantification for generalized Langevin dynamics

Here, we present efficient finite difference estimators for goal-oriented sensitivity indices with applications to the generalized Langevin equation (GLE). In particular, we apply these estimators to analyze an extended variable formulation of the GLE where other well known sensitivity analysis techniques such as the likelihood ratio method are not applicable to key parameters of interest. These easily implemented estimators are formed by coupling the nominal and perturbed dynamics appearing in the finite difference through a common driving noise or common random path. After developing a general framework for variance reduction via coupling, we demonstrate the optimality of the common random path coupling in the sense that it produces a minimal variance surrogate for the difference estimator relative to sampling dynamics driven by independent paths. In order to build intuition for the common random path coupling, we evaluate the efficiency of the proposed estimators for a comprehensive set of examples of interest in particle dynamics. These reduced variance difference estimators are also a useful tool for performing global sensitivity analysis and for investigating non-local perturbations of parameters, such as increasing the number of Prony modes active in an extended variable GLE.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Modeling and Simulating Dislocation Dynamics Near Sound Speeds in Cubic Crystals

The motion of crystal dislocations near sound speeds is investigated using atomic-scale molecular dynamics (MD) simulations and continuum equations of motion. A material’s deformation rate is influenced by whether dislocations can propagate in the high-velocity high-stress regime. In the femto-second resolved MD simulations, a shear stress is applied to a pair of partial edge dislocations in copper at T = 25 K, and the pair’s average position is tracked with time. Analysis for the steady-state motion reveals that the dislocations can move faster than at least two transverse sound speeds depending on the applied stress. On the modeling side, the dislocation glide equations govern the evolution of the local displacement field according to material elastic constants. A second-order accurate finite difference scheme is applied to solve the screw dislocation equation in the continuum limit for aluminum. This method proved difficult to implement due to numerical instabilities surrounding the dislocation core, so alternative solving schemes may be needed. Combining an external acceleration with the equation of motion should reveal whether the strain energy density remains finite at transonic and supersonic glide speeds.

36 MATERIALS SCIENCE↗
NASA NTRS1979

Micromechanics of intraply hybrid composites: Elastic and thermal properties

Composite micromechanics are used to derive equations for predicting the elastic and thermal properties of unidirectional intraply hybrid composites. The results predicted using these equations are compared with those predicted using approximate equations based on the rule of mixtures, linear laminate theory, finite element analysis and limited experimental data. The comparisons for three different intraply hybrids indicate that all four methods predict approximately the same elastic properties and are in good agreement with measured data. The micromechanics equations and linear laminate theory predict about the same values for thermal expansion coefficients. The micromechanics equations predict through-the-thickness properties which are in good agreement with the finite element results.

Chamis, C. C.↗
NASA NTRS1997

Solution of Aeroacoustic Problems by a Nonlinear, Hybrid Method

Category 1, problem 3 (scattering of sound by a sphere) and category 2, problem 1 (spherical source in a cylindrical duct subject to uniform flow) are solved in generalized coordinates using the nonlinear Euler equations together with nonreflecting boundary conditions. A temporally and spatially fourth-order accurate finite-difference, Runge-Kutta time-marching technique is employed for the near-field calculations and a Kirchhoff method is employed for the prediction of far-field sound. Computations are all performed on parallel processors using the data-parallel paradigm.

Oezyoeruek, Yusuf↗
DOE OSTI2015

Fourth-order partial differential equation noise removal on welding images

Partial differential equation (PDE) has become one of the important topics in mathematics and is widely used in various fields. It can be used for image denoising in the image analysis field. In this paper, a fourth-order PDE is discussed and implemented as a denoising method on digital images. The fourth-order PDE is solved computationally using finite difference approach and then implemented on a set of digital radiographic images with welding defects. The performance of the discretized model is evaluated using Peak Signal to Noise Ratio (PSNR). Simulation is carried out on the discretized model on different level of Gaussian noise in order to get the maximum PSNR value. The convergence criteria chosen to determine the number of iterations required is measured based on the highest PSNR value. Results obtained show that the fourth-order PDE model produced promising results as an image denoising tool compared with median filter.

42 ENGINEERING↗
DOE OSTI2024

Adaptive mesh refinement in binary black holes simulations

Abstract We discuss refinement criteria for the Berger–Rigoutsos (block-based) refinement algorithm in our numerical relativity code GR-Athena++ in the context of binary black hole (BBH) merger simulations. We compare three different strategies: the ‘box-in-box’ approach, the ‘sphere-in-sphere’ approach and a local criterion for refinement based on the estimation of truncation error of the finite difference scheme. We extract and compare gravitational waveforms using the three different mesh refinement methods and compare their accuracy against a calibration waveform and demonstrate that the sphere-in-sphere approach provides the best strategy overall when considering computational cost and the waveform accuracy. Ultimately, we demonstrate the capability of each mesh refinement method in accurately simulating gravitational waves from BBH systems—a crucial aspect for their application in next-generation detectors. We quantify the mismatch achievable with the different strategies by extrapolating the gravitational wave mismatch to higher resolution.

Astronomy & Astrophysics↗
DOE OSTI2016

TH-CD-BRA-09: Towards Absolute Dose Measurement in MRI-Linac and Gamma-Knife: Design and Construction of An MR-Compatible Water Calorimeter

Purpose: The purpose of this work is to design and optimize a portable Water Calorimeter (WC) for use in a commercial MRI-linac and Gamma-knife in addition to conventional radiotherapy linacs. Water calorimeters determine absorbed dose to water at a point by measuring radiation-induced temperature rise of the volume (the two are related by the medium specific heat capacity). In this formalism, one important correction factor is heat transfer correction k-ht. It compensates for heat gain/loss due to conductive and convective effects, and is numerically calculated as ratio of temperature rise in the absence of heat loss to that in the presence of heat loss. Operating at 4°C ensures convection is minimal. Methods: A commercial finite element software was used to evaluate several WC designs with different insulation materials and thicknesses; channels allowing coolant to travel around WC (to sustain WC at 4°C) were modeled, and worst-case scenario variation in the temperature of the coolant was simulated for optimization purposes (2.6 mK/s). Additionally, several calorimeter vessel design parameters (front/back glass thickness/separation, diameter) were also simulated and optimized. Optimization is based on minimizing long term calorimeter drift (24h) as well as variation and magnitude of k-ht. Results: The final selected WC design reached a modest drift of 11µK/s after 15h for the worst-case coolant temperature variation. This design consists of coolant channels being encompassed on both sides by cryogel insulation. For the MRI-linac beam, glass thickness plays the largest effect on k-ht with variation of upto 0.6% in the first run for thicknesses ranging between 0.5–1.7mm. Subsequent runs vary only within 0.1% with glass thickness. Other factors such as vessel radius and top/bottom glass separation have sub 0.1% effects on k-ht. Conclusion: An MR-safe 4°C stagnant WC appropriate for dosimetry in MRI-linac and Gamma-Knife was designed, optimized, and construction is nearly completed. NSERC Discovery Grant RGPIN-435608.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2018

Modeling underwater noise propagation from marine hydrokinetic power devices through a time-domain, velocity-pressure solution

Marine hydrokinetic (MHK) devices generate electricity from the motion of tidal and ocean currents, as well as ocean waves, to provide an additional source of renewable energy available to the United States. These devices are a source of anthropogenic noise in the marine ecosystem and must meet regulatory guidelines that mandate a maximum amount of noise that may be generated. In the absence of measured levels from in situ deployments, a model for predicting the propagation of sound from an array of MHK sources in a real environment is essential. A set of coupled, linearized velocity-pressure equations in the time-domain are derived and presented in this paper, which are an alternative solution to the Helmholtz and wave equation methods traditionally employed. Discretizing these equations on a three-dimensional (3D), finite-difference grid ultimately permits a finite number of complex sources and spatially varying sound speeds, bathymetry, and bed composition. The solution to this system of equations has been parallelized in an acoustic-wave propagation package developed at Sandia National Labs, called Paracousti. This work presents the broadband sound pressure levels from a single source in two-dimensional (2D) ideal and Pekeris wave-guides and in a 3D domain with a sloping boundary. Furthermore, the paper concludes with demonstration of Paracousti for an array of MHK sources in a simple wave-guide.

13 HYDRO ENERGY↗
DOE OSTI2017

Thermally Effective and Efficient Cooling Technologies for Advanced Gas Turbine Systems

This report presents the research which has been undertaken under DOE grant DE-FE-0011875 as part of NETL’s UTSR (University Turbine System Research) Program. The experimental portion of this work was undertaken at the University of North Dakota (UND) under the direction of Dr. Forrest Ames. The computational work was undertaken at LSU (Louisiana State University), UM (University of Memphis), and at IIT (Illinois Institute of Technology) under the direction of Dr. Sumanta Acharya. The research program was organized into a four phase program which included: 1) development of a FDA/FEA (finite difference analysis/finite element analysis) cooling model, 2) experimental/computational investigation of appropriate internal cooling methods, 3) experimental/computational investigation of external heat transfer and film cooling, 4) development of an enhanced FDA/FEA model along with the development of a warm cascade test for comparison purposes. The experimental work at UND included research on both internal cooling methods and external heat transfer. The internal cooling work included the design and testing of multiple variable hole-size “incremental impingement” configurations along with testing of two converging rounded diamond pedestal arrays for the trailing edge region. The external heat transfer work included the building of a large scale cascade facility to test a vane design with a large leading edge and an aft loaded suction surface. UND has acquired both midspan and full suction surface heat transfer measurements as well as prepared for film cooling measurements under the grant. LSU/UM/IIT has developed and conducted wall resolved LES computations for the incremental impingement geometries and provided rational understanding of the effect of hole sizes and their distributions. Additionally they have provided design correlations for the incremental impingement heat transfer coefficients. They have also developed LES methods to predict flow through film cooling plenums, which have pedestals along with the development of slot film cooling downstream from the plenums. Additionally, they have developed wall resolved LES to predict vane surface heat transfer including bypass transition without the use of a transition model. Through these simulations they provide understanding of the aerodynamic and surface heat transfer behavior under a range of influences such as freestream turbulence. The present study has produced some very useful research outcomes.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2017

Magnus: A New Resistive MHD Code with Heat Flow Terms

We present a new magnetohydrodynamic (MHD) code for the simulation of wave propagation in the solar atmosphere, under the effects of electrical resistivity—but not dominant—and heat transference in a uniform 3D grid. The code is based on the finite-volume method combined with the HLLE and HLLC approximate Riemann solvers, which use different slope limiters like MINMOD, MC, and WENO5. In order to control the growth of the divergence of the magnetic field, due to numerical errors, we apply the Flux Constrained Transport method, which is described in detail to understand how the resistive terms are included in the algorithm. In our results, it is verified that this method preserves the divergence of the magnetic fields within the machine round-off error (∼1×10{sup −12}). For the validation of the accuracy and efficiency of the schemes implemented in the code, we present some numerical tests in 1D and 2D for the ideal MHD. Later, we show one test for the resistivity in a magnetic reconnection process and one for the thermal conduction, where the temperature is advected by the magnetic field lines. Moreover, we display two numerical problems associated with the MHD wave propagation. The first one corresponds to a 3D evolution of a vertical velocity pulse at the photosphere–transition–corona region, while the second one consists of a 2D simulation of a transverse velocity pulse in a coronal loop.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Development of a three dimensional numerical water quality model for continental shelf applications

A model to predict the distribution of water quality parameters in three dimensions was developed. The mass transport equation was solved using a non-dimensional vertical axis and an alternating-direction-implicit finite difference technique. The reaction kinetics of the constituents were incorporated into a matrix method which permits computation of the interactions of multiple constituents. Methods for the computation of dispersion coefficients and coliform bacteria decay rates were determined. Numerical investigations of dispersive and dissipative effects showed that the three-dimensional model performs as predicted by analysis of simpler cases. The model was then applied to a two dimensional vertically averaged tidal dynamics model for the Providence River. It was also extended to a steady state application by replacing the time step with an iteration sequence. This modification was verified by comparison to analytical solutions and applied to a river confluence situation.

Spaulding, M.↗