Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Laminar and turbulent nozzle-jet flows and their acoustic near-field

We investigate numerically the effects of nozzle-exit flow conditions on the jet-flow development and the near-field sound at a diameter-based Reynolds number of Re{sub D} = 18 100 and Mach number Ma = 0.9. Our computational setup features the inclusion of a cylindrical nozzle which allows to establish a physical nozzle-exit flow and therefore well-defined initial jet-flow conditions. Within the nozzle, the flow is modeled by a potential flow core and a laminar, transitional, or developing turbulent boundary layer. The goal is to document and to compare the effects of the different jet inflows on the jet flow development and the sound radiation. For laminar and transitional boundary layers, transition to turbulence in the jet shear layer is governed by the development of Kelvin-Helmholtz instabilities. With the turbulent nozzle boundary layer, the jet flow development is characterized by a rapid changeover to a turbulent free shear layer within about one nozzle diameter. Sound pressure levels are strongly enhanced for laminar and transitional exit conditions compared to the turbulent case. However, a frequency and frequency-wavenumber analysis of the near-field pressure indicates that the dominant sound radiation characteristics remain largely unaffected. By applying a recently developed scaling procedure, we obtain a close match of the scaled near-field sound spectra for all nozzle-exit turbulence levels and also a reasonable agreement with experimental far-field data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Spectral Upscaling for Graph Laplacian Problems with Application to Reservoir Simulation

Here, we consider coarsening procedures for graph Laplacian problems written in a mixed saddle-point form. In that form, in addition to the original (vertex) degrees of freedom (dofs), we also have edge degrees of freedom. We extend previously developed aggregation-based coarsening procedures applied to both sets of dofs to now allow more than one coarse vertex dof per aggregate. Those dofs are selected as certain eigenvectors of local graph Laplacians associated with each aggregate. Additionally, we coarsen the edge dofs by using traces of the discrete gradients of the already constructed coarse vertex dofs. These traces are defined on the interface edges that connect any two adjacent aggregates. The overall procedure is a modification of the spectral upscaling procedure developed in for the mixed finite element discretization of diffusion type PDEs which has the important property of maintaining inf-sup stability on coarse levels and having provable approximation properties. We consider applications to partitioning a general graph and to a finite volume discretization interpreted as a graph Laplacian, developing consistent and accurate coarse-scale models of a fine-scale problem.

58 GEOSCIENCES↗
DOE OSTI2012

Sensitivity analysis of coupled criticality calculations

Perturbation theory based sensitivity analysis is a vital part of todays' nuclear reactor design. This paper presents an extension of standard techniques to examine coupled criticality problems with mutual feedback between neutronics and an augmenting system (for example thermal-hydraulics). The proposed procedure uses a neutronic and an augmenting adjoint function to efficiently calculate the first order change in responses of interest due to variations of the parameters describing the coupled problem. The effect of the perturbations is considered in two different ways in our study: either a change is allowed in the power level while maintaining criticality (power perturbation) or a change is allowed in the eigenvalue while the power is constrained (eigenvalue perturbation). The calculated response can be the change in the power level, the reactivity worth of the perturbation, or the change in any functional of the flux, the augmenting dependent variables and the input parameters. To obtain power- and criticality-constrained sensitivities power- and k-reset procedures can be applied yielding identical results. Both the theoretical background and an application to a one dimensional slab problem are presented, along with an iterative procedure to compute the necessary adjoint functions using the neutronics and the augmenting codes separately, thus eliminating the need of developing new programs to solve the coupled adjoint problem. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Virtual Design and Analysis of the Balancing Element of an External Gear Machine Considering Cavitation and Mixed Lubrication Effects

This report presents a simulation based approach for sizing the axial balance elements of a pressure compensated external gear machine. The reference hydraulic unit must be able to guarantee operation in two quadrants, as either a pump or a motor. It also needs to operate in a wide operating range of speed and pressure, as it is to be used as primary hydraulic unit for an electro-hydraulic actuator (EHA). The design procedure builds upon past work at the author’s center, but it extends it to the case of multiple quadrant units. Also, a unique method for modeling of both mixed lubrication and cavitation that might occur in the lubricating film is used. After describing the design procedure, the paper discusses the features of the balancing of a reference unit, along with the effects of both cavitation and mixed lubrication.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

A harmonic polynomial cell (HPC) method for 3D Laplace equation with application in marine hydrodynamics

We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods, e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

SU-F-T-347: An Absolute Dose-Volume Constraint Based Deterministic Optimization Framework for Multi-Co60 Source Focused Radiotherapy

Purpose: Treatment plan optimization in multi-Co60 source focused radiotherapy with multiple isocenters is challenging, because dose distribution is normalized to maximum dose during optimization and evaluation. The objective functions are traditionally defined based on relative dosimetric distribution. This study presents an alternative absolute dose-volume constraint (ADC) based deterministic optimization framework (ADC-DOF). Methods: The initial isocenters are placed on the eroded target surface. Collimator size is chosen based on the area of 2D contour on corresponding axial slice. The isocenter spacing is determined by adjacent collimator sizes. The weights are optimized by minimizing the deviation from ADCs using the steepest descent technique. An iterative procedure is developed to reduce the number of isocenters, where the isocenter with lowest weight is removed without affecting plan quality. The ADC-DOF is compared with the genetic algorithm (GA) using the same arbitrary shaped target (254cc), with a 15mm margin ring structure representing normal tissues. Results: For ADC-DOF, the ADCs imposed on target and ring are (D100>10Gy, D50,10, 0<12Gy, 15Gy and 20Gy) and (D40<10Gy). The resulting D100, 50, 10, 0 and D40 are (9.9Gy, 12.0Gy, 14.1Gy and 16.2Gy) and (10.2Gy). The objectives of GA are to maximize 50% isodose target coverage (TC) while minimize the dose delivered to the ring structure, which results in 97% TC and 47.2% average dose in ring structure. For ADC-DOF (GA) techniques, 20 out of 38 (10 out of 12) initial isocenters are used in the final plan, and the computation time is 8.7s (412.2s) on an i5 computer. Conclusion: We have developed a new optimization technique using ADC and deterministic optimization. Compared with GA, ADC-DOF uses more isocenters but is faster and more robust, and achieves a better conformity. For future work, we will focus on developing a more effective mechanism for initial isocenter determination.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗
DOE OSTI2023

Predeployment progress of the Canister Deposition Field Demonstration

This report updates the high-level test plan for evaluating surface deposition on three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUTECH Horizontal Modular Storage (NUHOMS) Advanced Horizontal Storage Modules (AHSMs) from Orano (formerly Transnuclear Inc.) and provides a summary of the surface sampling activities that have been conducted to date. The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or requirements. One goal of this testing is to collect defensible and detailed dust deposition measurements from the surface of dry storage canisters in a marine coastal environment to guide chloride-induced stress corrosion cracking (CISCC) research. Another goal is to provide data for the validation of computational fluid dynamics (CFD) based deposition modeling. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the decay heat and thermal hydraulic environment. Test and heater design is supported by detailed CFD modeling. Instrumentation throughout the canister, storage module, and environment will provide extensive information about the thermal-hydraulic behavior of horizontal dry cask storage systems. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer detailed quantification and composition of the deposited particulates from a realistic storage environment. Discussions of a potential host site for the Canister Deposition Field Demonstration (CDFD) are ongoing. Until a host site is chosen, testing of key CDFD hardware components including the heater assemblies, power skid, and remote data acquisition system will continue. Functional testing of the finalized heater assemblies and test apparatus started this fiscal year. These initial heater tests have shown the assemblies are performing within design specifications. Staged surface sampling of a mockup of a canister outside the AHSM on a transfer skid was also performed. Refinements to the sampling procedures and techniques were captured from observation of these activities and lessons-learned debriefs. These updated sampling procedures and techniques are planned to be tested again in the field using the mockup in order to assure personnel are using the most accurate and repeatable methods possible prior to deployment for actual CDFD testing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

RECONSTRUCTION OF INDIVIDUAL DOSES DUE TO MEDICAL EXPOSURES FOR MEMBERS OF THE TECHA RIVER COHORT

Purpose: To describe a methodology for reconstruction of doses due to medical exposures for members of the Techa River Cohort (TRC) who received diagnostic radiation at the clinic of the Urals Research Center for Radiation Medicine (URCRM) in 1952–2005. To calculate doses of medical exposure for the TRC members and compare with the doses that resulted from radioactive contamination of the Techa River. Material and Methods: Reconstruction of individual medical doses is based on data on x-ray diagnostic procedures available for each person examined at the URCRM clinics and values of absorbed dose in 12 organs per typical x-ray procedure calculated with the use of a mathematical phantom. Personal data on x-ray diagnostic examinations have been complied in the computerized “Registry of x-ray diagnostic procedures.” Sources of information are archival registry books from the URCRM x-ray room (available since 1956) and records on x-ray diagnostic procedures in patient-case histories (since 1952). The absorbed doses for 12 organs of interest have been evaluated per unit typical x-ray procedure with account taken of the x-ray examination parameters characteristic for the diagnostic machines used at the URCRM clinics. These parameters have been evaluated from published data on technical characteristics of the x-ray diagnostic machines used at the URCRM clinics in 1952–1988 and taken from the x-ray room for machines used at the URCRM in 1989–2005. Absorbed doses in the 12 organs per unit typical x-ray procedure have been calculated with use of a special computer code, EDEREX, developed at the Saint-Petersburg Research Institute of Radiation Hygiene after Professor P.V. Ramzaev. Individual accumulated doses of medical exposure have been calculated with a computer code, MEDS (Medical Exposure Dosimetry System), specifically developed at the URCRM. Results: At present, the “Registry of x-ray diagnostic procedures” contains information on individual x-ray examinations for over 9,500 persons including 6,415 TRC members. Statistical analysis of the Registry data showed that the more frequent types of examinations were fluoroscopy and radiography of the chest and fluoroscopy of the stomach and the esophagus. Average absorbed doses accumulated by year 2005 calculated for the 12 organs varied from 4 mGy for testes to 40 mGy for bone surfaces. Maximum individual medical doses could reach 500–650 mGy and in some cases exceeded doses from exposure at the Techa River. Conclusions: For the first time the doses of medical exposure were calculated and analyzed for members of the Techa River Cohort who received diagnostic radiation at the URCRM clinics. These results are being used in radiation-risk analysis to adjust for this source of confounding exposure in the TRC.

dose reconstruction↗
DOE OSTI2026

Demonstration of Algorithm for Sensor Placement Optimization using Simulation Data

This deliverable reports FY26 progress in advancing a neural-network-based Green’s-function framework for reconstructing neutron-flux distributions from ex-core measurements and for translating reconstruction requirements into a practical detector-layout strategy. Building on the FY25 formulation, the present work had two main objectives: (1) refine and re-evaluate the reconstruction methodology on an updated Purdue University Reactor Number One (PUR-1) model, and (2) develop a systematic, Green’s-function-guided procedure for boundary detector placement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Computational Design of a Novel Medium-Carbon, Low-Alloy Steel Microalloyed with Niobium

The design of a new steel with specific properties is always challenging owing to the complex interactions of many variables. In this work, this challenge is dealt with by combining metallurgical principles with computational thermodynamics and kinetics to design a novel steel composition suitable for thermomechanical processing and induction heat treatment to achieve a hardness level in excess of 600 HV with the potential for good fracture toughness. CALPHAD-based packages for the thermodynamics and kinetics of phase transformations and diffusion, namely Thermo-Calc{sup ®} and JMatPro{sup ®}, have been combined with an interdendritic segregation tool (IDS) to optimize the contents of chromium, molybdenum and niobium in a proposed medium-carbon low-manganese steel composition. Important factors taken into account in the modeling and optimization were hardenability and as-quenched hardness, grain refinement and alloying cost. For further investigations and verification, the designed composition, i.e., in wt.% 0.40C, 0.20Si, 0.25Mn, 0.90Cr, 0.50Mo, was cast with two nominal levels of Nb: 0 and 0.012 wt.%. The results showed that an addition of Nb decreases the austenite grain size during casting and after slab reheating prior to hot rolling. Validation experiments showed that the predicted properties, i.e., hardness, hardenability and level of segregation, for the designed composition were realistic. It is also demonstrated that the applied procedure could be useful in reducing the number of experiments required for developing compositions for other new steels.

36 MATERIALS SCIENCE↗
DOE OSTI2018

ReaxFF molecular dynamics simulation of intermolecular structure formation in acetic acid-water mixtures at elevated temperatures and pressures

The intermolecular structure formation in liquid and supercritical acetic acid-water mixtures was explored using ReaxFF-based molecular dynamics simulations. The microscopic structures of acetic acid-water mixtures with different acetic acid mole fractions (1.0 ≥ x HAc ≥ 0.2) at ambient and critical conditions were examined. The potential energy surface associated with the dissociation of acetic acid molecules was calculated using a metadynamics procedure to optimize the dissociation energy of ReaxFF potential. At ambient conditions, depending on the acetic acid concentration, either acetic acid clusters or water clusters are dominant in the liquid mixture. When acetic acid is dominant (0.4 ≤ x HAc ), cyclic dimers and chain structures between acetic acid molecules are present in the mixture. Both structures disappear at increased water content of the mixture. It was found by simulations that the acetic acid molecules released from these dimer and chain structures tend to stay in a dipole-dipole interaction. These structural changes are in agreement with the experimental findings. When switched to critical conditions, the long-range interactions (e.g., second or fourth neighbor) disappear and the water-water and acetic acid-acetic acid structural formations become disordered. The simulated radial distribution function for water-water interactions is in agreement with experimental and computational studies. The first neighbor interactions between acetic acid and water molecules are preserved at relatively lower temperatures of the critical region. As higher temperatures are reached in the critical region, these interactions were observed to weaken. These simulations suggest that ReaxFF molecular dynamics simulations are an appropriate tool for studying supercritical water/organic acid mixtures.

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

Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems

The objective of this paper is to present a local bounds preserving stabilized finite element scheme for hyperbolic systems on unstructured meshes based on continuous Galerkin (CG) discretization in space. A CG semi-discrete scheme with low order artificial dissipation that satisfies the local extremum diminishing (LED) condition for systems is used to discretize a system of conservation equations in space. The low order artificial diffusion is based on approximate Riemann solvers for hyperbolic conservation laws. In this case we consider both Rusanov and Roe artificial diffusion operators. In the Rusanov case, two designs are considered, a nodal based diffusion operator and a local projection stabilization operator. The result is a discretization that is LED and has first order convergence behavior. To achieve high resolution, limited antidiffusion is added back to the semi-discrete form where the limiter is constructed from a linearity preserving local projection stabilization operator. The procedure follows the algebraic flux correction procedure usually used in flux corrected transport algorithms. To further deal with phase errors (or terracing) common in FCT type methods, high order background dissipation is added to the antidiffusive correction. The resulting stabilized semi-discrete scheme can be discretized in time using a wide variety of time integrators. Numerical examples involving nonlinear scalar Burgers equation, and several shock hydrodynamics simulations for the Euler system are considered to demonstrate the performance of the method. In conclusion, for time discretization, Crank–Nicolson scheme and backward Euler scheme are utilized.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Measurement of joint kinematics using a conventional clinical single-perspective flat-panel radiography system

Purpose: The ability to accurately measure joint kinematics is an important tool in studying both normal joint function and pathologies associated with injury and disease. The purpose of this study is to evaluate the efficacy, accuracy, precision, and clinical safety of measuring 3D joint motion using a conventional flat-panel radiography system prior to its application in an in vivo study. Methods: An automated, image-based tracking algorithm was implemented to measure the three-dimensional pose of a sparse object from a two-dimensional radiographic projection. The algorithm was tested to determine its efficiency and failure rate, defined as the number of image frames where automated tracking failed, or required user intervention. The accuracy and precision of measuring three-dimensional motion were assessed using a robotic controlled, tibiofemoral knee phantom programmed to mimic a subject with a total knee replacement performing a stair ascent activity. Accuracy was assessed by comparing the measurements of the single-plane radiographic tracking technique to those of an optical tracking system, and quantified by the measurement discrepancy between the two systems using the Bland-Altman technique. Precision was assessed through a series of repeated measurements of the tibiofemoral kinematics, and was quantified using the across-trial deviations of the repeated kinematic measurements. The safety of the imaging procedure was assessed by measuring the effective dose of ionizing radiation associated with the x-ray exposures, and analyzing its relative risk to a human subject. Results: The automated tracking algorithm displayed a failure rate of 2% and achieved an average computational throughput of 8 image frames/s. Mean differences between the radiographic and optical measurements for translations and rotations were less than 0.08 mm and 0.07 Degree-Sign in-plane, and 0.24 mm and 0.6 Degree-Sign out-of-plane. The repeatability of kinematics measurements performed using the radiographic tracking technique was better than {+-}0.09 mm and 0.12 Degree-Sign in-plane, and {+-}0.70 mm and {+-}0.07 Degree-Sign out-of-plane. The effective dose associated with the imaging protocol used was 15 {mu}Sv for 10 s of radiographic cine acquisition. Conclusions: This study demonstrates the ability to accurately measure knee-joint kinematics using a single-plane radiographic measurement technique. The measurement technique can be easily implemented at most clinical centers equipped with a modern-day radiographic x-ray system. The dose of ionizing radiation associated with the image acquisition represents a minimal risk to any subjects undergoing the examination.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Solving The Long-Standing Problem Of Nuclear Reactions At The Highest Microscopic Level: Annual Continuation And Progress Report

The aim of this project is to develop a comprehensive framework that will lead to a fundamental description of both structural properties and reactions of light nuclei in terms of constituent protons and neutrons interacting through nucleon-nucleon (NN) and three-­nucleon (NNN) forces. This project will provide the research community with the theoretical and computational tools that will enable: (1) an accurate prediction for fusion reactions that power stars and Earth-­based fusion facilities; (2) an improved description of the spectroscopy of exotic nuclei, including light Borromean systems; and (3) a fundamental understanding of the three-­nucleon force in nuclear reactions and nuclei at the drip line. To achieve this goal, we build upon a promising technique emerged recently as a candidate to reach a fundamental description of low-­energy binary reactions between light ions, that is the ab initio no-­core shell model combined with the resonating-­group method (NCSM/RGM). This approach has demonstrated the capability to describe binary reactions below the three-­body breakup threshold based, up to now, on similarity-­renormalization-­group (SRG) 5 evolved NN only potentials. To advance the understanding of nuclear reactions at low energies and light exotic nuclei, this project aims at extending the NCSM/RGM approach to include the full range of NNN interactions as well as the treatment of three-­cluster bound and continuum states. Three-­nucleon interactions are unavoidable components of a fundamental nuclear Hamiltonian obtained in a low-­energy effective theory. In addition, three-­nucleon force terms are induced by the SRG procedure and have to be taken into account for such a transformation to be unitary in many-­body calculations. At the same time, the introduction of three-­body cluster states is key to achieve a microscopic description of Borromean systems as well as three-­body breakup reactions. This project will both enhance the fundamentality and enlarge the scope of our microscopic description of nuclear properties. A successful completion of this project will result in improved accuracy of the 3 He( 3 He,2p) 4 He and 3 He(α,γ) 7 Be reaction rates and consequently, in enhancement of the predictive capability of the standard solar model. In addition, we will study also the mirror reactions 3 H( 3 H,2n) 4 He and 3 H(α,γ) 7 Li (a key reaction for the production of 7 Li in the standard big-­bang nucleosynthesis), and the spectroscopy of the 6 He, 6 Be, and 11 Li nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A semiparametric latent factor model for large scale temporal data with heteroscedasticity

Large scale temporal data have flourished in a vast array of applications, and their sophisticated structures, especially the heteroscedasticity among subjects with inter- and intra-temporal dependence, have fueled a great demand for new statistical models. In this paper, with covariate information, we consider a flexible model for large scale temporal data with subject-specific heteroscedasticity. Formally, the model employs latent semiparametric factors to simultaneously account for the subject-specific heteroscedasticity and the contemporaneous and/or serial correlations. The subject-specific heteroscedasticity is modeled as the product of the unobserved factor process and subject’s covariate effect, which is further characterized via additive models. For estimation, we propose a two-step procedure. First, the latent factor process and nonparametric loading are recovered through projection-based methods, and following, we estimate the regression components by approaches motivated from the generalized least squares. By scrupulously examining the non-asymptotic rates for recovering the factor process and its loading, we show the consistency and efficiency of estimated regression coefficients in the absence of prior knowledge of latent factor process and subject’s covariate effect. Here, the statistical guarantees remain valid even for finite time points that makes our method particularly appealing when the subjects significantly outnumber the observation time points. Using comprehensive simulations, we demonstrate the finite sample performance of our method, which corroborates the theoretical findings. Finally, we apply our method to a data set of air quality and energy consumption collected at 129 monitoring sites in the United States in 2015.

97 MATHEMATICS AND COMPUTING↗