Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced basis method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

Scattering by a groove in an impedance plane

An analysis of two-dimensional scattering from a narrow groove in an impedance plane is presented. The groove is represented by a impedance surface and the problem reduces to that of scattering from an impedance strip in an otherwise uniform impedance plane. On the basis of this model, appropriate integral equations are constructed using a form of the impedance plane Green's functions involving rapidly convergent integrals. The integral equations are solved by introducing a single basis representation of the equivalent current on the narrow impedance insert. Both transverse electric (TE) and transverse magnetic (TM) polarizations are treated. The resulting solution is validated by comparison with results from the standard boundary integral method (BIM) and a high frequency solution. It is found that the presented solution for narrow impedance inserts can be used in conjunction with the high frequency solution for the characterization of impedance inserts of any given width.

Bindiganavale, Sunil↗

Correlation Function Approach for Diffusion in Confined Geometries

This paper describes a formalism for extracting spatially varying transport coefficients from simulations of a molecular fluid in a nano channel. This approach is applied to self-diffusion of a Lennard-Jones fluid confined between two parallel surfaces. A numerical grid is laid over the domain confining the fluid, and fluid properties are projected onto the grid cells. The time correlation functions between properties in different grid cells are calculated and can be used as the basis for a fitting procedure for extracting spatially varying diffusion coefficients from the simulation. Results for the Lennard-Jones system show that transport behavior varies sharply near the liquid-solid boundary and that the changes depend on the details of the liquid-solid interaction. A quantitative difference between the reduced and detailed models is discussed. It is found that the difference could be associated with assumptions about the form of the transport equations at molecular scales in lieu of problems with the method itself. The study suggests that this approach to fitting molecular simulations to continuum equations may guide the development of appropriate coarse-grained equations to model transport phenomena at nanometer scales.

Nanoscale flow, Transport, molecular simulation↗

Parallel, adaptive finite element methods for conservation laws

We construct parallel finite element methods for the solution of hyperbolic conservation laws in one and two dimensions. Spatial discretization is performed by a discontinuous Galerkin finite element method using a basis of piecewise Legendre polynomials. Temporal discretization utilizes a Runge-Kutta method. Dissipative fluxes and projection limiting prevent oscillations near solution discontinuities. A posteriori estimates of spatial errors are obtained by a p-refinement technique using superconvergence at Radau points. The resulting method is of high order and may be parallelized efficiently on MIMD computers. We compare results using different limiting schemes and demonstrate parallel efficiency through computations on an NCUBE/2 hypercube. We also present results using adaptive h- and p-refinement to reduce the computational cost of the method.

Biswas, Rupak↗

Medical Data Architecture Project Status

The Medical Data Architecture (MDA) project supports the Exploration Medical Capability (ExMC) risk to minimize or reduce the risk of adverse health outcomes and decrements in performance due to in-flight medical capabilities on human exploration missions. To mitigate this risk, the ExMC MDA project addresses the technical limitations identified in ExMC Gap Med 07: We do not have the capability to comprehensively process medically-relevant information to support medical operations during exploration missions. This gap identifies that the current International Space Station (ISS) medical data management includes a combination of data collection and distribution methods that are minimally integrated with on-board medical devices and systems. Furthermore, there are variety of data sources and methods of data collection. For an exploration mission, the seamless management of such data will enable an increasingly autonomous crew than the current ISS paradigm. The MDA will develop capabilities that support automated data collection, and the necessary functionality and challenges in executing a self-contained medical system that approaches crew health care delivery without assistance from ground support.To attain this goal, the first year of the MDA project focused on reducing technical risk, developing documentation and instituting iterative development processes that established the basis for the first version of MDA software (or Test Bed 1). Test Bed 1 is based on a nominal operations scenario authored by the ExMC Element Scientist. This narrative was decomposed into a Concept of Operations that formed the basis for Test Bed 1 requirements. These requirements were successfully vetted through the MDA Test Bed 1 System Requirements Review, which permitted the MDA project to begin software code development and component integration. This paper highlights the MDA objectives, development processes, and accomplishments, and identifies the fiscal year 2017 milestones and deliverables in the upcoming year.

Exploration medical system↗

Benchmark study of a new simplified DFN model for shearing of intersecting fractures and faults

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al. 46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.

58 GEOSCIENCES↗

A Two-Step Control Approach for Torque Ripple and Vibration Reduction in Switched Reluctance Motor Drives

The objective of this work is to reduce the torque ripple and vibration in Switched Reluctance Motor (SRM) drives. To achieve such a goal, a two-step approach is proposed. In the first step, the reference current is profiled using an optimization method, aimed to reduce the torque pulsation. In the second, the optimum current profile is applied to the SRM through an adaptive hysteresis band controller to actively cancel the radial vibration. The proposed technique uses an optimization procedure based on the Field Reconstruction Method (FRM) and Non-Derivative Optimization Method to find the optimal current profile that mitigates the torque pulsation, however, the investigation has shown that the proposed excitation leads to a higher radial vibration, which has been mitigated by using an adaptive hysteresis band controller. The theoretical basis for the proposed approach is initially presented, followed by system modeling and simulation. In addition, experimental verification has been conducted using an 8/6 SRM to show the feasibility of the proposed approach. The proposed method alleviates the acoustic noise and torque pulsation in SRM and offers an effective step towards the prevalent use of SRM in a wide range of applications.

42 ENGINEERING↗

Status of noise technology for advanced supersonic cruise aircraft

Developments in acoustic technology applicable to advanced supersonic cruise aircraft, particularly those which relate to jet noise and its suppression are reviewed. The noise reducing potential of high radius ratio, inverted velocity profile coannular jets is demonstrated by model scale results from a wide range of nozzle geometries, including some simulated flight cases. These results were verified statistically at large scale on a variable cycle engine (VCE) testbed. A preliminary assessment of potential VCE noise sources such as fan and core noise is made, based on the testbed data. Recent advances in the understanding of flight effects are reviewed. The status of component noise prediction methods is assessed on the basis of recent test data, and the remaining problem areas are outlined.

Stone, J. R.↗

Medical Data Architecture Project Status

The Medical Data Architecture (MDA) project supports the Exploration Medical Capability (ExMC) risk to minimize or reduce the risk of adverse health outcomes and decrements in performance due to in-flight medical capabilities on human exploration missions. To mitigate this risk, the ExMC MDA project addresses the technical limitations identified in ExMC Gap Med 07: We do not have the capability to comprehensively process medically-relevant information to support medical operations during exploration missions. This gap identifies that the current International Space Station (ISS) medical data management includes a combination of data collection and distribution methods that are minimally integrated with on-board medical devices and systems. Furthermore, there are variety of data sources and methods of data collection. For an exploration mission, the seamless management of such data will enable an increasingly autonomous crew than the current ISS paradigm. The MDA will develop capabilities that support automated data collection, and the necessary functionality and challenges in executing a self-contained medical system that approaches crew health care delivery without assistance from ground support. To attain this goal, the first year of the MDA project focused on reducing technical risk, developing documentation and instituting iterative development processes that established the basis for the first version of MDA software (or Test Bed 1). Test Bed 1 is based on a nominal operations scenario authored by the ExMC Element Scientist. This narrative was decomposed into a Concept of Operations that formed the basis for Test Bed 1 requirements. These requirements were successfully vetted through the MDA Test Bed 1 System Requirements Review, which permitted the MDA project to begin software code development and component integration. This paper highlights the MDA objectives, development processes, and accomplishments, and identifies the fiscal year 2017 milestones and deliverables in the upcoming year.

Data management↗

Butterfly, Vinylidene-Like, Monobridged and Trans Structures of Si 2 H 2 + : Comparison to the Well-Characterized Neutral Si 2 H 2

This paper investigates four of the constitutional isomers of Si 2 H 2 + , namely the butterfly, vinylidene-like, monobridged, and trans structures. These isomer geometries were all studied using the CCSD(T) method with basis sets as large as cc-pV5Z. Higher level methods CCSDT and CCSDT(Q) were used for final energetics. It is found that the butterfly isomer has the lowest energy, followed by vinylidene-like and monobridged structures; the trans structure has the highest energy of the four. All structures were compared with their neutral counterparts. Partial charges are computed and it is found that all isomers share the positive charge between the silicon atoms. NBO analysis shows that the cation Si-Si bond order is reduced by 0.5 relative to the neutral versions of the same isomer in the case of the butterfly, vinylidene-like, and monobridged; and reduced by 0.6 for the trans isomer. Vibrational frequencies, infrared intensities, and dipole moments are predicted to encourage the spectroscopic identification of Si 2 H 2 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduction of weather effects in the calculation of sea ice concentration from microwave radiances

A technique is presented which improves existing methods of calculating sea ice concentrations from microwave radiances by reducing weather-related effects over open ocean areas and in the vicinity of marginal sea ice zones. Winds, atmospheric water vapor, cloud liquid water, and rain increase the microwave emission over these regions and thus result in erroneous values of computed sea ice concentration. The method described is based on the microwave spectral properties of sea ice and ice-free ocean and utilizes ratios of the polarized radiances at the 0.81-cm (37 GHz) and 1.7-cm (18 GHz) wavelengths. Following a discussion of the physical basis for this technique, examples are provided which demonstrate its utility. While the technique was developed for use with the Nimbus 7 scanning multichannel microwave radiometer data, it is applicable also to data from other microwave radiometers operating in a similar wavelength range.

Gloersen, P.↗

Reanalysis procedure for large structural systems

Global-basis-vector approximate reanalysis techniques for use in automated structural optimization schemes are developed and demonstrated. The vibrational response of a modified structure (MS) is estimated by lumping the design variables into a single tracing parameter, applying an operator splitting procedure to express the FEM equations of the MS in terms of the original-structure equations plus correction terms, and reducing the MS equations via a classical Bubnov-Galerkin scheme. The sensitivity of the vibrational response to structural modifications is evaluated, and the relationship between this method and the preconditioned conjugate-gradient method (Noor and Peters, 1988) is explored. Numerical results for linear static and free vibration problems involving beamlike lattices, double-layered hexahedral grids, and structural gridworks are presented in tables and graphs and briefly characterized.

Noor, Ahmed K.↗

K-Means Cluster Study for Radiofrequency Propagation Characterization

The objective of this study is to design a simple method for mining radio frequency (RF) propagation data. The study explored the characteristics of a large dataset of propagation experiments conducted over the span of years and using several ground stations around the world. Furthermore, this study developed simple predictive models that can be used for link characterization and overall propagation behavior description, without the need for physical measurements on-site. It is understood that such statistical learning has several drawbacks in terms of accuracy and precision. K-means clustering was used to characterize the data set in a way never explored before in an attempt to create useful tools that reduce cost, time and risk. K-means clustering was used to characterize the data set. Cosine distance was used as a method to determine the optimal number for clustering each feature. Dependence and independence analysis was performed to explore intra and inter-sensitivity between the presented features, with respect to each other and time. Several predicative models were generated and evaluated with respect to a test set to assess a measure of prediction accuracy and precision. A simple method for data analysis was developed and tested as the basis for further studies and future refinement to produce optimal performing models.

Cognitive↗

Reactor and fuel cycle performance of light water reactor fuel with 235 U enrichments above 5%

Recent advances in nuclear fuel materials research, particularly on the topic of accident-tolerant fuels, have brought up potential opportunities for expanding the operating envelope of existing light water reactors. As many of the performance improvements offered by these technologies may be most fully realized by increasing fuel enrichment beyond the standard 5% limit, this paper examines the potential reactor performance and fuel cycle performance of low-enriched uranium oxide fueled light water reactors by generically considering pressurized water reactors with 235 U enrichment from 5 to 7%. Advanced cladding, including accident-tolerant cladding, has the potential to increase fuel burnup limits related to hydrogen in the cladding that coincide with those limits associated with end-of-life reactivity. Therefore, higher enrichment will be necessary in order to realize the higher fuel burnups. This work includes evaluation of the fuel cycle length, discharge burnup, reactivity coefficients, and fuel cycle performance, including radioactive waste and environmental impact metrics per unit energy generated. The analysis was performed using the evaluation metrics from the US Department of Energy Office of Nuclear Energy Fuel Cycle Evaluation and Screening Study. The reactor performance and safety analysis show that enrichments between 5 and 7% would have similar fuel temperature and moderator temperature coefficients. However, the soluble boron coefficient would decrease in magnitude, requiring more corrosive boric acid in the coolant or other methods of reactivity control during the fuel cycle. At these higher enrichments the maximum burnup at the rim of the fuel pellet would increase by almost a factor of two, which is expected to impact the formation of high-burnup structure in the fuel and the corresponding thermo-mechanical fuel properties. The fuel cycle performance assessment shows that increasing enrichment reduces the quantity of high-level waste disposed per unit energy generated, but it increases the natural resource requirements normalized to a gigawatt-electricity-per-year basis. Another impact is the slightly higher discharge burnup, resulting in somewhat different activity levels of the spent nuclear fuel and high-level waste radioactivity at 100 and 100,000 years after fuel discharge. The environmental impacts—including land use, water use, carbon emission, and radiological exposure—are of the same magnitude per unit energy generated. Finally, however, the impacts are distributed differently. Less than 5% enrichment has marginally more impact on the back-end of the fuel cycle, and greater than 5% enrichment has marginally more impact on the front-end of the fuel cycle. Ultimately, no neutronic or reactor safety hindrances to employing light water reactor fuel with enrichments greater than 5% are identified; given the achievable reactor performance benefits with advanced fuels, further practical exploration of increased enrichment fuel is recommended.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗

Electronic structure theory with molecular point group symmetries on quantum annealers

Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan–Wigner encoding based on the full Boolean symmetry group $\mathbb{Z}$$^{k}_{2}$ with our new implementation of the Xia–Bian–Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H 2 , LiH, He 2 , H 2 O, O 2 , N 2 , Li 2 , F 2 , CO, BH 3 , NH 3 , and CH 4 , with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method’s Hamiltonian embedding, or for certain molecules such as He 2 and H 2 O. Here, we provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Approach to the origin of turbulence on the basis of two-point kinetic theory

Equations for the fluctuation correlation in an incompressible shear flow are derived on the basis of kinetic theory, utilizing the two-point distribution function which obeys the BBGKY hierarchy equation truncated with the hypothesis of 'ternary' molecular chaos. The step from the molecular to the hydrodynamic description is accomplished by a moment expansion which is a two-point version of the thirteen-moment method, and which leads to a series of correlation equations, viz., the two-point counterparts of the continuity equation, the Navier-Stokes equation, etc. For almost parallel shearing flows the two-point equation is separable and reduces to two Orr-Sommerfeld equations with different physical implications.

Tsuge, S.↗