Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Domain decomposition”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Memory-efficient emulation of physical tabular data using quadtree decomposition

Computationally expensive functions are sometimes replaced in simulations with an emulator that approx-imates the true function (e.g., equations of state, wavelength-dependent opacity, or composition-dependent materials properties). For functions that have a constrained domain of interest, this can be done by discretizing the domain and performing a local interpolation on the tabulated function values of each local domain. For these so-called tabular data methods, the method of discretizing the domain and mapping the input space to each subdomain can drastically influence the memory and computational costs of the emulator. This is especially true for functions that vary drastically in different regions. We present a method for domain discretization and mapping that utilizes quadtrees, which results in significant reductions in the size of the emulator with minimal increases to computational costs or loss of global accuracy. We apply our method to the electron-positron Helmholtz free energy equation of state and show over an order of magnitude reduction in memory costs for reasonable levels of numerical accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Tensor factorization for elucidating mechanisms of piezoresponse relaxation via dynamic Piezoresponse Force Spectroscopy

Spatially resolved time and voltage-dependent polarization dynamics in PbTiO 3 thin films is explored using dynamic piezoresponse force microscopy (D-PFM) in conjunction with interferometric displacement sensing. This approach gives rise to 4D data sets containing information on bias-dependent relaxation dynamics at each spatial location without long-range electrostatic artifacts. To interpret these data sets in the absence of defined physical models, we employ a non-negative tensor factorization method which clearly presents the data as a product of simple behaviors allowing for direct physics interpretation. Correspondingly, we perform phase-field modeling finding the existence of ‘hard’ and ‘soft’ domain wall edges. This approach can be extended to other multidimensional spectroscopies for which even exploratory data analysis leads to unsatisfactory results due to many components in the decomposition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Rheology and dynamics of a solvent segregation driven gel (SeedGel)

Bicontinuous structures promise applications in a broad range of research fields, such as energy storage, membrane science, and biomaterials. Kinetically arrested spinodal decomposition is found responsible for stabilizing such structures in different types of materials. A recently developed solvent segregation driven gel (SeedGel) is demonstrated to realize bicontinuous channels thermoreversibly with tunable domain sizes by trapping nanoparticles in a particle domain. As the mechanical properties of SeedGel are very important for its future applications, a model system is characterized by temperature-dependent rheology. The storage modulus shows excellent thermo-reproducibility and interesting temperature dependence with the maximum storage modulus observed at an intermediate temperature range (around 28 °C). SANS measurements are conducted at different temperatures to identify the macroscopic solvent phase separation during the gelation transition, and solvent exchange between solvent and particle domains that is responsible for this behavior. Further, the long-time dynamics of the gel is further studied by X-ray Photon Correlation Spectroscopy (XPCS). The results indicate that particles in the particle domain are in a glassy state and their long-time dynamics are strongly correlated with the temperature dependence of the storage modulus.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Randomized Projection for Rank-Revealing Matrix Factorizations and Low-Rank Approximations

Rank-revealing matrix decompositions provide an essential tool in spectral analysis of matrices, including the Singular Value Decomposition (SVD) and related low-rank approximation techniques. QR with Column Pivoting (QRCP) is usually suitable for these purposes, but it can be much slower than the unpivoted QR algorithm. For large matrices, the difference in performance is due to increased communication between the processor and slow memory, which QRCP needs in order to choose pivots during decomposition. Our main algorithm, Randomized QR with Column Pivoting (RQRCP), uses randomized projection to make pivot decisions from a much smaller sample matrix, which we can construct to reside in a faster level of memory than the original matrix. This technique may be understood as trading vastly reduced communication for a controlled increase in uncertainty during the decision process. Furthermore, for rank-revealing purposes, the selection mechanism in RQRCP produces results that are the same quality as the standard algorithm, but with performance near that of unpivoted QR (often an order of magnitude faster for large matrices). Additionally, we also propose two formulas that facilitate further performance improvements. The first efficiently updates sample matrices to avoid computing new randomized projections. The second avoids large trailing updates during the decomposition in truncated low-rank approximations. Our truncated version of RQRCP also provides a key initial step in our truncated SVD approximation, TUXV. These advances open up a new performance domain for large matrix factorizations that will support efficient problem-solving techniques for challenging applications in science, engineering, and data analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Tensor decompositions for count data that leverage stochastic and deterministic optimization

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the maximum likelihood estimator (MLE) of the Poisson CPD model. Here, this work presents two new algorithms that extend state-of-the-art local methods for Poisson CPD. Hybrid GCP-CPAPR combines Generalized Canonical Decomposition (GCP) with stochastic optimization and CP Alternating Poisson Regression (CPAPR), a deterministic algorithm, to increase the probability of converging to the MLE over either method used alone. Restarted CPAPR with SVDrop uses a heuristic based on the singular values of the CPD model unfoldings to identify convergence toward optimizers that are not the MLE and restarts within the feasible domain of the optimization problem, thus reducing overall computational cost when using a multi-start strategy. We provide empirical evidence that indicates our approaches outperform existing methods with respect to converging to the Poisson CPD MLE.

CPAPR↗
DOE OSTI2024

Unified wavefront singularity characterization of three-dimensional elastodynamic time-domain half-space Green's function under impulsive boundary and internal loads

Founded on a novel analytical formulation that led to a rigorous yet compact path-integral representation of the time-domain elastodynamic half-space Green's function, a unified analysis of the possible occurrence of different singular wavefront behaviour in the response under arbitrary impulsive internal or surface point loads at arbitrary source-receiver locations is presented. With the decomposition of the general solution into distinct initiating and reflected wave group integrals that share a common factored format and simple contour definitions, the mathematical framework is shown to allow a straightforward identification of the specific conditions and the particular wave groups that are responsible for the singular wavefront phenomena without resorting to advanced analytic function theories or asymptotic methods. Analytic characterizations of the nature, strength and direction of all intrinsic singular wavefront behaviours of the three-dimensional Green's function in three canonical cases of source-receiver configurations are given in a dual integral-closed form format to facilitate their theoretical understanding as well as computational applications. Graphical illustrations of their variation with the source-receiver configuration and the medium's Poisson's ratio together with relevant comparison and clarifications of some classical treatments are included.

Green’s functions↗
DOE OSTI2026

BrickQA: Bridging the Semantic Gap in Building Operations with Dynamic Graph Exploration

While standardized ontologies like the Brick schema address data heterogeneity in Building Automation Systems (BAS), accessing this semantic data remains a challenge as domain experts often lack the expertise to formulate complex SPARQL queries. To bridge this gap, we present BrickQA, a Large Language Model (LLM)-based framework that translates natural language into executable SPARQL queries through structured query decomposition, dynamic schema exploration, and inline validation. BrickQA utilizes an iterative reasoning agent to actively navigate graph topology through dynamic exploration actions without requiring exhaustive context injection or model fine-tuning. This approach effectively mitigates hallucinations, particularly in large-scale building knowledge graphs. Empirical evaluation on BuildingQA, a standardized benchmark, demonstrates that BrickQA significantly outperforms ReAct baselines, delivering a 0.291–0.355 absolute F1 improvement while achieving 3 × –12.7 × higher token cost-efficiency. Beyond these metrics, the framework maintains structural fidelity across heterogeneous buildings and remains resilient to ambiguous queries without requiring site-specific fine-tuning. Furthermore, a case study on operational analytics validates the framework’s capability to handle temporal and aggregation constraints, effectively transforming abstract semantic models into actionable facility management insights.1

Ko, Yun-Dam↗
DOE OSTI2022

Distributed Optimization Approaches with Discrete Variables in the Power Distribution Systems

Traditionally, centralized approaches have predominantly been used for the power system operation and control. With increasing penetration of small-scale distributed energy resources (DERs) in the distribution network, especially independently owned renewable resources, distributed algorithms can serve as a potential alternative for improving scalability, resiliency and addressing privacy concerns. However, the complexity of distributed algorithms significantly increases with the integration of the legacy devices, the operation of which depend on discrete control variables. This paper aims to provide a review of the distributed optimization algorithms incorporating discrete control variables for the power distribution system. While the research in this domain is still at its nascence, an extensive comparison of the approaches in the literature for applying quadratic penalty, branch and bound,ordinal optimization and proximal operator to handle discrete variables in the framework of ADMM and dual decomposition have been addressed. Future research direction in this field have been also provided.

Adan, Jannatul↗
DOE OSTI2022

Cyber-Attack Identification of Synchrophasor Data Via VMD and Multifusion SVM

A large amount of synchrophasor data in the wide area measurement system (WAMS) needs to be collected and transmitted to the phasor data concentrator, thereby increasing the possibility of being attacked by hackers. The attacked data are therefore hidden into the normal synchrophasor data so that the synchrophasor data based application will be affected. To remedy this problem, an identification framework is proposed to detect the data cyber-attack in WAMS utilizing variational mode decomposition (VMD) and multifusion support vector machine (MSVM). First, VMD is used to transform the attacked data into multiple modal components. Thereafter, a novel MSVM is employed to classify the deterministic features using the proposed linear combined multikernel (LCM). Further, this LCM can fuse multiple types of features, including the time, frequency, and statistical domains of the synchrophasor data. Utilizing the actual data from FNET/GridEye, different experiments are conducted under multiple attack strengths and types. The results demonstrate that the identification framework has higher precision and robustness compared with other conventional classifiers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Comprehensive Kinetics on the C 7 H 7 Potential Energy Surface under Combustion Conditions

The automated kinetics workflow code, KinBot, was used to explore and characterize the regions of the C 7 H 7 potential energy surface that are relevant to combustion environments and especially soot inception. We first explored the lowest-energy region, which includes the benzyl, fulvenallene + H, and cyclopentadienyl + acetylene entry points. We then expanded the model to include two higher-energy entry points, vinylpropargyl + acetylene and vinylacetylene + propargyl. The automated search was able to uncover the pathways from the literature. In addition, three important new routes were discovered: a lower-energy pathway connecting benzyl with vinylcyclopentadienyl, a decomposition mechanism from benzyl that results in side-chain hydrogen atom loss to produce fulvenallene + H, and shorter and lower energy routes to the dimethylene-cyclopentenyl intermediates. We systematically reduced the extended model to a chemically relevant domain composed of 63 wells, 10 bimolecular products, 87 barriers, and 1 barrierless channel and constructed a master equation using the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level of theory to provide rate coefficients for chemical modeling. Our calculated rate coefficients show excellent agreement with measured ones. We also simulated concentration profiles and calculated branching fractions from the important entry points to provide an interpretation of this important chemical landscape.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Flow structures of the cross-flow over a five-layer helically coiled steam generator geometry

The helically-coiled heat exchanger (HCSG) offers advantages over straight tubes such as compactness in geometry, increased heat transfer coefficients, and capability to absorb thermal expansion. In this study, the time-resolved velocity field data of shell-side crossflow over a five-layer helically coiled steam generator geometry was obtained using Particle Image Velocimetry (PIV) between adjacent rods in four regions with up to 5000 frames per second (5 kHz) at approximately Re u = 3600. The slant rod bundle created different flow patterns in different regions. Proper Orthogonal Decomposition (POD) revealed high-energy-mode flow structures. Due to the confined flow channel geometry, the large flow structures were not advected resulting in no strong POD mode paring. Continuous Wavelet Transform (CWT) visualized the characteristics of flow fluctuation in time and frequency domain simultaneously, which ranged in Strouhal number, S u , between 0.03 and 1.17. In the wake regions of the rods, flows changed in both pattern and magnitude over time. Once a pattern formed, it continued for relatively shorter than straight tube arrays, before changing the pattern again. The histograms of the POD time coefficients presented this multimodal flow characteristics. The current study focuses more on whole flow field analyses in the selected regions rather than local point-wise analyses. Finally, this study revealed that various Strouhal numbers can exist in the helically coiled heat exchanger geometry, which was not observed in straight tube arrays. In particular, Strouhal numbers smaller than 0.1 do not seem to be caused by vortex shedding but seem to be the result of the combination of complex geometry and possible multimodal trends. Based on the current observation, further works would involve intensive instantaneous point-wise analyses on the interactions between important locations such as separation points, stagnation points, shear layers, etc.

42 ENGINEERING↗
DOE OSTI2023

Spatio–Temporal Coarse–Graining Decomposition of the Global Ocean Geostrophic Kinetic Energy

We expand on a recent determination of the first global energy spectrum of the ocean’s surface geostrophic circulation using a coarse-graining (CG) method. We compare spectra from CG to those from spherical harmonics by treating land in a manner consistent with the boundary conditions. While the two methods yield qualitatively consistent domain-averaged results, spherical harmonics spectra are too noisy at gyre-scales (> 1000 km). More importantly, spherical harmonics are inherently global and cannot provide local information connecting scales with currents geographically. CG shows that the extra-tropics mesoscales (100–500 km) have a root-mean-square (rms) velocity of ~15 cm/s, which increases to ~30–40 cm/s locally in the Gulf Stream and Kuroshio and to ~16–28 cm/s in the ACC. There is notable hemispheric asymmetry in mesoscale energy-per-area, which is higher in the north due to continental boundaries. We estimate that ≈25–50% of total geostrophic energy is at scales smaller than 100 km, and is un(der)-resolved by pre-SWOT satellite products. Spectra of the time-mean circulation show that most of its energy (up to 70%) resides in stationary eddies with characteristic scales smaller than (< 500 km). This highlights the preponderance of ‘standing’ small-scale structures in the global ocean due to the temporally coherent forcing by boundaries. By coarse-graining in space and time, we compute the first spatio-temporal global spectrum of geostrophic circulation from AVISO and NEMO. These spectra show that every length-scale evolves over a wide range of time-scales with a consistent peak at ≈200 km and ≈2–3 weeks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Discrete Empirical Interpolation Method Based Dynamic Load Model Reduction

Dynamic load models add significant complexity to bulk power system time-domain simulations. The complexity is due to the large number of ordinary differential equations (ODEs) introduced by the dynamic load components such as induction motors. It is challenging to derive reduced-order models (ROMs) for dynamic loads due to the nonlinear functions in their governing equations. This paper applies the discrete empirical interpolation method enhanced proper orthogonal decomposition (DEIM-POD) to approximate the full dynamic load model with the ROM that minimizes the projection error of the nonlinear functions in dynamic load ODEs onto their dominant modes. This approach only requires evaluation of nonlinear functions at selected observation points. The observation points selected by DEIM also provide information for screening critical load buses where dynamic load model parameters contribute the most to the accuracy of ROM across multiple contingencies. The proposed approach is validated on IEEE 9-bus, WECC 179-bus and 2384-bus Polish systems.

bulk power system↗
DOE OSTI2020

Comparison of time-frequency-analysis techniques applied in building energy data noise cancellation for building load forecasting: A real-building case study

Time-frequency analysis that disaggregates a signal in both time and frequency domain is an important supporting technique for building energy analysis such as noise cancellation in data-driven building load forecasting. There is a gap in the literature related to comparing various time–frequency-analysis techniques, especially discrete wavelet transform (DWT) and empirical mode decomposition (EMD), to guide the selection and tuning of time–frequency-analysis techniques in data-driven building load forecasting. This article provides a framework to conduct a comprehensive comparison among thirteen DWT/EMD techniques with various parameters in a load forecasting modeling task. A real campus building is used as a case study for illustration. The DWT and EMD techniques are also compared under various data-driven modeling algorithms for building load forecasting. The results in the case study show that the load forecasting models trained with noise-cancelled energy data have increased their accuracy to 9.6% on average tested under unseen data. This study also shows that the effectiveness of DWT/EMD techniques depends on the data-driven algorithms used for load forecasting modeling and the training data. Hence, DWT/EMD-based noise cancellation needs customized selection and tuning to optimize their performance for data-driven building load forecasting modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A reduced-order shifted boundary method for parametrized incompressible Navier–Stokes equations

Here, we investigate a projection-based reduced order model of the steady incompressible Navier–Stokes equations for moderate Reynolds numbers. In particular, we construct an “embedded” reduced basis space, by applying proper orthogonal decomposition to the Shifted Boundary Method, a high-fidelity embedded method recently developed. We focus on the geometrical parametrization through level-set geometries, using a fixed Cartesian background geometry and the associated mesh. This approach avoids both remeshing and the development of a reference domain formulation, as typically done in fitted mesh finite element formulations. Two-dimensional computational examples for one and three parameter dimensions are presented to validate the convergence and the efficacy of the proposed approach.

"embedded" basis↗
DOE OSTI2022

Variable Active Galactic Nuclei in the Galaxy Evolution Explorer Time Domain Survey

We searched the northern hemisphere fields of the Galaxy Evolution Explorer Time Domain Survey for galaxies with UV variability indicative of active galactic nuclei (AGNs). We identified 48 high-probability candidate AGNs from a parent sample of 1819 galaxies in the NASA Sloan Atlas catalog. We further characterized these systems using optical spectroscopic diagnostics, Wide-field Infrared Survey Explorer IR color selection criteria, and spectral energy distribution modeling. Of the 48 candidates, eight were identified as AGNs from optical emission lines, two were identified by their IR colors, and 28 were identified through spectral energy decomposition. Observational biases of each selection method are discussed in connecting these AGN subsamples to another. By selecting AGNs based on UV variability, we also identified six low-mass AGN candidates, all of which would have been missed by spectroscopic selection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Time and Frequency Analysis of Load Profile Data

Technology advancements and integration of modern advanced metering systems can monitor, forecast, inform, control, and operate the building's mechanical, electrical, and plumbing (MEP) systems. They offer a higher level of information, which can contribute to making smart buildings more energy efficient and to making them closer to becoming grid-interactive energy efficient buildings (GEB). This paper builds on the ongoing research on variability analysis of a case study building with a 1-minute load profile and examines the Discrete Wavelet Transform (DWT) process in the frequency domain to quantify the signal's energy in each bandwidth, with respect to each end-use category. Moreover, the amount of variability in the total variability is not similar among the end-use categories. This information is needed to understand the behavior of the variability in the frequency domain for future applications, such as generating synthetic load profiles with a similar frequency spectrum as the measured signal.

decomposition↗