Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Space radiation incident on SATS missions

A special orbital radiation study was conducted in order to evaluate mission encountered energetic particle fluxes. This information is to be supplied to the project subsystem engineers for their guidance in designing flight hardware to withstand the expected radiation levels. Flux calculations were performed for a set of 20 nominal trajectories placed at several altitudes and inclinations. Temporal variations in the ambient electron environment were considered and partially accounted for. Magnetic field calculations were performed with a current field model, extrapolated to the tentative SATS launch epoch with linear time terms. Orbital flux integrations ware performed with the latest proton and electron environment models, using new computational methods. The results are presented in graphical and tabular form. Estimates of energetic solar proton fluxes are given for a one year mission at selected integral energies ranging from 10 to 100 Mev, calculated for a year of maximum solar activity during the next solar cycle.

Stassinopoulos, E. G.↗
NASA NTRS1992

Simulated families: A test for different methods of family identification

A set of families generated in fictitious impact events (leading to a wide range of 'structure' in the orbital element space have been superimposed to various backgrounds of different densities in order to investigate the efficiency and the limitations of the methods used by Zappala et al. (1990) and by Bendjoya et al. (1990) for identifying asteroid families. In addition, an evaluation of the expected interlopers at different significance levels and the possibility of improving the definition of the level of maximum significant of a given family were analyzed.

Bendjoya, Philippe↗
NASA NTRS1993

Flight calibration assessment of HiRAP accelerometer data

A flight derived method of calibrating the High Resolution Accelerometer Package (HiRAP) flight data has been developed and is discussed for Shuttle Orbiter missions STS-35 and STS-40. These two mission data sets have been analyzed using ground calibration factors and flight derived calibration factors. This flight technique evolved early in the flight program when it was recognized that ground calibration factors are insufficient to determine absolute low-acceleration levels. The application of flight calibration factors to the data sets from these missions produced calibrated acceleration levels within an accuracy of less than +/- 1.5 microgravity of zero during a time in the flight when the acceleration level was known to be less than 1.0 microgravity. This analysis further confirms the theory that flight calibrations are required in order to obtain the absolute measurement of low-frequency, low-acceleration flight signals.

Blanchard, Robert C.↗
NASA NTRS1993

Stochastic dynamics of baroclinic waves

The maintenance of variance and attendant heat flux in linear, forced dissipative baroclinic shear flows subject to stochastic excitation is examined. The baroclinic problem, is intrinsically nonnormal and its stochastic dynamics is found to differ significantly from the more familiar stochastic dynamics of normal systems. When the shear is sufficiently great in comparison to dissipative effects, stochastic excitation supports highly enhanced variance levels in these nonnormal systems compared to variance levels supproted by the same forcing and dissipation in related normal systems. The eddy variance and associated heat flux are found to arise in response to transient amplification of a subset of forcing functions that obtain energy from the mean flow and project this energy on a distinct subset of response functions (EOFs) that are in turn distinct from the set of normal modes of the system. A method for obtaining the dominant forcing and response functions as well as the distribution of heat flux for a given flow is described.

Farrell, Brian F.↗
NASA NTRS2018

Trim Flight Conditions for a Low-Boom Aircraft Design Under Uncertainty

The purpose of this paper is to investigate the noise generation of a low-boom aircraft design in flight trim conditions under uncertainty. The deflection of control surfaces to maintain a trimmed flight state has the potential to change the perceived loudness at the ground. Furthermore, the uncertainties associated with the control surface deflections can complicate the overall uncertainty quantification. Incorporation of the uncertainties in the prediction of perceived sound levels during the design phase can lead to improved robustness. In this paper, a brief review of low-boom flight trim research is presented. Realistic flight trim conditions requiring control surface deflection are integrated into the current research efforts for uncertainty quantification and vehicle design. In addition, a generalized set of procedures for the characterization of uncertainties in flight trim conditions are introduced. In a case study of the application of these procedures, a 5 decibel average difference in the perceived level of loudness was found between clean (no deflections) and trimmed configurations. Also, uncertainties attributable to control surface deflection were found to account for, on average, over 50% of the total near field uncertainty. Uncertainty discretization methods implemented were able to give more insight into the overall global variances.

Phillips, Ben D↗
DOE OSTI2012

New approach to creation of geometrical module for nuclear reactor neutron transport computer simulation analysis

This paper presents the new approach to creation of geometrical module for nuclear reactor neutron transport computer simulation analysis so called the differential cross method. It is elaborated for detecting boards between physical zones. It is proposed to use GMSH open source mesh editor extended by some features: a special option and a special kind of mesh (cubic background mesh).This method is aimed into Monte Carlo Method as well as for deterministic neutron transport methods. Special attention is attended for reactor core composed of a set of material zones with complicate geometrical boundaries. The idea of this approach is described. In general case method works for 3-D space. Algorithm of creation of the geometrical module is given. 2-D neutron transport benchmark-test for RBMK reactor cluster cell is described. It demonstrates the ability of this approach to provide flexible definition of geometrical meshing with preservation of curved surface or any level of heterogeneity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

On the accuracy of Hipparcos using binary stars as a calibration tool

Stellar binary systems, specifically those that present the most accurate available orbital elements, are a reliable tool to test the accuracy of astrometric observations. We selected all 35 binaries with these characteristics. Our objective is to provide standard uncertainties for the positions and parallaxes measured by Hipparcos relative to this trustworthy set, as well as to check supposed correlations between several parameters (measurement residuals, positions, magnitudes, and parallaxes). In addition, using the high-confidence subset of visual–spectroscopic binaries, we implemented a validation test of the Hipparcos trigonometric parallaxes of binary systems that allowed the evaluation of their reliability. Standard and non-standard statistical analysis techniques were applied in order to achieve well-founded conclusions. In particular, errors-in-variables models such as the total least-squares method were used to validate Hipparcos parallaxes by comparison with those obtained directly from the orbital elements. Previously, we executed Thompson's τ technique in order to detect suspected outliers in the data. Furthermore, several statistical hypothesis tests were carried out to verify if our results were statistically significant. A statistically significant trend indicating larger Hipparcos angular separations with respect to the reference values in 5.2 ± 1.4 mas was found at the 10{sup −8} significance level. Uncertainties in the polar coordinates θ and ρ of 1.°8 and 6.3 mas, respectively, were estimated for the Hipparcos observations of binary systems. We also verified that the parallaxes of binary systems measured in this mission are absolutely compatible with the set of orbital parallaxes obtained from the most accurate orbits at least at the 95% confidence level. This methodology allows us to better estimate the accuracy of Hipparcos observations of binary systems. Indeed, further application to the data collected by Gaia should yield a standard procedure to compare both data sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Linearization Techniques for Controlled Piecewise Deterministic Markov Processes; Application to Zubov's Method

We aim at characterizing domains of attraction for controlled piecewise deterministic processes using an occupational measure formulation and Zubov's approach. Firstly, we provide linear programming (primal and dual) formulations of discounted, infinite horizon control problems for PDMPs. These formulations involve an infinite-dimensional set of probability measures and are obtained using viscosity solutions theory. Secondly, these tools allow to construct stabilizing measures and to avoid the assumption of stability under concatenation for controls. The domain of controllability is then characterized as some level set of a convenient solution of the associated Hamilton-Jacobi integral-differential equation. The theoretical results are applied to PDMPs associated to stochastic gene networks. Explicit computations are given for Cook's model for gene expression.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Water interactions with hydrophobic groups: Assessment and recalibration of semiempirical molecular orbital methods

In this work, we present a study of the ability of different semiempirical methods to describe intermolecular interactions in water solution. In particular, we focus on methods based on the Neglect of Diatomic Differential Overlap approximation. Significant improvements of these methods have been reported in the literature in the past years regarding the description of non-covalent interactions. In particular, a broad range of methodologies has been developed to deal with the properties of hydrogen-bonded systems, with varying degrees of success. In contrast, the interactions between water and a molecule containing hydrophobic groups have been little analyzed. Indeed, by considering the potential energy surfaces obtained using different semiempirical Hamiltonians for the intermolecular interactions of model systems, we found that none of the available methods provides an entirely satisfactory description of both hydrophobic and hydrophilic interactions in water. In addition, a vibrational analysis carried out in a model system for these interactions, a methane clathrate cluster, showed that some recent methods cannot be used to carry out studies of vibrational properties. Following a procedure established in our group [M. I. Bernal-Uruchurtu, M. T. C. Martins-Costa, C. Millot, and M. F. Ruiz-López, J. Comput. Chem. 21, 572 (2000); W. Harb, M. I. Bernal-Uruchurtu, and M. F. Ruiz-López, Theor. Chem. Acc. 112, 204 (2004)], we developed new parameters for the core-core interaction terms based on fitting potential energy curves obtained at the MP2 level for our model system. We investigated the transferability of the new parameters to describe a system, having both hydrophilic and hydrophobic groups, interacting with water. We found that only by introducing two different sets of parameters for hydrophilic and hydrophobic hydrogen atom types we are able to match the features of the ab initio calculated properties. Once this assumption is made, a good agreement with the MP2 reference is achieved. The results reported in this work provide therefore a direction for future developments of semiempirical approaches that are still required to investigate chemical processes in biomolecules and in large disordered systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Two-Level Least Squares Methods in Krylov Subspaces

Two-level least squares acceleration approaches are applied to the Chebyshev acceleration method and the restarted conjugate residual method in solving systems of linear algebraic equations with sparse unsymmetric coefficient matrices arising from finite volume or finite element approximations of boundary-value problems on irregular grids. Application of the proposed idea to other iterative restarted processes also is considered. The efficiency of the algorithms suggested is investigated numerically on a set of model Dirichlet problems for the convection-diffusion equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: A Comparison Between Energy - and Variance-based Approaches

An efficient implementation of geometrical derivatives at the Hartree–Fock (HF) and current-density functional theory (CDFT) levels is presented for the study of molecular structure in strong magnetic fields. The required integral derivatives are constructed using a hybrid McMurchie–Davidson and Rys quadrature approach, which combines the amenability of the former to the evaluation of derivative integrals with the efficiency of the latter for basis sets with high angular momentum. In addition to its application to evaluating derivatives of four-center integrals, this approach is also applied to gradients using the resolution-of-the-identity approximation, enabling efficient optimization of molecular structure for many-electron systems under a strong magnetic field. The CDFT contributions have been implemented for a wide range of density functionals up to and including the meta-GGA level with current-density dependent contributions and (range-separated) hybrids for the first time. Illustrative applications are presented to the OH and benzene molecules, revealing the rich and complex chemistry induced by the presence of an external magnetic field. Challenges for geometry optimization in strong fields are highlighted, along with the requirement for careful analysis of the resulting electronic structure at each stationary point. The importance of correlation effects is examined by comparison of results at the HF and CDFT levels. The present implementation of molecular gradients at the CDFT level provides a cost-effective approach to the study of molecular structure under strong magnetic fields, opening up many new possibilities for the study of chemistry in this regime.

Comparison↗
DOE OSTI2016

REPEATABILITY OF SPITZER/IRAC EXOPLANETARY ECLIPSES WITH INDEPENDENT COMPONENT ANALYSIS

The research of effective and reliable detrending methods for Spitzer data is of paramount importance for the characterization of exoplanetary atmospheres. To date, the totality of exoplanetary observations in the mid- and far-infrared, at wavelengths >3 μm, have been taken with Spitzer. In some cases, in past years, repeated observations and multiple reanalyses of the same data sets led to discrepant results, raising questions about the accuracy and reproducibility of such measurements. Morello et al. (2014, 2015) proposed a blind-source separation method based on the Independent Component Analysis of pixel time series (pixel-ICA) to analyze InfraRed Array Camera (IRAC) data, obtaining coherent results when applied to repeated transit observations previously debated in the literature. Here we introduce a variant to the pixel-ICA through the use of wavelet transform, wavelet pixel-ICA, which extends its applicability to low-signal-to-noise-ratio cases. We describe the method and discuss the results obtained over 12 eclipses of the exoplanet XO3b observed during the “Warm Spitzer” era in the 4.5 μm band. The final results are reported, in part, also in Ingalls et al. (2016), together with results obtained with other detrending methods, and over 10 synthetic eclipses that were analyzed for the “IRAC Data Challenge 2015.” Our results are consistent within 1σ with the ones reported in Wong et al. (2014) and with most of the results reported in Ingalls et al. (2016), which appeared on arXiv while this paper was under review. Based on many statistical tests discussed in Ingalls et al. (2016), the wavelet pixel-ICA method performs as well as or better than other state-of-art methods recently developed by other teams to analyze Spitzer/IRAC data, and, in particular, it appears to be the most repeatable and the most reliable, while reaching the photon noise limit, at least for the particular data set analyzed. Another strength of the ICA approach is its highest objectivity, as it does not use prior information about the instrument systematics, making it a promising method to analyze data from other observatories. The self-consistency of individual measurements of eclipse depth and phase curve slope over a span of more than three years proves the stability of Warm Spitzer/IRAC photometry within the error bars, at the level of 1 part in 10{sup 4} in stellar flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Offshore wind energy forecasting sensitivity to sea surface temperature input in the Mid-Atlantic

Abstract. As offshore wind farm development expands, accurate wind resource forecasting over the ocean is needed. One important yet relatively unexplored aspect of offshore wind resource assessment is the role of sea surface temperature (SST). Models are generally forced with reanalysis data sets, which employ daily SST products. Compared with observations, significant variations in SSTs that occur on finer timescales are often not captured. Consequently, shorter-lived events such as sea breezes and low-level jets (among others), which are influenced by SSTs, may not be correctly represented in model results. The use of hourly SST products may improve the forecasting of these events. In this study, we examine the sensitivity of model output from the Weather Research and Forecasting model (WRF) 4.2.1 to different SST products. We first evaluate three different data sets: the Multiscale Ultrahigh Resolution (MUR25) SST analysis, a daily, 0.25∘ × 0.25∘ resolution product; the Operational Sea Surface Temperature and Ice Analysis (OSTIA), a daily, 0.054∘ × 0.054∘ resolution product; and SSTs from the Geostationary Operational Environmental Satellite 16 (GOES-16), an hourly, 0.02∘ × 0.02∘ resolution product. GOES-16 is not processed at the same level as OSTIA and MUR25; therefore, the product requires gap-filling using an interpolation method to create a complete map with no missing data points. OSTIA and GOES-16 SSTs validate markedly better against buoy observations than MUR25, so these two products are selected for use with model simulations, while MUR25 is at this point removed from consideration. We run the model for June and July of 2020 and find that for this time period, in the Mid-Atlantic, although OSTIA SSTs overall validate better against in situ observations taken via a buoy array in the area, the two products result in comparable hub-height (140 m) wind characterization performance on monthly timescales. Additionally, during hours-long flagged events (< 30 h each) that show statistically significant wind speed deviations between the two simulations, both simulations once again demonstrate similar validation performance (differences in bias, earth mover's distance, correlation, and root mean square error on the order of 10−1 or less), with GOES-16 winds validating nominally better than OSTIA winds. With a more refined GOES-16 product, which has been not only gap-filled but also assimilated with in situ SST measurements in the region, it is likely that hub-height winds characterized by GOES-16-informed simulations would definitively validate better than those informed by OSTIA SSTs.

17 WIND ENERGY↗
NASA NTRS1980

Synthetic aperture radar range - Azimuth ambiguity design and constraints

Problems concerning the design of a system for mapping a planetary surface with a synthetic aperture radar (SAR) are considered. Given an ambiguity level, resolution, and swath width, the problems are related to the determination of optimum antenna apertures and the most suitable pulse repetition frequency (PRF). From the set of normalized azimuth ambiguity ratio curves, the designer can arrive at the azimuth antenna length, and from the sets of normalized range ambiguity ratio curves, he can arrive at the range aperture length or pulse repetition frequency. A procedure based on this design method is shown in an example. The normalized curves provide results for a SAR using a uniformly or cosine weighted rectangular antenna aperture.

Mehlis, J. G.↗
NASA NTRS1988

Extraction of topography from side-looking satellite systems - A case study with SPOT simulation data

A simulated data set was used to evaluate techniques for extracting topography from side-looking satellite systems for an area of northwest Washington state. A negative transparency orthophotoquad was digitized at a spacing of 85 microns, resulting in an equivalent ground distance of 9.86 m between pixels and a radiometric resolution of 256 levels. A bilinear interpolation was performed on digital elevation model data to generate elevation data at a 9.86-m resolution. The nominal orbital characteristics and geometry of the SPOT satellite were convoluted with the data to produce simulated panchromatic HRV digital stereo imagery for three different orbital paths and techniques for reconstructing topographic data were developed. Analyses with the simulated HRV data and other data sets show that the method is effective.

Ungar, Stephen G.↗
NASA NTRS1994

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simuulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rarefaction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the after body sting are emphasizes in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗
DOE OSTI2012

Modern energy density functional and the current status of the equation of state of nuclear matter

We first describe a method, based on the simulated annealing approach, for determining a modern energy density functional within the Skyrme Hartree-Fock (HF) theory by carrying out a fit to extensive set of experimental data with additional constraints on the Skyrme parameters. Next, we review the HF-based random phase approximation (RPA) approach for calculating properties of giant resonances. We then present results of calculations for the centroid energies of giant resonances within the HF-based RPA and discuss the current status of the equation of state of nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗