Search NASA⌕ Search

SEARCH · Search NASA

Results for “methods - statistical”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠− 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK DG, collocation, and CG methods.

Numerical Methods for Ordinary Differential Equati↗
DOE OSTI2017

Simulating Bubble Deformations Using the Intersection Marker (ISM) Front-Tracking method

Surface and volume tracking methods have been under continuous development over the last 50+ years since the first MAC (Marker and Cell) and SLIC (Simple Line Interface Calculation) methods appeared in the 1960's. The sustained interest in simulating multiphase flow stem from their many practical applications in the physical world; from biological systems, to chemical systems, to boiling dynamics and atmospheric heat and mass transfer. The ability to capture the details of interface evolution as well as the physics of interfacial heat and mass transfer has thus been a tempting goal for many numerical researchers. Numerical analysts, apart from the well-established and popular VOF PLIC (Volume of Fluid, Piecewise Linear Interface Construction) method, can now call upon even more computationally intensive options such as the LS (Level Set) method, the FT (Front Tracking) method, the MoF (Moment of Fluid) method and hybrid methods.

42 ENGINEERING↗
DOE OSTI2021

A minimally invasive, efficient method for propagation of full-field uncertainty in solid dynamics

In this work, we present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

42 ENGINEERING↗
DOE OSTI2022

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2023

Yet another parameter-free shape optimization method

The use of node coordinates as design variables in shape optimization offers a larger design space than computer-aided design (CAD)-based shape parameterizations. It also allows for the optimization of legacy designs, i.e., a finite element mesh from an existing design can be readily optimized to meet new performance requirements without involving a CAD model. However, it is well known that the node coordinate parameterization method is fraught with numerical difficulties, which makes it impractical to use. This has led to several of “parameter-free” shape optimization methods that seek the advantages and avoid the pitfalls of the naïve node coordinate parameterization method. These methods come in two main varieties: sensitivity filtering (or gradient smoothing) and consistent filtering. The latter is analogous to the density filter method used in topology optimization (TO). In this work, we use the PDE filter from TO and energy-based filters to implement consistent shape optimization filtering schemes easily and efficiently. Numerical experiments demonstrate that consistent methods are more robust than sensitivity filtering methods.

42 ENGINEERING↗
DOE OSTI2017

Comprehensive Evaluation of Fast-Response, Reynolds-Averaged Navier–Stokes, and Large-Eddy Simulation Methods Against High-Spatial-Resolution Wind-Tunnel Data in Step-Down Street Canyons

Three computational fluid dynamics (CFD) methods with different levels of flow-physics modelling are comprehensively evaluated against high-spatial-resolution wind-tunnel velocity data from step-down street canyons (i.e., a short building downwind of a tall building). The first method is a semi-empirical fast-response approach using the Quick Urban Industrial Complex (QUIC-URB) model. The second method solves the Reynolds-averaged Navier–Stokes (RANS) equations, and the third one utilizes a fully-coupled fluid-structure interaction large-eddy simulation (LES) model with a grid-turbulence inflow generator. Unlike typical point-by-point evaluation comparisons, here the entire two-dimensional wind-tunnel dataset is used to evaluate the dynamics of dominant flow topological features in the street canyon. Each CFD method is scrutinized for several geometric configurations by varying the downwind-to-upwind building-height ratio (H d /H u ) and street canyon-width to building-width aspect ratio (S / W) for inflow winds perpendicular to the upwind building front face. Disparities between the numerical results and experimental data are quantified in terms of their ability to capture flow topological features for different geometric configurations. Ultimately, all three methods qualitatively predict the primary flow topological features, including a saddle point and a primary vortex. But, the secondary flow topological features, namely an in-canyon separation point and secondary vortices, are only well represented by the LES method despite its failure for taller downwind building cases. Misrepresentation of flow-regime transitions, exaggeration of the coherence of recirculation zones and wake fields, and overestimation of downwards vertical velocity into the canyon are the main defects in QUIC-URB, RANS and LES results, respectively. All three methods underestimate the updrafts and, surprisingly, QUIC-URB outperforms RANS for the streamwise velocity component, while RANS is superior to QUIC-URB for the vertical velocity component in the street canyon.

58 GEOSCIENCES↗
DOE OSTI2017

Quantifying and reducing the differences in forest CO 2 -fluxes estimated by eddy covariance, biometric and chamber methods: A global synthesis

Carbon dioxide (CO 2 ) fluxes between terrestrial ecosystems and the atmosphere are primarily measured with eddy covariance (EC), biometric, and chamber methods. However, it is unclear why the estimates of CO 2 -fluxes, when measured using these different methods, converge at some sites but diverge at others. We synthesized a novel global dataset of forest CO 2 -fluxes to evaluate the consistency between EC and biometric or chamber methods for quantifying CO 2 budget in forests. The EC approach, comparing with the other two methods, tended to produce 25% higher estimate of net ecosystem production (NEP, 0.52Mg C ha-1 yr-1), mainly resulting from lower EC-estimated Re; 10% lower ecosystem respiration (Re, 1.39Mg C ha-1 yr-1); and 3% lower gross primary production (0.48 Mg C ha-1 yr-1) The discrepancies between EC and the other methods were higher at sites with complex topography and dense canopies versus those with flat topography and open canopies. Forest age also influenced the discrepancy through the change of leaf area index. The open-path EC system induced >50% of the discrepancy in NEP, presumably due to its surface heating effect. These results provided strong evidence that EC produces biased estimates of NEP and Re in forest ecosystems. A global extrapolation suggested that the discrepancies in CO 2 fluxes between methods were consistent with a global underestimation of Re, and overestimation of NEP, by the EC method. Accounting for these discrepancies would substantially improve the our estimates of the terrestrial carbon budget .

Wang, Xingchang↗
DOE OSTI2019

The effect of functional unit and co-product handling methods on life cycle assessment of an algal biorefinery

Microalgae are a promising feedstock for sustainable fuel and nutritional products. To be economically viable, algal biorefineries will need to produce multiple products, such as nutraceuticals and animal feed, in addition to fuel. Life cycle assessment (LCA) has been applied widely to algal biorefineries, but results vary across studies. Here we address two methodological challenges that can lead to variability in LCAs of multi-product algal biorefineries: the choice of functional unit and the method used to handle co-products. LCA was conducted for an algal biorefinery that produces two representative products—fuel and animal feed—using three functional units: (1) 1 MJ fuel, (2) 1 kg animal feed, and (3) 1 ha of production area, coupled with the system expansion method. For the fuel-based functional unit, the system expansion method was compared to impact allocation based on mass, energy, and market value of the biorefinery products. The choice of functional unit affects the results for several environmental indicators. For ecosystem quality, the choice of functional unit can change the overall balance of environmental impacts from harmful to beneficial. The effect of functional unit on climate change depends on the boundary considered in the model—i.e., “well to gate” or “well to wheels”. When a fuel-based functional unit is used, the choice of co-product handling method (i.e., allocation versus the system expansion method) affects environmental impacts, such as ecosystem quality and climate change, differently. For LCAs that use a fuel-based functional unit to compare algal biofuel to other fuels, multiple co-product handling methods should be analyzed and presented to stakeholders. A land based functional unit coupled with the system expansion method can be used to improve LCAs that compare among multi-product algal biorefineries.

09 BIOMASS FUELS↗
DOE OSTI2019

Spatially dependent embedded self-shielding method for nonuniform temperature distribution

The embedded self-shielding method (ESSM) is simple and robust, and applicable to any geometrical configuration to account for resonance self-shielding by performing fixed source transport calculations. ESSM is very similar to the subgroup method, which has been successful in the MPACT neutronics code developed by CASL. The subgroup method introduces errors in a least-squares fitting process that is used to generate subgroup quadrature sets and cannot be applied to thermal energy groups, as ESSM introduces error for spatially dependent self-shielded cross sections inside the fuel. This study focused on integrating the conventional ESSM with an improved Stoker-Weiss method to obtain spatially dependent resonance self-shielded cross sections for a nonuniform fuel temperature distribution. Here, the proposed method would serve as a substitute for the subgroup method in high fidelity simulations using CASL’s VERA. Benchmark results show that this new, spatially dependent ESSM (SDESSM) is an appropriate candidate to replace the subgroup method.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A systematic comparison of machine learning methods for modeling of dynamic processes applied to combustion emission rate modeling

Ten established, data-driven dynamic algorithms are surveyed and a practical guide for understanding these methods generated. Existing Python programming packages for implementing each algorithm are acknowledged, and the model equations necessary for prediction are presented. A case study on a coal-fired power plant’s NO x emission rates is performed, directly comparing each modeling method’s performance on a mutual system. Each model is evaluated by its root mean squared error (RMSE) on out-of-sample future horizon predictions. Optimal hyperparameters are identified using either an exhaustive search or genetic algorithm. The top five model structures of each method are used to recursively predict future NO x emission rates over a 60-step time horizon. The RMSE at each future timestep is determined, and the recursive output prediction trends compared against measurements in time. The GRU neural network is identified as the best candidate for representing the system, demonstrating accurate and stable predictions across the future horizon by all considered models, while satisfactory performance was observed in several of the ARX/NARX formulations. Finally, these efforts have contributed 1) a concise resource of multiple proven dynamic machine learning methods, 2) a practical guide explaining the use of these methods, effectively lowering the “barrier-to-entry” of deploying such models in control systems, 3) a comparison study evaluating each method’s performance on a mutual system, 4) demonstration of accurate multi-timestep emissions modeling suitable for systems-level control, and 5) generalizable results demonstrating the suitability of each method for prediction over a multi-step future horizon to other complex dynamic systems.

42 ENGINEERING↗
DOE OSTI2023

Generalized moving least squares vs. radial basis function finite difference methods for approximating surface derivatives

Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. Further, the second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The Adjoint Petrov–Galerkin method for non-linear model reduction

Here, we formulate a new projection-based reduced-order modeling technique for non-linear dynamical systems. The proposed technique, which we refer to as the Adjoint Petrov–Galerkin (APG) method, is derived by decomposing the generalized coordinates of a dynamical system into a resolved coarse-scale set and an unresolved fine-scale set. A Markovian finite memory assumption within the Mori–Zwanzig formalism is then used to develop a reduced-order representation of the coarse scales. This procedure leads to a closed reduced-order model that displays commonalities with the adjoint stabilization method used in finite elements. The formulation is shown to be equivalent to a Petrov–Galerkin method with a non-linear, time-varying test basis, thus sharing some similarities with the Least-Squares Petrov–Galerkin method. Theoretical analysis examining a priori error bounds and computational cost is presented. Numerical experiments on the compressible Navier–Stokes equations demonstrate that the proposed method can lead to improvements in numerical accuracy, robustness, and computational efficiency over the Galerkin method on problems of practical interest. Improvements in numerical accuracy and computational efficiency over the Least-Squares Petrov–Galerkin method are observed in most cases.

42 ENGINEERING↗
DOE OSTI2021

Sparse invariant domain preserving discontinuous Galerkin methods with subcell convex limiting

In this paper, we develop high-order nodal discontinuous Galerkin (DG) methods for hyperbolic conservation laws that satisfy invariant domain preserving properties using subcell flux corrections and convex limiting. These methods are based on a subcell flux corrected transport (FCT) methodology that involves blending a high-order target scheme with a robust, low-order invariant domain preserving method that is obtained using a graph viscosity technique. Furthermore, the new low-order discretizations are based on sparse stencils which do not increase with the polynomial degree of the high-order DG method. As a result, the accuracy of the low-order method does not degrade when used with high-order target methods. The method is applied to both scalar conservation laws, for which the discrete maximum principle is naturally enforced, and to systems of conservation laws such as the Euler equations, for which positivity of density and a minimum principle for specific entropy are enforced. Numerical results are presented on a number of benchmark test cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Direct Discontinuous Galerkin methods for the reacting multi-component flow equations

The Direct Discontinuous Galerkin (DDG (Liu and Yan, 2008)) method and a counterpart with Interface Correction (DDGIC (Danis and Yan, 2022)) are extended to compute diffusion terms that arise when solving the compressible multi-component flow equations in thermochemical nonequilibrium. Thermodynamic properties, transport properties, chemical reaction rates, and energy exchange terms are computed using Mutation++ (Scoggins et al., 2020). The DG method is applied on unstructured grids, where the accuracy and convergence rates can be sensitive to the numerical method chosen for parabolic terms. A method for determining the homogeneity tensor of the flow equations required for DDGIC is shown. The convergence properties of the DDG methods are studied and compared to the Interior Penalty (IP) method. A number of numerical experiments are conducted to assess the accuracy and performance of the method. The numerical results and convergence studies indicate that DDG and DDGIC provide accurate solutions and perform well for general flows in thermochemical nonequilibrium.

Diffusion↗
DOE OSTI2017

Method for estimating the social sustainability of infrastructure projects

Nowadays, sustainability assessments tend to focus on the biophysical and economic considerations of the built environment. Social facets are generally underestimated when investment in infrastructure projects is appraised. This paper proposes a method to estimate the contribution of infrastructure projects to social sustainability. This method takes into account the interactions of an infrastructure with its environment in terms of the potential for short and long-term social improvement. The method is structured in five stages: (1) social improvement criteria and goals to be taken into account are identified and weighed; (2) an exploratory study is conducted to determine transfer functions; (3) each criterion is homogenized through value functions; (4) the short and long-term social improvement indices are established; and finally, (5) social improvement indices are contrasted to identify the socially selected alternatives and to assign an order of priority. The method was implemented in six alternatives for road infrastructure improvement. The results of the analysis show that the method can distinguish the contribution to social sustainability of different infrastructure projects and location contexts, according to early benefits and potential long-term equitable improvement. This method can be applied prior to the implementation of a project and can complement environmental and economic sustainability assessments. - Highlights: •Method to select suitable infrastructure projects from the social sustainability point of view •Emphasizes social interactions of the infrastructure in the short and long term •Distinguishes the social sustainability of infrastructure projects in different locations •Efficiency of a social contribution in terms of early social benefits and a long-term distribution •Supports early decision-making of public agencies regarding infrastructure projects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Image phase shift invariance based cloud motion displacement vector calculation method for ultra-short-term solar PV power forecasting

Irradiance received on the earth's surface is the main factor that affects the output power of solar PV plants, and is chiefly determined by the cloud distribution seen in a ground-based sky image at the corresponding moment in time. It is the foundation for those linear extrapolation-based ultra-short-term solar PV power forecasting approaches to obtain the cloud distribution in future sky images from the accurate calculation of cloud motion displacement vectors (CMDVs) by using historical sky images. Theoretically, the CMDV can be obtained from the coordinate of the peak pulse calculated from a Fourier phase correlation theory (FPCT) method through the frequency domain information of sky images. The peak pulse is significant and unique only when the cloud deformation between two consecutive sky images is slight enough, which is likely possible for a very short time interval (such as 1?min or shorter) with common changes in the speed of cloud. Sometimes, there will be more than one pulse with similar values when the deformation of the clouds between two consecutive sky images is comparatively obvious under fast changing cloud speeds. This would probably lead to significant errors if the CMDVs were still only obtained from the single coordinate of the peak value pulse. However, the deformation estimation of clouds between two images and its influence on FPCT-based CMDV calculations are terrifically complex and difficult because the motion of clouds is complicated to describe and model. Therefore, to improve the accuracy and reliability under these circumstances in a simple manner, an image-phase-shift-invariance (IPSI) based CMDV calculation method using FPCT is proposed for minute time scale solar power forecasting. First, multiple different CMDVs are calculated from the corresponding consecutive images pairs obtained through different synchronous rotation angles compared to the original images by using the FPCT method. Second, the final CMDV is generated from all of the calculated CMDVs through a centroid iteration strategy based on its density and distance distribution. Third, the influence of different rotation angle resolution on the final CMDV is analyzed as a means of parameter estimation. Simulations under various scenarios including both thick and thin clouds conditions indicated that the proposed IPSI-based CMDV calculation method using FPCT is more accurate and reliable than the original FPCT method, optimal flow (OF) method, and particle image velocimetry (PIV) method.

14 SOLAR ENERGY↗
DOE OSTI2021

An Experimental Study of Uncertainty Considerations Associated with Predicting Auto-ignition Timing using Livengood-Wu Integral Method

The application of the Livengood-Wu (LW) integral method as a tool to estimate knock onset in spark ignited (SI) engines and combustion phasing in advanced compression ignition (ACI) engines has been demonstrated through simulations several times. In this study, the effect of uncertainties associated with parameters required for the LW integral method, when used as a tool for model based control of ignition timing in an ACI engine, were experimentally studied using five full boiling range gasoline fuels. As a first step, the method was applied to experimental data from a rapid compression machine and it was found that the ability of the LW integral method to predict ignition timing was very sensitive to the performance of the chemical kinetic model of each fuel. The method was subsequently applied to experimental data from a single-cylinder gasoline engine with simple approximations for the LW integral input parameters, and it was found that the predicted time of ignition was significantly different from the actual start of combustion. Systematic evaluation of various parametric uncertainties conducted thereafter showed that the uncertainty in cylinder charge temperature has the greatest influence. Improved methods of estimating cylinder charge temperature are proposed to account for the previously determined corrections, to enable the use of the LW integral method for model based control of ignition timing.

Compression Ignition↗