Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational reflection”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,117 records · Page 62

Knowledge-guided learning with curated prior genetic biomarkers for robust model interpretation

Abstract Motivation Knowledge-guided learning offers effective and robust model training strategies in data-scarce settings by incorporating established domain knowledge, thereby enhancing generalization, robustness, and interpretability. By contrast, conventional deep learning approaches rely purely on data-driven learning, which can limit robust model interpretability, particularly in high-dimensional settings with limited size samples. In computational biology, knowledge-guided learning has primarily leveraged network- and structural-based knowledge, leading to biologically interpretable representations and enhanced predictive performance compared to conventional approaches. However, curated biomarkers, one of the most accessible forms of biological knowledge, remain largely unexplored within knowledge-guided paradigms. Results In this study, we propose a model-agnostic training paradigm, Biomarker-driven Explainable Prior-guided Learning (BioExPL), that can be applied to any neural networks that incorporates curated prior knowledge. BioExPL enforces neural networks to reflect curated biomarker priors in their latent representations through a novel knowledge-alignment loss. BioExPL consistently demonstrated significantly improved predictive performance and enhanced model interpretability with minimized computational overhead in simulation studies and intensive experiments on multiple cancer datasets. BioExPL not only integrates prior curated knowledge into the model but also accurately identifies unknown associated signals additionally. BioExPL is model-agnostic and domain-independent, enabling its integration into diverse neural network architectures. Availability and implementation The open-source is publicly available at: https://github.com/datax-lab/BioExPL.

Baek, Beomsu [Department of Computer Science, Univ↗

Two regimes in the regularity of sunspot number

Sunspot numbers WN display quasi-periodical variations that undergo regime changes. These irregularities could indicate a chaotic system and be measured by Lyapunov exponents. We define a functional λ (an 'irregularity index') that is close to the (maximal) Lyapunov exponent for dynamical systems and well defined for series with a random component: this allows one to work with sunspot numbers. We compute λ for the daily WN from 1850 to 2012 within 4 yr sliding windows: λ exhibit sharp maxima at solar minima and secondary maxima at solar maxima. This pattern is reflected in the ratio R of the amplitudes of the main versus secondary peaks. Two regimes have alternated in the past 150 yr, R1 from 1850 to 1915 (large λ and R values) and R2 from 1935 to 2005 (shrinking difference between main and secondary maxima, R values between 1 and 2). We build an autoregressive model consisting of Poisson noise plus an 11 yr cycle and compute its irregularity index. The transition from R1 to R2 can be reproduced by strengthening the autocorrelation a of the model series. The features of the two regimes are stable for model and WN with respect to embedding dimension and delay. Near the time of the last solar minimum (∼2008), the irregularity index exhibits a peak similar to the peaks observed before 1915. This might signal a regime change back from R2 to R1 and the onset of a significant decrease of solar activity.

79 ASTRONOMY AND ASTROPHYSICS↗

An efficient sparse matrix multiplication scheme for the CYBER 205 computer

This paper describes the development of an efficient algorithm for computing the product of a matrix and vector on a CYBER 205 vector computer. The desire to provide software which allows the user to choose between the often conflicting goals of minimizing central processing unit (CPU) time or storage requirements has led to a diagonal-based algorithm in which one of four types of storage is selected for each diagonal. The candidate storage types employed were chosen to be efficient on the CYBER 205 for diagonals which have nonzero structure which is dense, moderately sparse, very sparse and short, or very sparse and long; however, for many densities, no diagonal type is most efficient with respect to both resource requirements, and a trade-off must be made. For each diagonal, an initialization subroutine estimates the CPU time and storage required for each storage type based on results from previously performed numerical experimentation. These requirements are adjusted by weights provided by the user which reflect the relative importance the user places on the two resources. The adjusted resource requirements are then compared to select the most efficient storage and computational scheme.

Lambiotte, Jules J., Jr.↗

WE-EF-207-06: Dedicated Cone-Beam Breast CT with Laterally-Shifted Detector: Monte Carlo Evaluation of X-Ray Scatter Distribution and Scatter-To-Primary Ratio

Purpose: To determine the spatial distribution of x-ray scatter and scatter-to-primary ratio (SPR) in projections during cone-beam breast CT (CBBCT) with laterally-shifted detector that results in coronal (fan-angle) truncation. Methods: We hypothesized that CBBCT with coronal truncation would lower SPR due to reduction in irradiated breast volume, and that the location of maximum x-ray scatter fluence (scatter-peak) in the detector plane can be determined from the ratio of irradiated-to-total breast volume, breast dimensions and system geometry. Monte Carlo simulations (GEANT4) reflecting a prototype CBBCT system were used to record the position-dependent primary and scatter x-ray photon fluence incident on the detector without coronal truncation (full fan-angle, 2f=24-degrees) and with coronal truncation (fan-angle, f+ f=12+2.7-degrees). Semi-ellipsoidal breasts (10/14/18-cm diameter, chest-wall to nipple length: 0.75xdiameter, 2%/14%/100% fibroglandular content) aligned with the axis-of-rotation (AOR) were modeled. Mono-energy photons were simulated and weighted for 2 spectra (49kVp, 1.4-mm Al HVL; 60kVp, 3.76-mm Al HVL). In addition to SPR, the scatter maps were analyzed to identify the location of the scatter-peak. Results: For CBBCT without fan-angle truncation, the scatter-peaks were aligned with the projection of the AOR onto the detector for all breasts. With truncated fan-beam, the scatter-peaks were laterally-shifted from the projection of the AOR along the fan-angle direction by 14/38/70-pixels for 10/14/18-cm diameter breasts. The corresponding theoretical shifts were 14.8/39.7/68-pixels (p=0.47, 2-tailed paired-ratio t-test). Along the cone-angle, the shift in scatter-peaks between truncated and full-fan angle CBBCT were 2/2/4 -pixels for 10/14/18-cm diameter breasts. CBBCT with fan-angle truncation reduced SPR by 14/22/28% for 10/14/18-cm diameter breasts. 60kVp reduced SPR by 21–25% compared to 49kVp. Peak SPR for CBBCT with fan-angle truncation (60kVp) were 0.09/0.25/0.73 for 10/14/18-cm diameter breasts. Conclusion: CBBCT with laterally-shifted detector geometry and with appropriate kVp/beam quality reduces SPR. If residual scatter needs correction, the location corresponding to scatter-peak can be analytically computed. This work was supported in part by NIH R01 CA128906. The contents are solely the responsibility of the authors and do not reflect the official views of the NIH or NCI.

60 APPLIED LIFE SCIENCES↗

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗

Fundamental analysis of the failure of polymer-based fiber reinforced composites

A mathematical model is described which will permit predictions of the strength of fiber reinforced composites containing known flaws to be made from the basic properties of their constituents. The approach was to embed a local heterogeneous region (LHR) surrounding the crack tip into an anisotropic elastic continuum. The model should (1) permit an explicit analysis of the micromechanical processes involved in the fracture process, and (2) remain simple enough to be useful in practical computations. Computations for arbitrary flaw size and orientation under arbitrary applied load combinations were performed from unidirectional composites with linear elastic-brittle constituent behavior. The mechanical properties were nominally those of graphite epoxy. With the rupture properties arbitrarily varied to test the capability of the model to reflect real fracture modes in fiber composites, it was shown that fiber breakage, matrix crazing, crack bridging, matrix-fiber debonding, and axial splitting can all occur during a period of (gradually) increasing load prior to catastrophic fracture. The computations reveal qualitatively the sequential nature of the stable crack process that precedes fracture.

Kanninen, M. F.↗

Response to reflected-force feedback to fingers in teleoperations

Reflected-force feedback is an important aspect of teleoperations. The objective is to determine the ability of the human operator to respond to that force. Telerobotics operation is simulated by computer control of a motor-driven device with capabilities for programmable force feedback and force measurement. A computer-controlled motor drive is developed that provides forces against the fingers as well as (angular) position control. A load cell moves in a circular arc as it is pushed by a finger and measures reaction forces on the finger. The force exerted by the finger on the load cell and the angular position are digitized and recorded as a function of time by the computer. Flexure forces of the index, long and ring fingers of the human hand in opposition to the motor driven load cell are investigated. Results of the following experiments are presented: (1) Exertion of maximum finger force as a function of angle; (2) Exertion of target finger force against a computer controlled force; and (3) Test of the ability to move to a target force against a force that is a function of position. Averaged over ten individuals, the maximum force that could be exerted by the index or long finger is about 50 Newtons, while that of the ring finger is about 40 Newtons. From the tests of the ability of a subject to exert a target force, it was concluded that reflected-force feedback can be achieved with the direct kinesthetic perception of force without the use of tactile or visual clues.

Sutter, P. H.↗

Integral geometry and holography

We present a mathematical framework which underlies the connection between information theory and the bulk spacetime in the AdS 3 /CFT 2 correspondence. A key concept is kinematic space: an auxiliary Lorentzian geometry whose metric is defined in terms of conditional mutual informations and which organizes the entanglement pattern of a CFT state. When the field theory has a holographic dual obeying the Ryu-Takayanagi proposal, kinematic space has a direct geometric meaning: it is the space of bulk geodesics studied in integral geometry. Lengths of bulk curves are computed by kinematic volumes, giving a precise entropic interpretation of the length of any bulk curve. We explain how basic geometric concepts -- points, distances and angles -- are reflected in kinematic space, allowing one to reconstruct a large class of spatial bulk geometries from boundary entanglement entropies. In this way, kinematic space translates between information theoretic and geometric descriptions of a CFT state. As an example, we discuss in detail the static slice of AdS 3 whose kinematic space is two-dimensional de Sitter space.

97 MATHEMATICS AND COMPUTING↗

Reflection and extinction of light by self-assembled monolayers of a quinque-thiophene derivative: A coherent scattering approach

Scattering matrix theory is used to describe resonant optical properties of molecular monolayers. Three types of coupling are included: exciton-exciton, exciton-photon, and exciton-phonon coupling. We use the K-matrix formalism, developed originally to describe neutron scattering spectra in nuclear physics to compute the scattering of polaritons by phonons. This perturbation approach takes into account the three couplings and allows one to go beyond molecular exciton theory without the need of introducing additional boundary conditions for the polariton. We demonstrate that reflection, absorption, and extinction of light by 2D self-assembled monolayers of molecules containing quinque-thiophene chromophoric groups can be calculated. The extracted coherence length of the Frenkel exciton is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Operational MODIS Cloud Optical and Microphysical Property Product: Overview of the Collection 6 Algorithm and Preliminary Results

Operational Moderate Resolution Imaging Spectroradiometer (MODIS) retrievals of cloud optical and microphysical properties (part of the archived products MOD06 and MYD06, for MODIS Terra and Aqua, respectively) are currently being reprocessed along with other MODIS Atmosphere Team products. The latest "Collection 6" processing stream, which is expected to begin production by summer 2012, includes updates to the previous cloud retrieval algorithm along with new capabilities. The 1 km retrievals, based on well-known solar reflectance techniques, include cloud optical thickness, effective particle radius, and water path, as well as thermodynamic phase derived from a combination of solar and infrared tests. Being both global and of high spatial resolution requires an algorithm that is computationally efficient and can perform over all surface types. Collection 6 additions and enhancements include: (i) absolute effective particle radius retrievals derived separately from the 1.6 and 3.7 !-lm bands (instead of differences relative to the standard 2.1 !-lm retrieval), (ii) comprehensive look-up tables for cloud reflectance and emissivity (no asymptotic theory) with a wind-speed interpolated Cox-Munk BRDF for ocean surfaces, (iii) retrievals for both liquid water and ice phases for each pixel, and a subsequent determination of the phase based, in part, on effective radius retrieval outcomes for the two phases, (iv) new ice cloud radiative models using roughened particles with a specified habit, (v) updated spatially-complete global spectral surface albedo maps derived from MODIS Collection 5, (vi) enhanced pixel-level uncertainty calculations incorporating additional radiative error sources including the MODIS L1 B uncertainty index for assessing band and scene-dependent radiometric uncertainties, (v) and use of a new 1 km cloud top pressure/temperature algorithm (also part of MOD06) for atmospheric corrections and low cloud non-unity emissivity temperature adjustments.

Platnick, Steven↗

Helioseismic Constraints on New Solar Models from the MoSEC Code

Evolutionary solar models are computed using a new stellar evolution code, MOSEC (Modular Stellar Evolution Code). This code has been designed with carefully controlled truncation errors in order to achieve a precision which reflects the increasingly accurate determination of solar interior structure by helioseismology. A series of models is constructed to investigate the effects of the choice of equation of state (OPAL or MHD-E, the latter being a version of the MHD equation of state recalculated by the author), the inclusion of helium and heavy-element settling and diffusion, and the inclusion of a simple model of mixing associated with the solar tachocline. The neutrino flux predictions are discussed, while the sound speed of the computed models is compared to that of the sun via the latest inversion of SOI-NMI p-mode frequency data. The comparison between models calculated with the OPAL and MHD-E equations of state is particularly interesting because the MHD-E equation of state includes relativistic effects for the electrons, whereas neither MHD nor OPAL do. This has a significant effect on the sound speed of the computed model, worsening the agreement with the solar sound speed. Using the OPAL equation of state and including the settling and diffusion of helium and heavy elements produces agreement in sound speed with the helioseismic results to within about +.-0.2%; the inclusion of mixing slightly improves the agreement.

Elliott, J. R.↗

Non-Reflecting Regions for Finite Difference Methods in Modeling of Elastic Wave Propagation in Plates

Solution of the wave equation using techniques such as finite difference or finite element methods can model elastic wave propagation in solids. This requires mapping the physical geometry into a computational domain whose size is governed by the size of the physical domain of interest and by the required resolution. This computational domain, in turn, dictates the computer memory requirements as well as the calculation time. Quite often, the physical region of interest is only a part of the whole physical body, and does not necessarily include all the physical boundaries. Reduction of the calculation domain requires positioning an artificial boundary or region where a physical boundary does not exist. It is important however that such a boundary, or region, will not affect the internal domain, i.e., it should not cause reflections that propagate back into the material. This paper concentrates on the issue of constructing such a boundary region.

Kishoni, Doron↗

International Conference on Numerical Methods in Fluid Dynamics, 7th, Stanford University, Stanford and Moffett Field, CA, June 23-27, 1980, Proceedings

Topics discussed include polygon transformations in fluid mechanics, computation of three-dimensional horseshoe vortex flow using the Navier-Stokes equations, an improved surface velocity method for transonic finite-volume solutions, transonic flow calculations with higher order finite elements, the numerical calculation of transonic axial turbomachinery flows, and the simultaneous solutions of inviscid flow and boundary layer at transonic speeds. Also considered are analytical solutions for the reflection of unsteady shock waves and relevant numerical tests, reformulation of the method of characteristics for multidimensional flows, direct numerical simulations of turbulent shear flows, the stability and separation of freely interacting boundary layers, computational models of convective motions at fluid interfaces, viscous transonic flow over airfoils, and mixed spectral/finite difference approximations for slightly viscous flows.

Reynolds, W. C.↗

Multispectral imaging and analysis system

Arrays of charge coupled devices or linear detector arrays simultaneously obtain spectral reflectance data of different wavelengths for a target area. Several accommodating a particular bandwidth, are individually associated with each array. Data from the arrays are read out in parallel and applied to a computer or microprocessor for processing. The microprocessor serves to analyze the data in real time and if possible, in accordance with hard-wired algorithms. The data are then displayed as an image on an appropriate display unit and also recorded for further use. The display system may be operationally connected to receive a terrain image such that the target area and the analyzed spectral reflectance data are superimposed and simultaneously displayed.

Goetz, A. F. H.↗

Surface reflectance measurements in the UV from an airborne platform. I

The reflectance of naturally occurring surfaces is measured from a twin engine Cessna 402B aircraft both spectrally (0.29-0.40 micron) with a compact double monochromator and broadband (0.29-1.2 micron) with a UV enhanced photodiode. The measurement system, which is computer controlled and electrically isolated from the aircraft, consists of upward and downward looking hemispheric diffusers, filters, a rotating 90 deg mirror, a focusing lens, and a double monochromator/PMT or the UV photodiode. Measurements are taken at several altitudes enabling the empirical determination of backscatter and attenuation effects on the reflectance. The results are presented for pine forest canopy, green farmland, open ocean, and brown farmland as a function of wavelength and altitude.

Doda, D. D.↗

Quantification of Pediatric Abdominal Organ Motion With a 4-Dimensional Magnetic Resonance Imaging Method

Purpose: To characterize respiration-induced abdominal organ motion in children receiving radiation treatment with a 4-dimensional (4D) magnetic resonance imaging (MRI) method. Methods and Materials: We analyzed free-breathing coronal 4D MRI datasets acquired from 35 patients (aged 1-20 years) with abdominal tumors. A deformable image registration of the 4D MRI datasets was performed to derive motion trajectories of selected anatomic landmarks, from which organ motions were quantified. The association between organ motion and patient characteristics was investigated and compared with previous studies. The relation between patient height and organ motion was further investigated to predict organ motion in prospective patients. Results: Organ motion and its individual variation were reduced in younger patients (eg, kidney peak-to-peak motion <5 mm for all but 1 patient aged ≤8 years), although special motion management may be warranted in some adolescents. The liver and spleen exhibited greater motion than did the kidneys, while intraorgan variation was present. The motions in the liver and kidneys agreed with those reported by the previous 4D computed tomography studies. Individual variations of organ motion in younger patients were due, in part, to changes in respiration rate, which ostensibly reflected the effect of anesthesia. The prediction of organ motion was limited by large individual variations, particularly for older patients. Conclusions: The 4D MRI acquisition method and motion analysis described in this study provide a nonionizing approach to understand age-associated organ motion, which aids in the planning of abdominal radiation therapy for pediatric patients. Use of 4D MRI facilitates monitoring of changes in target motion patterns during treatment courses and in various studies of the effect of organ motion on radiation treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗