Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set methods”

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,009 records · Page 56

Conotoxin Prediction: New Features to Increase Prediction Accuracy

Conotoxins are toxic, disulfide-bond-rich peptides from cone snail venom that target a wide range of receptors and ion channels with multiple pathophysiological effects. Conotoxins have extraordinary potential for medical therapeutics that include cancer, microbial infections, epilepsy, autoimmune diseases, neurological conditions, and cardiovascular disorders. Despite the potential for these compounds in novel therapeutic treatment development, the process of identifying and characterizing the toxicities of conotoxins is difficult, costly, and time-consuming. This challenge requires a series of diverse, complex, and labor-intensive biological, toxicological, and analytical techniques for effective characterization. While recent attempts, using machine learning based solely on primary amino acid sequences to predict biological toxins (e.g., conotoxins and animal venoms), have improved toxin identification, these methods are limited due to peptide conformational flexibility and the high frequency of cysteines present in toxin sequences. This results in an enumerable set of disulfide-bridged foldamers with different conformations of the same primary amino acid sequence that affect function and toxicity levels. Consequently, a given peptide may be toxic when its cysteine residues form a particular disulfide-bond pattern, while alternative bonding patterns (isoforms) or its reduced form (free cysteines with no disulfide bridges) may have little or no toxicological effects. Similarly, the same disulfide-bond pattern may be possible for other peptide sequences and result in different conformations that all exhibit varying toxicities to the same receptor or to different receptors. We present here new features, when combined with primary sequence features to train machine learning algorithms to predict conotoxins, that significantly increase prediction accuracy.

collisional cross section↗

Comparison of wind-farm control strategies under realistic offshore wind conditions: wake quantities of interest

Wind-farm control strategies aim to increase the efficiency, and therefore lower the levelized cost of energy, of wind farms. This is done by using turbine settings such as the yaw angle, blade pitch angles, or generator torque to manipulate the wake that negatively affects downstream turbines in the farm. Two inherently different wind-farm control methods have been identified in the literature: wake steering (WS) and active wake mixing (AWM). As one of two companion papers focused on understanding practical aspects of these two wind-farm control strategies using large-eddy simulation (LES), we below analyze the wake quantities of interest for a single wind turbine performing WS and AWM, while the companion article (Frederik et al., 2025) focuses on turbine quantities of interest including power and structural loads for the same computational setup and also includes two-turbine arrays with full and partial wake overlap. The simulations, which are based in the LES solver AMR-Wind, are tailored to have inflow conditions representative of measurements from a site off the East Coast of the US, including with strong veer and low turbulence. The turbine, which is modeled in OpenFAST and coupled to the LES, is the IEA 15 MW, an open-source offshore design. After presenting an overview of the wake recovery for the different wake-control cases, the analysis probes the fluid-dynamic causes for the different performance of the arrays reported in the companion article by examining control volumes around the wakes and the budget of the mean-flow kinetic energy (MKE) within these volumes. In the high veer environment considered, the MKE recovery is dominated by mean convection, and this is shown to especially benefit the WS strategy when a neighboring turbine is directly downstream: there is ≈65 % more available power for a downstream turbine than in the baseline case, and this power is gained primarily through mean convection on the left-tip and top-tip faces of the control volume. However, the case with imperfect knowledge of the exact wind direction favors the pulse-type AWM strategy, largely because of ≈9 % increased turbulent entrainment from aloft versus the baseline that could be related to an apparent resistance to skewing in the pulsed wake. The general reduced effectiveness of helix-type and other individual-pitch-based AWM strategies for inflow with high veer and low turbulence as reported in the companion paper is due, in part, to low magnitudes of phase-averaged turbulent entrainment. Two main findings of this study are thus that veer has a significant impact on the effectiveness of different wake-control strategies and that pulse-type AWM may be a useful strategy when the objective is power maximization in realistic, offshore flow environments with imperfect knowledge of the exact wake overlap position on the downstream turbine.

17 WIND ENERGY↗

Flow structure in continuous flow electrophoresis chambers

There are at least two ways that hydrodynamic processes can limit continiuous flow electrophoresis. One arises from the sensitivity of the flow to small temerature gradients, especially at low flow rates and power levels. This sensitivity can be suppressed, at least in principle, by providing a carefully tailored, stabilizing temperature gradient in the cooling system that surrounds the flow channel. At higher power levels another limitation arises due to a restructuring of the main flow. This restructuring is caused by buoyancy, which is in turn affected by the electro-osmotic crossflow. Approximate solutions to appropriate partial differential equations have been computed by finite difference methods. One set of results is described here to illustrate the strong coupling between the structure of the main (axial) flow and the electro-osmotic flow.

Deiber, J. A.↗

Fan Acoustic Flight Effects on the PAA & ASN Flight Test

Accurate and reliable prediction of noise emitted from turbofan aircraft engines is an ongoing challenge. For continued progress, it is necessary to assess the performance of the most recent prediction methods in the NASA Aircraft NOise Prediction Program (ANOPP), and to examine the effects of flight condition on fan noise. The assessment is based on an extensive acoustic dataset acquired as part of the Propulsion Airframe Aeroacoustics and Aircraft System Noise (PAA&ASN) Flight Test conducted on a Boeing 787-10 aircraft. Engine operating condition was varied to study the characteristics of fan noise in flight. In a subset of the flight test matrix, the influence of the aft duct acoustic liner was studied by covering the lined area with aluminum tape to simulate a hardwall duct. The combination of test data at different flight conditions, along with rigorous predictions of jet noise, is used in this paper to extract fan noise from the measured total aircraft noise levels. The measured fan noise levels are compared with those predicted using the latest methods for fan source noise, liner attenuation, and PAA effect prediction. Detailed analyses are undertaken to understand agreement between data and various aspects of the prediction models, including spectral shape, directivity, and response to throttle setting. An additional investigation was undertaken to study the effect of aircraft speed on fan noise, as current models do not directly account for this parameter.

Aircraft System Noise↗

Parameterizing the Variability and Uncertainty of Wind and Solar in CEMs

We present current and improved methods for estimating the capacity value and curtailment impacts from variable generation (VG) in capacity expansion models (CEMs). The ideal calculation of these variability metrics is through an explicit co-optimized investment-dispatch model using multiple years of VG and load data. Because of data and computational limitations, existing CEMs typically approximate these metrics using a subset of all hours from a single year and/or using statistical methods, which often do not capture the tail-event impacts or the broader set of interactions between VG, storage, and conventional generators. In our proposed new methods, we use hourly generation and load values across all hours of the year to characterize the (1) contribution of VG to system capacity during high load hours, (2) the curtailment level of VG, and (3) the reduction in VG curtailment due to storage and shutdown of select thermal generators. Using CEM model outputs from a preceding model solve period, we apply these methods to exogenously calculate capacity value and curtailment metrics for the subsequent model solve period. Preliminary results suggest that these hourly methods offer improved capacity value and curtailment representations of VG in the CEM from existing approximation methods without additional computational burdens.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Assimilation of scatterometer winds into surface pressure fields using a variational method

A variational formulation was used to assimilate Seasat-A scatterometer (SASS) surface wind measurements near and during a severe storm in the North Atlantic into conventional National Meteorological Center sea level pressure fields. An estimate of the relative vorticity at every point on a grid was calculated using each of these two data sets. A solution to a modified geostrophic stream function is found subject to the constraints that (1) the relative vorticities calculated from the data agree as closely as possible with the relative vorticities from the variational solution, and that (2) the average kinetic energy is a minimum. Results are obtained which support the idea that averaged satellite data can be treated as synoptic data. Direct substitution rather than a time-weighted insertion made from SASS winds generally resulted in more accurate pressure analyses. In addition, this relatively simple model provides surface pressure fields which agree extremely well with surface truth and the results of other investigators who required additional sources of input data into more complex models. It will be possible to obtain improved wind field maps from future scatterometer pressure fields in mid-latitudes.

Harlan, J., Jr.↗

Reduction of Convection in Closed Tube Vapor Transport Experiments

The primary objective of this effort was to develop a method for suppressing convective flows during the growth of mercurous chloride crystals by vapor transport in closed tubes to levels approaching those obtained in the microgravity environment. Mercurous chloride was chosen because it is a technologically interesting acoustical optical material whose optical properties are believed to be affected by convective flows. Since the Grashof number scales as the cube of the smallest dimension in the flow system, reduction of the size scale can be extremely effective in reducing unwanted convective flows. However, since materials of practical interest must be grown at least on the cm scale, reduction of the overall growth system is not feasible. But if the region just above the growing crystal could be restricted to a few mm, considerable reduction in flow velocity would result. By suspending an effusive barrier in the growth ampoule just above the growth interface, it should be possible to reduce the convective velocity in this vicinity to levels approaching flows in microgravity. If successful, this growth technique will offer a screening test for proposed space experiments that involve vapor transport to see if reduction of convection will result in improved material and will set a new standard against which the improvements obtained in microgravity may be judged. In addition, it may provide an improved method for preparing materials on Earth whose growth is affected adversely by convection. If the properties of this material can be improved there is a potential commercial interest from Brimrose Inc., who has agreed to fabricate and test devices from the crystals we have grown. This report describes the development of the growth facility, the purification processes developed for preparing the starting material, and the results from growth experiments with and without the effusive baffle. Mercurous chloride turned out to be a more difficult material to deal with than originally anticipated. At growth temperatures, it is extremely sensitive to practically any impurity which causes it to form oxychlorides and/or to decompose into elemental mercury and bichloride of mercury. We were unable to find a suitable method for protecting the magnetic material used to suspend the effusion barrier from the attack of mercurous chloride vapor. Although we were successful in growing single crystals of mercurous chloride without the effusion baffle, they exhibited severe microcracking which we attribute to wall-induced thermal stresses. This leads us to believe that uncontrolled convection may not be the most important problem in the development of this material and a new growth process was attempted that eliminates the wall-induced stress. Unfortunately, the grant ran out before this new method could be adequately tested.

Naumann, R. J.↗

Characteristic archaebacterial 16S rRNA oligonucleotides

A method of analyzing 16S rRNA catalog data has been developed in which groupings at various taxonomic levels can be characterized in terms of specific "signature" oligonucleotides. This approach provides an alternative means for evaluating higher order branching possibilities and can be used to assess the phylogenetic position of isolates that are poorly placed by the usual clustering procedures. This signature approach has been applied to forty archaebacterial catalogs and every oligonucleotide with significant signature value has been identified. Sets of specific oligonucleotides were identified for every major group on a dendrogram produced by cluster analysis procedures. Signatures that would establish between group relationships were also sought and found. In the case of the Methanobacteriaceae the clustering methods suggest a specific relationship to the Methanococcaceae. This inclusion is in fact supported by six strong signature oligonucleotides. However there are also significant numbers of signature oligonucleotides supporting a specific relationship of the Methanobacteriaceae to either the Halobacteriaceae or the Methanomicrobiaceae. Thus the placement of the Methanobacteriaceae is less certain than the usual dendrograms imply. The signature approach also was used to assess the phylogenetic position of Thermoplasma acidophilum which is found to be more closely related to the methanogen/halophile Division than to the sulfur dependent Division of the archaebacteria. This does not imply however that Thermoplasma acidophilum is properly regarded as being in the methanogen/halophile Division.

NASA Discipline Exobiology↗

R2 relaxometry based MR imaging for estimation of liver iron content: A comparison between two methods

Purpose: To compare the reproducibility and accuracy of R2-relaxometry MRI for estimation of liver iron concentration (LIC) between in-house analysis and FDA-approved commercially available third party results. Methods: All MR studies were performed on a 1.5T scanner. Multi-echo spin-echo scans with a fixed TR and increasing TE values of 6 ms, 9 ms, 12 ms, 15 ms, and 18 ms (spaced at 3 ms intervals) were used. Post-processing of the images to calculate mean relaxivity, R2, included drawing of regions of interest to include the whole liver on mid-slice. The relationship between liver R2 values and estimated LIC calculated with in-house analysis and values reported by an external company (FerriScan{sup ®}, Resonance Health, Australia) were assessed with correlation coefficients and Bland–Altman difference plots. Continuous variables are presented as mean ± standard deviation. Significance was set at p value < 0.05. Results: 474 studies from 175 patients were included in the study (mean age 10.4 ± 4.2 years (range 1–18 years); 254 studies from girls, 220 studies from boys). LIC ranged from 0.6 to 43 mg/g dry tissue, covering a broad range from normal levels to extremely high iron levels. Linearity between proprietary and in-house methods was excellent across the observed range for R2 (31.5 to 334.8 s{sup −1}); showing a correlation coefficient of r = 0.87, p < 0.001. Bland–Altman R2 difference plot between the two methods shows a mean bias of + 21.5 s{sup −1} (range − 47.0 to + 90.0 s{sup −1} between two standard deviations). LIC reported by FerriScan{sup ®} compared with LIC estimated in-house with R2 as reported by FerriScan{sup ®} agreed strongly, (r = 1.0, p < 0.001). Conclusion: R2 relaxometry MR imaging for liver iron concentration estimation is reproducible between proprietary FDA-approved commercial software and in-house analysis methods.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Accurate complex scaling of three dimensional numerical potentials

The complex scaling method, which consists in continuing spatial coordinates into the complex plane, is a well-established method that allows to compute resonant eigenfunctions of the time-independent Schroedinger operator. Whenever it is desirable to apply the complex scaling to investigate resonances in physical systems defined on numerical discrete grids, the most direct approach relies on the application of a similarity transformation to the original, unscaled Hamiltonian. We show that such an approach can be conveniently implemented in the Daubechies wavelet basis set, featuring a very promising level of generality, high accuracy, and no need for artificial convergence parameters. Complex scaling of three dimensional numerical potentials can be efficiently and accurately performed. By carrying out an illustrative resonant state computation in the case of a one-dimensional model potential, we then show that our wavelet-based approach may disclose new exciting opportunities in the field of computational non-Hermitian quantum mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Composite vibrational spectroscopy of the group 12 difluorides: ZnF{sub 2}, CdF{sub 2}, and HgF{sub 2}

The vibrational spectra of group 12 difluorides, MF{sub 2} (M = Zn, Cd, Hg), were investigated via coupled cluster singles, doubles, and perturbative triples, CCSD(T), including core correlation, with a series of correlation consistent basis sets ranging in size from triple-zeta through quintuple-zeta quality, which were then extrapolated to the complete basis set (CBS) limit using a variety of extrapolation procedures. The explicitly correlated coupled cluster method, CCSD(T)-F12b, was employed as well. Although exhibiting quite different convergence behavior, the F12b method yielded the CBS limit estimates closely matching more computationally expensive conventional CBS extrapolations. The convergence with respect to basis set size was examined for the contributions entering into composite vibrational spectroscopy, including those from higher-order correlation accounted for through the CCSDT(Q) level of theory, second-order spin-orbit coupling effects assessed within four-component and two-component relativistic formalisms, and vibrational anharmonicity evaluated via a perturbative treatment. Overall, the composite results are in excellent agreement with available experimental values, except for the CdF{sub 2} bond-stretching frequencies compared to spectral assignments proposed in a matrix isolation infrared and Raman study of cadmium difluoride vapor species [Loewenschuss et al., J. Chem. Phys. 50, 2502 (1969); Givan and Loewenschuss, J. Chem. Phys. 72, 3809 (1980)]. These assignments are called into question in the light of the composite results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Search for neutron decay into an antineutrino and a neutral kaon in 0.401 megaton-years exposure of Super-Kamiokande

We searched for bound neutron decay via 𝑛 → $\bar{𝜈}$ +𝐾 0 predicted by the grand unified theories in 0.401 Mton·years exposure of all pure water phases in the Super-Kamiokande detector. About 4.4 times more data than in the previous search have been analyzed by a new method including a spectrum fit to kaon invariant mass distributions. No significant data excess has been observed in the signal regions. As a result of this analysis, we set a lower limit of 7.8 × 10 32 years on the neutron lifetime at a 90% confidence level.

Cherenkov detectors↗

A Data-Driven Methodology for Probabilistic Wind Power Ramp Forecasting

With increasing wind penetration, wind power ramps (WPRs) are currently drawing great attention to balancing authorities, since these wind ramps largely affect power system operations. To help better manage and dispatch the wind power, this paper develops a data-driven probabilistic WPR forecasting (p-WPRF) method based on a large number of simulated scenarios. A machine learning technique is first adopted to forecast the basic wind power forecasting scenario and produce the historical forecasting errors. To accurately model the distribution of wind power forecasting errors, a generalized Gaussian mixture model is developed and the cumulative distribution function (CDF) is also analytically deduced. The inverse transform method based on the CDF is used to generate a large number of forecasting error scenarios. An optimized swinging door algorithm is adopted to extract all the WPRs from the complete set of wind power forecasting scenarios. The p-WPRF is generated based on all generated scenarios under different weather and time conditions. Numerical simulations on publicly available wind power data show that the developed p-WPRF method can predict WPRs with a high level of reliability and accuracy.

17 WIND ENERGY↗

A 2D ion chamber array audit of wedged and asymmetric fields in an inhomogeneous lung phantom

Purpose: The Australian Clinical Dosimetry Service (ACDS) has implemented a new method of a nonreference condition Level II type dosimetric audit of radiotherapy services to increase measurement accuracy and patient safety within Australia. The aim of this work is to describe the methodology, tolerances, and outcomes from the new audit. Methods: The ACDS Level II audit measures the dose delivered in 2D planes using an ionization chamber based array positioned at multiple depths. Measurements are made in rectilinear homogeneous and inhomogeneous phantoms composed of slabs of solid water and lung. Computer generated computed tomography data sets of the rectilinear phantoms are supplied to the facility prior to audit for planning of a range of cases including reference fields, asymmetric fields, and wedged fields. The audit assesses 3D planning with 6 MV photons with a static (zero degree) gantry. Scoring is performed using local dose differences between the planned and measured dose within 80% of the field width. The overall audit result is determined by the maximum dose difference over all scoring points, cases, and planes. Pass (Optimal Level) is defined as maximum dose difference ≤3.3%, Pass (Action Level) is ≤5.0%, and Fail (Out of Tolerance) is >5.0%. Results: At close of 2013, the ACDS had performed 24 Level II audits. 63% of the audits passed, 33% failed, and the remaining audit was not assessable. Of the 15 audits that passed, 3 were at Pass (Action Level). The high fail rate is largely due to a systemic issue with modeling asymmetric 60° wedges which caused a delivered overdose of 5%–8%. Conclusions: The ACDS has implemented a nonreference condition Level II type audit, based on ion chamber 2D array measurements in an inhomogeneous slab phantom. The powerful diagnostic ability of this audit has allowed the ACDS to rigorously test the treatment planning systems implemented in Australian radiotherapy facilities. Recommendations from audits have led to facilities modifying clinical practice and changing planning protocols.

60 APPLIED LIFE SCIENCES↗

SU-E-T-780: Use Robustness Optimization (RO) Method to Improve the Planning Efficiency for Pencil Beam Scanning Cranial Spinal Irradiation

Purpose: We evaluate the feasibility of using robustness optimization (RO) function to improve the planning efficiency of pencil beam scanning (PBS) craniospinal irradiation (CSI) with gradient matching technique. Methods: A CSI patient was planned with 2 lateral brain fields and 4 posterior fields to cover the entire spine to maximal field of 24 cm × 20 cm on a compact PBS gantry, ProteusONE. CSI plans were generated using traditional volumetric gradient dose optimization (VGDO) and robustness optimization (RO) method respectively. In traditional VGDO, besides the sectioned spine target volumes, gradient volume (GV) were generated as 4 equally spaced structures e.g. 80%, 60%, 40%, and 20% of prescription dose. In RO method, only sectioned spine target volumes with an overlap of 4cm were created. In the robustness optimization settings, 5mm uncertainty in superior and inferior direction was defined for auto gradient optimization. Dosimetric metrics of conformity number (CN), homogeneity index (HI), and maximal cord doses were compared in Raystation version 4.6.100.6. Results: In VGDO method, total 16 GV structures and five 100% dose level target structures were contoured compared to total 5 target structures in RO method which saves 30 min in contour. With the same PTV coverage (95% volume receive 30.6Gy prescription dose), maximum cord dose is 32.64Gy in VGDO and 31.94Gy in RO. HI is 1.03 and 1.04 for VGDO and RO respectively. CN is 0.93 and 0.94 for VGDO and RO respectively. Conclusions: The dosimetric comparison demonstrated both methods are equivalent in terms of plan quality. With robust optimization for CSI gradient matching, it efficiently reduces the amount of planning target contour structure by factor of 4 and thus improves the planning efficiency especially for 4 or more gradient junction area.

60 APPLIED LIFE SCIENCES↗

BDEM (Discrete-element-simulator for high-solids granular flows) [SWR-22-72]

BDEM is a discrete element method based simulation tool developed specifically for modeling high-solids granular flows that include polydispersity, heat-transfer, moving boundaries and chemistry. Our solver provides facilities for simulating spherical/non-spherical particles with modified contact and friction models in complex dynamic geometries defined using level-sets or triangulated files. The solver is developed on top of NREL's open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Simulations at the scale of several millions to billion particles have been performed using this software on large scale computing resources. This software can be applied to non-reacting solids dominant flows in silos, hoppers and screw conveyors as well as in high temperature reacting systems such as screw kilns and auger reactors.

Sitaraman, Hariswaran↗