Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization 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

Efficient Near-Field to Mid-Field Sonic Boom Propagation Using a High-Order Space Marching Method

An efficient strategy for propagating sonic boom signatures from a near-field Computational Fluid Dynamics (CFD) solution to the mid-field is presented. The method is based on a high-order accurate finite-difference discretization of the 3D Euler equations on a specially designed curvilinear grid and a single sweep space marching solution algorithm. The new approach leads to more than a factor of two reduction in overall computational resources compared to the current method used to propagate near-field sonic booms to the ground. Accuracy and efficiency of the near-field to mid-field process is demonstrated using a selection of test cases from the AIAA Sonic Boom Prediction Workshops. Azimuthal dependence of nonlinear wave propagation from the near-field to mid-field is analyzed along with its effects on the ground level noise.

Housman, Jeffrey A.↗

Ice Accretion with Varying Surface Tension

During an icing encounter of an aircraft in flight, super-cooled water droplets impinging on an airfoil may splash before freezing. This paper reports tests performed to determine if this effect is significant and uses the results to develop an improved scaling method for use in icing test facilities. Simple laboratory tests showed that drops splash on impact at the Reynolds and Weber numbers typical of icing encounters. Further confirmation of droplet splash came from icing tests performed in the NaSA Lewis Icing Research Tunnel (IRT) with a surfactant added to the spray water to reduce the surface tension. The resulting ice shapes were significantly different from those formed when no surfactant was added to the water. These results suggested that the droplet Weber number must be kept constant to properly scale icing test conditions. Finally, the paper presents a Weber-number-based scaling method and reports results from scaling tests in the IRT in which model size was reduced up to a factor of 3. Scale and reference ice shapes are shown which confirm the effectiveness of this new scaling method.

Bilanin, Alan J.↗

Newton like: Minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗

Newton-like minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessing the technical aspects of ocean-alkalinity-enhancement approaches

Ocean alkalinity enhancement (OAE) is an emerging strategy that aims to mitigate climate change by increasing the alkalinity of seawater. This approach involves increasing the alkalinity of the ocean to enhance its capacity to absorb and store carbon dioxide (CO 2 ) from the atmosphere. This chapter presents an overview of the technical aspects associated with the full range of OAE methods being pursued and discusses implications for undertaking research on these approaches. Various methods have been developed to implement OAE, including the direct injection of alkaline liquid into the surface ocean; dispersal of alkaline particles from ships, platforms, or pipes; the addition of minerals to coastal environments; and the electrochemical removal of acid from seawater. Each method has its advantages and challenges, such as scalability, cost effectiveness, and potential environmental impacts. The choice of technique may depend on factors such as regional oceanographic conditions, alkalinity source availability, and engineering feasibility. This chapter considers electrochemical methods, the accelerated weathering of limestone, ocean liming, the creation of hydrated carbonates, and the addition of minerals to coastal environments. In each case, the technical aspects of the technologies are considered, and implications for best-practice research are drawn. The environmental and social impacts of OAE will likely depend on the specific technology and the local context in which it is deployed. Therefore, it is essential that the technical feasibility of OAE is undertaken in parallel with, and informed by, wider impact assessments. While OAE shows promise as a potential climate change mitigation strategy, it is essential to acknowledge its limitations and uncertainties. Further research and development are needed to understand the long-term effects, optimize techniques, and address potential unintended consequences. OAE should be viewed as complementary to extensive emission reductions, and its feasibility may be improved if it is operated using energy and supply chains with minimal CO 2 emissions.

54 ENVIRONMENTAL SCIENCES↗

Design and Analysis Tools for Supersonic Inlets

Computational tools are being developed for the design and analysis of supersonic inlets. The objective is to update existing tools and provide design and low-order aerodynamic analysis capability for advanced inlet concepts. The Inlet Tools effort includes aspects of creating an electronic database of inlet design information, a document describing inlet design and analysis methods, a geometry model for describing the shape of inlets, and computer tools that implement the geometry model and methods. The geometry model has a set of basic inlet shapes that include pitot, two-dimensional, axisymmetric, and stream-traced inlet shapes. The inlet model divides the inlet flow field into parts that facilitate the design and analysis methods. The inlet geometry model constructs the inlet surfaces through the generation and transformation of planar entities based on key inlet design factors. Future efforts will focus on developing the inlet geometry model, the inlet design and analysis methods, a Fortran 95 code to implement the model and methods. Other computational platforms, such as Java, will also be explored.

Slater, John W.↗

Photometry on Structured Backgrounds: Local Pixel-wise Infilling by Regression

Photometric pipelines struggle to estimate both the flux and flux uncertainty for stars in the presence of structured backgrounds such as filaments or clouds. However, it is exactly stars in these complex regions that are critical to understanding star formation and the structure of the interstellar medium. We develop a method, similar to Gaussian process regression, which we term local pixel-wise infilling (LPI). Using a local covariance estimate, we predict the background behind each star and the uncertainty of that prediction in order to improve estimates of flux and flux uncertainty. We show the validity of our model on synthetic data and real dust fields. We further demonstrate that the method is stable even in the crowded field limit. While we focus on optical-IR photometry, this method is not restricted to those wavelengths. We apply this technique to the 34 billion detections in the second data release of the Dark Energy Camera Plane Survey. In addition to removing many >3σ outliers and improving uncertainty estimates by a factor of ~2–3 on nebulous fields, we also show that our method is well behaved on uncrowded fields. The entirely post-processing nature of our implementation of LPI photometry allows it to easily improve the flux and flux uncertainty estimates of past as well as future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

The Trojan Horse Method: A Nuclear Physics Tool for Astrophysics

The Trojan Horse Method (THM) represents an indirect path to determine the bare nucleus astrophysical S-factor for reactions among charged particles at astrophysical energies. This is achieved by measuring the quasi-free cross section of a suitable three-body process. The method is also suited to study neutron-induced reactions, especially in the case of radioactive ion beams. A comprehensive review of the theoretical as well as experimental features behind the THM is presented here. An overview is given of some recent applications to demonstrate the method's practical use for reactions that have a great impact on selected astrophysical scenarios.

Physics↗

Towards Integrating Data Quality Assessments and Radiometer Uncertainty for Determining the Expanded Uncertainty of Three-Component Solar Radiation Measurements

Accurate solar irradiance data are fundamental for determining the design and performance characteristics of photovoltaic systems. The uncertainty of solar irradiance measurements depends on many factors including radiometer design, calibration, installation, maintenance, and operational environment. The key contributors to this uncertainty can be classified as the measurement uncertainty of a particular radiometer and the operational uncertainty determined for the time of measurement. Radiometer measurement uncertainty estimates (U R ) can be based on well-established methods used as part of the radiometer calibration process. Estimates of operational uncertainties (U o ) require consideration of additional site-specific factors that affect data quality. A method is needed for establishing the accuracy of solar irradiance data by integrating an existing data quality process and measurement uncertainty estimates for specific radiometers. An algorithm has been developed to integrate data quality analyses and measurement uncertainty estimates for three-component solar irradiance data: global horizontal (total hemispheric) irradiance, direct normal (beam) irradiance, and diffuse horizontal (sky) irradiance collected at one- to 60-minute intervals. The algorithm has been tested using one-minute irradiance measurements. The goal of the project is to distribute a user-friendly software package based on the new algorithm.

data integrity↗

On Variable Geometric Factor Systems for Top-Hat Electrostatic Space Plasma Analyzers

Even in the relatively small region of space that is the Earth's magnetosphere, ion and electron fluxes can vary by several orders of magnitude. Top-hat electrostatic analyzers currently do not possess the dynamic range required to sample plasma under all conditions. The purpose of this study was to compare, through computer simulation, three new electrostatic methods that would allow the sensitivity of a sensor to be varied through control of its geometric factor (GF) (much like an aperture on a camera). The methods studied were inner filter plates, split hemispherical analyzer (SHA) and top-cap electrode. This is the first discussion of the filter plate concept and also the first study where all three systems are studied within a common analyzer design, so that their relative merits could be fairly compared. Filter plates were found to have the important advantage that they facilitate the reduction in instrument sensitivity whilst keeping all other instrument parameters constant. However, it was discovered that filter plates have numerous disadvantages that make such a system impracticable for a top-hat electrostatic analyzer. It was found that both the top-cap electrode and SHA are promising variable geometric factor system (VGFS) concepts for implementation into a top-hat electrostatic analyzer, each with distinct advantages over the other.

Kataria, Dhiren O.↗

Mapping Mangrove Extent and Change: A Globally Applicable Approach

This study demonstrates a globally applicable method for monitoring mangrove forest extent at high spatial resolution. A 2010 mangrove baseline was classified for 16 study areas using a combination of ALOS PALSAR and Landsat composite imagery within a random forests classifier. A novel map-to-image change method was used to detect annual and decadal changes in extent using ALOS PALSAR/JERS-1 imagery. The map-to-image method presented makes fewer assumptions of the data than existing methods, is less sensitive to variation between scenes due to environmental factors (e.g., tide or soil moisture) and is able to automatically identify a change threshold. Change maps were derived from the 2010 baseline to 1996 using JERS-1 SAR and to 2007, 2008 and 2009 using ALOS PALSAR. This study demonstrated results for 16 known hotspots of mangrove change distributed globally, with a total mangrove area of 2,529,760 ha. The method was demonstrated to have accuracies consistently in excess of 90% (overall accuracy: 92.2–93.3%, kappa: 0.86) for mapping baseline extent. The accuracies of the change maps were more variable and were dependent upon the time period between images and number of change features. Total change from 1996 to 2010 was 204,850 ha (127,990 ha gain, 76,860 ha loss), with the highest gains observed in French Guiana (15,570 ha) and the highest losses observed in East Kalimantan, Indonesia (23,003 ha). Changes in mangrove extent were the consequence of both natural and anthropogenic drivers, yielding net increases or decreases in extent dependent upon the study site. These updated maps are of importance to the mangrove research community, particularly as the continual updating of the baseline with currently available and anticipated spaceborne sensors. It is recommended that mangrove baselines are updated on at least a 5-year interval to suit the requirements of policy makers.

Thomas, Nathan↗

Disentangling error structures of precipitation datasets using decision trees

Characterizing error structures in precipitation products not only facilitates their proper applications for scientific and practical purposes but also helps improve their retrieval algorithms and processing methods. Despite the fact that multiple precipitation products have been assessed in the literature, factors that affect their error structures remain inadequately addressed. By interpreting 60 binary decision trees, this study disentangles the error characteristics of precipitation products in terms of their spatiotemporal patterns and geographical factors. Three independent precipitation products - two satellite-based and one reanalysis datasets: the Integrated Multi-satellitE Retrievals for GPM (Global Precipitation Measurement) late run (IMERG-L), Soil Moisture to Rain-Advanced SCATterometer (SM2RAIN-ASCAT), and the Modern-Era Retrospective analysis for Research and Applications, Version 2 uncorrected precipitation output (MERRA2-UC), are evaluated across the contiguous United States from 2010 to 2019. Here, the ground-based Stage IV precipitation dataset is used as the ground truth. Results indicate that the MERRA2-UC outperforms the IMERG-L and SM2RAIN-ASCAT with higher accuracy and more stable interannual patterns for the analysis period. Decision trees cross-assess three spatiotemporal factors and find that the underestimation of MERRA2-UC occurs in the east of the Rocky Mountains, and SM2RAIN-ASCAT underestimates precipitation over high latitudes, especially in winter. Additionally, the decision tree method ascribes system errors to nine different geographical characteristics, of which the distance to the coast, soil type, and DEM are the three dominant features. On the other hand, the land cover type, topography position index, and aspect are three relatively weak factors.

54 ENVIRONMENTAL SCIENCES↗

Multi-Area Model-Free State Estimation via Distributed Tensor Decomposition

This paper proposes a model-free method for distribution system state estimation based on tensor completion using canonical polyadic decomposition. In particular, we consider a setting where the network is divided into multiple areas. The measured physical quantities at buses located in the same area are processed by an area controller. A third-order tensor is constructed to collect these measured quantities. The measurements are analyzed locally to recover the full state information of the network. A closed-form iterative algorithm based on the alternating direction method of multipliers is developed to obtain the low-rank factors of the whole network state tensor where information exchange happens only between neighboring areas. To demonstrate the efficacy of the developed algorithm, numerical simulations are carried out using an IEEE test system.

alternating direction method of multipliers↗

Large-eddy simulation of non-vaporizing sprays using the spectral-element method

Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.

97 MATHEMATICS AND COMPUTING↗

Computing a Stability Spectrum by Use of the HHT

The Hilbert-Huang transform (HHT) is part of the mathematical basis of a method of calculating a stability spectrum. This method can be regarded as an extended and improved version of a prior HHT-based method of calculating a damping spectrum. In the prior method, information on positive damping (which leads to stability) and negative damping (which leads to instability) becomes mixed into a single squared damping loss factor. Hence, there is no way to distinguish between stability and instability by examining a damping spectrum. In contrast, in the present stabilityspectrum method, information on the mathematical sign of the damping is retained, making it possible to identify reThe Hilbert-Huang transform (HHT) is part of the mathematical basis of a method of calculating a stability spectrum. This method can be regarded as an extended and improved version of a prior HHT-based method of calculating a damping spectrum. In the prior method, information on positive damping (which leads to stability) and negative damping (which leads to instability) becomes mixed into a single squared damping loss factor. Hence, there is no way to distinguish between stability and instability by examining a damping spectrum. In contrast, in the present stabilityspectrum method, information on the mathematical sign of the damping is retained, making it possible to identify regions of instability in a spectrum. This method is expected to be especially useful for analyzing vibration-test data for the purpose of predicting vibrational instabilities in structures (e.g., flutter in airplane wings).

Huang, Norden E.↗

The HUNTER Dynamic Human Reliability Analysis Tool: Development of a Module for Performance Shaping Factors

Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) is a framework to support dynamic human reliability analysis (HRA) in communication with a variety of methods and tools. In this paper, how we have developed one of the HUNTER modules, the individual module for evaluating performance shaping factors (PSFs), is introduced. The PSF refers to any factor that influences human performance such as workload or complexity. It has been used for highlighting human errors and adjusting the error probabilities in the existing HRA. We consider the eight PSFs suggested in the Standardized Plant Analysis Risk-HRA (SPAR-H) method, which is the representative HRA method widely used in the nuclear field. To support our dynamic modeling using the eight SPAR-H PSFs, we reviewed human performance literature and developed data-based mathematical models to rate and quantify PSFs in the context of dynamic HRA. We also design the individual module composed of the two functions: 1) the PSF qualification function that automatically or manually evaluates a PSF level, and 2) the PSF quantification function that dynamically or statically determines the PSF multiplier values and integrate them to adjust human error probabilities (HEPs). How each function works with the SPAR-H PSFs and how the PSFs adjust the HEPs are investigated through literature and discussed in this paper.

99 GENERAL AND MISCELLANEOUS↗

Daily and Seasonal Variations of Shallow Seismic Velocities in Southern California From Joint Analysis of H/V Ratios and Autocorrelations of Seismic Waveforms

Establishing a baseline of ongoing secular velocity variations at the subsurface can improve the accuracy of detecting and interpreting short-term velocity changes, and advance the understanding of observed seismic motions and the behavior of subsurface materials. Toward these goals, we develop and apply a deconvolved autocorrelation (DA) method to estimate regional daily and seasonal changes of seismic velocities in southern California. The DA method combines advantages of traditional autocorrelation and Horizontal-to-Vertical Spectral Ratio, and is used to analyze over 10 years of data recorded by 50 stations. The results indicate widespread daily and seasonal changes of up to 10% and 4%, respectively, in the top tens of meters of the crust. The thickness of the surface layer, distance from the coast, and topographic variations are important factors controlling the amplitudes of the resolved velocity variations. The results suggest that changes of soil moisture and thermoelastic strain are likely dominant factors affecting the daily and seasonal variations, respectively. The developed DA method can improve the accuracy and robustness of estimated changes of subsurface materials at other locations.

58 GEOSCIENCES↗