Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 2”

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 901 records · Page 50

Neutrino emission from binary neutron star mergers: characterising light curves and mean energies

Abstract Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resulting r -process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several $$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$ 10 53 erg s - 1 ) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of Proton Activity Gaps between Electrodes on Open-Circuit Potential of H 2 /H 2 and H 2 /Air Cells

Polymer -electrolyte-fuel-cell open-circuit voltages (OCVs) are exactly defined by equation (1), where cathode and anode proton activities [(aH+)cathode and (aH+)anode, respectively] usually are identical, so the third term in the right-hand side of equation (2) is ignored. OCV=E0+RT/2F*ln(a1/2 O2*(a2 H+)cathode/aH2O)-RT/2F*ln((a2 H+)anode/aH2) (1) =E0+RT/2F*ln(a1/2 O2*aH2/aH2O)+RT/2F*ln((a2 H+)cathode/(a2 H+)anode) (2) Water vapor pressure is a colligative property that fundamentally correlates to electrolyte concentrations in aqueous solutions. Proton activity is a function of acid concentration, such as pH, when electrolytes are acids. In polymer-electrolyte membranes, water vapor pressure and acid concentration are understood as relative humidity (RH) and water uptake (λ), respectively, where λ represents number of water molecules per sulfonic acid molecule. Several investigations have reported the relationship between RH and λ, meaning that proton activities and associated water uptakes are intimately related to RH. In actual fuel cell operation, cathode RH is determined by ambient-atmosphere and/or humidifier RH(s), and anode RH depends on hydrogen-circulator RH. Therefore, RH is not always identical at both electrodes, and the difference between electrode RHs is considerable during dry operation of polymer electrolyte fuel cells. Therefore, the third term in the right-hand side of equation (2) may be significant for dry operation. We measured OCVs when hydrogen was supplied to both electrodes at 80°C. One electrode (A) was fixed at 30% RH, while RH at the other electrode (B) was varied (0, 5, 10, 20, and 30%). Measured OCVs varied from 0 to 75 mV. For fuel cell tests, electrode A was supplied with hydrogen at 30% RH; electrode B, oxygen at 0, 5, 10, 20, and 30% RH. OCVs deviated from that measured when RH at electrode B was 30%, increasing from 0 to 60 mV with decreasing RH at electrode B. Results are also shown in Figure 1. Proton activities of both electrodes were thermodynamically calculated. The Gibbs–Duhem relation was applied to obtain molar Gibbs free energies of water and sulfonic acid, and proton activity coefficient was calculated using the Gibbs free energy of sulfonic acid and the relationship between RH and λ1–4, assuming that protons and sulfonic anions show identical ionic-activity coefficients. OCVs were estimated using the third term in the right-hand side of equation (2). Results are shown in Figure 1. Fuel-cell current–voltage performance was poor when RHs at the anode and cathode were 30 and 20%, respectively. To determine kinetic current, we measured the oxygen-reduction reaction (ORR) using a rotating-disk electrode (RDE) in concentrated-acid aqueous solutions, which modeled catalyst-layer ionomers. Kinetic currents decreased with acid concentrations. References T. A. Zawodzinski, Jr., C. Derouin, S. Radzinski, R. J. Sherman, V. T. Smith, T. E. Springer and S. Gottesfeld , J. Electrochem. Soc., 140,1041 (1993) P. K. Das and A. Z. Weber, Proceedings of the ASME 2013 11th Fuel Cell Science, Engineering and Technology Conference, Fuel Cell 18010 (2013) V. A. Sethuraman, J. W. Weidner, A. T. Haug, S. Motupally,b and L. V. Protsailo, J. Electrochem. Soc., 155, B50 (2008) A. Kusoglu and A. Z. Weber, Chem. Rev., 117, 987 (2017) Figure 1

Yoshida, Toshihiko↗

Circulation of the atmosphere from the surface to 100 km

A reference model of the atmospheric circulation on Venus based on available observations is presented. The reference atmosphere has the following main features: (1) the entire atmosphere below 85 km moves predominantly from east to west in the planet's reference system (in the same direction as the rotation of the solid planet itself) with the possible exception of the lowest 10 km where velocities are low, (2) a jet is present near the cloud-top level at 45 deg latitude in both hemispheres, with a magnitude of approximately 100 m/s, (3) a weak meridional (north-south component) flow directed towards either pole is superimposed on the zonal (east-west motion at cloud-top level (about 68 km) altitudes, and, (4) eddies or wave motions are present in the atmosphere, with amplitudes of less than 16/ms in the upper atmosphere.

Kerzhanovich, V. V.↗

Experimental investigation of transonic oscillating cascade aerodynamics

Fundamental experiments are performed in the NASA Lewis Transonic Oscillating Cascade Facility to investigate the subsonic and transonic aerodynamics of cascaded airfoils executing torsion mode oscillations at realistic values of reduced frequency. In particular, an unsteady aerodynamic influence coefficient technique is developed and utilized. In this technique, only one airfoil in the cascade is oscillated at a time, with the resulting airfoil surface unsteady pressure distribution measured on one dynamically instrumented reference airfoil. The unsteady aerodynamics of an equivalent cascade with all airfoils oscillating at any specified interblade phase angle are then determined through a vector summation of these data. These influence coefficient determined oscillating cascade data were correlated with: (1) data obtained in this cascade with all airfoils oscillating at several interblade phase angle values; and (2) predictions from a classical linearized unsteady cascade model.

Buffum, Daniel H.↗

Novel Gamma Attenuation Methods Utilized for Waste Characterization of UO{sub 2}F{sub 2}-Bearing Materials at the Paducah GDP - 17327

The Paducah Gaseous Diffusion Plant (PGDP) is currently under deactivation with demolition in the decades to follow. UO{sub 2}F{sub 2}-bearing waste is generated at PGDP by UF{sub 6} hydrolysing to materials such as pipes, valves, PPE, etc. when UF{sub 6} is exposed to humid air. In addition, Al{sub 2}O{sub 3} and NaF trap materials contain significant amounts of Uranium. The UO{sub 2}F{sub 2} deposits within the generated waste can be of sub-gram levels to several hundreds of localized grams and are much denser than the surrounding interfering matrix materials and may be of very complex geometric and chemical forms such as with the Al{sub 2}O{sub 3} and NaF trap materials. The main objective is to characterize the fissile (U-235) content of the UO{sub 2}F{sub 2}-bearing waste so as to meet the Data Quality Objectives (DQO) for proper disposal. If a waste container cannot meet the DQO then it must be subjected to costly remediation, repackaging and re-characterization processes. We have developed, at PGDP, several novel model-based methods and techniques that vastly improve the standard gamma attenuation correction methods and do not require Working Reference Materials (WRM) for calibration. These new methods involve rotating and non-rotating 18.9 liter (5 U.S. gallon) waste containers with and without a transmission modality. A new model-based variant of differential attenuation has been developed that provides a high-quality self-absorption correction that does not require knowledge of the Uranium enrichment or isotopic ratios. For more challenging waste containers, a method using multi-view (non-rotating) 18.9 liter waste containers has been considered using a recently developed Advanced ISOCS{sup C} Uncertainty Estimator (AIUE) package. The new methods suffice to significantly improve accuracy and precision and to lower the Total Measurement Uncertainty (TMU) leading to satisfactory DQO's that were not previously achievable using the standard gamma attenuation correction methods. In this paper we present real 18.9 liter waste container and Working Reference Material (WRM) data subjected to the new methods and compared to standard methods depicting the significant improvement in the DQO that is realized. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A three-dimensional sectional representation of aerosol mixing state for simulating optical properties and cloud condensation nuclei

Light absorption by black carbon (BC) particles emitted from fossil fuel combustion depends on the how thickly they are coated with non-refractory species such as ammonium, sulfate, nitrate, organics, and water. The cloud condensation nuclei (CCN) activation property of a particle depends on its dry size and the hygroscopicities of all the individual species mixed together. It is therefore necessary to represent both size and mixing state of aerosols to reliably predict their climate-relevant properties in atmospheric models. Here we describe and evaluate a novel sectional framework in the Model for Simulating Aerosol Interactions and Chemistry, referred to as MOSAIC-mix, that represents the mixing state by resolving aerosol dry size (D dry ), BC dry mass fraction (w BC ), and hygroscopicity (κ). Using ten idealized urban plume scenarios in which different types of aerosols evolve over 24 hours under a range of atmospherically relevant environmental conditions, we examine errors in CCN concentrations and optical properties with respect to a more explicit aerosol mixing state representation. We find that only a small number of w BC and κ bins are needed to achieve significant reductions in the errors, and propose a configuration consisting of 24 D dry bins, 2 w BC bins, and 2 κ bins that gives 24-hour average errors of about 5% or less in CCN concentrations and optical properties, 3-4 times lower than those from size-only-resolved simulations. Here, these results show that MOSAIC-mix is suitable for use in regional and global models to examine the effects of evolving aerosol mixing states on aerosol-radiation-cloud feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Measurement of prompt D 0 , D + , D *+ , and DS+$ {\mathrm{D}}_{\mathrm{S}}^{+} $ production in p–Pb collisions at sNN$ \sqrt{{\mathrm{s}}_{\mathrm{NN}}} $ = 5.02 TeV

The measurement of the production of prompt D 0 , D + , D *+, and $D_S^+$ mesons in proton–lead (p–Pb) collisions at the centre-of-mass energy per nucleon pair of \mathrm{\sqrt{s_{NN}}}$ = 5.02 TeV, with an integrated luminosity of 292 ± 11 μb -1 , are reported. Differential production cross sections are measured at mid-rapidity (-96 < y cms < 0.04) as a function of transverse momentum (p T ) in the intervals 0 < p T < 36 GeV/c for D 0 , 1 < p T < 36 GeV/c for D + and D *+ , and 2 < p T < 24 GeV/c for D + mesons. For each species, the nuclear modification factor R pPb is calculated as a function of p T using a proton-proton (pp) reference measured at the same collision energy. The results are compatible with unity in the whole p T range. The average of the non-strange D mesons R pPb is compared with theoretical model predictions that include initial-state effects and parton transport model predictions. The p T dependence of the D 0 , D + , and D *+ nuclear modification factors is also reported in the interval 1 < p T < 36 GeV/c as a function of the collision centrality, and the central-to-peripheral ratios are computed from the D-meson yields measured in different centrality classes. The results are further compared with charged-particle measurements and a similar trend is observed in all the centrality classes. The ratios of the p T -differential cross sections of D 0 , D + , D *+ , and $D_S^+$ mesons are also reported. The $D_S^+$ and D + yields are compared as a function of the charged-particle multiplicity for several p T intervals. No modification in the relative abundances of the four species is observed with respect to pp collisions within the statistical and systematic uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Context. Analysis of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE) data suggests the existence of a clear distinction between two sequences of disc stars in the [α/Fe] versus [Fe/H] abundance ratio space, known as the high- and low-α sequence, respectively. This dichotomy also emerges from an analysis of the vertical distribution of the [α/Fe] abundance ratio. Aims. We aim to test whether the revised two-infall chemical evolution models designed to reproduce the low- and high-α sequences in the [α/Fe] versus [Fe/H] ratios in the solar neighbourhood are also capable of predicting the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Methods. Along with the chemical composition of the simple stellar populations born at different Galactic times predicted by our reference chemical evolution models in the solar vicinity, we provide their maximum vertical height above the Galactic plane |zmax| computed assuming the relation between the vertical action and stellar age in APOGEE thin-disc stars. Result. The vertical distribution of the [Mg/Fe] abundance ratio predicted by the reference chemical evolution models is in agreement with that observed when combining the APOGEE DR16 data (chemical abundances) with the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-α sequence stars). Including the high-α disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z|< 2 kpc. However, our model predicts an overly flat (almost constant) growth of the maximum vertical height |zmax| quantity as a function of [Mg/Fe] for high-α objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are that: (i) the APOGEE sample with |z|< 2 kpc is more likely than ours to be contaminated by halo stars, causing the median values to be kinematically hotter, and (ii) external perturbations – such as minor mergers – that the Milky Way experienced in the past could have heated up the disc, and the heating of the orbits cannot be modeled by only scattering processes. Assuming a disc dissection based on chemistry for APOGEE-DR16 stars (|z|< 2 kpc), the observed |zmax| distributions for high-α and low-α sequences are in good agreement with our model predictions if we consider the errors in the vertical action estimates in the calculation. Moreover, a better agreement between predicted and observed stellar distributions at different Galactic vertical heights is achieved if asteroseismic ages are included as a constraint in the best-fit model calculations. Conclusions. The signature of a delayed gas infall episode, which gives rise to a hiatus in the star formation history of the Galaxy, are imprinted both in the [Mg/Fe] versus [Fe/H] relation and in vertical distribution of [Mg/Fe] abundances in the solar vicinity.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scale translation yields insights into gas adsorption under nanoconfinement

This work describes a scale-translating simulation framework to investigate gas adsorption behavior in nanoconfined pores. The framework combines molecular simulations (MSs), equation of state (EoS), and lattice Boltzmann (LB) simulations. MSs reveal the physics of methane adsorption in nano-sized pores, where input values of fugacity coefficients are optimized based on EoS predictions. Then, an LB free-energy model, which incorporates a viral EoS, upscales intermolecular forces and estimates adsorption behavior via a proposed fluid–wall interaction model. Armed with the values of the LB interaction parameter as a function of pressure, the LB model is used to predict fluid behavior in irregular nanopores, and the results are validated against reference MS data. The LB model is then used to study adsorption behavior at a continuum scale in representative organic shale nanopores based on finely characterized Vaca Muerta shale samples. Furthermore, the results show that methane adsorption could significantly increase contained fluids by 10%–25% in pores smaller than 20 nm. However, in larger pores (40 nm to 90 nm), adsorption's impact diminishes to 2%–3%, suggesting sorption's negligible role beyond a 40 nm pore size.

74 ATOMIC AND MOLECULAR PHYSICS↗

In Situ Liquid Electron Microscope Cells Strongly Attenuate Electrochemical Behavior

In situ electrochemical Scanning Transmission Electron Microscopy (ec-S/TEM) probes the dynamics of electrochemical processes in real time and at high spatial resolution. This method utilizes microfabricated electrochemical cells, with enclosed dimensions of 800 μm wide, 800 μm long, and 1–2 μm tall, with electron-transparent silicon nitride windows. Working, counter, and reference electrodes in the configuration studied are deposited Pt with a planar surface area of about 1 μm 2 and a height of 50 nm. The cell confines the electrolyte within a 1 μm tall channel, and the microfabricated electrodes have a non-standard geometry, are coplanar, and are not front-facing. As such, standard assumptions of electrochemical experiments do not apply. COMSOL modeling was used to compare and contrast differences in electrochemical behavior between conventional experimental setups and the in situ ec-S/TEM cells. Cell height strongly affected voltammogram measurements: peak magnitudes and shifts correlate directly to cell height. Reactant concentration above the electrode quickly decreases during deposition, leading to a planar diffusion-dominated regime. The concentration gradient prioritizes particle growth along the thin edges of the electrode and not on the planar part of the electrode. Experiments done in the ec-S/TEM ex situ of a TEM verify the model’s accuracy.

Wittman, Reed M. [Oak Ridge National Laboratory (O↗

Advanced Hydrogen Turbine Development

Under the sponsorship of the U.S. Department of Energy (DOE) National Energy Technology Laboratories, Siemens has completed the Advanced Hydrogen Turbine Development Program to develop an advanced gas turbine for incorporation into future coal-based Integrated Gasification Combined Cycle (IGCC) plants. All the scheduled DOE Milestones were completed and significant technical progress was made in the development of new technologies and concepts. Advanced computer simulations and modeling, as well as subscale, full scale laboratory, rig and engine testing were utilized to evaluate and select concepts for further development. Program Requirements of: A 3 to 5 percentage point improvement in overall plant combined cycle efficiency when compared to the reference baseline plant; 20 to 30 percent reduction in overall plant capital cost when compared to the reference baseline plant; and NOx emissions of 2 PPM out of the stack. were all met. The program was completed on schedule and within the allotted budget

42 ENGINEERING↗

Robustness of adaptive control algorithms in the presence of unmodeled dynamics

This paper reports the outcome of an exhaustive analytical and numerical investigation of stability and robustness properties of a wide class of adaptive control algorithms in the presence of unmodeled dynamics and output disturbances. The class of adaptive algorithms considered are those commonly referred to as model-reference adaptive control algorithms, self-tuning controllers, and dead-beat adaptive controllers; they have been developed for both continuous-time systems and discrete-time systems. The existing adaptive control algorithms have been proven to be globally asymptotically stable under certain assumptions, the key ones being (1) that the number of poles and zeroes of the unknown plant are known, and (2) that the primary performance criterion is related to good command following. These theoretical assumptions are too restrictive from an engineering point of view. Real plants always contain unmodeled high-frequency dynamics and small delays, and hence no upper bound on the number of the plant poles and zeroes exists. Also real plants are always subject to unmeasurable output additive disturbances, although these may be guide small. Hence, it is important to critically examine the stability robustness properties of the existing adaptive algorithms when some of the theoretical assumptions are removed; in particular, their stability and performance properties in the presence of unmodeled dynamics and output disturbances. Previously announced in STAR as N83-16061

Rohrs, C. E.↗

Computation of turbulent boundary layers employing the defect wall-function method

In order to decrease overall computational time requirements of spatially-marching parabolized Navier-Stokes finite-difference computer code when applied to turbulent fluid flow, a wall-function methodology, originally proposed by R. Barnwell, was implemented. This numerical effort increases computational speed and calculates reasonably accurate wall shear stress spatial distributions and boundary-layer profiles. Since the wall shear stress is analytically determined from the wall-function model, the computational grid near the wall is not required to spatially resolve the laminar-viscous sublayer. Consequently, a substantially increased computational integration step size is achieved resulting in a considerable decrease in net computational time. This wall-function technique is demonstrated for adiabatic flat plate test cases from Mach 2 to Mach 8. These test cases are analytically verified employing: (1) Eckert reference method solutions, (2) experimental turbulent boundary-layer data of Mabey, and (3) finite-difference computational code solutions with fully resolved laminar-viscous sublayers. Additionally, results have been obtained for two pressure-gradient cases: (1) an adiabatic expansion corner and (2) an adiabatic compression corner.

Brown, Douglas L.↗

An Integrated Tool for the Coupled Thermal and Mechanical Analysis of Pyrolyzing Heatshield Materials

Materials that pyrolyze at elevated temperature have been commonly used as thermal protection materials in hypersonic flight, and advanced pyrolyzing materials for this purpose continue to be developed. Because of the large temperature gradients that can arise in thermal protection materials, significant thermal stresses can develop. Advanced applications of pyrolytic materials are calling for more complex heatshield configurations, making accurate thermal stress analysis more important, and more challenging. For non-pyrolyzing materials, many finite element codes are available and capable of performing coupled thermal-mechanical analyses. These codes do not, however, have a built-in capability to perform analyses that include pyrolysis effects. When a pyrolyzing material is heated, one or more components of the original virgin material pyrolyze and create a gas. This gas flows away from the pyrolysis zone to the surface, resulting in a reduction in surface heating. A porous residue, referred to as char, remains in place of the virgin material. While the processes involved can be complex, it has been found that a simple physical model in which virgin material reacts to form char and pyrolysis gas, will yield satisfactory analytical results. Specifically, the effects that must be modeled include: (1) Variation of thermal properties (density, specific heat, thermal conductivity) as the material composition changes; (2) Energy released or absorbed by the pyrolysis reactions; (3) Energy convected by the flow of pyrolysis gas from the interior to the surface; (4) The reduction in surface heating due to surface blowing; and (5) Chemical and mass diffusion effects at the surface between the pyrolysis gas and edge gas Computational tools for the one-dimensional thermal analysis these materials exist and have proven to be reliable design tools. The objective of the present work is to extend the analysis capabilities of pyrolyzing materials to axisymmetric configurations, and to couple thermal and mechanical analyses so that thermal stresses may be efficiently and accurately calculated.

Pronchick, Stephen W.↗

Coupled Attitude and Orbit Dynamics and Control in Formation Flying Systems

Formation flying systems can range from global constellations offering extended service coverage to clusters of highly coordinated vehicles that perform distributed sensing. Recently, the use of groups of micro-satellites in the areas of near Earth explorations, deep space explorations, and military applications has received considerable attention by researchers and practitioners. To date, most proposed control strategies are based on linear models (e.g., Hill-Clohessy-Wiltshire equations) or nonlinear models that are restricted to circular reference orbits. Also, all models in the literature are uncoupled between relative position and relative attitude. In this paper, a generalized dynamic model is proposed. The reference orbit is not restricted to the circular case. In this formulation, the leader or follower satellite can be in either a circular or an elliptic orbit. In addition to maintaining a specified relative position, the satellites are also required to maintain specified relative attitudes. Thus the model presented couples vehicle attitude and orbit requirements. Orbit perturbations are also included. In particular, the J(sub 2) effects are accounted in the model. Finally, a sliding mode controller is developed and used to control the relative attitude of the formation and the simulation results are presented.

Yunjun Xu↗

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗

Satellite Power Systems (SPS) concept definition study. Volume 2, part 2: System engineering

The latest technical and programmatic developments are considered as well as expansions of the Rockwell SPS cost model covering each phase of the program through the year 2030. Comparative cost/economic analyses cover elements of the satellite, construction system, space transportation vehicles and operations, and the ground receiving station. System plans to define time phased costs and planning requirements that support major milestones through the year 2000. A special analysis is included on natural resources required to build the SPS reference configuration. An appendix contains the SPS Work Breakdown Structure and dictionary along with detail cost data sheet on each system and main element of the program. Over 200 line items address DDT&E, theoretical first unit, investment cost per satellite, and operations charges for replacement capital and normal operations and maintenance costs.

Hanley, G. M.↗