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.

SourceResearch connectionFrom yearThrough year

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

At least 1,729 records · Page 96

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Unstable Air-Sea Interaction in the Extratropical North Atlantic

The possibility of coupled modes in the extratropical North Atlantic has fascinated the climate community since 1960's. A significant aspect of such modes is an unstable air-sea interaction, also called positive feedback, where disturbances between the atmosphere and ocean grow unbound. If a delayed response exists before the negative feedback takes effect, an oscillatory behaviour will develop. Here we explore the relationship between heat flux (positive upward) and sea surface temperature (SST). Positive feedback is characterized by a cross-correlation between the two where correlation maintains a negative sign whether SST or heat flux leads. We use model results and observations to argue that in the North Atlantic there exist regions with positive feedback. The two main locations coincide with the well-known north-south SST dipole where anomalies of opposite sign occupy areas east of Florida and north-east of Newfoundland. We show that oceanic dynamics, wave propagation and advection, give rise to oceanic anomalies in these regions. Subsequently these anomalies are amplified by atmosphere- ocean interaction: thus a positive feedback.

Hakkinen, Sirpa↗

NASA NTRS2011

Ablator Modeling: Why Not Much Has Changed Over the Past 45 plus Years

In the late 1950s, the earliest models describing the thermal response of ablative materials were based on the heat of ablation concept, which is an empirical approach that was reasonable for the types of materials of interest at that time. In the early-mid 60s the models were expanded to include pyrolysis since organic resin composites became the TPS materials of interest. However, surface recession was still predominantly modeled via empirical correlation. The development of the 1-D CMA finite difference code in the mid-late 60s introduced the thermochemical ablation approach for gas/surface interactions. Since that time investigators have developed finite volume and finite element codes, in 1-D, 2-D and 3-D, but the basic modeling has not evolved significantly. Models describing internal gas pressure due to pyrolysis, particle impact erosion, in-depth radiant transport, etc., have been added to address specific problems, but the fundamental modeling has not evolved. The reasons for this stagnation, as viewed by the author, will be described.

Laub, B.↗

NASA NTRS2005

Comprehensive Micromechanics-Analysis Code - Version 4.0

Version 4.0 of the Micromechanics Analysis Code With Generalized Method of Cells (MAC/GMC) has been developed as an improved means of computational simulation of advanced composite materials. The previous version of MAC/GMC was described in "Comprehensive Micromechanics-Analysis Code" (LEW-16870), NASA Tech Briefs, Vol. 24, No. 6 (June 2000), page 38. To recapitulate: MAC/GMC is a computer program that predicts the elastic and inelastic thermomechanical responses of continuous and discontinuous composite materials with arbitrary internal microstructures and reinforcement shapes. The predictive capability of MAC/GMC rests on a model known as the generalized method of cells (GMC) - a continuum-based model of micromechanics that provides closed-form expressions for the macroscopic response of a composite material in terms of the properties, sizes, shapes, and responses of the individual constituents or phases that make up the material. Enhancements in version 4.0 include a capability for modeling thermomechanically and electromagnetically coupled ("smart") materials; a more-accurate (high-fidelity) version of the GMC; a capability to simulate discontinuous plies within a laminate; additional constitutive models of materials; expanded yield-surface-analysis capabilities; and expanded failure-analysis and life-prediction capabilities on both the microscopic and macroscopic scales.

Arnold, S. M.↗

NASA NTRS2012

The Aerosol-Monsoon Climate System of Asia

In Asian monsoon countries such as China and India, human health and safety problems caused by air-pollution are worsening due to the increased loading of atmospheric pollutants stemming from rising energy demand associated with the rapid pace of industrialization and modernization. Meanwhile, uneven distribution of monsoon rain associated with flash flood or prolonged drought, has caused major loss of human lives, and damages in crop and properties with devastating societal impacts on Asian countries. Historically, air-pollution and monsoon research are treated as separate problems. However a growing number of recent studies have suggested that the two problems may be intrinsically intertwined and need to be studied jointly. Because of complexity of the dynamics of the monsoon systems, aerosol impacts on monsoons and vice versa must be studied and understood in the context of aerosol forcing in relationship to changes in fundamental driving forces of the monsoon climate system (e.g. sea surface temperature, land-sea contrast etc.) on time scales from intraseasonal variability (~weeks) to climate change (~ multi-decades). Indeed, because of the large contributions of aerosols to the global and regional energy balance of the atmosphere and earth surface, and possible effects of the microphysics of clouds and precipitation, a better understanding of the response to climate change in Asian monsoon regions requires that aerosols be considered as an integral component of a fully coupled aerosol-monsoon system on all time scales. In this paper, using observations and results from climate modeling, we will discuss the coherent variability of the coupled aerosol-monsoon climate system in South Asia and East Asia, including aerosol distribution and types, with respect to rainfall, moisture, winds, land-sea thermal contrast, heat sources and sink distributions in the atmosphere in seasonal, interannual to climate change time scales. We will show examples of how elevated absorbing aerosols (dust and black carbon) may interact with monsoon dynamics to produce feedback effects on the atmospheric water cycle, leading to in accelerated melting of snowpacks over the Himalayas and Tibetan Plateau, and subsequent changes in evolution of the pre-monsoon and peak monsoon rainfall, moisture and wind distributions in South Asia and East Asia.

Lau, William K. M.↗

NASA NTRS2013

Quantifying Effects of Voids in Woven Ceramic Matrix Composites

Randomness in woven ceramic matrix composite architecture has been found to cause large variability in stiffness and strength. The inherent voids are an aspect of the architecture that may cause a significant portion of the variability. A study is undertaken to investigate the effects of many voids of random sizes and distributions. Response surface approximations were formulated based on void parameters such as area and length fractions to provide an estimate of the effective stiffness. Obtaining quantitative relationships between the properties of the voids and their effects on stiffness of ceramic matrix composites are of ultimate interest, but the exploratory study presented here starts by first modeling the effects of voids on an isotropic material. Several cases with varying void parameters were modeled which resulted in a large amount of variability of the transverse stiffness and out-of-plane shear stiffness. An investigation into a physical explanation for the stiffness degradation led to the observation that the voids need to be treated as an entity that reduces load bearing capabilities in a space larger than what the void directly occupies through a corrected length fraction or area fraction. This provides explanation as to why void volume fraction is not the only important factor to consider when computing loss of stiffness.

Goldsmith, Marlana B.↗

NASA NTRSDate not supplied

TPSAS-NF1676L-17931-DND

Many geophysical variables including temperature, clouds, and precipitation exhibit robust diurnal cycles in response to the daily cycle of solar insolation. Due to the fundamental nature of this variability, it is critical that weather and climate models accurately represent the diurnal cycle. Numerical model, however, have difficulty reproducing this observed diurnal cycle behavior, which leads to systematic errors in model representation of earth energy budget terms: including TOA and surface radiation, precipitation, and surface latent and sensible heat fluxes. This study evaluates the regional diurnal cycle in the Tropics within reanalysis models (ERA-Interim and MERRA) and quantifies systematic errors in the simulated TOA flux and precipitation due to biases in the diurnal cycle representation. The focus of this study is to quantify the importance of diurnal cycle simulation to systematic bias in reanalysis data set climatologies.

Patrick C Taylor↗

DOE OSTI2020

Toroidal modeling of runaway electron loss due to 3-D fields in DIII-D and COMPASS

The 3-D field induced relativistic runaway electron (RE) loss has been simulated for DIII-D, COMPASS and ITER plasmas, utilizing the MARS-F code incorporated with the recently developed and updated RE orbit module (REORBIT). Modeling shows effectively 100% loss of a post-disruption, high-current runaway beam in DIII-D, due to the 1 kG level of magnetic field perturbation produced by a fast growing n = 1 resistive kink instability. The RE loss is shown to be independent of the particle energy or the initial location of particles in the configuration space. Applied resonant magnetic perturbation (RMP) fields from in-vessel coils are not effective for RE beam mitigation in DIII-D, but do produce finite (>10%) RE loss in COMPASS, consistent with experimental observations in above two devices. The major reasons for this difference in RE control by RMP between these two devices are (i) the coil proximity to the RE beam and (ii) the effective coil current scaling versus the machine size and the toroidal magnetic field. About 10% RE loss fraction is also predicted for an ITER 15 MA scenario with pre-disruption plasma, highlighting the role of the plasma response. Up to 30% loss is computed, however, by artificially scaling the equilibrium pressure to zero. This is due to the more resistive plasma response and stronger resulting field line stochasticity. Distributions of the lost REs to the limiting surface show poloidally peaked profile near the high-field-side in both DIII-D and COMPASS, covering about 100o poloidal angle. Higher perturbation field level and/or higher particle energy also result in REs being lost to the low-field-side of the limiting surface of these two devices, increasing the effective wetted area. Finally, for the modeled ITER plasmas, the applied RMP field, with optimal poloidal spectrum that maximizes the edge localized mode control, leads to REs being lost to the lower divertor region of the limiting surface within a narrow poloidal band.

3-D field↗

NASA NTRS2001

Modelling Greenland Outlet Glaciers

The objective of this project was to develop simple yet realistic models of Greenland outlet glaciers to better understand ongoing changes and to identify possible causes for these changes. Several approaches can be taken to evaluate the interaction between climate forcing and ice dynamics, and the consequent ice-sheet response, which may involve changes in flow style. To evaluate the icesheet response to mass-balance forcing, Van der Veen (Journal of Geophysical Research, in press) makes the assumption that this response can be considered a perturbation on the reference state and may be evaluated separately from how this reference state evolves over time. Mass-balance forcing has an immediate effect on the ice sheet. Initially, the rate of thickness change as compared to the reference state equals the perturbation in snowfall or ablation. If the forcing persists, the ice sheet responds dynamically, adjusting the rate at which ice is evacuated from the interior to the margins, to achieve a new equilibrium. For large ice sheets, this dynamic adjustment may last for thousands of years, with the magnitude of change decreasing steadily over time as a new equilibrium is approached. This response can be described using kinematic wave theory. This theory, modified to pertain to Greenland drainage basins, was used to evaluate possible ice-sheet responses to perturbations in surface mass balance. The reference state is defined based on measurements along the central flowline of Petermann Glacier in north-west Greenland, and perturbations on this state considered. The advantage of this approach is that the particulars of the dynamical flow regime need not be explicitly known but are incorporated through the parameterization of the reference ice flux or longitudinal velocity profile. The results of the kinematic wave model indicate that significant rates of thickness change can occur immediately after the prescribed change in surface mass balance but adjustments in flow rapidly diminish these rates to a few cm/yr at most. The time scale for adjustment is of the order of a thousand years or so.

vanderVeen, Cornelis↗

DOE OSTI2021

Real-time estimation and control of divertor surface heat flux on the DⅢ-D tokamak

We report Future tokamaks will require robust technologies for the mitigation of heat exhaust onto the plasma-facing components. As a first step towards this development, a system has been developed at DIII-D that estimates and controls in real-time the heat flux to the PFCs. Real-time estimation of the peak power flux from this model-based approach is validated with off-line infra-red measurements for various DIII-D plasma discharges. A nonlinear free-boundary simulation code (GSevolve) is used for simulating the closed loop response and for the off-line determination of the control parameters. The implementation and first experimental results of the application during the DIII-D plasma campaign are reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1995

Three-dimensional position and shape of the bow shock and their variation with Alfvenic, sonic and magnetosonic Mach numbers and interplanetary magnetic field orientation

A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.

Peredo, M.↗

NASA NTRS1993

GCM simulations of volcanic aerosol forcing. I - Climate changes induced by steady-state perturbations

The response of the climate system to a temporally and spatially constant amount of volcanic particles is simulated using a general circulation model (GCM). The optical depth of the aerosols is chosen so as to produce approximately the same amount of forcing as results from doubling the present CO2 content of the atmosphere and from the boundary conditions associated with the peak of the last ice age. The climate changes produced by long-term volcanic aerosol forcing are obtained by differencing this simulation and one made for the present climate with no volcanic aerosol forcing. The simulations indicate that a significant cooling of the troposphere and surface can occur at times of closely spaced multiple sulfur-rich volcanic explosions that span time scales of decades to centuries. The steady-state climate response to volcanic forcing includes a large expansion of sea ice, especially in the Southern Hemisphere; a resultant large increase in surface and planetary albedo at high latitudes; and sizable changes in the annually and zonally averaged air temperature.

Pollack, James B.↗

DOE OSTI2025

Temperature and Strain Sensing Characteristics of a 128° YX-Cut LiNbO 3 Rayleigh-Mode SAW Sensor From Room to Cryogenic Temperatures

Accurate, passive, and wireless monitoring of cryogenic hardware is essential for high-energy physics, space propulsion, and biomedical instrumentation. Here, this study quantifies the coupled temperature-strain behavior of Rayleigh-mode surface acoustic-wave (SAW) delay-line sensors fabricated on 128° YX-cut LiNbO 3 . A nonlinear finite element (FE) model incorporating Varshni-based elastic constants, higher-order thermal expansion, and temperature-dependent piezo- and dielectric coefficients was developed and validated experimentally between 280K and 80K. Free-standing (first test condition) and bonded/wired (second test condition) devices exhibited indistinguishable thermal responses; the average temperature coefficient of delay (TCD) in the critical cryogenic range from 130K down to 80K differed by only 0.15 ppm/K (0.32%), confirming that bonding-induced stress is negligible. Over 280-80K the measured TCD was 61.77 ppm/K, while the FE model predicted an equivalent temperature coefficient of frequency (TCF) of −62.74 ppm/K with an overall coefficient of determination R2 = 0.998. In the critical cryogenic interval 130-80K the TCD fell to 47.66 ppm/K, indicating improved thermal stability at low temperature. Controlled loading (0-300 με) revealed a strain coefficient of delay (SCD) that rises from 0.53 ± 0.02 ppm/με at 300K to 1.05 ± 0.02 ppm/με at 80K. This modest sensitivity confirms that, for temperature sensing, strain is a second-order perturbation above 135K but must be compensated at deeper cryogenic levels. Overall, this work establishes a predictive multiphysics model together with repeatable wired measurements that confirm the suitability of SAW sensors for temperature and strain monitoring in extreme cryogenic environments, while also providing a baseline for future wireless implementations.

25 ENERGY STORAGE↗

DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗

DOE OSTI2022

Plasma-Catalyst Reactivity Control of Surface Nitrogen Species through Plasma-Temperature-Programmed Hydrogenation to Ammonia

Nonthermal plasma activation of N 2 can facilitate nitrogen adsorption on metal catalysts at low bulk temperatures and atmospheric pressure. Here, we apply a plasma-assisted temperature-programmed reaction (plasma-TPRxn) for ammonia (NH 3 ) synthesis using sequential exposure of a silica-supported metal catalyst to N 2 plasma followed by thermal hydrogen treatment while ramping the temperature to decouple the plasma activation of N 2 from surface catalyzed hydrogenation steps. This approach eliminates the effects from bulk plasma phase reactions, thereby allowing for direct interrogation of plasma activated nitrogen on the active metal surfaces. We confirm previously reported spectroscopic observations that show plasma-generated surface nitrogen can be converted to NH 3 through surface catalyzed pathways. Further, we demonstrate that the ammonia desorption peak temperature is sensitive to metal, with Pt desorbing NH 3 at the lowest temperature. Unsteady state microkinetic models of desorption kinetics as a function of initial N coverage and metal recover observed trends in NH 3 desorption temperatures and confirm that observed results reflect hydrogenation of plasma-induced N accommodation at each surface. In total, we show that the hydrogenation ability of the catalyst after plasma activation of N 2 is responsible for the reactivity trends observed in plasma-assisted NH 3 synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Intensification of future subsurface marine heatwaves in an eddy-resolving model

A shift in depth range enables marine organisms to adapt to marine heatwaves (MHWs). Subsurface MHWs could limit this pathway, yet their response to climate warming remains unclear. Here, using an eddy-resolving Earth system model forced under a high emission scenario, we project a robust global increase in subsurface MHWs driven by rising subsurface mean temperatures and enhanced temperature variability. Historically, maximum MHW intensity occurs around 100m depth, which shifts to the faster-warming surface under greenhouse warming. However, removing the long-term warming trend yields an increase in subsurface MHW intensity and annual days greater than that at the surface, especially in large marine ecosystem regions, primarily due to increased variability. Additionally, days of the surface and subsurface concurrent event increase ten times more than those of individual events. Our study highlights a heightened threat to marine organisms under global warming, as the increased subsurface heatwaves reduce their refuge options.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

High frequency beam oscillation keyhole dynamics in laser melting revealed by in-situ x-ray imaging

The metal additive manufacturing industry is actively developing instruments and strategies to enable higher productivity, optimal build quality, and controllable as-built microstructure. A beam controlling technique, laser oscillation has shown potential in all these aspects in laser welding; however, few attempts have been made to understand the underlying physics of the oscillating keyholes/melt pools which are the prerequisites for these strategies to become a useful tool for laser-based additive manufacturing processes. Here, to address this gap, we utilized a synchrotron-based X-ray operando technique to image the dynamic keyhole oscillation in Ti-6Al-4V using a miniature powder bed fusion setup. We found good agreement between the experimental observations and simulations performed with a validated Lattice Boltzmann multiphysics model. The study revealed the continuous and periodic fluctuations in the characteristic keyhole parameters that are unique to the oscillating laser beam processing and responsible for the chevron pattern formation at solidification. In particular, despite the intrinsic longer-range fluctuation, the oscillating technique displayed potential for reducing keyhole instability, mitigating porosity formation, and altering surface topology. These insights on the oscillating keyhole dynamics can be useful for the future development and application of this technique.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Utilization of coupled eigenmodes in Akiyama atomic force microscopy probes for bimodal multifrequency sensing

Akiyama atomic force microscopy probes represent a unique means of combining several of the desirable properties of tuning fork and cantilever probe designs. As a hybridized mechanical resonator, the vibrational characteristics of Akiyama probes result from a complex coupling between the intrinsic vibrational eigenmodes of its constituent tuning fork and bridging cantilever components. Here, through a combination of finite element analysis modeling and experimental measurements of the thermal vibrations of Akiyama probes, we identify a complex series of vibrational eigenmodes and measure their frequencies, quality factors, and spring constants. We then demonstrate the viability of Akiyama probes to perform bimodal multi-frequency force sensing by performing a multimodal measurement of a surface's nanoscale photothermal response using photo-induced force microscopy imaging techniques. Further performing a parametric search over alternative Akiyama probe geometries, we propose two modified probe designs to enhance the capability of Akiyama probes to perform sensitive bimodal multifrequency force sensing measurements.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2024

Radiological Safety Analysis Computer Program V 7.2

Users generates a fission product inventory from either reactor operating history or a nuclear criticality event. RSAC-7 models the effects of high-efficiency particulate air filters or other cleanup systems and calculates the decay and ingrowth during transport through processes, facilities, and the environment. Doses are calculated for inhalation, air immersion, ground surface, ingestion, and cloud gamma pathways. RSAC-7 is used as a tool to evaluate accident conditions in emergency response scenarios, radiological sabotage events, and safety basis accident consequences.

Jordan, RobertaA. [Idaho National Laboratory (INL)↗