Search NASA⌕ Search

SEARCH · Search NASA

Results for “discretization methods”

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,657 records · Page 92

Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper

A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.

36 MATERIALS SCIENCE↗

Measurements of Heat Transfer, Flow, and Pressures in a Simulated Turbine Blade Internal Cooling Passage

An experimental study was made to obtain quantitative information on heat transfer, flow, and pressure distribution in a branched duct test section that had several significant features of an internal cooling passage of a turbine blade. The objective of this study was to generate a set of experimental data that could be used for validation of computer codes that would be used to model internal cooling. Surface heat transfer coefficients and entrance flow conditions were measured at nominal entrance Reynolds numbers of 45,000, 335,000, and 726,000. Heat transfer data were obtained by using a steady-state technique in which an Inconel heater sheet is attached to the surface and coated with liquid crystals. Visual and quantitative flow-field data from particle image velocimetry measurements for a plane at midchannel height for a Reynolds number of 45,000 were also obtained. The flow was seeded with polystyrene particles and illuminated by a laser light sheet. Pressure distribution measurements were made both on the surface with discrete holes and in the flow field with a total pressure probe. The flow-field measurements yielded flow-field velocities at selected locations. A relatively new method, pressure sensitive paint, was also used to measure surface pressure distribution. The pressure paint data obtained at Reynolds numbers of 335,000 and 726,000 compared well with the more standard method of measuring pressures by using discrete holes.

Russell, Louis M.↗

Robust Multigrid Smoothers for Three Dimensional Elliptic Equations with Strong Anisotropies

We discuss the behavior of several plane relaxation methods as multigrid smoothers for the solution of a discrete anisotropic elliptic model problem on cell-centered grids. The methods compared are plane Jacobi with damping, plane Jacobi with partial damping, plane Gauss-Seidel, plane zebra Gauss-Seidel, and line Gauss-Seidel. Based on numerical experiments and local mode analysis, we compare the smoothing factor of the different methods in the presence of strong anisotropies. A four-color Gauss-Seidel method is found to have the best numerical and architectural properties of the methods considered in the present work. Although alternating direction plane relaxation schemes are simpler and more robust than other approaches, they are not currently used in industrial and production codes because they require the solution of a two-dimensional problem for each plane in each direction. We verify the theoretical predictions of Thole and Trottenberg that an exact solution of each plane is not necessary and that a single two-dimensional multigrid cycle gives the same result as an exact solution, in much less execution time. Parallelization of the two-dimensional multigrid cycles, the kernel of the three-dimensional implicit solver, is also discussed. Alternating-plane smoothers are found to be highly efficient multigrid smoothers for anisotropic elliptic problems.

Llorente, Ignacio M.↗

Structure-preserving numerical discretizations for domains with boundaries

This SAND report documents Exploratory Express LDRD Project 223790, "Structure-preserving numerical discretizations for domains with boundaries", which developed a method to incorporate consistent treatment of domain boundaries and arbitrary boundary conditions in discrete exterior calculus (DEC) for arbitrary polygonal (2D) and tensor-product structure prism (3D) grids. The new DEC required the development of novel discrete exterior derivatives, boundary operators, wedge products and Hodge stars. This was accomplished through the use of boundary extension and the blending of known 2D operators on the interior with 1D operators on the boundary. The Hodge star was based on the Voronoi Hodge star, and retained the limitation of a triangular circumcentric primal or dual grid along with low-order accuracy. In addition to the new DEC, two related software packages were written: one for the study of DEC operators on arbitrary polygonal and polyhedral grids using both symbolic and numerical approaches and one for a (thermal) shallow water testbed using TRiSK-type numerics. Immediately relevant (already funded, through CANGA) followup work is the development of a high-order, geometrically flexible Hodge star and structure-preserving, high-order, oscillation-limiting transport operators (using WENO) for n-forms on arbitrary 2D and 3D grids. This will provide all of the machinery required for a high-order version of TRiSK with boundaries on arbitrary 2D and tensor-product 3D grids, which is applicable to both the atmospheric (CRM in E3SM-MMF) and oceanic (MPAS-O) components of E3SM.

97 MATHEMATICS AND COMPUTING↗

Trajectory Planning by Preserving Flexibility: Metrics and Analysis

In order to support traffic management functions, such as mitigating traffic complexity, ground and airborne systems may benefit from preserving or optimizing trajectory flexibility. To help support this hypothesis trajectory flexibility metrics have been defined in previous work to represent the trajectory robustness and adaptability to the risk of violating safety and traffic management constraints. In this paper these metrics are instantiated in the case of planning a trajectory with the heading degree of freedom. A metric estimation method is presented based on simplifying assumptions, namely discrete time and heading maneuvers. A case is analyzed to demonstrate the estimation method and its use in trajectory planning in a situation involving meeting a time constraint and avoiding loss of separation with nearby traffic. The case involves comparing path-stretch trajectories, in terms of adaptability and robustness along each, deduced from a map of estimated flexibility metrics over the solution space. The case demonstrated anecdotally that preserving flexibility may result in enhancing certain factors that contribute to traffic complexity, namely reducing proximity and confrontation.

Idris, Husni R.↗

A sixth order Mehrstellen scheme with an application to the Method of Local Corrections for the 3D Poisson equation

We present a sixth order finite difference scheme for Poisson’s equation when discretized with the compact 27-point stencil based on Mehrstellen corrections of the forcing function term f. Our approach results in a sixth order accurate solution error as opposed to a fourth-order error imposed by the classical Mehrstellen correction for the 19-point and 27-point stencils. The present study is a continuation of former work of Spotz and Carey (1996) on compact finite difference schemes for Poisson’s equation where sixth order convergence may be obtained under the assumption that the fourth order derivatives of f are determined analytically. Specifically, we show that sixth order convergence can still be attained when only values of f at grid points are available. The sixth order Mehrstellen scheme is further coupled with a Method of Local Corrections (MLC) 3D Poisson solver improving to sixth order accuracy the results reported in Kavouklis and Colella (2019). The MLC test case considered involves an adaptive grid that comprises 7.5 billion cells.

97 MATHEMATICS AND COMPUTING↗

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

Discontinuous Galerkin↗

Domain decomposition for the SPN solver MINOS

In this article we present a domain decomposition method for the mixed SPN equations, discretized with Raviart-Thomas-Nedelec finite elements. This domain decomposition is based on the iterative Schwarz algorithm with Robin interface conditions to handle communications. After having described this method, we give details on how to optimize the convergence. Finally, we give some numerical results computed in a realistic 3D domain. The computations are done with the MINOS solver of the APOLLO3 (R) code. (authors)

97 MATHEMATICS AND COMPUTING↗

The Activity-tracking paradigm in discrete-event modeling and simulation: The case of spatially continuous distributed systems

From a modeling and simulation perspective, studying dynamic systems consists of focusing on changes in states. According to the precision of state changes, generic algorithms can be developed to track the activity of sub-systems. This paper aims at describing and applying this more natural and intuitive way to describe and implement dynamic systems. Activity is defined mathematically. A generic application case of diffusion is experimented with to compare the efficiency of quantized state methods using this new approach with traditional methods which do not focus computations on active areas. Our goal is to demonstrate that the concept of activity can estimate the computational effort required by a quantized state method. Specifically, when properly designed, a discrete-event simulator for such a method achieves a reduction in the number of state transitions that more than compensates for the overhead it imposes.

99 GENERAL AND MISCELLANEOUS↗

Perfectly Matched Layers versus discrete transparent boundary conditions in quantum device simulations

Discrete transparent boundary conditions (DTBC) and the Perfectly Matched Layers (PML) method for the realization of open boundary conditions in quantum device simulations are compared, based on the stationary and time-dependent Schrödinger equation. The comparison includes scattering state, wave packet, and transient scattering state simulations in one and two space dimensions. The Schrödinger equation is discretized by a second-order Crank–Nicolson method in case of DTBC. For the discretization with PML, symmetric second-, fourth-, and sixth-order spatial approximations as well as Crank–Nicolson and classical Runge–Kutta time-integration methods are employed. In two space dimensions, a ring-shaped quantum waveguide device is simulated in the stationary and transient regime. As an application, a simulation of the Aharonov–Bohm effect in this device is performed, showing the excitation of bound states localized in the ring region. The numerical experiments show that the results obtained from PML are comparable to those obtained using DTBC, while keeping the high numerical efficiency and flexibility as well as the ease of implementation of the former method. -- Highlights: •In-depth comparison between discrete transparent boundary conditions (DTBC) and PML. •First 2-D transient scattering state simulations using DTBC. •First 2-D transient scattering state simulations of the Aharonov–Bohm effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical Schemes for Rough Parabolic Equations

This paper is devoted to the study of numerical approximation schemes for a class of parabolic equations on (0,1) perturbed by a non-linear rough signal. It is the continuation of Deya (Electron. J. Probab. 16:1489-1518, 2011) and Deya et al. (Probab. Theory Relat. Fields, to appear), where the existence and uniqueness of a solution has been established. The approach combines rough paths methods with standard considerations on discretizing stochastic PDEs. The results apply to a geometric 2-rough path, which covers the case of the multidimensional fractional Brownian motion with Hurst index H>1/3.

97 MATHEMATICS AND COMPUTING↗

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗

Large Deviations for Level Sets of a Branching Brownian Motion and Gaussian Free Fields

We study deviation probabilities for the number of high positioned particles in branching Brownian motion and confirm a conjecture of Derrida and Shi. We also solve the corresponding problem for the two-dimensional discrete Gaussian free field. Our method relies on an elementary inequality for inhomogeneous Galton–Watson processes.

97 MATHEMATICS AND COMPUTING↗

In-Situ SEM High Strain Rate Testing of Large Diameter Micropillars Followed by TEM and EBSD Postmortem Analysis

Background Dislocation dynamic simulations are intended as a tool to understand and predict the mechanical behavior of metallic materials, but its prediction has never been directly verified by experiments due to differences in specimen strain rate and size. Objective In this work, a comprehensive experimental framework is proposed to attempt direct comparison between experiments and discrete dislocation dynamics (DDD) modelling. Methods By integrating high-throughput sample fabrication and a customized testing apparatus, the sample size and strain rate typically employed in DDD simulations are explored experimentally. Constitutive properties such as stress-strain response are measured, and microstructural information is obtained from transmission electron microscopy (TEM) imaging, electron backscatter diffraction (EBSD), and TEM-based orientation mapping. Results Magnesium and copper were selected, as case studies, to demonstrate the newly developed experimental procedure. Measured stress-strain responses for Mg are consistent with those obtained with a miniaturized Hopkison bar experiments. Furthermore, by exploiting the validated workflow, the effect of strain rate on micropillar heterogeneous deformation and associated dislocation plasticity were revealed. Conclusion The work establishes a methodology for the systematic study of not only metals but also other materials and structures at the microscale and high strain rates.

42 ENGINEERING↗

Eulerian-Eulerian large eddy simulation of two-phase dilute bubbly flows

An efficient mathematical model is presented for conducting large eddy simulation (LES) of turbulent bubbly flows at low gas hold-ups. The mathematical model is based on the filtered multi-fluid equations derived by application of the component-weighted volume-averaging process of each phase. A two-fluid model is applied to the filtered set of equations, yielding a more computationally tractable set of filtered equations. The filtered equations are formulated in generalized curvilinear coordinates to simulate turbulent gas-liquid bubbly flows in complex geometries. The Germano model is used to account for sub-grid turbulent stresses. The gas and liquid momentum equations are discretized with the hybrid staggered/non-staggered method. The performance of the developed mathematical model is demonstrated by simulating the benchmark test case of a turbulent gas-liquid bubbly flow in a rectangular bubble column. The accuracy of the proposed method is assessed via comparison with the experimental data and other numerical results in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A low-rank solver for the stochastic unsteady Navier–Stokes problem

Here we study a low-rank iterative solver for the unsteady Navier–Stokes equations for incompressible flows with a stochastic viscosity. The equations are discretized using the stochastic Galerkin method, and we consider an all-at-once formulation where the algebraic systems at all the time steps are collected and solved simultaneously. The problem is linearized with Picard’s method. To efficiently solve the linear systems at each step, we use low-rank tensor representations within the Krylov subspace method, which leads to significant reductions in storage requirements and computational costs. Combined with effective mean-based preconditioners and the idea of inexact solve, we show that only a small number of linear iterations are needed at each Picard step. The proposed algorithm is tested with a model of flow in a two-dimensional symmetric step domain with different settings to demonstrate the computational efficiency.

97 MATHEMATICS AND COMPUTING↗

Selection of suitable alternatives to reduce the environmental impact of road traffic noise using a fuzzy multi-criteria decision model

Road traffic noise is one of the most significant environmental impacts generated by transport systems. To this regard, the recent implementation of the European Environmental Noise Directive by Public Administrations of the European Union member countries has led to various noise action plans (NAPs) for reducing the noise exposure of EU inhabitants. Every country or administration is responsible for applying criteria based on their own experience or expert knowledge, but there is no regulated process for the prioritization of technical measures within these plans. This paper proposes a multi-criteria decision methodology for the selection of suitable alternatives against traffic noise in each of the road stretches included in the NAPs. The methodology first defines the main criteria and alternatives to be considered. Secondly, it determines the relative weights for the criteria and sub-criteria using the fuzzy extended analytical hierarchy process as applied to the results from an expert panel, thereby allowing expert knowledge to be captured in an automated way. A final step comprises the use of discrete multi-criteria analysis methods such as weighted sum, ELECTRE and TOPSIS, to rank the alternatives by suitability. To illustrate an application of the proposed methodology, this paper describes its implementation in a complex real case study: the selection of optimal technical solutions against traffic noise in the top priority road stretch included in the revision of the NAP of the regional road network in the province of Almeria (Spain).

54 ENVIRONMENTAL SCIENCES↗