Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Form”

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 127 records · Page 7

Thermal Modelling of the In-Situ Consolidation of Automated Fiber Placement Of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Christopher J. Stelter↗

Elastic Response of Crimped Collagen Fibrils

A physiologic constitutive expression is presented in algorithmic format for the elastic response of wavy collagen fibrils found in soft connective tissues. The model is based on the observation that crimped fibrils have a three-dimensional structure at the micrometer scale that we approximate as a helical spring. The symmetry of this waveform allows the force/displacement relationship derived from Castigliano's theorem to be solved in closed form. Model predictions are in good agreement with experimental observations for mitral-valve chordae tendineae

Freed, Alan D.↗

Linear Multivariable Regression Models for Prediction of Eddy Dissipation Rate from Available Meteorological Data

Linear multivariable regression models for predicting day and night Eddy Dissipation Rate (EDR) from available meteorological data sources are defined and validated. Model definition is based on a combination of 1997-2000 Dallas/Fort Worth (DFW) data sources, EDR from Aircraft Vortex Spacing System (AVOSS) deployment data, and regression variables primarily from corresponding Automated Surface Observation System (ASOS) data. Model validation is accomplished through EDR predictions on a similar combination of 1994-1995 Memphis (MEM) AVOSS and ASOS data. Model forms include an intercept plus a single term of fixed optimal power for each of these regression variables; 30-minute forward averaged mean and variance of near-surface wind speed and temperature, variance of wind direction, and a discrete cloud cover metric. Distinct day and night models, regressing on EDR and the natural log of EDR respectively, yield best performance and avoid model discontinuity over day/night data boundaries.

MCKissick, Burnell T.↗

Framework for Estimating Performance and Associated Uncertainty for Modified Aircraft Configurations

Flight testing has been the historical standard for determining aircraft airworthiness - however, increases in the cost of flight testing and the accuracy of inexpensive CFD promote certification by analysis to reduce or replace flight testing. A framework is introduced to predict the performance in the special case of a modification to an existing, previously certified aircraft. This framework uses a combination of existing flight test or high fidelity data of the original aircraft as well as lower fidelity data of the original and modified configurations. Two methods are presented which estimate the model form uncertainty of the modified configuration, which is then used to conduct non-deterministic simulations. The framework is applied to an example aircraft system with simulated flight test data to demonstrate the ability to predict the performance and associated uncertainty of modified aircraft configurations. However, it is important that the models and methods used are applicable and accurate throughout the intended use domain. The factors and limitations of the framework are explored to determine the range of applicability of the framework. The effects of these factors on the performance and uncertainty results are demonstrated using the example aircraft system. The framework is then applied to NASA’s X-57 Maxwell and each of its modifications. The estimated performance and associated uncertainties are then compared to the airworthiness criteria to evaluate the potential of the framework as a component to the certification by analysis process.

Uncertainty Quantification↗

Acoustic Investigation of Jet Mixing Noise in Dual Stream Nozzles

In an earlier study, a prediction model for jet noise in dual stream jets was proposed that is founded on velocity scaling laws in single stream jets and similarity features of the mean velocity and turbulent kinetic energy in dual stream flows. The model forms a composite spectrum from four component single-stream jets each believed to represent noise-generation from a distinct region in the actual flow. While the methodology worked effectively at conditions considered earlier, recent examination of acoustic data at some unconventional conditions indicate that further improvements are necessary in order to expand the range of applicability of the model. The present work demonstrates how these predictions compare with experimental data gathered by NASA and industry for the purpose of examining the aerodynamic and acoustic performance of such nozzles for a wide range of core and fan stream conditions. Of particular interest are jets with inverted velocity and temperature profiles and the appearance of a second spectral peak at small aft angles to the jet under such conditions. It is shown that a four-component spectrum succeeds in modeling the second peak when the aft angle refraction effects are properly incorporated into the model. A tradeoff of noise emission takes place between two turbulent regions identified as transition and fully mixed regions as the fan stream velocity exceeds that of the core stream. The effect of nozzle discharge coefficients will also be discussed.

Khavaran, Abbas↗

Bio-Optical Properties and Ocean Color Algorithms for Coastal Waters Influenced by the Mississippi River During a Cold Front Passage

During the passage of a cold front in March 2002, bio-optical properties examined in coastal waters impacted by the Mississippi River indicated westward advective flows and increasing river discharge containing a larger nonalgal particle content contributed significantly to surface optical variability. A comparison of seasonal data from three cruises indicated spectral models of absorption and scattering to be generally consistent with other coastal environments, while their parameterization in terms of chlorophyll a concentration (Chl) showed seasonal variability. The exponential slope of the colored dissolved organic matter (CDOM) averaged 0.0161 plus or minus 0.00054 per nanometer, and for nonalgal absorption it averaged 0.011 per nanometer with deviations from general trends observed due to anomalous water properties. Although the phytoplankton specific absorption coefficients varied over a wide range (0.02 to 0.1 square meters (mg Chl) sup -1)) being higher in offshore surface waters, values of phytoplankton absorption spectra at the SeaWiFS wavebands were highly correlated to modeled values. The normalized scattering spectral shapes and the mean spectrum were in agreement to observations in other coastal waters, while the backscattering ratios were on average lower in phytoplankton dominated surface waters (0.0101 plus or minus 0.002) and higher in near-bottom waters (0.0191 plus or minus 0.0045) with low Chl. Average percent differences in remote sensing reflectance R (sub rs) derived form modeled and in-eater radiometric measurements were highest in the blue wavebands (52%) and at sampling stations with a ore stratified water column. Estimates of Chl and CDOM absorption derived from SeaWiFS images generated using regional empirical algorithms were highly correlated to in situ data.

D'Sa Eurico J.↗

Determining the Mass Density Along Magnetic Field Lines from Toroida Eigenfrequencies

Toroidal eigenfrequencies can be used to remotely sense the equatorial mass density rho(sub eq) and the density dependence along a magnetic field line. Here we present improvements to the method of Schulz [1996], which allows rho(sub eq) and the power law index alpha (for mass density along a field line proportional to R(sup -alpha), where R is the radial distance from the center of the Earth) to be determined from the y intercept and slope of a plot of toroidal frequency versus toroidal harmonic number n. Our modifications include a model form for eigenfrequencies with a fractional precision of 0.0005 for -6 less than or = alpha less than or = 6 and 2 less than or = L less than or = 8 (accuracy is doubtful beyond L = 5) and an iterative procedure for getting more accurate results than those found using Schulz's method. In addition, we do an analysis of the effect of random measurement errors. Observed frequencies need to be accurate to approx. 6% (3%) of the fundamental frequency in order to determine rho(sub eq) (alpha) to a precision of 30% (unity). We then apply our method to data generated using the Global Core Plasma Model for plasmaspheric mass density; our analysis demonstrates clearly bow the alpha index represents the mass density dependence on the outer part of the field line (R/(LR(sub E)) greater than or approx. 2/3).

Denton, R.↗

A preliminary tidal analysis of TOPEX/POSEIDON altimetry

Approximately 12 months of data from the TOPEX/POSEIDON satellite altimeter mission are analyzed for the major short-period oceanic tides. A harmonic analysis is performed on data captured within bins defined on a deep-ocean grid, which, owing to tidal aliasing considerations, must have a relatively coarse spatial resolution. Our analysis is in terms of corrections to the Schwiderski and Cartwright-Ray models, and it confirms many of the Schwiderski differences previously reported by Cartwright and Ray. Our differences with respect to the Geosat-based Cartwright and Ray model form a sectorial pattern in M(sub 2) with high/low differences separated roughly 180 deg in longitude. We suggest that these sectorial errors were most likely induced by Geosat's relatively large orbit error. Comparisons to independent data validate the improved TOPEX/POSEIDON solutions; in-situ 'ground truth' shows M(sub 2) root mean square (RMS) differences of 4.10 cm (Schwiderski), 3.86 cm (cartwright and Ray), 2.63 cm (this paper). Global rates of energy dissipation confirm earlier estimates for M(sub 2), and show improved agreement with satellite tracking studies for K(sub 1) and S(sub 2). These preliminary exercises confirm that TOPEX/POSEIDON should results in a new generation of improved global tidal models.

Schrama, E. J. O.↗

Modeling of Stress and Temperature Effects on Creep of Reduced Activation Ferritic-Martensitic Steel Alloy F82H (Tertiary Creep Modeling of RAFM Steel)

A Bayesian optimization procedure is presented for calibrating a multi-mechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600 °C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework [1]. Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. Validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel at 600 °C. The model is further extended using a temperature dependent scaling law approach to simulate creep response between 550 °C and 650 °C. The efficacy of this extension is compared with the previously used scaling law approach for Grade 91 steel.

36 MATERIALS SCIENCE↗

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗

Calibration of RAFM Micromechanical Model for Creep Using Bayesian Optimization for Functional Output

A Bayesian optimization procedure is presented for calibrating a multimechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on Fe(8–9)%Cr are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600°C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework (Huang et al., 2021, “Bayesian optimization of functional output in inverse problems,” Optim. Eng., 22, pp. 2553–2574). Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. In conclusion, validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel.

42 ENGINEERING↗

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗

Subject-specific multi-scale modeling of the fate of inhaled aerosols

Determining the fate of inhaled aerosols in the respiratory system is essential in assessing the potential toxicity of inhaled airborne materials, responses to airborne pathogens, or in improving inhaled drug delivery. The availability of high-resolution clinical lung imaging and advances in the reconstruction of lung airways from CT images have led to the development of subject-specific in-silico 3D models of aerosol dosimetry, often referred to as computational fluid-particle-dynamics (CFPD) models. As CFPD models require extensive computing resources, they are typically confined to the upper and large airways. These models can be combined with lower-dimensional models to form multiscale models that predict the transport and deposition of inhaled aerosols in the entire respiratory tract. Understanding where aerosols deposit is only the first of potentially several key events necessary to predict an outcome, being a detrimental health effect or a therapeutic response. To that end, multiscale approaches that combine CFPD with physiologically-based pharmacokinetics (PBPK) models have been developed to evaluate the absorption, distribution, metabolism, and excretion (ADME) of toxic or medicinal chemicals in one or more compartments of the human body. CFPD models can also be combined with host cell dynamics (HCD) models to assess regional immune system responses. Here, this paper reviews the state of the art of these different multiscale approaches and discusses the potential role of personalized or subject-specific modeling in respiratory health.

60 APPLIED LIFE SCIENCES↗

Creation of Self-Semi-Interpenetrating Network Structures in PIM-1 Membranes for Enhanced Physical Aging Resistance

A series of self-semi-interpenetrating network (ssIPN) thin films based on PIM-1 structure were developed by end-cross-linking telechelic PIM-1 oligomers end-capped with curable carboxylic acid groups to form model networks, which are penetrated by linear high-molecular-weight PIM-1 chains. ssIPN films with systematically varied network content ranging from 10 to 30 wt % were comprehensively examined on their microstructure and gas permeation properties. Fresh PIM-1 ssIPN films exhibited gas separation performances close to those of as-cast linear PIM-1, where the films closely followed the upper bound trade-off line, gaining as much selectivity as they lose gas permeability as the network content increases. This indicates that gas permeability is largely preserved despite cross-linking. Wide-angle X-ray scattering supported this observation, with peaks shifting toward lower d-spacing as oligomer content increased, suggesting tighter chain packing with higher cross-linkable oligomer loading. In physical aging studies over two months, all ssIPN films outperform linear PIM-1 for H 2 /CH 4 and O 2 /N 2 separations. Notably, the PIM-1 ssIPN with 10 wt % network content showed the best physical aging resistance, with negligible permeability loss even after one month. Finally, these findings highlight a promising macromolecular strategy for enhancing the physical aging resistance of microporous polymer membranes for gas separation.

PIM-1↗

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗

Summary of UMI Chemical Work

Our initial work for the Urban Modeling Integration (UMI) project was to validate the QUIC model dense gas dispersion of chlorine advection up a ramp. The validation task was to use QUIC to model an experiment done in a wind tunnel. Naturally the QUIC model uses a slumping cylindrical model where inside the cylinder the gas is treated as dense but outside ideal. In trying to validate using the wind tunnel experiment, due to the scale of the wind tunnel experiment and ramp, The model formed a very clear, discrete, concentration barrier at the edge of the cylinder that likely doesn’t line up well with the real dense gas advection/diffusion. This effect is much less noticeable at a larger scale such as an instantaneous release of a 25 ton tank.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrical conductivity of condensed molecular hydrogen in the giant planets

Theoretical interpretation of several phenomena concerning Jupiter and Saturn depends upon the electrical conductivity of molecular hydrogen which, according to present models, forms the outermost layer of both planets. The layer starts at the transition pressure between the metallic and the molecular form of hydrogen, that is around 1 Mbar, and extends to the outside limits of the atmosphere. Whether at the highest pressures (and temperatures) this layer is a solid or a dense fluid is not certain. In any case, the fluid is in supercritical condition so that there is only a gradual transition from dense liquid to a gaseous form. The two theories which require specific values of the conductivity of the condensed molecular hydrogen are those pertaining to the generation of a magnetic field in the liquid hydrogen rather than in the deep metallic interior (HIDE, 1967), and those concerned with the electromagnetic coupling and exchange of angular momentum between the liquid core and the solid molecular hydrogen mantle.

Smoluchowski, R.↗

The atmosphere of Titan

Spectroscopic and radiometric observations are used to form models for the Titan atmospheric composition and circulation.

Hunten, D. M.↗