Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 487 records · Page 27

Creation of a Rapid High-Fidelity Aerodynamics Module for a Multidisciplinary Design Environment

In the traditional aerospace vehicle design process, each successive design phase is accompanied by an increment in the modeling fidelity of the disciplinary analyses being performed. This trend follows a corresponding shrinking of the design space as more and more design decisions are locked in. The correlated increase in knowledge about the design and decrease in design freedom occurs partly because increases in modeling fidelity are usually accompanied by significant increases in the computational expense of performing the analyses. When running high fidelity analyses, it is not usually feasible to explore a large number of variations, and so design space exploration is reserved for conceptual design, and higher fidelity analyses are run only once a specific point design has been selected to carry forward. The designs produced by this traditional process have been recognized as being limited by the uncertainty that is present early on due to the use of lower fidelity analyses. For example, uncertainty in aerodynamics predictions produces uncertainty in trajectory optimization, which can impact overall vehicle sizing. This effect can become more significant when trajectories are being shaped by active constraints. For example, if an optimal trajectory is running up against a normal load factor constraint, inaccuracies in the aerodynamic coefficient predictions can cause a feasible trajectory to be considered infeasible, or vice versa. For this reason, a trade must always be performed between the desired fidelity and the resources available. Apart from this trade between fidelity and computational expense, it is very desirable to use higher fidelity analyses earlier in the design process. A large body of work has been performed to this end, led by efforts in the area of surrogate modeling. In surrogate modeling, an up-front investment is made by running a high fidelity code over a Design of Experiments (DOE); once completed, the DOE data is used to create a surrogate model, which captures the relationships between input variables and responses into regression equations. Depending on the dimensionality of the problem and the fidelity of the code for which a surrogate model is being created, the initial DOE can itself be computationally prohibitive to run. Cokriging, a modeling approach from the field of geostatistics, provides a desirable compromise between computational expense and fidelity. To do this, cokriging leverages a large body of data generated by a low fidelity analysis, combines it with a smaller set of data from a higher fidelity analysis, and creates a kriging surrogate model with prediction fidelity approaching that of the higher fidelity analysis. When integrated into a multidisciplinary environment, a disciplinary analysis module employing cokriging can raise the analysis fidelity without drastically impacting the expense of design iterations. This is demonstrated through the creation of an aerodynamics analysis module in NASA’s OpenMDAO framework. Aerodynamic analyses including Missile DATCOM, APAS, and USM3D are leveraged to create high fidelity aerodynamics decks for parametric vehicle geometries, which are created in NASA’s Vehicle Sketch Pad (VSP). Several trade studies are performed to examine the achieved level of model fidelity, and the overall impact to vehicle design is quantified.

Aerodynamic Properties↗

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced multimaterial shape optimization methods as applied to advanced manufacturing of wind turbine generators

Abstract Currently, many utility‐scale wind turbine generator original equipment manufacturers are dependent on imported rare earth permanent magnets, which are susceptible to market risks from cost instability. To lower the production costs of these generators and stay competitive in the market, several small wind manufacturers are pursuing continuous improvements to both generator design and manufacturing. However, traditional design and manufacturing methods have yielded marginal improvements in wind power performance. This work presents novel methods to redesign a baseline 15‐kW wind turbine generator with reduced rare‐earth permanent magnets by leveraging cutting‐edge three‐dimensional (3D) printed polymer‐bonded permanent magnets and steel. Symmetric, asymmetric, and multimaterial‐magnet parametrization methods are introduced for shape optimization. We extend the symmetric and asymmetric methods to the back iron in the stator to further investigate the impact and opportunities for performance improvements with lesser active materials. We employ a design‐of‐experiments approach with parametric computer‐aided design for shape generation and evaluate different designs by magneto‐thermal modeling and finite‐element analysis. We use adaptive sampling technique to identify better performing designs with lesser magnet mass, higher efficiency, and lower cogging torque when compared with the baseline generator. Asymmetric pole designs resulted in a magnet mass in the range of 4.77–5.37 kg, which was 27%–35% lighter than the baseline generator, suggesting that a new design freedom exists that can be enabled by advanced manufacturing, such as 3D printing. Shaping the back iron in the stator resulted in material savings in electrical steel of up to 14.62 kg, which was 20% lighter than the baseline stator. We conducted a structural analysis to evaluate an optimized asymmetric rotor design from the point of view of mechanical integrity and air‐gap stiffness. The magnetically optimal shape profile was shown as having a positive impact on the radial stiffness, and an optimal solution was discovered to reduce the structural mass by nearly 30 kg, which was 29% lighter than the baseline.

17 WIND ENERGY↗

Physics-Informed Active Learning With Simultaneous Weak-Form Latent Space Dynamics Identification

The parametric greedy latent space dynamics identification (gLaSDI) framework has demonstrated promising potential for accurate and efficient modeling of high-dimensional nonlinear physical systems. However, it remains challenging to handle noisy data. Here, to enhance robustness against noise, we incorporate the weak-form estimation of nonlinear dynamics (WENDy) into gLaSDI. In the proposed weak-form gLaSDI (WgLaSDI) framework, an autoencoder and WENDy are trained simultaneously to discover intrinsic nonlinear latent-space dynamics of high-dimensional data. Compared with the standard sparse identification of nonlinear dynamics (SINDy) employed in gLaSDI, WENDy enables variance reduction and robust latent space discovery, therefore leading to more accurate and efficient reduced-order modeling. Furthermore, the greedy physics-informed active learning in WgLaSDI enables adaptive sampling of optimal training data on the fly for enhanced modeling accuracy. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including viscous and inviscid Burgers' equations, time-dependent radial advection, and the Vlasov equation for plasma physics. With data that contains 5%–10% Gaussian white noise, WgLaSDI outperforms gLaSDI by orders of magnitude, achieving 1%–7% relative errors. Compared with the high-fidelity models, WgLaSDI achieves 121 to 1779x speed-up.

97 MATHEMATICS AND COMPUTING↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

GridTool: A surface modeling and grid generation tool

GridTool is designed around the concept that the surface grids are generated on a set of bi-linear patches. This type of grid generation is quite easy to implement, and it avoids the problems associated with complex CAD surface representations and associated surface parameterizations. However, the resulting surface grids are close to but not on the original CAD surfaces. This problem can be alleviated by projecting the resulting surface grids onto the original CAD surfaces. GridTool is designed primary for unstructured grid generation systems. Currently, GridTool supports VGRID and FELISA systems, and it can be easily extended to support other unstructured grid generation systems. The data in GridTool is stored parametrically so that once the problem is set up, one can modify the surfaces and the entire set of points, curves and patches will be updated automatically. This is very useful in a multidisciplinary design and optimization process. GridTool is written entirely in ANSI 'C', the interface is based on the FORMS library, and the graphics is based on the GL library. The code has been tested successfully on IRIS workstations running IRIX4.0 and above. The memory is allocated dynamically, therefore, memory size will depend on the complexity of geometry/grid. GridTool data structure is based on a link-list structure which allows the required memory to expand and contract dynamically according to the user's data size and action. Data structure contains several types of objects such as points, curves, patches, sources and surfaces. At any given time, there is always an active object which is drawn in magenta, or in their highlighted colors as defined by the resource file which will be discussed later.

Samareh-Abolhassani, Jamshid↗

Reliable photometric membership (RPM) of galaxies in clusters – I. A machine learning method and its performance in the local universe

ABSTRACT We introduce a new method to determine galaxy cluster membership based solely on photometric properties. We adopt a machine learning approach to recover a cluster membership probability from galaxy photometric parameters and finally derive a membership classification. After testing several machine learning techniques (such as stochastic gradient boosting, model averaged neural network and k-nearest neighbours), we found the support vector machine algorithm to perform better when applied to our data. Our training and validation data are from the Sloan Digital Sky Survey main sample. Hence, to be complete to $M_r^* + 3$, we limit our work to 30 clusters with $z$phot-cl ≤ 0.045. Masses (M200) are larger than $\sim 0.6\times 10^{14} \, \mathrm{M}_{\odot }$ (most above $3\times 10^{14} \, \mathrm{M}_{\odot }$). Our results are derived taking in account all galaxies in the line of sight of each cluster, with no photometric redshift cuts or background corrections. Our method is non-parametric, making no assumptions on the number density or luminosity profiles of galaxies in clusters. Our approach delivers extremely accurate results (completeness, C $\sim 92{\rm{ per\ cent}}$ and purity, P $\sim 87{\rm{ per\ cent}}$) within R200, so that we named our code reliable photometric membership. We discuss possible dependencies on magnitude, colour, and cluster mass. Finally, we present some applications of our method, stressing its impact to galaxy evolution and cosmological studies based on future large-scale surveys, such as eROSITA, EUCLID, and LSST.

Lopes, Paulo A. A.↗

Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms

This numerical study compares the wave field generated by the spectral wave action balance code, SNL-SWAN, to the linear-wave boundary-element method (BEM) code, WAMIT. The objective of this study is to assess the performance of SNL-SWAN for modeling wave field effects produced by individual wave energy converters (WECs) and wave farms comprising multiple WECs by comparing results from SNL-SWAN with those produced by the BEM code WAMIT. BEM codes better model the physics of wave-body interactions and thus simulate a more accurate near-field wave field than spectral codes. In SNL-SWAN, the wave field’s energy extraction is modeled parametrically based on the WEC’s power curve. The comparison between SNL-SWAN and WAMIT is made over a range of incident wave conditions, including short-, medium-, and long-wavelength waves with various amounts of directional spreading, and for three WEC archetypes: a point absorber (PA), a pitching flap (PF) terminator, and a hinged raft (HR) attenuator. Individual WECs and wave farms of five WECs in various configuration were studied with qualitative comparisons made of wave height and spectra at specific locations, and quantitative comparisons of the wave fields over circular arcs around the WECs as a function of radial distance. Results from this numerical study demonstrate that in the near-field, the difference between SNL-SWAN and WAMIT is relatively large (between 20% and 50%), but in the far-field from the array the differences are minimal (between 1% and 5%). The resultant wave field generated by the two different numerical approaches is highly dependent on parameters such as: directional wave spreading, wave reflection or scattering, and the WEC’s power curve.

16 TIDAL AND WAVE POWER↗

Research on Orbital Plasma-Electrodynamics (ROPE)

Since the development of probe theory by Langmuir and Blodgett, the problem of current collection by a charged spherically or cylindrically symmetric body has been investigated by a number of authors. This paper overviews the development of a fully three-dimensional particle simulation code which can be used to understand the physics of current collection in three dimensions and can be used to analyze data resulting from the future tethered satellite system (TSS). According to the TSS configurations, two types of particle simulation models were constructed: a simple particle simulation (SIPS) and a super particle simulation (SUPS). The models study the electron transient response and its asymptotic behavior around a three dimensional, highly biased satellite. The potential distribution surrounding the satellite is determined by solving Laplace's equation in the SIPS model and by solving Poisson's equation in the SUPS model. Thus, the potential distribution in space is independent of the density distribution of the particles in the SUPS model but it does depend on the density distribution of the particles in the SUPS model. The evolution of the potential distribution in the SUPS model is described. When the spherical satellite is charged to a highly positive potential and immersed in a plasma with a uniform magnetic field, the formation of an electron torus in the equatorial plane (the plane in perpendicular to the magnetic field) and elongation of the torus along the magnetic field are found in both the SIPS and the SUPS models but the shape of the torus is different. The areas of high potential that exist in the polar regions in the SUPS model exaggerate the elongation of the electron torus along the magnetic field. The current collected by the satellite for different magentic field strengths is investigated in both models. Due to the nonlinear effects present in SUPS, the oscillating phenomenon of the current collection curve during the first 10 plasma periods can be seen (this does not appear in SIPS). From the parametric studies, it appears that the oscillating phenomenon of the current collection curve occurs only when the magnetic field strength is less than 0.2 gauss for the present model.

Wu, S. T.↗

Radiative Forcing by Contrails

A parametric study of the instantaneous radiative impact of contrails is presented using three different radiative transfer models for a series of model atmospheres and cloud parameters. Contrails are treated as geometrically and optically thin plane parallel homogeneous cirrus layers in a static atmospheres The ice water content is varied as a function of ambient temperature. The model atmospheres include tropical, mid-latitude, and subarctic summer and winter atmospheres Optically thin contrails cause a positive net forcing at top of the atmosphere. At the surface the radiative forcing is negative during daytime. The forcing increases with the optical depth and the amount of contrail cover. At the top of the atmosphere a mean contrail cover of 0.1% with average optical depth of 0.2 to 0.5 causes about 0.01 to 0.03 W/m(exp 2)a daily mean instantaneous radiative forcing. Contrails cool the surface during the day and heat the surface during the night, and hence reduce the daily temperature amplitude The net effect depends strongly on the daily variation of contrail cloud cover. The indirect radiative forcing due to particle changes in natural cirrus clouds may be of the same magnitude as the direct one due to additional cover.

Meerkoetter, R.↗

The Implications of Collisions on the Spatial Profile of Electric Potential and the Space-Charge-Limited Current

The space-charge-limited current (SCLC) in a vacuum diode is given by the Child-Langmuir law (CLL), whose electric potential ϕ(x) ∝ (x/D) 4/3 , where x is the spatial coordinate across the gap and D is the gap separation distance. For a collisional diode, SCLC is given by the Mott-Gurney law (MGL) and ϕ(x) ∝ (x/D) 3/2 . Here, we apply a capacitance argument for SCLC and use the transit time from a recent exact solution for collisional SCLC to show that ϕ(x) ∝ (x/D)ξ for a general collisional gap, where 4/3 ≤ ξ ≤ 3/2 . Furthermore, ξ is strictly a function of νT, where ν is the collision frequency and T is the electron transit time. Using this definition of ξ, we estimate the spatial dependence of the electron velocity and use the gap capacitance to derive an analytic equation for collisional SCLC that agrees within ~4.5% of the exact solution that requires solving parametrically through T. This analytic equation for general ξ asymptotically recovers the CLL as ν → 0 and the MGL as ν → ∞. As a result, matching these limits shows that ξ ≈ 1.40 and V ∝ D 2 ν 2 at the transition from a vacuum to a collisional diode for any device condition.

Anodes↗

Gaining Perspective on Unconventional Well Design Choices through Play-level Application of Machine Learning Modeling

The recent development of unconventional oil and gas (O&G) reservoirs has led to an abundant hydrocarbon supply, both domestically and globally. However, there is a continued push to develop new and innovative approaches to improve exploration and extraction efficiencies and overall well productivity moving forward. Substantial improvements in unconventional O&G development are expected through optimized well completion and stimulation strategies aimed at maximizing well productivity. Optimizing well designs will require tailoring to the distinctive geologic conditions present for any newly placed well. To better evaluate the impact of well design attributes and their associated interactions on productivity in a major unconventional play, multivariate machine learning-based models that use empirical datasets were developed. A gradient boosted regression tree (GBRT) algorithm was applied. GBRT has been narrowly investigated for O&G applications but enables straightforward parametric importance and influence evaluation, as well as assessment of parameter interaction effects. Models were trained on well design and locational parameters that serve as a proxy for variable geologic conditions to estimate two types of productivity indicator response variables strongly correlated to estimated ultimate recovery (EUR). The dataset utilized consists of over 7,000 well observations that cover the majority of the productive region of the Marcellus Shale. Model performance was evaluated and algorithm parameters tuned by analyzing the goodness-of-fit for simulated results against observed data in a cross-validation approach. Models were found capable of 73–79 percent prediction accuracy on held out testing data of gas equivalent production and can be used to inform future well design and placement decisions for increasing EUR per well and improving overall field-level recovery. Study results indicate that Marcellus well performance improves most with upscaling perforated interval lengths and water and proppant volumes per foot; but relative productivity improvements are spatially dependent across the play. Finally, optimal combinations of water and proppant on well performance were found to vary depending on well location, emphasizing the utility of data-driven models capable of broad application across a play of interest for informing tailored well design approaches prior to their field deployment.

04 OIL SHALES AND TAR SANDS↗

Constraining the scatter in the galaxy–halo connection at Milky Way masses

In this work, we develop and implement two new methods for constraining the scatter in the relationship between galaxies and dark matter haloes. These new techniques are sensitive to the scatter at low halo masses, making them complementary to previous constraints that are dependent on clustering amplitudes or rich galaxy groups, both of which are only sensitive to more massive haloes. In both of our methods, we use a galaxy group finder to locate central galaxies in the Sloan Digital Sky Survey main galaxy sample. Our first technique uses the small-scale cross-correlation of central galaxies with all lower mass galaxies. This quantity is sensitive to the satellite fraction of low-mass galaxies, which is in turn driven by the scatter between haloes and galaxies. The second technique uses the kurtosis of the distribution of line-of-sight velocities between central galaxies and neighbouring galaxies. This quantity is sensitive to the distribution of halo masses that contain the central galaxies at fixed stellar mass. Theoretical models are constructed using peak halo circular velocity, Vpeak, as our property to connect galaxies to haloes, and all comparisons between theory and observation are made after first passing the model through the group-finding algorithm. We parametrize scatter as a lognormal distribution in M* at fixed Vpeak, σ[M*|Vpeak]. The cross-correlation technique yields a constraint of σ[M*|Vpeak] = 0.27 ± 0.05 dex at a mean Vpeak of 168 km s-1, corresponding to a scatter in log M* at fixed Mh of σ[M*|Mh] = 0.38 ± 0.06 dex at $M_\text{h}=10^{11.8}\, \text{M}_\odot$. The kurtosis technique yields σ[M*|Vpeak] = 0.30 ± 0.03 at Vpeak = 209 km s-1, corresponding to σ[M*|Mh] = 0.34 ± 0.04 at $M_\text{h}=10^{12.2}\, \text{M}_\odot$. The values of σ[M*|Mh] are significantly larger than the constraints at higher masses, in agreement with the results of hydrodynamic simulations. This increase is only partly due to the scatter between Vpeak and Mh, and it represents an increase of nearly a factor of 2 relative to the values inferred from clustering and group studies at high masses.

59 BASIC BIOLOGICAL SCIENCES↗

Improved Parametric Models for Explosion Pressure Signals Derived From Large Datasets

Accurate recording and characterization of explosion-induced pressure signals are key components of the forensic analysis of explosion events in the atmosphere. Parametric overpressure models based on several key waveform features (peak overpressure, positive pulse duration, and impulse) are widely used to estimate explosion energy in terms of trinitrotoluene equivalent yield. However, those models are often developed by a limited dataset, including only a few events or recordings at relatively short propagation distances. Here, we develop empirical waveform-parameter models based on a regression analysis of a large set of data curated from four chemical explosion experiments including 16 detonations. We measured peak overpressure and impulse for positive and negative phases from 1000 pressure signals recorded at local ranges (<20 km) with scaled distance up to 8000 m/kg 1/3 . Additionally, the measured waveform parameters showed large variation with respect to observing distances indicating the effects of atmospheric propagation. In this study, a second-order polynomial model was used in a least-squares regression to account for those propagation effects and to improve data fitting. In addition to model parameters for waveform features, we also determined range-dependent model uncertainties based on data variance. The model uncertainty represents the prediction error of our models and can be critical to evaluating the uncertainty of yield estimate.

58 GEOSCIENCES↗

Influence of Cooling Channel Geometry on the Thermal Response in Silicon Nitride Plates Studied

Engine manufacturers are continually attempting to improve the performance and efficiency of internal combustion engines. Usually they raise the operating temperature or reduce the cooling air requirement for the hot section turbine components. However, the success of these attempts depends on finding materials that are lightweight, are strong, and can withstand high temperatures. Ceramics are among the top candidate materials considered for such harsh applications. They hold low-density, high-temperature strength, and thermal conductivity, and they are undergoing investigation as potential materials for replacing nickel-base alloys and superalloys that are currently used for engine hot-section components. Ceramic structures can withstand higher operating temperatures and a harsh combustion environment. In addition, their low densities relative to metals help reduce component mass. The long-term objectives of the High Temperature Propulsion Components (HOTPC) Project are to develop manufacturing technology, thermal and environmental barrier coatings (TBC/EBC), and the analytical modeling capability to predict thermomechanical stresses in minimally cooled silicon nitride turbine nozzle vanes under simulated engine conditions. Two- and three-dimensional finite element analyses with TBC were conducted at the NASA Glenn Research Center. Nondestructive evaluation was used to determine processing defects. The study included conducting preliminary parametric analytical runs of heat transfer and stress analyses under steady-state conditions to demonstrate the feasibility of using cooled Si3N4 parts for turbine applications. The influence of cooling-channel shapes (such as circular, square, and ascending-order cooling channels) on cooling efficiency and thermal stresses was investigated. Temperature distributions were generated for all cases considered under both cooling and no-cooling conditions, with air being the cooling medium. The table shows the magnitude of the maximum and minimum temperature obtained for the plates under air cooling. Each channel's cross-sectional shape delivered a different temperature; however, the two-dimensional analyses for circular and square or equal-side rectangular holes produced close results. Moreover, the model of the panel with ascending order cooling channels experienced the lowest temperature. A difference of near 260 C was found among the three cooling-hole configurations investigated. The ascending-order cooling channels arrangement showed superior performance by attaining the lowest temperature (1077 C) in comparison to the circular (1379 C) and square (1343 C) channels for the same cooling-hole size. This indicates that the panel with ascending-order cooling channels is the most suitable configuration regardless of the complexity involved in its manufacture. More details pertaining to this study are reported.

Abdul-Aziz, Ali↗

Systematic exploration of heavy element nucleosynthesis in protomagnetar outflows

ABSTRACT We study the nucleosynthesis products in neutrino-driven winds from rapidly rotating, highly magnetized and misaligned protomagnetars using the nuclear reaction network SkyNet. We adopt a semi-analytic parametrized model for the protomagnetar and systematically study the capabilities of its neutrino-driven wind for synthesizing nuclei and eventually producing ultra-high energy cosmic rays (UHECRs). We find that for neutron-rich outflows (Ye < 0.5), synthesis of heavy elements ($\overline{A}\sim 20-65$) is possible during the first $\sim 10\, {\rm s}$ of the outflow, but these nuclei are subjected to composition-altering photodisintegration during the epoch of particle acceleration at the dissipation radii. However, after the first $\sim 10\, {\rm s}$ of the outflow, nucleosynthesis reaches lighter elements ($\overline{A}\sim 10-50$) that are not subjected to subsequent photodisintegration. For proton-rich (Ye ≥ 0.5) outflows, synthesis is more limited ($\overline{A}\sim 4-15$). These suggest that while protomagnetars typically do not synthesize nuclei heavier than second r-process peak elements, they are intriguing sources of intermediate/heavy mass UHECRs. For all configurations, the most rapidly rotating protomagnetars are more conducive for nucleosynthesis with a weaker dependence on the magnetic field strength.

79 ASTRONOMY AND ASTROPHYSICS↗

Energy Saving Quantification on Ductless Heat Pump (DHP) in Existing Homes

In residential retrofit applications, ductless mini-split heat pumps (DHP) are often reported to have high-energy savings potential, depending on the system they are supplementing or replacing. However, recently, there have been a number of utility studies and analyses indicating these energy savings are not being achieved when the existing system is left in place for backup heating or air conditioning. The Pacific Northwest National Laboratory conducted a three-phase project to help determine which control strategies would have the most energy savings impact in various climate zones around the United States. The first phase developed the standalone simulation model for the PNNL Lab Homes and investigated the energy-saving potential from different control strategies and HVAC system configurations. The second phase focused on conducting experiments in the PNNL Lab Homes, which tested the most promising solutions that were modeled in the first phase. The third phase used the field data to calibrate the simulation model and then extrapolated the results to different climate locations and different building sizes. This report focuses on the third phase of the study, including five major parts: model calibration, parametric model setup, results, a sensitivity analysis of air leakage rate, and conclusions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Polar Heat Flow on Io

Recently, Galileo spacecraft data have revealed Io's polar regions to be much warmer than previously expected. This unexpected development came from Photo-Polarimeter Radiometer (PPR) data which show that the minimum night temperatures are in the range of 90-95 K virtually everywhere on Io. The minimum night temperatures show no dependence upon latitude and, when away from the sunset terminator, they show no dependence upon time of night. This is indeed bizarre behavior for surface units which generally had been assumed to be passive with respect to Io's pervasive volcanism. Night temperatures of 90-95 K at high, polar latitudes are particularly hard to explain. Even assuming infinite thermal inertia, at these latitudes there is insufficient sunlight to support these warm night temperatures. Thus, through the process of elimination of other possibilities, we come to the conclusion that these surfaces are volcanically heated. Taking previously passive units and turning them into new sources of heat flow is a radical departure from previous thermophysical model paradigms. However, the geological interpretation is straight forward. We are simply seeing the effect of old, cool lava flows which cover most of the surface of Io but yet have some heat to radiate. Under these new constraints, we have taken on the challenge of formulating a physical model which quantitatively reproduces all of the observations of Io's thermal emission. In the following we introduce a new parametric model which suffices to identify a previously unrecognized polar component of Io's heat flow.

Veeder, G. J.↗