Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Response surface model”

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

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

At least 1,171 records · Page 65

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗

Quantifying CanESM5 and EAMv1 sensitivities to Mt. Pinatubo volcanic forcing for the CMIP6 historical experiment

Large volcanic eruptions reaching the stratosphere have caused marked perturbations to the global climate including cooling at the Earth's surface, changes in large-scale circulation and precipitation patterns and marked temporary reductions in global ocean heat content. Many studies have investigated these effects using climate models; however, uncertainties remain in the modelled response to these eruptions. This is due in part to the diversity of forcing datasets that are used to prescribe the distribution of stratospheric aerosols resulting from these volcanic eruptions, as well as uncertainties in optical property derivations from these datasets. To improve this situation for the sixth phase of the Coupled Model Intercomparison Project (CMIP6), a two-step process was undertaken. First, a combined stratospheric aerosol dataset, the Global Space-based Stratospheric Aerosol Climatology (GloSSAC; 1979–2016), was constructed. Next, GloSSAC, along with information from ice cores and Sun photometers, was used to generate aerosol distributions, characteristics and optical properties to construct a more consistent stratospheric aerosol forcing dataset for models participating in CMIP6. This “version 3” of the stratospheric aerosol forcing has been endorsed for use in all contributing CMIP6 simulations. Recent updates to the underlying GloSSAC from version 1 to version 1.1 affected the 1991–1994 period and necessitated an update to the stratospheric aerosol forcing from version 3 to version 4. As version 3 remains the official CMIP6 input, quantification of the impact on radiative forcing and climate is both relevant and timely for interpreting results from experiments such as the CMIP6 historical simulations. This study uses two models, the Canadian Earth System Model version 5 (CanESM5) and the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1), to estimate the difference in instantaneous radiative forcing in simulated post-Pinatubo climate response when using version 4 instead of version 3. Differences in temperature, precipitation and radiative forcings are generally found to be small compared to internal variability. An exception to this is differences in monthly temperature anomalies near 24 km altitude in the tropics, which can be as large as 3 °C following the eruption of Mt. Pinatubo.

58 GEOSCIENCES↗

Results of 1 MeV proton irradiation of front and back surfaces of silicon solar cells

Several silicon solar cells with and without back surface fields (BSF), having thicknesses of 200 microns and 63 microns were irradiated with 1 MeV protons having fluences between 1 times 10 to the 10th power and 1 times 10 to the 12th power p/square cm. The irradiation was performed using both normal and isotropic incidence on the front as well as back surfaces of the solar cells. The results of the back surface irradiations are analyzed using a model in which irradiation induced defects across the high-low (BSF) junction are considered. It is concluded that degradation of the high-low junction is responsible for the severe performance loss in thinner cells when irradiated from the rear.

Anspaugh, B. E.↗

More accurate predictions with transonic Navier-Stokes methods through improved turbulence modeling

Significant improvements in predictive accuracies for off-design conditions are achievable through better turbulence modeling; and, without necessarily adding any significant complication to the numerics. One well established fact about turbulence is it is slow to respond to changes in the mean strain field. With the 'equilibrium' algebraic turbulence models no attempt is made to model this characteristic and as a consequence these turbulence models exaggerate the turbulent boundary layer's ability to produce turbulent Reynolds shear stresses in regions of adverse pressure gradient. As a consequence, too little momentum loss within the boundary layer is predicted in the region of the shock wave and along the aft part of the airfoil where the surface pressure undergoes further increases. Recently, a 'nonequilibrium' algebraic turbulence model was formulated which attempts to capture this important characteristic of turbulence. This 'nonequilibrium' algebraic model employs an ordinary differential equation to model the slow response of the turbulence to changes in local flow conditions. In its original form, there was some question as to whether this 'nonequilibrium' model performed as well as the 'equilibrium' models for weak interaction cases. However, this turbulence model has since been further improved wherein it now appears that this turbulence model performs at least as well as the 'equilibrium' models for weak interaction cases and for strong interaction cases represents a very significant improvement. The performance of this turbulence model relative to popular 'equilibrium' models is illustrated for three airfoil test cases of the 1987 AIAA Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 'nonequilibrium' turbulence model is currently being applied to wing flows for which similar improvements in predictive accuracy are being realized.

Johnson, Dennis A.↗

Origin of orbital debris impacts on LDEF's trailing surfaces

A model was developed to determine the origin of orbital impacts measured on the training surfaces of LDEF. The model calculates the expected debris impact crater distribution around LDEF as a function of debris orbital parameters. The results show that only highly elliptical, low inclination orbits could be responsible for these impacts. The most common objects left in this type of orbit are orbital transfer stages used by the U.S. and ESA to place payloads into geosynchronous orbit. Objects in this type of orbit are difficult to catalog by the U.S. Space Command; consequently there are independent reasons to believe that the catalog does not adequately represent this population. This analysis concludes that the relative number of cataloged objects with highly elliptical, low inclination orbits must be increased by a factor of 20 to be consistent with the LDEF data.

Kessler, Donald J.↗

Infrared radiation parameterizations for the minor CO2 bands and for several CFC bands in the window region

Fast and accurate parameterizations have been developed for the transmission functions of the CO2 9.4- and 10.4-micron bands, as well as the CFC-11, CFC-12, and CFC-22 bands located in the 8-12-micron region. The parameterizations are based on line-by-line calculations of transmission functions for the CO2 bands and on high spectral resolution laboratory measurements of the absorption coefficients for the CFC bands. Also developed are the parameterizations for the H2O transmission functions for the corresponding spectral bands. Compared to the high-resolution calculations, fluxes at the tropopause computed with the parameterizations are accurate to within 10 percent when overlapping of gas absorptions within a band is taken into account. For individual gas absorption, the accuracy is of order 0-2 percent. The climatic effects of these trace gases have been studied using a zonally averaged multilayer energy balance model, which includes seasonal cycles and a simplified deep ocean. With the trace gas abundances taken to follow the Intergovernmental Panel on Climate Change Low Emissions 'B' scenario, the transient response of the surface temperature is simulated for the period 1900-2060.

Kratz, David P.↗

The rate of chemical weathering of pyrite on the surface of Venus

This abstract reports results of an experimental study of the chemical weathering of pyrite (FeS2) under Venus-like conditions. This work, which extends the earlier study by Fegley and Treiman, is part of a long range research program to experimentally measure the rates of thermochemical gas-solid reactions important in the atmospheric-lithospheric sulfur cycle on Venus. The objectives of this research are (1) to measure the kinetics of thermochemical gas-solid reactions responsible for both the production (e.g., anhydrite formation) and destruction (e.g., pyrrhotite oxidation) of sulfur-bearing minerals on the surface of Venus and (2) to incorporate these and other constraints into holistic models of the chemical interactions between the atmosphere and surface of Venus. Experiments were done with single crystal cubes of natural pyrite (Navajun, Logrono, Spain) that were cut and polished into slices of known weight and surface area. The slices were isothermally heated at atmospheric pressure in 99.99 percent CO2 (Coleman Instrument Grade) at either 412 C (685 K) or 465 C (738 K) for time periods up to 10 days. These two isotherms correspond to temperatures at about 6 km and 0 km altitude, respectively, on Venus. The reaction rate was determined by measuring the weight loss of the reacted slices after removal from the furnace. The reaction products were characterized by X-ray diffraction, scanning electron microscopy, and energy dispersive spectroscopy on the SEM.

Fegley, B., Jr.↗

The Martian valley networks: Origin by niveo-fluvial processes

The valley networks may hold the key to unlocking the paleoclimatic history of Mars. These enigmatic landforms may be regarded as the Martian equivalent of the Rosetta Stone. Therefore, a more thorough understanding of their origin and evolution is required. However, there is still no consensus among investigators regarding the formation (runoff vs. sapping) of these features. Recent climatic modeling precludes warm (0 degrees C) globally averaged surface temperatures prior to 2 b.y. when solar luminosity was 25-30 percent less than present levels. This paper advocates snowmelt as the dominant process responsible for the formation of the dendritic valley networks. Evidence for Martian snowfall and subsequent melt has been discussed in previous studies.

Rice, J. W., Jr.↗

Transient Catalytic Combustor Model With Detailed Gas and Surface Chemistry

In this work, we numerically investigate the transient combustion of a premixed gas mixture in a narrow, perfectly-insulated, catalytic channel which can represent an interior channel of a catalytic monolith. The model assumes a quasi-steady gas-phase and a transient, thermally thin solid phase. The gas phase is one-dimensional, but it does account for heat and mass transfer in a direction perpendicular to the flow via appropriate heat and mass transfer coefficients. The model neglects axial conduction in both the gas and in the solid. The model includes both detailed gas-phase reactions and catalytic surface reactions. The reactants modeled so far include lean mixtures of dry CO and CO/H2 mixtures, with pure oxygen as the oxidizer. The results include transient computations of light-off and system response to inlet condition variations. In some cases, the model predicts two different steady-state solutions depending on whether the channel is initially hot or cold. Additionally, the model suggests that the catalytic ignition of CO/O2 mixtures is extremely sensitive to small variations of inlet equivalence ratios and parts per million levels of H2.

Struk, Peter M.↗

Nonlinear Aeroelastic Analysis Using a Time-Accurate Navier-Stokes Equations Solver

A method to simulate limit cycle oscillation (LCO) due to control surface freeplay using a modified CFL3D, a time-accurate Navier-Stokes computational fluid dynamics (CFD) analysis code with structural modeling capability, is presented. This approach can be used to analyze aeroelastic response of aircraft with structural behavior characterized by nonlinearity in the force verses displacement curve. A limited validation of the method, using very low Mach number experimental data for a three-degrees-of-freedom (pitch/plunge/flap deflection) airfoil model with flap freeplay, is also presented.

Kuruvila, Geojoe↗

The Crucial Role of Error Correlation for Uncertainty Modeling of CFD-Based Aerodynamics Increments

The Ares I ascent aerodynamics database for Design Cycle 3 (DAC-3) was built from wind-tunnel test results and CFD solutions. The wind tunnel results were used to build the baseline response surfaces for wind-tunnel Reynolds numbers at power-off conditions. The CFD solutions were used to build increments to account for Reynolds number effects. We calculate the validation errors for the primary CFD code results at wind tunnel Reynolds number power-off conditions and would like to be able to use those errors to predict the validation errors for the CFD increments. However, the validation errors are large compared to the increments. We suggest a way forward that is consistent with common practice in wind tunnel testing which is to assume that systematic errors in the measurement process and/or the environment will subtract out when increments are calculated, thus making increments more reliable with smaller uncertainty than absolute values of the aerodynamic coefficients. A similar practice has arisen for the use of CFD to generate aerodynamic database increments. The basis of this practice is the assumption of strong correlation of the systematic errors inherent in each of the results used to generate an increment. The assumption of strong correlation is the inferential link between the observed validation uncertainties at wind-tunnel Reynolds numbers and the uncertainties to be predicted for flight. In this paper, we suggest a way to estimate the correlation coefficient and demonstrate the approach using code-to-code differences that were obtained for quality control purposes during the Ares I CFD campaign. Finally, since we can expect the increments to be relatively small compared to the baseline response surface and to be typically of the order of the baseline uncertainty, we find that it is necessary to be able to show that the correlation coefficients are close to unity to avoid overinflating the overall database uncertainty with the addition of the increments.

Hemsch, Michael J.↗

Active control of wake/blade-row interaction noise through the use of blade surface actuators

A combined analytical/computational approach for controlling of the noise generated by wake/blade-row interaction through the use of anti-sound actuators on the blade surfaces is described. A representative two-dimensional section of a fan stage, composed of an upstream fan rotor and a downstream fan exit guide vane (FEGV), is examined. An existing model for the wakes generated by the rotor is analyzed to provide realistic magnitudes for the vortical excitations imposed at the inlet to the FEGV. The acoustic response of the FEGV is determined at multiples of the blade passing frequency (BPF) by using the linearized unsteady flow analysis, LINFLO. Acoustic field contours are presented at each multiple of BPF illustrating the generated acoustic response disturbances. Anti-sound is then provided by placing oscillating control surfaces, whose lengths and locations are specified arbitrarily, on the blades. An analysis is then conducted to determine the complex amplitudes required for the control surface motions to best reduce the noise. It is demonstrated that if the number of acoustic response modes to be controlled is equal to the number of available independent control surfaces, complete noise cancellation can be achieved. A weighted least squares minimization procedure for the control equations is given for cases in which the number of acoustic modes exceeds the number of available control surfaces. The effectiveness of the control is measured by the magnitude of a propagating acoustic response vector, which is related to the circumferentially averaged sound pressure level (SPL), and is minimized by a standard least-squares minimization procedure.

Kousen, Kenneth A.↗

Bayesian calibration of a physics-based crystal plasticity and damage model

In this work, we present a model parameter calibration procedure for a physics-based crystal plasticity model. The calibration process utilizes a powerful statistics-based Bayesian calibration method. Calibration of the crystal plasticity parameters makes use of experimentally-measured data, i.e. compressive stress–strain response, from $\langle$100$\rangle$ and $\langle$123$\rangle$ single crystal copper dynamically loaded via Kolsky bar tests. The calibration of damage parameters is achieved using experimentally-measured free-surface velocity history data from plate impact test on $\langle$100$\rangle$ and $\langle$110$\rangle$ single crystal copper, which generates shock compression followed by dynamic tensile failure. A validation assessment is then carried out by comparing the calibrated model predictions and experimental measurements of the dynamic tensile damage generated in an impacted bicrystal copper plate. Lastly, a model-informed rationale for the experimentally-observed dependence of the spatial distribution of ductile damage (porosity) on crystallography is provided.

42 ENGINEERING↗

Evolution of oxygen and stratification and their relationship in the North Pacific Ocean in CMIP6 Earth system models

Abstract. This study examines the linkages between the upper-ocean (0–200 m) oxygen (O2) content and stratification in the North Pacific Ocean using four Earth system models (ESMs), an ocean hindcast simulation, and an ocean reanalysis. The trends and variability in oceanic O2 content are driven by the imbalance between physical supply and biological demand. Physical supply is primarily controlled by ocean ventilation, which is responsible for the transport of O2-rich surface waters to the subsurface. Isopycnic potential vorticity (IPV), a quasi-conservative tracer proportional to density stratification that can be evaluated from temperature and salinity measurements, is used herein as a dynamical proxy for ocean ventilation. The predictability potential of the IPV field is evaluated through its information entropy. The results highlight a strong O2–IPV connection and somewhat higher (as compared to the rest of the basin) predictability potential for IPV across the tropical Pacific, where the El Niño–Southern Oscillation occurs. This pattern of higher predictability and strong anticorrelation between O2 and stratification is robust across multiple models and datasets. In contrast, IPV at mid-latitudes has low predictability potential and its center of action differs from that of O2. In addition, the locations of extreme events or hotspots may or may not differ between the two fields, with a strong model dependency, which persists in future projections. On the one hand, these results suggest that it may be possible to monitor ocean O2 in the tropical Pacific based on a few observational sites co-located with the more abundant IPV measurements; on the other, they lead us to question the robustness of the IPV–O2 relationship in the extratropics. The proposed framework helps to characterize and interpret O2 variability in relation to physical variability and may be especially useful in the analysis of new observation-based data products derived from the BGC-Argo float array in combination with the traditional but far more abundant Argo data.

Novi, Lyuba↗

Global Air Quality and Health Co-benefits of Mitigating Near-term Climate Change Through Methane and Black Carbon Emission Controls

Tropospheric ozone and black carbon (BC), a component of fine particulate matter (PM < or = 2.5 microns in aerodynamic diameter; PM2.5), are associated with premature mortality and they disrupt global and regional climate. Objectives: We examined the air quality and health benefits of 14 specific emission control measures targeting BC and methane, an ozone precursor, that were selected because of their potential to reduce the rate of climate change over the next 20-40 years. Methods: We simulated the impacts of mitigation measures on outdoor concentrations of PM2.5 and ozone using two composition-climate models, and calculated associated changes in premature PM2.5‑ and ozone-related deaths using epidemiologically derived concentration-response functions. Results: We estimated that, for PM2.5 and ozone, respectively, fully implementing these measures could reduce global population-weighted average surface concentrations by 23-34% and 7-17% and avoid 0.6-4.4 and 0.04-0.52 million annual premature deaths globally in 2030. More than 80% of the health benefits are estimated to occur in Asia. We estimated that BC mitigation measures would achieve approximately 98% of the deaths that would be avoided if all BC and methane mitigation measures were implemented, due to reduced BC and associated reductions of nonmethane ozone precursor and organic carbon emissions as well as stronger mortality relationships for PM2.5 relative to ozone. Although subject to large uncertainty, these estimates and conclusions are not strongly dependent on assumptions for the concentration-response function. Conclusions: In addition to climate benefits, our findings indicate that the methane and BC emission control measures would have substantial co-benefits for air quality and public health worldwide, potentially reversing trends of increasing air pollution concentrations and mortality in Africa and South, West, and Central Asia. These projected benefits are independent of carbon dioxide mitigation measures. Benefits of BC measures are underestimated because we did not account for benefits from reduced indoor exposures and because outdoor exposure estimates were limited by model spatial resolution.

air quality↗

Mass spectral imaging showing the plant growth-promoting rhizobacteria's effect on the Brachypodium awn

The plant growth-promoting rhizobacteria (PGPR) on the host plant surface play a key role in biological control and pathogenic response in plant functions and growth. However, it is difficult to elucidate the PGPR effect on plants. Such information is important in biomass production and conversion. Brachypodium distachyon (Brachypodium), a genomics model for bioenergy and native grasses, was selected as a C3 plant model; and the Gram-negative Pseudomonas fluorescens SBW25 ( P.) and Gram-positive Arthrobacter chlorophenolicus A6 ( A.) were chosen as representative PGPR strains. The PGPRs were introduced to the Brachypodium seed's awn prior to germination, and their possible effects on the seeding and growth were studied using different modes of time-of-flight secondary ion mass spectrometry (ToF-SIMS) measurements, including a high mass-resolution spectral collection and delayed image extraction. We observed key plant metabolic products and biomarkers, such as flavonoids, phenolic compounds, fatty acids, and auxin indole-3-acetic acid in the Brachypodium awns. Furthermore, principal component analysis and two-dimensional imaging analysis reveal that the Brachypodium awns are sensitive to the PGPR, leading to chemical composition and morphology changes on the awn surface. Our results show that ToF-SIMS can be an effective tool to probe cell-to-cell interactions at the biointerface. This work provides a new approach to studying the PGPR effects on awn and shows its potential for the research of plant growth in the future.

metabolite↗

Soil Moisture Initialization Error and Subgrid Variability of Precipitation in Seasonal Streamflow Forecasting

Offline simulations over the conterminous United States (CONUS) with a land surface model are used to address two issues relevant to the forecasting of large-scale seasonal streamflow: (i) the extent to which errors in soil moisture initialization degrade streamflow forecasts, and (ii) the extent to which a realistic increase in the spatial resolution of forecasted precipitation would improve streamflow forecasts. The addition of error to a soil moisture initialization field is found to lead to a nearly proportional reduction in streamflow forecast skill. The linearity of the response allows the determination of a lower bound for the increase in streamflow forecast skill achievable through improved soil moisture estimation, e.g., through satellite-based soil moisture measurements. An increase in the resolution of precipitation is found to have an impact on large-scale streamflow forecasts only when evaporation variance is significant relative to the precipitation variance. This condition is met only in the western half of the CONUS domain. Taken together, the two studies demonstrate the utility of a continental-scale land surface modeling system as a tool for addressing the science of hydrological prediction.

Streamflow Prediction↗

Towards Improved Decision Making for the Smarter Grid: Geomagnetic Disturbance Implicit Modeling Uncertainty Quantification

Geomagnetic disturbances can greatly disrupt the power grid operation. This is caused by low-frequency changes in the earth’s magnetic field, which upon interaction with the deep earth conductivity induces a changing electric field at the earth’s surface. This in turn induces quasi-direct currents which have the potential to cause protective devices to operate, widespread damage to high-voltage transformers, or voltage collapse due to induced reactive power losses. While these impacts have been known since the 1940s, only recently have standards for GMD response been implemented. Because these standards are heavily based on models, it is important to understand when it is appropriate to make certain modeling assumptions. In this work, implicit modeling uncertainty is quantified for the purpose of enhanced decision making for disturbance impact mitigation. Considering electromagnetic models, modeling assumptions of non-uniform grounding, imbalanced operation, transformer saturation, and geomagnetically induced currents, voltage estimation values are analyzed. Statistical characteristics of modeling uncertainty are obtained by analysis of synthetically generated data. Results statistically characterize geomagnetic disturbances modeling uncertainty, providing an important tool for impact mitigation.

grid analysis↗