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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2023

An ICA-Based HVAC Load Disaggregation Method Using Smart Meter Data

This paper presents an independent component analysis (ICA) based unsupervised-learning method for heat, ventilation, and air-conditioning (HVAC) load disaggregation using row-resolution (i.e., 15 minutes) smart meter data. We first demonstrate that the electricity consumption profiles on mild-temperature days can be used to approximate the base load on hot days. A residual load profile can then be calculated by subtracting the mild-day load profile from the hot-day load profile. The residual load profiles are processed using ICA for HVAC load extraction. An optimization-based algorithm is proposed for post-adjustment of the ICA results, considering two bounding factors for enhancing the robustness of the ICA algorithm. First, we use the hourly HVAC energy bounds computed from the relationship between HVAC load and temperature to remove unrealistic HVAC load spikes. Second, we exploit the dependency between the daily nocturnal and diurnal loads extracted from historical meter data to smooth the base load profile. Pecan Street data with sub-metered HVAC data were used to test and verify the proposed methods. Simulation results demonstrated that the proposed method is computationally efficient and robust across multiple customers.

Kim, Hyeonjin↗

DOE OSTI2018

Numerical Analysis on Hydraulic Power Take-Off for Wave Energy Converter and Power Smoothing Methods: Preprint

One of the primary challenges for wave energy converter (WEC) systems is the fluctuating nature of wave resources, which require the WEC components to be designed to handle loads (i.e., torques, forces, and powers) that are many times greater than the average load. This approach requires a much greater power take-off (PTO) capacity than the average power output and indicates a higher cost for the PTO. Moreover, additional design requirements, such as battery storage, are needed, particularly for practical electrical grid connection, and can be a problem for sensitive equipment (e.g., radar, computing devices, and sensors). Therefore, it is essential to investigate potential methodologies to reduce the overall power fluctuation while trying to optimize the power output from WECs. In this study, a detailed hydraulic PTO model was developed and coupled with a time-domain hydrodynamics model (WEC-Sim) to evaluate the PTO efficiency for WECs and the trade-off between power output and fluctuation using different power smoothing methods, including energy storage, pressure relief mechanism, and a powerbased setpoint control method. The study also revealed that the maximum power fluctuation for WECs can be significantly reduced by one order of magnitude when these power smoothing methods are applied.

16 TIDAL AND WAVE POWER↗

DOE OSTI2026

Data for Propagation Method and Planting Density Influence Canopy Developmental Transition and Biomass Productivity in Miscanthus × giganteus

Understanding how establishment practices influence the mechanisms underlying Miscanthus × giganteus (miscanthus) productivity and canopy development is critical for optimizing management. Data was collected during the juvenile (2011–2013) and mature (2024) phases of a long-term field experiment established in Urbana, Illinois, to evaluate the effects of propagation method (plug propagation [PP] and rhizome propagation [RP]), planting density (1.0, 0.75, and 0.25 plants m⁻²), and nitrogen application (0 and 67 kg N ha⁻¹) on end-of-season biomass yield, tiller mass, tiller density, and tiller height. Linear regression models identified the dominant predictors of yield across stand ages and management regimes. Planting density, nitrogen (N) application, and propagation method significantly influenced early yield and canopy development. During the juvenile phase, biomass yield was driven by tiller density due to canopy expansion; in the mature phase, yield became driven by tiller mass. The PP plots produced higher tiller density than the RP plots, resulting in faster canopy closure and higher juvenile-phase yields. Rhizome-propagated (RP) plots produced lower tiller density, but individual tillers were 3.3–6.4 g tiller−1 heavier than PP tillers. After the canopy reached equilibrium, the PP and RP yields were similar because greater RP tiller mass compensated for its lower tiller density. Higher planting density resulted in greater yield and tiller density during the second year (2012), but this effect was absent from the third year (2013) onward. In the juvenile phase, N fertilization enhanced yield by 1.6–3.4 Mg ha−1. Initiating fertilization in 2013 on unfertilized plots produced biomass similar to that in fertilized plots, suggesting yield recovery in the mature phase. These findings revealed that establishment strategies, including propagation method and planting density, influence juvenile miscanthus canopy development and productivity, transitioning from tiller-density- to mass-dominated yields, but not mature phase productivity.

Miscanthus↗

DOE OSTI2018

Innovative thermal distillation method using solar heat localization

The thermal distillation uses energy to evaporate water and subsequently condense the vapor to liquid form. Waste energy from industrial processes, electricity, or solar energy are the examples of an input energy for thermal distillation systems. In this paper, we have designed, fabricated, and tested an eco-friendly and low- cost distillation system which uses solar heat localization for evaporation of water and later condense the vapor back to liquid as liquid water. The unique feature of this method is the high rate of evaporation when compared to other thermal solar distillation/desalination systems. In this method, the incoming heat is confined to the porous media where bottom layer acts a porous insulator. The capillary force in the porous structure pumps up the liquid water to the hot surface, within the porous medium, where liquid heats up and changes to vapor phase. Since there is no moving part in this method, and due to its independency to electricity, it a great candidate for deployment in remote areas that access to grid power is limited and maintenance is not feasible.

Jaladi, Divya↗

DOE OSTI2023

Universal Passive Synchronization Method for Grid-Forming Inverters Without Mode Transition

Power systems are transforming with increasing levels of inverter-based resources (IBRs). This transformation requires critical roles of grid-forming (GFM) inverters replacing synchronous generators for bulk power system stabilization and ancillary services, also allowing flexible power system operation, such as microgrid that is operated by multiple GFM IBRs to achieve system resilience against contingencies. To realize the resilient power systems allowing flexible in-and-out operation of GFM IBRs potentially programmed with different primary controls, a synchronization method universally applicable, i.e., independent of control types, would be beneficial to ease the integration process, but it has not been actively studied. To fill the gap, this paper proposes a universal synchronization method that achieves a passive synchronization to enable a smooth transition in a grid with off-nominal system parameters, i.e., voltage and frequency. The logic proposed requires no modification on the primary control, thus applicable to any type of GFMs with a voltage reference input. To validate the concept, a simulation of an IEEE 13-bus benchmark system modified with 3 GFM inverters is presented. It simulates an inverter-driven black start scenario in which GFM inverters autonomously turn on and connect to the grid under heavy loading, using the synchronization logic. The case study demonstrates that GFM inverters can tune their voltage reference to smoothly synchronize without severe transients, and contribute to a seamless black start of the grid under unbalanced load conditions. Two GFM methods - Droop and dispatchable virtual oscillator control - are used for the demo to validate feasibility and interoperability of the passive synchronization.

black start↗

NASA NTRS1974

Fracture control methods for space vehicles. Volume 1: Fracture control design methods

A systematic approach for applying methods for fracture control in the structural components of space vehicles consists of four major steps. The first step is to define the primary load-carrying structural elements and the type of load, environment, and design stress levels acting upon them. The second step is to identify the potential fracture-critical parts by means of a selection logic flow diagram. The third step is to evaluate the safe-life and fail-safe capabilities of the specified part. The last step in the sequence is to apply the control procedures that will prevent damage to the fracture-critical parts. The fracture control methods discussed include fatigue design and analysis methods, methods for preventing crack-like defects, fracture mechanics analysis methods, and nondestructive evaluation methods. An example problem is presented for evaluation of the safe-crack-growth capability of the space shuttle crew compartment skin structure.

Liu, A. F.↗

NASA NTRS1982

Equivalence of the generalized Lie-Hori method and the method of averaging

In this investigation, a comparison is made of two methods for developing perturbation theories for non-canonical dynamical systems. The methods compared are the generalized Lie-Hori method and the method of averaging. In the comparison presented here, the equivalence of the methods up to the second order in the small parameter is shown. However, the approach used can be extended to demonstrate the equivalence for higher orders. To illustrate the equivalence both Duffing's equation and the van der Pol equation are solved using each method.

Ahmed, A.↗

NASA NTRS1984

Equivalence of the generalized Lie-Hori method and the method of averaging

In this investigation, a comparison is made of two methods for developing perturbation theories for non-canonical dynamical systems. The methods compared are the generalized Lie-Hori method and the method of averaging. In the comparison presented here, the equivalence of the methods up to the second order in the small parameter is shown. However, the approach used can be extended to demonstrate the equivalence for higher orders. To illustrate the equivalence Duffing's equation, the van der Pol equation and the oscillator with quadratic damping problem are solved using each method.

Ahmed, A. H.↗

NASA NTRS1993

A new method of time difference measurement: The time difference method by dual phase coincidence points detection

In the high accurate measurement of periodic signals, the greatest common factor frequency and its characteristics have special functions. A method of time difference measurement - the time difference method by dual 'phase coincidence points' detection is described. This method utilizes the characteristics of the greatest common factor frequency to measure time or phase difference between periodic signals. It can suit a very wide frequency range. Measurement precision and potential accuracy of several picoseconds were demonstrated with this new method. The instrument based on this method is very simple, and the demand for the common oscillator is low. This method and instrument can be used widely.

Zhou, Wei↗

NASA NTRS2002

A Multifunctional Interface Method for Coupling Finite Element and Finite Difference Methods: Two-Dimensional Scalar-Field Problems

A multifunctional interface method with capabilities for variable-fidelity modeling and multiple method analysis is presented. The methodology provides an effective capability by which domains with diverse idealizations can be modeled independently to exploit the advantages of one approach over another. The multifunctional method is used to couple independently discretized subdomains, and it is used to couple the finite element and the finite difference methods. The method is based on a weighted residual variational method and is presented for two-dimensional scalar-field problems. A verification test problem and a benchmark application are presented, and the computational implications are discussed.

Ransom, Jonathan B.↗

NASA NTRS2000

New digital measurement methods for left ventricular volume using real-time three-dimensional echocardiography: comparison with electromagnetic flow method and magnetic resonance imaging

AIM: The aim of this study was to investigate the feasibility and accuracy of using symmetrically rotated apical long axis planes for the determination of left ventricular (LV) volumes with real-time three-dimensional echocardiography (3DE). METHODS AND RESULTS: Real-time 3DE was performed in six sheep during 24 haemodynamic conditions with electromagnetic flow measurements (EM), and in 29 patients with magnetic resonance imaging measurements (MRI). LV volumes were calculated by Simpson's rule with five 3DE methods (i.e. apical biplane, four-plane, six-plane, nine-plane (in which the angle between each long axis plane was 90 degrees, 45 degrees, 30 degrees or 20 degrees, respectively) and standard short axis views (SAX)). Real-time 3DE correlated well with EM for LV stroke volumes in animals (r=0.68-0.95) and with MRI for absolute volumes in patients (r-values=0.93-0.98). However, agreement between MRI and apical nine-plane, six-plane, and SAX methods in patients was better than those with apical four-plane and bi-plane methods (mean difference = -15, -18, -13, vs. -31 and -48 ml for end-diastolic volume, respectively, P<0.05). CONCLUSION: Apically rotated measurement methods of real-time 3DE correlated well with reference standards for calculating LV volumes. Balancing accuracy and required time for these LV volume measurements, the apical six-plane method is recommended for clinical use.

Evaluation Studies↗

NASA NTRS1987

Atmospheric H2O2 measurement: comparison of cold trap method with impinger bubbling method

Collection of atmospheric H2O2 was performed by a cold trap method using dry ice-acetone as the refrigerant. The air was drawn by a pump into a glass gas trap immersed in the dry ice-acetone slush in a dewar flask at a flow rate of 2.5 l min-1 for approximately 2 h. Collection efficiency was > 99% and negligible interferences by O3, SO2 or organic matter with the collected H2O2 in the trap were observed. This method was compared with the air impinger bubbling method which has been previously described (Kok et al., 1978a, b, Envir. Sci. Technol. 12, 1072-1080). The measured total peroxide (H2O2 + organic peroxide) values in a series of aim samples collected by the impinger bubbling method (0.06-3.7 ppb) were always higher than those obtained by the cold trap method (0.02-1.2 ppb). Laboratory experiments suggest that the difference in values between the two methods probably results from the aqueous phase generation of H2O2 and organic peroxide in the impinger solution by a reaction of atmospheric O3 with olefinic and aromatic compounds. If these O3-organic compound reactions which occur in the impinger also occur in aqueous droplets in the atmosphere, the process could be very important for aqueous phase generation of H2O2 in clouds and rainwater.

NASA Discipline Exobiology↗

DOE OSTI2022

Comparison of spherical harmonics method and discrete ordinates method for radiative transfer in a turbulent jet flame

Here, in this study, we systematically compared the accuracy and computational cost of two popular solution methods for the radiative transfer equation (RTE): the spherical harmonics method (P N ) and the discrete ordinates method (DOM). We first investigated convergence characteristics of different orders of P N and DOM in a series of 1D homogeneous configurations with varying optical thicknesses. Both solvers perform better for optically thicker cases. The accuracy of P N methods increases with its order, , but the gain in accuracy reduces with the increase in , i.e., improvement of P 7 over P 5 is less than that of P 3 over P 1 . This decreasing trend becomes more prominent as the optical thickness decreases. On the other hand, DOM’s accuracy increases almost linearly with the increase in the number of ordinates (or polar angles in this study) in all cases. While comparing the directional profile of radiative intensity, both solvers perform better when the radiative intensity is more isotropic. These solvers were then connected with a full spectrum k-distribution (FSK) spectral model and used to perform radiation-coupled simulations of a turbulent jet flame in an axi-symmetric cylindrical domain. Results obtained from P 1 to P 7 approximations for P N , and 2 x 4, 4 x 4, 4 x 8, 8 x 8 finite angles for DOM are compared with that from an optically thin model, and a reference solution from line-by-line (LBL) photon Monte Carlo (PMC) method. The choice of radiation solver shows a noticeable impact on the temperature distribution of the flame. The P N solvers lead to slightly higher radiant fractions and the DOM solvers lead to slightly lower radiant fractions than the PMC benchmark solution. Finally, the computational costs of each of these solvers are also reported and an intermittent evaluation / time blending scheme to improve the computational efficiency of radiation solvers in radiation-coupled simulations are also demonstrated.

42 ENGINEERING↗

DOE OSTI2023

The Random Ray Method Versus Multigroup Monte Carlo: The Method of Characteristics in OpenMC and SCONE

The Random Ray Method (TRRM) is a recently developed approach to solving neutral particle transport problems based on the Method of Characteristics. While the method previously has been implemented only in closed-source or limited-functionality codes, this work describes its implementation in two open-source Monte Carlo codes: OpenMC and SCONE. The random ray implementations required small modifications to the existing Multigroup Monte Carlo (MGMC) solvers, offering a rare venue for redundant, fine-grained, "apples-to-apples" speed and accuracy comparisons between transport methods. To this end, TRRM and MGMC solvers are evaluated against each other using each code's native capabilities on reactor eigenvalue problems with different degrees of energy discretization. On the C5G7 benchmark (featuring only seven energy groups), TRRM achieves a maximum pin power error comparable to or lower than that of MGMC for a given run time. On a problem with 69 energy groups, MGMC is found to scale more efficiently, obtaining a lower pin power error for a given run time. However, the defining difference between the two transport methods is found to be their vastly different uncertainty distributions. Specifically, TRRM is found to maintain similar levels of accuracy and uncertainty throughout the simulation domain whereas MGMC can exhibit orders-of-magnitude greater errors in areas of the problem that feature low neutron flux. For instance, TRRM provided an up to 373 times speed advantage compared with MGMC for computing the flux in low-flux regions in the moderator surrounding the C5G7 core.

42 ENGINEERING↗

DOE OSTI2018

An Adaptive Reduced Basis Collocation Method Based on PCM ANOVA Decomposition for Anisotropic Stochastic PDEs

The combination of reduced basis and collocation methods enables efficient and accurate evaluation of the solutions to parameterized PDEs. Here, we study the stochastic collocation methods that can be combined with reduced basis methods to solve high-dimensional parameterized stochastic PDEs. We also propose an adaptive algorithm using a probabilistic collocation method (PCM) and ANOVA decomposition. This procedure involves two stages. First, the method employs an ANOVA decomposition to identify the effective dimensions, i.e., subspaces of the parameter space in which the contributions to the solution are larger, and sort the reduced basis solution in a descending order of error.Then, the adaptive search refines the parametric space by increasing the order of polynomials until the algorithm is terminated by a saturation constraint. We demonstrate the effectiveness of the proposed algorithm for solving a stationary stochastic convection-diffusion equation, a benchmark problem chosen because solutions contain steep boundary layers and anisotropic features. Lastly, we show that two stages of adaptivity are critical in a benchmark problem with anisotropic stochasticity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

On the implicit immersed boundary method in coupled discrete element and lattice Boltzmann method

The coupled discrete element method and lattice Boltzmann method (DEMLBM) has increasingly drawn attention of researchers in geomechanics due to its mesoscopic nature since 2000. Immersed boundary method (IBM) and immersed moving boundary (IMB) are two popular schemes for coupling fluid particle in DEMLBM. This work aims at coupling DEM and LBM using the latest IBM algorithm and investigating its accuracy, computational efficiency, and applicability. Two benchmark tests, interstitial fluid flow in an ideal packing and single particle sedimentation in viscous fluid, are carried out to demonstrate the accuracy of IBM through semi-empirical Ergun equation, finite element method (FEM), and IMB. Then, simulations of particle migration with relatively large velocity in Poiseuille flow are utilized to address limitations of IBM in DEMLBM modeling. In addition, advantages and deficiencies of IBM are discussed and compared with IMB. It is found that the accuracy of IBM can be only guaranteed when sufficient boundary points are used and it is not suitable for geomechanical problems involving large fluid or particle velocity.

58 GEOSCIENCES↗

DOE OSTI2017

Stability analysis of nonlinear two-grid method for multigroup neutron diffusion problems

We present theoretical analysis of a nonlinear acceleration method for solving multigroup neutron diffusion problems. This method is formulated with two energy grids that are defined by (i) fine-energy groups structure and (ii) coarse grid with just a single energy group. The coarse-grid equations are derived by averaging of the multigroup diffusion equations over energy. The method uses a nonlinear prolongation operator. We perform stability analysis of iteration algorithms for inhomogeneous (fixed-source) and eigenvalue neutron diffusion problems. To apply Fourier analysis the equations of the method are linearized about solutions of infinite-medium problems. The developed analysis enables us to predict convergence properties of this two-grid method in different types of problems. Numerical results of problems in 2D Cartesian geometry are presented to confirm theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗