SEARCH · Engineering Papers

Results for “Level set method”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗

NASA NTRS2012

Elements of Nondestructive Examination for the Visual Inspection of Composite Structures

Visual inspection (VI) of composite structures provides an effective, wide field survey to ensure design and material compliance is maintained for the entire service life of the manufactured component. Visual inspection is one of the most commonly used nondestructive inspection methods to assess surface defects on composites. By applying visual inspection elements, mechanical damage that could affect component residual strength can be identified and a disposition reached before a potentially catastrophic event occurs. This technique is non-contact and applied at all stages of the composite structure s processing and use. Additionally, VI is required to be performed in all service environments until decommission. By following Level I and Level II damage accept/reject criteria set forth in this document, the material review (MR) process can be initiated. This nondestructive evaluation (NDE) method should be complemented with additional NDE methods to best understand the nature of the observed indications. When VI identification is followed by interrogation using effective complimentary NDE methods, the final Material Review Board (MRB) disposition of the component can be effectively achieved. Although elements of the visual inspection method are discussed and expressly required by the various range standards (KNPR 8715.3 and AFSPCMAN 91-710), emphasis should be placed on supporting the interpretation of recorded visual observations. Sound engineering practices and procedures should be applied with the interpretation of nondestructive testing (NDT) results when the residual strength data specific to a vessel design is incomplete or absent.

Tommy B. Yoder↗

NASA NTRS2003

Detecting Edges in Images by Use of Fuzzy Reasoning

A method of processing digital image data to detect edges includes the use of fuzzy reasoning. The method is completely adaptive and does not require any advance knowledge of an image. During initial processing of image data at a low level of abstraction, the nature of the data is indeterminate. Fuzzy reasoning is used in the present method because it affords an ability to construct useful abstractions from approximate, incomplete, and otherwise imperfect sets of data. Humans are able to make some sense of even unfamiliar objects that have imperfect high-level representations. It appears that to perceive unfamiliar objects or to perceive familiar objects in imperfect images, humans apply heuristic algorithms to understand the images

Dominguez, Jesus A.↗

DOE OSTI2019

Multi-scale topology optimization of multi-material structures with controllable geometric complexity — Applications to heat transfer problems

This study presents a topology optimization method to design assemblies of periodic cellular materials with controllable geometric complexity. The framework is based on a novel multi-scale and multi-material design model in which the structure, the layout of the material subdomains, and their micro-structures are optimized concurrently. To allow a tight control over its geometrical complexity, the layout at the macro-scale is described by level-set fields, parameterized by geometric primitives. The micro-scale geometry is represented through a density approach. A nonlinear programming algorithm drives the optimization process using design sensitivities computed by the discrete adjoint method. The proposed design framework is studied with heat transfer problems. Practical design problems, such as a heat sink and a thermal storage unit with phase change are discussed. The macro-scale analysis model relies on a generalized transient diffusion equation. At the micro-scale, homogenization is used to compute equivalent material properties. The numerical examples in this study show that the optimized multi-scale multi-material layouts outperform the corresponding mono-scale structures while maintaining geometric simplicity at the macro-scale level.

42 ENGINEERING↗

NASA NTRSDate not supplied

A Spectral Fingerprinting Method for Deriving Consistent Climate Data Records from Multiple Satellite IR Sounders

Deriving Climate Data Records (CDRs) from multiple IR sounders such AIRS on Aqua, CrIS on S-NPP, NOAA20 and JPSS-2, IASI on Metop A, B, and C requires from current operational products are challenging due to different radiative transfer models and retireavl algorithms used for processing level 2 data. We developed a Climate Fingerprinting Sounder Product (ClimFiSP) algorithm, which uses a single set of radiative kernels a robust spectral fingerprinting method to performs retrievals using spatiotemporally averaged L1 hyperspectral radiances directly. The ClimFiSP algorithm provides accurate data fusion CDR products from multiple satellite sensors. We have applied this method to both AIRS and CrIS (on SNPP and on NOAA 20) data and generated two decades climate data records for atmospheric temperature, water vapor, cloud, trace gases, and surface skin temperature. We plan to add IASI to the CDR data set in the future. The ClimFiSP is being transitioned to NASA GES DISC data center for public access.

Xu Liu↗

DOE OSTI2019

Efficient propagation of systematic uncertainties from calibration to analysis with the SnowStorm method in IceCube

Efficient treatment of systematic uncertainties that depend on a large number of nuisance parameters is a persistent difficulty in particle physics and astrophysics experiments. Where low-level effects are not amenable to simple parameterization or re-weighting, analyses often rely on discrete simulation sets to quantify the effects of nuisance parameters on key analysis observables. Such methods may become computationally untenable for analyses requiring high statistics Monte Carlo with a large number of nuisance degrees of freedom, especially in cases where these degrees of freedom parameterize the shape of a continuous distribution. In this paper we present a method for treating systematic uncertainties in a computationally efficient and comprehensive manner using a single simulation set with multiple and continuously varied nuisance parameters. Finally, this method is demonstrated for the case of the depth-dependent effective dust distribution within the IceCube Neutrino Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2000

Fixed-Order Mixed Norm Designs for Building Vibration Control

This study investigates the use of H2, mu-synthesis, and mixed H2/mu methods to construct full order controllers and optimized controllers of fixed dimensions. The benchmark problem definition is first extended to include uncertainty within the controller bandwidth in the form of parametric uncertainty representative of uncertainty in the natural frequencies of the design model. The sensitivity of H2 design to unmodeled dynamics and parametric uncertainty is evaluated for a range of controller levels of authority. Next, mu-synthesis methods are applied to design full order compensators that are robust to both unmodeled dynamics and to parametric uncertainty. Finally, a set of mixed H2/mu compensators are designed which are optimized for a fixed compensator dimension. These mixed norm designs recover the H2 design performance levels while providing the same levels of robust stability as the mu designs. It is shown that designing with the mixed norm approach permits higher levels of controller authority for which the H2 designs are destabilizing. The benchmark problem is that of an active tendon system. The controller designs are all based on the use of acceleration feedback.

Whorton, Mark S.↗

DOE OSTI2020

Description of Sensor Assignment Optimization Method as Deployed on a Multi-Node Cluster

Data analytic methods are being developed to address the problem of how to assign a sensor set in a nuclear facility such that a requisite level of process monitoring capability is realized and that the sensor set is sufficiently rich to determine the status of the individual sensors with respect to need for calibration. There is an awareness in the nuclear industry that data analytics combined with rich sensor sets represent a means to improve operations and reduce costs. In the industry the calibration problem has been previously approached as an empirical data-driven problem with several methods having been developed. However, the experience of the utilities over the past ten years with these methods indicates that the absence of physics-based information renders the data-driven approach less reliable. Complicating factors such as the inherent variability of operation (both equipment alignment and operating condition) can confound a pure data-driven approach while there are no rigorous guidelines for determining what constitutes an adequate sensor set. The solution under development to overcome these shortcomings supplements the data analytic method with process information in a so-called process-constrained data-analytic approach. Simple balance equations are written for generic components (e.g., mechanical pump, valve, and heat exchanger). These do not require a priori knowledge of process parameters, such as heat transfer coefficients or friction factors. All that is needed on the part of the utility user is to identify the components and how they are connected. This report describes the development of a parallel computing capability for determining the optimal sensor set. The optimal sensor set problem suffers from the curse of dimensionality. Computation time increases exponentially as the size of the system grows. To overcome this difficulty a pre-conditioner algorithm is developed to find an approximate solution close the actual solution. This serves as a seed for the full-blown algorithm and acts to constrain the space that must searched. The optimization algorithms are described and the implementation on a parallel computing platform is described. The application of the method to a use case we are solving in collaboration with our utility partner served to illustrate how the default sensor set in a nuclear plant may not provide sufficient coverage to infer sensor calibration status.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2018

Fundamental Vibrational Analyses of the HCN Monomer, Dimer and Associated Isotopologues

In this work we provide high level ab initio treatments of the structures, vibrational frequencies, and electronic energies of the HCN monomer and dimer systems along with several isotopologues. Additionally, the plethora of information related to this system within the literature is summarized and serves as a basis for comparison with the results of this paper. The geometry of the dimer and monomer are reported at the all electroncoupled-cluster singles, doubles, and perturbative triples level of theory [AE-CCSD(T)] with the correlation consistent quadruple-zeta quality basis sets with extra core functions (cc-pCVQZ) from Dunning. The theoretical geometries and electronic structures are further analyzed through the use of the Natural Bond Orbital (NBO) method and Natural Resonance Theory (NRT). At the AE-CCSD(T)/cc-pCVQZ level of theory, the full cubic with semi-diagonal quartic force field for nine dimer and four monomer isotopologues (the parent isotopologue along with 15N, 13C, and D derivatives) were obtained to treat the anharmonicity of the vibrations via second order vibrational perturbation theory (VPT2). Lastly, the enthalpy change associated with the formation of the dimer from two monomer units was determined using the focal point analysis. Computations including coupled-cluster through perturbative quadruples as well as basis sets up to six-zeta quality, including core functions (cc-pCVXZ, X=D,T,Q,5,6) were used to extrapolate to the AE-CCSDT(Q)/CBS energy associated with this hydrogen-bond forming process. After appending anharmonic zero-point vibrational, relativistic, and diagonal Born-Oppenheimer corrections, we report a value of -3.93 kcal mol -1 for the enthalpy of formation. To our knowledge, each set of results (geometries, vibrational frequencies, and energetics) reported in this study represents the highest-level and most reliable theoretical predictions reported for this system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Studies on Workability of Geopolymer Solidification - 18171

In Japan, intermediate-level radioactive waste, such as spent resin, is solidified and set in sealed containers without venting. To reduce the amount of hydrogen generated due to water radiolysis in stored intermediate-level waste, geopolymer solidification has been studied. Geopolymer solidification is an excellent method for conditioning intermediate-level waste from the viewpoint of its high thermal stability and low processing cost. However, it has a problem that the apparent viscosity increases drastically a few minutes after the onset of setting. This may cause some trouble in the process of pouring the geopolymer into drums and mixing it with waste. In this study, in order to control the workability of a geopolymer mixture, we investigated the effects of the compounding ratio and the kinds of chemical agents used. We confirmed that the geopolymer reaction was slowed down by reducing the amount of calcium and adding boron. The uniaxial compressive strength of the geopolymer was 3-80 MPa. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1976

Assessment of NDE Reliability Data

Twenty sets of relevant Nondestructive Evaluation (NDE) reliability data have been identified, collected, compiled, and categorized. A criterion for the selection of data for statistical analysis considerations has been formulated. A model to grade the quality and validity of the data sets has been developed. Data input formats, which record the pertinent parameters of the defect/specimen and inspection procedures, have been formulated for each NDE method. A comprehensive computer program has been written to calculate the probability of flaw detection at several confidence levels by the binomial distribution. This program also selects the desired data sets for pooling and tests the statistical pooling criteria before calculating the composite detection reliability. Probability of detection curves at 95 and 50 percent confidence levels have been plotted for individual sets of relevant data as well as for several sets of merged data with common sets of NDE parameters.

Yee, B. G. W.↗

DOE OSTI2024

Formally Verified ZTA Requirements for OT/ICS Environments with Isabelle/HOL

The clean energy transformation includes the integration of distributed energy resources with the power grid, which has led to a substantial increase in the complexity of power grids infrastructure and the underlying operational technology environment. Power grids infrastructure represents an operational technology environment that has become a system of systems, integrating heterogeneous devices which are both software-and hardware-intensive; as a result, there are increasing demands to exploit advances in the commodity of software-hardware infrastructures to improve energy systems requirements such as cybersecurity and resilience. In such a setting, system requirements at different levels mix, which leads to vulnerabilities and undesirable outcomes. The use of formal methods to characterize and prove system requirements removes ambiguity, increases automation, and provides high levels of assurance and reliability. In this paper, we contribute a methodology and a framework for the system-level verification of zero trust architecture requirements in operational technology environments. We define a formal specification for the core functionalities of operational technology environments, the corresponding invariants, and security proofs. Of particular note is our modular approach for the formal verification of asynchronous interactions in operational technology environments. The formal specification and the proofs have been mechanized using the interactive theorem proving environment Isabelle/HOL.

formal methods↗

DOE OSTI2024

On the Impact of Bus Dwelling on Macroscopic Fundamental Diagrams

Network macroscopic fundamental diagrams (MFDs) have recently been shown to exist in real-world urban traffic networks. When present, MFDs can be used to model traffic dynamics within an urban network by dividing the network into a set of spatially compact homogeneous regions and tracking the average level of congestion in each region. Existing analytical methods to estimate MFD mostly focus on the behavior of a single type of vehicle and do not capture the patterns of mixed traffic (e.g., cars and buses). The existence of buses matters since a bus will block the movements of other vehicles when it dwells at the bus stop. This paper proposes an analytical method to estimate the impact of bus dwelling on a network’s MFD based on the network’s geometric features, traffic control strategies, and bus operation parameters, and validates the performance of the proposed method using simulations based on microscopic traffic models. Comparisons of the analytical and simulation results show that the proposed analytical method can generally provide a good estimate of the lower bound and upper bound of the network’s MFD.

Xu, Guanhao↗

DOE OSTI2026

North Slope of Alaska XSAPR b1 Data Processing Report: May 2025-April 2026

The U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) User Facility operates an X-band Scanning ARM Precipitation Radar (XSAPR) at the North Slope of Alaska (NSA) atmospheric observatory in Utqiaġvik. This radar provides year-round dual-polarization measurements of snow and ice in an arctic environment. Over the past few years, there has been greater focus on the NSA XSAPR, and recent efforts by ARM staff have resulted in quality-controlled, b-level data (Rocque et al. 2026b). This report outlines the processing of the second set of b-level data released for this radar from May 2025 through April 2026 following similar methods from Rocque et al. (2026b).

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1990

Compiling global name-space programs for distributed execution

Distributed memory machines do not provide hardware support for a global address space. Thus programmers are forced to partition the data across the memories of the architecture and use explicit message passing to communicate data between processors. The compiler support required to allow programmers to express their algorithms using a global name-space is examined. A general method is presented for analysis of a high level source program and along with its translation to a set of independently executing tasks communicating via messages. If the compiler has enough information, this translation can be carried out at compile-time. Otherwise run-time code is generated to implement the required data movement. The analysis required in both situations is described and the performance of the generated code on the Intel iPSC/2 is presented.

Koelbel, Charles↗

NASA NTRS2010

Method of Obtaining High Resolution Intrinsic Wire Boom Damping Parameters for Multi-Body Dynamics Simulations

The goal of NASA's Magnetospheric MultiScale (MMS) mission is to understand magnetic reconnection with sensor measurements from four spinning satellites flown in a tight tetrahedron formation. Four of the six electric field sensors on each satellite are located at the end of 60- meter wire booms to increase measurement sensitivity in the spin plane and to minimize motion coupling from perturbations on the main body. A propulsion burn however, might induce boom oscillations that could impact science measurements if oscillations do not damp to values on the order of 0.1 degree in a timely fashion. Large damping time constants could also adversely affect flight dynamics and attitude control performance. In this paper, we will discuss the implementation of a high resolution method for calculating the boom's intrinsic damping, which was used in multi-body dynamics simulations. In summary, experimental data was obtained with a scaled-down boom, which was suspended as a pendulum in vacuum. Optical techniques were designed to accurately measure the natural decay of angular position and subsequently, data processing algorithms resulted in excellent spatial and temporal resolutions. This method was repeated in a parametric study for various lengths, root tensions and vacuum levels. For all data sets, regression models for damping were applied, including: nonlinear viscous, frequency-independent hysteretic, coulomb and some combination of them. Our data analysis and dynamics models have shown that the intrinsic damping for the baseline boom is insufficient, thereby forcing project management to explore mitigation strategies.

Yew, Alvin G.↗

NASA NTRS1994

Relativistic and Correlation Effects in CuH, AgH and AuH: Comparison of Various Relativistic Methods

The effects of relativity on the bond lengths, dissociation energies, and harmonic vibrational frequencies of the 1Epsilon(+) electronic ground states of the group IB hydrides CuH, AgH and AuH have been evaluated with a variety of ab initio methods. These properties were investigated with moderately-sized basis sets at the self-consistent field Hartree Fock (SCF HF) level and with second-order Moller-Plesset (MP2) perturbation theory for electron correlation. Comparisons were made between all-electron results using the nonrelativistic Hamiltonian, perturbation theory (PT) at first-order with only the one-electron non-fine structure terms of the Breit-Pauli Hamiltonian, the spin-free Douglas-Kroll (DK) transformed Dirac Hamiltonian and the untransformed Dirac Hamiltonian, and results using two sets of relativistic effective core potentials (RECPs). The expected trends of bond length decrease, dissociation energy increase and harmonic frequency increase with both relativity and correlation are found. Both sets of RECPs are shown to give good results, if accompanied by a reasonable basis set. The DK method is demonstrated to be an inexpensive, reliable approximation to the DHF method.

Collins, Charlene L.↗

NASA NTRS1994

Accuracy assessment of the large-scale dynamic ocean topography from TOPEX/POSEIDON altimetry

The quality of TOPEX/POSEIDON determinations of the global scale dynamic ocean topography have been assessed by determining mean topography solutions for successive 10-day repeat cycles and by examining the temporal changes in the sea surface topography to identify known features. The assessment is based on the analysis of TOPEX altimeter data cycles 1 through 36. Important errors in the tide model used to correct the altimeter data have been identified. The errors were reduced significantly by use of a new tide model derived with the TOPEX/POSEIDON measurements. Maps of the global 1-year mean topography, produced using four of the most accurate of the marine geoid, show that the largest error in the dynamic ocean topography show expected features, such as the known annual hemispherical sea surface rise and fall and the seasonal variability due to monsoon influence in the Indian Ocean. Changes in the sequence of 10-day topography maps show the development and propagation of an equatorial Kelvin wave in the Pacific beginning in December 1992 with a propagation velocity of approximately 3 m/s. The observations are consistent with observed changes in the equatorial trade winds, and with tide gauge and other in situ observations of the strengthening of the El Nino. Comparison of TOPEX-determine sea surface height at points near oceanic tide gauges shows agreement at the 4 cm root-mean-square (RMS) level over the tropical Pacific. The results show that the TOPEX altimeter data set can be used to map the ocean surface with a temporal resolution of 10 days and an accuracy which is insonsistent with traditional in situ methods for the determination of sea level variations.

Tapley, B. D.↗