SEARCH · Engineering Papers

Results for “Sensitivity formulations”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1996

An Electrothermodynamic Model with Distributed Properties for Effective Permittivities of Sea Ice

This paper presents a model to calculate the temperature dependence of effective permittivities for sea ice, a heterogeneous medium containing multiphase scatterers. With the strong permittivity fluctuation approach the model accounts for the electrodynamic scattering effect together with the quasi-static characteristics of multiple species and subspecies of inhomogeneities with distributed orientations, sizes, and shapes. Because of a preferential direction in the orientation distribution, the medium is effectively anisotropic. The size distribution is described with a probability density function in terms of normalized volumetric sizes. Scatterer shapes are nonuniform and have a general ellipsoidal form characterized by arbitrary axial ratios of correlation lengths which are related to physical geometries of the scatterers. In this formulation, sea ice consisting of solid ice, liquid brine, and gaseous inclusions is modeled to derive effective permittivities with thermodynamic phase redistribution and structural metamorphism. Theoretical results are in good agreement with experimental data at the C band frequency of 4.8 GHz for saline ice undergoing warming and cooling cycles. A competitive effect between the increase of liquid brine and the shape rounding of ellipsoidal scatterers at increasing temperatures explains the trend observed in measured data. Sensitivities of effective permittivities to structural and physical parameters characterizing sea ice are also studied.

Nghiem, S. V.↗

DOE OSTI2026

Instabilities and phase transitions in architected metamaterials: a gradient-enhanced continuum approach

Architected metamaterials such as foams and lattices exhibit a wide range of properties governed by microstructural instabilities and emerging phase transitions. Their macroscopic response–including energy dissipation during impact, large recoverable deformations, morphing between configurations, and auxetic behavior–remains difficult to capture with conventional continuum models, which often rely on discrete approaches that limit scalability. In this work, we propose a nonlocal continuum formulation that captures both stable and unstable responses of elastic architected metamaterials. The framework extends anisotropic hyperelasticity by introducing nonlocal variables and internal length scales reflective of microstructural features. Local polyconvex free-energy models are systematically augmented with two families of non-(poly)convex energies, enabling both metastable and bistable responses. Implementation in a finite element framework enables solution using a hybrid monolithic–staggered strategy. Simulations capture densification fronts, forward and reverse transitions, hysteresis loops, imperfection sensitivity, and globally coordinated auxetic modes. Overall, this framework provides a robust foundation for accelerated modeling of instability-driven phenomena in architected metamaterials, while enabling extensions to anisotropic, dissipative, and active systems as well as integration with data-driven and machine learning approaches.

42 ENGINEERING↗

DOE OSTI2024

Improving Stratocumulus Cloud Amounts in a 200‐m Resolution Multi‐Scale Modeling Framework Through Tuning of Its Interior Physics

Abstract High‐Resolution Multi‐scale Modeling Frameworks (HR)—global climate models that embed separate, convection‐resolving models with high enough resolution to resolve boundary layer eddies—have exciting potential for investigating low cloud feedback dynamics due to reduced parameterization and ability for multidecadal throughput on modern computing hardware. However low clouds in past HR have suffered a stubborn problem of over‐entrainment due to an uncontrolled source of mixing across the marine subtropical inversion manifesting as stratocumulus dim biases in present‐day climate, limiting their scientific utility. We report new results showing that this over‐entrainment can be partly offset by using hyperviscosity and cloud droplet sedimentation. Hyperviscosity damps small‐scale momentum fluctuations associated with the formulation of the momentum solver of the embedded large eddy simulation. By considering the sedimentation process adjacent to default one‐moment microphysics in HR, condensed phase particles can be removed from the entrainment zone, which further reduces entrainment efficiency. The result is an HR that can produce more low clouds with a higher liquid water path and a reduced stratocumulus dim bias. Associated improvements in the explicitly simulated sub‐cloud eddy spectrum are observed. We report these sensitivities in multi‐week tests and then explore their operational potential alongside microphysical retuning in decadal simulations at operational 1.5° exterior resolution. The result is a new HR having desired improvements in the baseline present‐day low cloud climatology, and a reduced global mean bias and root mean squared error of absorbed shortwave radiation. We suggest it should be promising for examining low cloud feedbacks with minimal approximation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Exploring the impacts of unprecedented climate extremes on forest ecosystems: hypotheses to guide modeling and experimental studies

Abstract. Climatic extreme events are expected to occur more frequently in the future, increasing the likelihood of unprecedented climate extremes (UCEs) or record-breaking events. UCEs, such as extreme heatwaves and droughts, substantially affect ecosystem stability and carbon cycling by increasing plant mortality and delaying ecosystem recovery. Quantitative knowledge of such effects is limited due to the paucity of experiments focusing on extreme climatic events beyond the range of historical experience. Here, we present a road map of how dynamic vegetation demographic models (VDMs) can be used to investigate hypotheses surrounding ecosystem responses to one type of UCE: unprecedented droughts. As a result of nonlinear ecosystem responses to UCEs that are qualitatively different from responses to milder extremes, we consider both biomass loss and recovery rates over time by reporting a time-integrated carbon loss as a result of UCE, relative to the absence of drought. Additionally, we explore how unprecedented droughts in combination with increasing atmospheric CO2 and/or temperature may affect ecosystem stability and carbon cycling. We explored these questions using simulations of pre-drought and post-drought conditions at well-studied forest sites using well-tested models (ED2 and LPJ-GUESS). The severity and patterns of biomass losses differed substantially between models. For example, biomass loss could be sensitive to either drought duration or drought intensity depending on the model approach. This is due to the models having different, but also plausible, representations of processes and interactions, highlighting the complicated variability of UCE impacts that still need to be narrowed down in models. Elevated atmospheric CO2 concentrations (eCO2) alone did not completely buffer the ecosystems from carbon losses during UCEs in the majority of our simulations. Our findings highlight the consequences of differences in process formulations and uncertainties in models, most notably related to availability in plant carbohydrate storage and the diversity of plant hydraulic schemes, in projecting potential ecosystem responses to UCEs. We provide a summary of the current state and role of many model processes that give way to different underlying hypotheses of plant responses to UCEs, reflecting knowledge gaps which in future studies could be tested with targeted field experiments and an iterative modeling–experimental conceptual framework.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1979

The performance of a seasonal global climatic model

A seasonal global climatic model developed by Sellers is studied. Several changes are made to correct errors in the original model and to allow its use in a time dependent simulation mode. The major changes are in the infrared radiation formulation and in the size of the time steps. The seasonal cycles of surface temperature, wind, ice, snow, albedo, horizontal heat transport, and vertical flux of solar and infrared radiation are compared to recent observations. The seasonal cycle of temperature is fairly well reproduced, but has too small an amplitude in the Northern Hemisphere high latitudes. An overestimate of albedo at the poles is found, which is related to errors in some of the other variables. The sensitivity of the model to CO2 and solar constant (O) changes is summarized. Doubling the CO2 amount causes a 2.45 K increase in global average surface temperature (T). Increasing Q by 1% increases T by 3.14K. Decreasing 0 by 1% decreases T by 4.65K. An annual average version of the model is more sensitive to changes in Q than the model with the seasonal cycle.

Robock, A.↗

NASA NTRS1982

Sensitivity of constrained linear inversions to the selection of the Lagrange multiplier

The influence of the choice of the Lagrange multiplier on constrained linear inversions is explored, with reference made to applications in inferring the columnar aerosol size distributions from spectral aerosol optical depth measurements. A range of the Lagrange multiplier is examined to find all positive solutions for the solution vector, which represents modifying factors to the assumed form of the size distribution. An iterative method is devised to constrain the calculations to consideration of only positive quantities and a requirement that the regression fit to data be consistent with measurement errors. The determination of the variances and covariances is formulated and applied to existing data sets for optical depth. Variances in the solution are found to be large for particle radii when the information content of the data is small.

King, M. D.↗

NASA NTRS1984

Shape design sensitivity analysis of built-up structures

Selection of the best shape of a fillet in a tension bar such that no yielding occurs has long attracted the attention of engineers. Dimensions and notations for the bar and fillet are shown. With symmetry, only the upper half of the bar is considered. The optimal design problem is to find a boundary shape to minimize the total area of the fillet such that no yielding occurs. The classical boundary value problem is reduced to a variational or energy related problem which not only has excellent properties of existence and uniqueness but also provides the mathematical foundation for finite element analysis. The variational formulation may be viewed as the principle of virtual work and the finite element method as a application of the Galerkin method to the variational equation for approximate solution of the boundary value problem.

Choi, K. K.↗

DOE OSTI2020

Quantitative comparisons of electron-scale turbulence measurements in NSTX via synthetic diagnostics for high-k scattering.

Two synthetic diagnostics are implemented for the high-k scattering system in NSTX [D. R. Smith, E. Mazzucato, W. Lee, H. K. Park, C. W. Domier, and N. C. Luhmann Jr., Rev. Sci. Instrum., 79, 123501 (2008)] allowing direct comparisons between synthetic and experimentally detected frequency and wavenumber spectra of electron-scale turbulence fluctuations. Synthetic diagnostics were formulated in real-space and in wavenumber space, and were deployed to realistic electron-scale simulations carried out with the GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186 (2003) 545]. A highly unstable ETG regime (electron temperature gradient mode) in a modest beta NSTX NBI-heated H-mode discharge was chosen for the analysis. Mapping the measured wavenumbers to field aligned coordinates showed that the high-k system is sensitive to fluctuations that are closer to the spectral peak in the density fluctuation wavenumber spectrum (streamers) than originally predicted. Analysis of synthetic spectra show that the frequency response of the detected fluctuations is dominated by Doppler shift and insensitive to the turbulence drive. The shape of the high-k density fluctuation wavenumber spectrum is sensitive to the ETG turbulence drive conditions, and can be reproduced in a sensitivity scan of the most pertinent turbulent drive terms in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Hypersonic flutter of a curved shallow panel with aerodynamic heating

The general equations describing the nonlinear fluttering oscillations of shallow, curved, heated orthotropic panels have been derived. The formulation takes into account the location of the panel on the surface of a generic hypersonic vehicle, when calculating the aerodynamic loads. It is also shown that third order piston theory produces unsteady aerodynamic loading which is in close agreement with that based upon direct solution of the Euler equations. Results, for simply supported panels, are obtained using Galerkin's method combined with direct numerical integration in time to compute stable limit cycle amplitudes. These results illustrate the sensitivity of the aeroelastic behavior to the unsteady aerodynamic assumptions, temperature, orthotropicity and flow orientation.

Bein, T.↗

DOE OSTI2021

Developing Technology Performance Level Assessments for Early-Stage Wave Energy Converter Technologies: Preprint

The advantage of using Technology Performance Level (TPL) in conjunction with Technology Readiness Level (TRL) assessments in guiding technology development trajectories to successful outcomes in less time, at less overall cost, and with less encountered risk has been well articulated in the literature. In partnership with industry and international collaborators, a TPL assessment methodology for grid-connected applications has been developed through the application of the systems engineering approach. Metrics under seven different categories have been developed, weighted based on their relative relevance, and combined to yield a composite score. The methodology has been implemented in a spreadsheet tool plus a web application specifically aimed at assessing early stage (TRL 1-3) concepts. The target use cases are (a) technology developers improving their design, to find fatal flaws early, to get feedback on current design, to identify areas of improvement that will yield the highest return on investment, (b) reviewers assessing technologies in competitions or for making funding decisions, (c) investor or project developer doing due diligence, (d) policy makers landscaping the technology domain for formulating R&D strategy. The methodology and the tools are undergoing continuous improvement based on the experience and lessons learnt from applying it to internal and external marine energy technology development projects. The methodology is also being adapted for assessing WECs servicing markets outside the continental grid - broadly categorized as Powering the Blue Economy (PBE) applications. Such applications have vastly different functional requirements entailing a modification of the methodology to account for their higher risk tolerance, reduced price sensitivities, lower power needs, different permitting protocols, etc. This paper presents the latest status of the TPL assessment methodology and tools, describes its adaptation to select PBE markets, and explores its extension to other domains where it could provide a comprehensive and holistic measure of a nascent or disruptive technology's technoeconomic performance potential.

metrics↗

DOE OSTI2021

Developing Technology Performance Level Assessments for Early-Stage Wave Energy Converter Technologies

metrics↗

DOE OSTI2022

High Volume Production of Nanoscale Elements

The utility of nanoscale elements spans numerous potential markets including healthcare, communications, energy, optics, and structures. The energy market in particular will require large volumes and is price sensitive. In this talk we will discuss strategies to scale nanoscale elements, including nanoparticle inks and nanoscale thin films, for the large area deposition (hundreds of square meters) using high through put processes. Historically, roll-to-roll production of thin films has been carried out at very high speeds over large areas using solution phase inks. Nanoparticles are produced in a bottom-up approach in solution phase that simplifies the processing and is amenable to ink formulations. Sintering and annealing of the thin films can be increased by using radiation thermal processing that integrates easily into roll-to-roll manufacturing. Results of printed films on the order of tens of nanometers to hundreds of nanometers are presented including studies of morphologies and uniformity over large areas. Solar cells devices using the roll-to-roll setup are demonstrated and a techno-economic analysis of the manufacturing is proposed. These results are presented with a comparison to the current photovoltaic technologies in production.

Druffel, Thad↗

NASA NTRS1989

Non-linear three dimensional spectral model of the Venusian thermosphere with super-rotation. I - Formulation and numerical technique. II - Temperature, composition and winds

A three-dimensional spectral model for the Venusion thermosphere is presented which uses spherical harmonics to represent the horizontal variations in longitude and latitude and which uses Fourier harmonics to represent the LT variations due to atmospheric rotation. A differencing scheme with tridiagonal block elimination is used to perform the height integration. Quadratic nonlinearities are taken into account. In the second part, numerical results obtained with the model are shown to reproduce the observed broad daytime maxima in CO2 and CO and the significantly larger values at dawn than at dusk. It is found that the diurnal variations in He are most sensitive to thermospheric superrotation, and that, given a globally uniform atmosphere as input, larger heating rates yield a larger temperature contrast between day and night.

Stevens-Rayburn, D. R.↗

DOE OSTI2026

PDE-constrained high-order mesh optimization

Here, we present a novel framework for PDE-constrained r-adaptivity of high-order meshes. The proposed method formulates mesh movement as an optimization problem, with an objective function defined as a convex combination of a mesh quality metric and a measure of the accuracy of the PDE solution obtained via finite element discretization. The proposed formulation achieves optimized, well-defined high-order meshes by integrating mesh quality control, PDE solution accuracy, and robust gradient regularization. We adopt the Target-Matrix Optimization Paradigm to control geometric properties across the mesh, independent of the PDE of interest. To incorporate the accuracy of the PDE solution, we introduce error measures that control the finite element discretization error. The implicit dependence of these error measures on the mesh nodal positions is accurately captured by adjoint sensitivity analysis. Additionally, a convolution-based gradient regularization strategy is used to ensure stable and effective adaptation of high-order meshes. We demonstrate that the proposed framework can improve mesh quality and reduce the error by up to 10 times for the solution of Poisson and linear elasto-static problems. The approach is general with respect to the dimensionality, the order of the mesh, the types of mesh elements, and can be applied to any PDE that admits well-defined adjoint operators.

Computer science↗

DOE OSTI2025

Formation trajectories of solution-processed perovskite thin films from mixed solvents

The engineering of mixed-solvent formulations and their evaporation conditions are key to reproducible perovskite coatings for high-performance photovoltaics. Here, we report a lumped-parameter evaporation model to predict the evolution of a perovskite ink liquid film over time (solvent ratio, solute concentration, and film thickness). The drying-rate model is validated via in situ film-thickness measurements, and the predicted transient liquid film state is mapped as a process path. These methods allow for the prediction of process sensitivity to local environmental factors and the understanding and visualization of a broader processing parameter space enabled through the coupling of process and ink engineering. Process maps are applied to create a new framework for scalable perovskite coating development with a goal of improving the reproducibility and transferability of perovskite fabrication. This approach is demonstrated with blade-coated FA 0.83 Cs 0.17 PbI 3 photovoltaic devices, improving the photovoltaic conversion efficiency from 17.5% ± 1.7% to 20.3% ± 0.6%.

14 SOLAR ENERGY↗

DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1975

Numerical techniques for the linear, nonadiabatic stellar pulsation problem

The linear, nonadiabatic eigenvalue problem is formulated using Castor's method for calculating radial pulsations of stellar models. Both left and right eigenvectors are calculated. Initial eigenvalues for the linear, nonadiabatic solutions are obtained from the adiabatic eigenvalues and left and right eigenvectors. The orthogonality relation is obtained. Simple formulas for the Newton method are given. The iteration procedure is constrained to improve convergence. The Newton method is less satisfactory than the secant method for difficult cases. The linear, nonadiabatic solutions are shown to be sensitive to the number of zones with tau smaller than 2/3, and the value of (P/P sub r) surface or tau (surface). Optimum values can be determined for the number of zones and tau (surface). Application of the method to Population II Cepheids is briefly presented.

Bednarek, T. A.↗

NASA NTRS2008

Rightsizing LISA

The LISA science requirements and conceptual design have been fairly stable for over a decade. In the interest of reducing costs, the LISA Project at NASA has looked for simplifications of the architecture, at downsizing of subsystems, and at descopes of the entire mission. This is a natural activity of the formulation phase, and one that is particularly timely in the current NASA budgetary context. There is, and will continue to be, enormous pressure for cost reduction from both ESA and NASA, reviewers and the broader research community. Here, we review the rationale for the baseline architecture. and report recent efforts to find simplifications and other reductions that might lead to savings. A few possible simplifications have been found in the LISA baseline architecture. In the interest of exploring cost sensitivity, one moderate and one aggressive descope have been evaluated; the cost savings are modest and the loss of science is not.

Stebbins, Robin T.↗