Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation 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,513 records · Page 84

Moments-based interface reconstruction, remap and advection

Here, we present a new moment-of-fluid (MOF 2 ) interface reconstruction method. It uses the zeroth, first, and second moments of the fragment of material inside a cell of the mesh to reconstruct a convex material polygon or a union of convex polygons that approximate the respective material fragment. The new method requires information about the material moments only for the cell under consideration. The MOF 2 method allows to exactly reproduce several convex shapes: corners, filaments, and some concave shapes: cell-complements to corners and filaments. Interface reconstruction is formulated as a local (for each cell), non-linear, equality constrained optimization problem, which does not require additional communication and allows for an efficient parallel implementation. We present an extensive set of test problems, both for interface reconstruction on a single cell, and for reconstruction of a variety of shapes on a variety of meshes. We describe how to perform two-material advection using the MOF 2 method and present the results for the classical advection tests. We also show the examples of material interface remapping needed in the framework of multi-material arbitrary Lagrangian-Eulerian methods, and give a brief description of a procedure that can be used to update the material moments on the Lagrangian stage of those methods.

97 MATHEMATICS AND COMPUTING↗

A Global, Multi-Resolution Approach to Regional Ocean Modeling

In this collaborative research project between Pennsylvania State University, Colorado State University and Florida State University, we mainly focused on developing multi-resolution algorithms which are suitable to regional ocean modeling. We developed hybrid implicit and explicit adaptive multirate time integration method to solve systems of time-dependent equations that present two ignorantly different scales. We studied the e ects of spatial simplicial meshes on the stability and the conditioning of fully discrete approximations. We also studies adaptive nite element method (AFEM) based upon the Centroidal Voronoi Tessellation (CVT) and superconvergent gradient recovery. Some of these techniques are now being used by geoscientists(such as those at LANL).

58 GEOSCIENCES↗

Moments-based interface reconstruction, remap and advection

We present a new moment-of-fluid (MOF 2 ) interface reconstruction method. It uses the zeroth, first, and second moments of the fragment of material inside a cell of the mesh to reconstruct a convex material polygon or a union of convex polygons that approximate the respective material fragment. The new method requires information about the material moments only for the cell under consideration. Furthermore, the MOF 2 method allows to exactly reproduce several convex shapes: corners, filaments, and some concave shapes: cell-complements to corners and filaments.

97 MATHEMATICS AND COMPUTING↗

Fast Automatic Knot Placement Method for Accurate B-spline Curve Fitting

The choice of knot vector has immense influence on the resulting accuracy of a B-spline approximation of a curve. However, despite the significance of this problem and the various solutions that were proposed in the literature, optimizing the number and placement of knots remains a difficult task. This paper presents a novel method for the approximation of a curve by a B-spline of arbitrary order, which automatically determines a knot vector that achieves high approximation quality. Additionally, at the core of our approach is a feature function that characterizes the amount and spatial distribution of geometric details in the input curve by estimating its derivatives. Knots are then selected in such a way as to evenly distribute the feature contents across their intervals. A comparison to the state of the art for a wide variety of curves shows that our method is faster and achieves more accurate reconstruction results, while typically reducing the number of necessary knots.

97 MATHEMATICS AND COMPUTING↗

Hartree-Fock calculation of the differential photoionization cross sections of small Li clusters

Cross sections and angular distribution parameters for the single-photon ionization of all electron orbitals of Li{sub 2−8} are systematically computed in a broad interval of the photoelectron kinetic energies for the energetically most stable geometry of each cluster. Calculations of the partial photoelectron continuum waves in clusters are carried out by the single center method within the Hartree-Fock approximation. We study photoionization cross sections per one electron and analyze in some details general trends in the photoionization of inner and outer shells with respect to the size and geometry of a cluster. The present differential cross sections computed for Li{sub 2} are in a good agreement with the available theoretical data, whereas those computed for Li{sub 3−8} clusters can be considered as theoretical predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Varying constants quantum cosmology

We discuss minisuperspace models within the framework of varying physical constants theories including Λ-term. In particular, we consider the varying speed of light (VSL) theory and varying gravitational constant theory (VG) using the specific ansätze for the variability of constants: c(a) = c{sub 0} a{sup n} and G(a)=G{sub 0} a{sup q}. We find that most of the varying c and G minisuperspace potentials are of the tunneling type which allows to use WKB approximation of quantum mechanics. Using this method we show that the probability of tunneling of the universe ''from nothing'' (a=0) to a Friedmann geometry with the scale factor a{sub t} is large for growing c models and is strongly suppressed for diminishing c models. As for G varying, the probability of tunneling is large for G diminishing, while it is small for G increasing. In general, both varying c and G change the probability of tunneling in comparison to the standard matter content (cosmological term, dust, radiation) universe models.

79 ASTRONOMY AND ASTROPHYSICS↗

Industrializing Offshore Wind Power with Serial Assembly and Lower-cost Deployment - Final Report

A team of engineers and contractors has developed a method to move offshore wind installation toward lower cost, faster deployment, and lower environmental impact. A combination of methods, some incremental and some breaks from past practice, interact to yield multiple improvements. Three designs were evaluated based on detailed engineering: 1) a 5 MW turbine on a jacket with pin piles (base case), 2) a 10 MW turbine on a conventional jacket with pin piles, assembled at sea, and 3) a 10 MW turbine on tripod jacket with suction buckets (caissons) and with complete turbine assembly on-shore. The larger turbine, assembly ashore, and the use of suction buckets together substantially reduce capital cost of offshore wind projects. Notable capital cost reductions are: changing from 5 MW to 10 MW turbine, a 31% capital cost reduction, and assembly on land then single-piece install at sea an additional 9% capital cost reduction. An estimated Design 4) estimates further cost reduction when equipment and processes of Design 3) are optimized, rather than adapted to existing equipment and process. Cost of energy for each of the four Designs are also calculated, yielding approximately the same percentage reductions. The methods of Design 3) analyzed here include accepted structures such as suction buckets used in new ways, innovations conceived but previously without engineering and economic validation, combined with new methods not previously proposed. Analysis of Designs 2) and 3) are based on extensive engineering calculations and detailed cost estimates. All design methods can be done with existing equipment, including lift equipment, ports and ships (except that design 4 assumes a more optimized ship). The design team consists of experienced offshore structure designers, heavy lift engineers, wind turbine designers, vessel operators, and marine construction contractors. Comparing the methods based on criteria of cost and deployment speed, the study selected the third design. That design is, in brief: a conventional turbine and tubular tower is mounted on a tripod jacket, in turn atop three suction buckets. Blades are mounted on the tower, not on the hub. The entire structure is built in port, from the bottom up, then assembled structures are queued in the port for deployment. During weather windows, the fully-assembled structures are lifted off the quay, lashed to the vessel, and transported to the deployment site. The vessel analyzed is a shear leg crane vessel with dynamic positioning like the existing Gulliver, or it could be a US-built crane barge. On site, the entire structure is lowered to the bottom by the crane vessel, then pumping of the suction buckets is managed by smaller service vessels. Blades are lifted into place by small winches operated by workers in the nacelle without lift vessel support. Advantages of the selected design include: cost and time at sea of the expensive lift vessel are significantly reduced; no jack up vessel is required; the weather window required for each installation is shorter; turbine structure construction is continuous with a queue feeding the weather-dependent installation process; pre-installation geotechnical work is faster and less expensive; there are no sound impacts on marine mammals, thus minimal spotting and no work stoppage Industrializing Offshore Wind Power 6 of 96 9 for mammal passage; the entire structure can be removed for decommissioning or major repairs; the method has been validated for current turbines up to 10 MW, and a calculation using simple scaling shows it usable up to 20 MW turbines.

17 WIND ENERGY↗

Whole-remnant and maximum-voxel SPECT/CT dosimetry in {sup 131}I-NaI treatments of differentiated thyroid cancer

Purpose: To investigate the possible differences between SPECT/CT based whole-remnant and maximum-voxel dosimetry in patients receiving radio-iodine ablation treatment of differentiated thyroid cancer (DTC). Methods: Eighteen DTC patients were administered 1.11 GBq of {sup 131}I-NaI after near-total thyroidectomy and rhTSH stimulation. Two patients had two remnants, so in total dosimetry was performed for 20 sites. Three SPECT/CT scans were performed for each patient at 1, 2, and 3–7 days after administration. The activity, the remnant mass, and the maximum-voxel activity were determined from these images and from a recovery-coefficient curve derived from experimental phantom measurements. The cumulated activity was estimated using trapezoidal-exponential integration. Finally, the absorbed dose was calculated using S-values for unit-density spheres in whole-remnant dosimetry and S-values for voxels in maximum-voxel dosimetry. Results: The mean absorbed dose obtained from whole-remnant dosimetry was 40 Gy (range 2–176 Gy) and from maximum-voxel dosimetry 34 Gy (range 2–145 Gy). For any given patient, the activity concentrations for each of the three time-points were approximately the same for the two methods. The effective half-lives varied (R = 0.865), mainly due to discrepancies in estimation of the longer effective half-lives. On average, absorbed doses obtained from whole-remnant dosimetry were 1.2 ± 0.2 (1 SD) higher than for maximum-voxel dosimetry, mainly due to differences in the S-values. The method-related differences were however small in comparison to the wide range of absorbed doses obtained in patients. Conclusions: Simple and consistent procedures for SPECT/CT based whole-volume and maximum-voxel dosimetry have been described, both based on experimentally determined recovery coefficients. Generally the results from the two approaches are consistent, although there is a small, systematic difference in the absorbed dose due to differences in the S-values, and some variability due to differences in the estimated effective half-lives, especially when the effective half-life is long. Irrespective of the method used, the patient absorbed doses obtained span over two orders of magnitude.

60 APPLIED LIFE SCIENCES↗

Bayesian Monte-Carlo Evaluation Framework of Differential and Integral Data [Slides]

This presentation discusses the Bayesian Monte Carlo evaluation framework of differential and integral data, specifically in relation to two methods that are based on Bayes theorem, which assumes that all PDFs are normal (Gaussian) and that they use a linear approximation for IBEs. The presentation states that basic components of the MC evaluation framework include the computation of random MC ensemble from ENDF File 32, the simulation of IBEs and R-matrix cross section compared to experimental data, and the computation of weighted averages. It also states that the application to U-233 indicates deviation from the conventional linear approximation. In conclusion, evaluation framework will require MCMC method e.g., Metropolis-Hastings (M.-H.).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Capacitance–voltage modeling of mid-wavelength infrared nBn detectors

Capacitance–voltage measurements are a powerful technique to determine doping profiles of semiconductor pn junctions and Schottky barrier diodes. The measurements were recently extended to III-V-based mid-wavelength nBn infrared detectors, and absorber doping densities have been extracted using the widely accepted Schottky approximation, where the potential drop across the device is assumed to be across the depleting absorber layer. However, this approach is limited to when the absorber region of the nBn is under high reverse bias and thus is only able to extract the absorber region doping profile. Here, we introduce a semi-analytical model that is capable of extracting barrier dopant polarity, doping concentration, and thickness, as well as contact and absorber layer doping concentrations, all from a capacitance–voltage measurement. Rather than solely considering the potential drop across the depleting layers, it considers the potential drop across the accumulating layer as well. This negative charge accumulation occurs for the contact and absorber layers in the case of reverse and forward biases, respectively. This allows for a single model to be applied to a capacitance–voltage curve at forward and reverse biases and it can provide regions of bias where the absorber transitions from depletion to accumulation. We compare the agreement of the semianalytical model with modeling results from commercially available finite element method software and experimental capacitance–voltage data. Finally, we show that the method is consistent with the Schottky approximation of extracting absorber doping densities at high reverse bias and discuss the model's limitations.

Physics↗

The backward phase flow method for the Eulerian finite time Lyapunov exponent computations

We propose a simple Eulerian approach to compute the moderate to long time flow map for approximating the Lyapunov exponent of a (periodic or aperiodic) dynamical system. The idea is to generalize a recently proposed backward phase flow method which is specially designed for long time level set propagation. Unlike the original phase flow method or the backward phase flow method, which is applicable only to autonomous systems, the current approach can also be applied to any time-dependent (periodic or aperiodic) flow. We will discuss the stability of the proposed method. Numerical examples will be given to demonstrate the effectiveness of the algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the evaporation of solar dark matter: spin-independent effective operators

As a part of the effort to investigate the implications of dark matter (DM)-nucleon effective interactions on the solar DM detection, in this paper we focus on the evaporation of the solar DM for a set of the DM-nucleon spin-independent (SI) effective operators. In order to put the evaluation of the evaporation rate on a more reliable ground, we calculate the non-thermal distribution of the solar DM using the Monte Carlo methods, rather than adopting the Maxwellian approximation. We then specify relevant signal parameter spaces for the solar DM detection for various SI effective operators. Based on the analysis, we determine the minimum DM masses for which the DM-nucleon coupling strengths can be probed from the solar neutrino observations. As an interesting application, our investigation also shows that evaporation effect can not be neglectd in a recent proposal aiming to solve the solar abundance problem by invoking the momentum-dependent asymmetric DM in the Sun.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A fresh look at the photoelectron spectrum of bromobenzene: A third-order non-Dyson electron propagator study

The valence-shell ionization spectrum of bromobenzene, as a representative halogen substituted aromatic, was studied using the non-Dyson third-order algebraic-diagrammatic construction [nD-ADC(3)] approximation for the electron propagator. This method, also referred to as IP-ADC(3), was implemented as a part of the Q-Chem program and enables large-scale calculations of the ionization spectra, where the computational effort scales as n{sup 5} with respect to the number of molecular orbitals n. The IP-ADC(3) scheme is ideally suited for investigating low-lying ionization transitions, so fresh insight could be gained into the cationic state manifold of bromobenzene. In particular, the present IP-ADC(3) calculations with the cc-pVTZ basis reveal a whole class of low-lying low-intensity two-hole-one-particle (2h-1p) doublet and quartet states, which are relevant to various photoionization processes. The good qualitative agreement between the theoretical spectral profile for the valence-shell ionization transitions generated with the smaller cc-pVDZ basis set and the experimental photoelectron spectrum measured at a photon energy of 80 eV on the PLÉIADES beamline at the Soleil synchrotron radiation source allowed all the main features to be assigned. Some theoretical aspects of the ionization energy calculations concerning the use of various approximation schemes and basis sets are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic Properties of KPb{sub 2}Cl{sub 5} and RbPb{sub 2}Br{sub 5} Doped with Er{sup 3+} and Yb{sup 3+}

Alkali halide crystals doped with Er{sup 3+} exhibit several properties that make them promising for practical use in optical fibers and converters of the emission spectrum. In this work, we studied the luminescence spectra of crystals as a function of the excitation radiation wavelength and the temperature dependences of the luminescence intensity in the up-conversion mechanism of excitation of the most intense lines (530, 550, and 650 nm) and upon excitation at 980 nm. In KPb{sub 2}Cl{sub 5}:Er{sup 3+}, RbPb{sub 2}Br{sub 5}:Er{sup 3+}, and KPb{sub 2}Cl{sub 5}:Er{sup 3+} + Yb{sup 3+} crystals, the luminescence kinetics were measured upon excitation at 313 and 365 nm, and the theoretical estimates of lifetimes were made using the Judd–Ofelt theory. The results can be used in the development of laser-based systems and for estimating the degree of approximation of calculations by the Judd–Ofelt method for these crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Monte Carlo methods and their analysis for Coulomb collisions in multicomponent plasmas

Highlights: •A general approach to Monte Carlo methods for multicomponent plasmas is proposed. •We show numerical tests for the two-component (electrons and ions) case. •An optimal choice of parameters for speeding up the computations is discussed. •A rigorous estimate of the error of approximation is proved. -- Abstract: A general approach to Monte Carlo methods for Coulomb collisions is proposed. Its key idea is an approximation of Landau–Fokker–Planck equations by Boltzmann equations of quasi-Maxwellian kind. It means that the total collision frequency for the corresponding Boltzmann equation does not depend on the velocities. This allows to make the simulation process very simple since the collision pairs can be chosen arbitrarily, without restriction. It is shown that this approach includes the well-known methods of Takizuka and Abe (1977) [12] and Nanbu (1997) as particular cases, and generalizes the approach of Bobylev and Nanbu (2000). The numerical scheme of this paper is simpler than the schemes by Takizuka and Abe [12] and by Nanbu. We derive it for the general case of multicomponent plasmas and show some numerical tests for the two-component (electrons and ions) case. An optimal choice of parameters for speeding up the computations is also discussed. It is also proved that the order of approximation is not worse than O(√(ε)), where ε is a parameter of approximation being equivalent to the time step Δt in earlier methods. A similar estimate is obtained for the methods of Takizuka and Abe and Nanbu.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hydrogen adsorption in metal-organic frameworks: The role of nuclear quantum effects

The role of nuclear quantum effects on the adsorption of molecular hydrogen in metal-organic frameworks (MOFs) has been investigated on grounds of Grand-Canonical Quantized Liquid Density-Functional Theory (GC-QLDFT) calculations. For this purpose, we have carefully validated classical H{sub 2}-host interaction potentials that are obtained by fitting Born-Oppenheimer ab initio reference data. The hydrogen adsorption has first been assessed classically using Liquid Density-Functional Theory and the Grand-Canonical Monte Carlo methods. The results have been compared against the semi-classical treatment of quantum effects by applying the Feynman-Hibbs correction to the Born-Oppenheimer-derived potentials, and by explicit treatment within the GC-QLDFT. The results are compared with experimental data and indicate pronounced quantum and possibly many-particle effects. After validation calculations have been carried out for IRMOF-1 (MOF-5), GC-QLDFT is applied to study the adsorption of H{sub 2} in a series of MOFs, including IRMOF-4, -6, -8, -9, -10, -12, -14, -16, -18, and MOF-177. Finally, we discuss the evolution of the H{sub 2} quantum fluid with increasing pressure and lowering temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vibrational and Thermal Properties of Oxyanionic Crystals

The vibrational and thermal properties of dolomite and alkali chlorates and perchlorates were studied in the gradient approximation of density functional theory using the method of a linear combination of atomic orbitals (LCAO). Long-wave vibration frequencies, IR and Raman spectra, and mode Gruneisen parameters were calculated. Equation-of-state parameters, thermodynamic potentials, entropy, heat capacity, and thermal expansion coefficient were also determined. The thermal expansion coefficient of dolomite was established to be much lower than for chlorates and perchlorates. The temperature dependence of the heat capacity at T > 200 K was shown to be generally governed by intramolecular vibrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗