Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average 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,351 records · Page 75

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL SCIENCES↗

A Model of Thermal Conductivity for Planetary Soils: 1. Theory for Unconsolidated Soils

We present a model of heat conduction for mono-sized spherical particulate media under stagnant gases based on the kinetic theory of gases, numerical modeling of Fourier s law of heat conduction, theoretical constraints on the gas thermal conductivity at various Knudsen regimes, and laboratory measurements. Incorporating the effect of the temperature allows for the derivation of the pore-filling gas conductivity and bulk thermal conductivity of samples using additional parameters (pressure, gas composition, grain size, and porosity). The radiative and solid-to-solid conductivities are also accounted for. Our thermal model reproduces the well-established bulk thermal conductivity dependency of a sample with the grain size and pressure and also confirms laboratory measurements finding that higher porosities generally lead to lower conductivities. It predicts the existence of the plateau conductivity at high pressure, where the bulk conductivity does not depend on the grain size. The good agreement between the model predictions and published laboratory measurements under a variety of pressures, temperatures, gas compositions, and grain sizes provides additional confidence in our results. On Venus, Earth, and Titan, the pressure and temperature combinations are too high to observe a soil thermal conductivity dependency on the grain size, but each planet has a unique thermal inertia due to their different surface temperatures. On Mars, the temperature and pressure combination is ideal to observe the soil thermal conductivity dependency on the average grain size. Thermal conductivity models that do not take the temperature and the pore-filling gas composition into account may yield significant errors.

Piqueux, S.↗

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗

Physics-Informed Deep Learning for Reconstruction of Spatial Missing Climate Information in the Antarctic

Understanding the influence of the Antarctic on the global climate is crucial for the prediction of global warming. However, due to very few observation sites, it is difficult to reconstruct the rational spatial pattern by filling in the missing values from the limited site observations. To tackle this challenge, regional spatial gap-filling methods, such as Kriging and inverse distance weighted (IDW), are regularly used in geoscience. Nevertheless, the reconstructing credibility of these methods is undesirable when the spatial structure has massive missing pieces. Inspired by image inpainting, we propose a novel deep learning method that demonstrates a good effect by embedding the physics-aware initialization of deep learning methods for rapid learning and capturing the spatial dependence for the high-fidelity imputation of missing areas. We create the benchmark dataset that artificially masks the Antarctic region with ratios of 30%, 50% and 70%. The reconstructing monthly mean surface temperature using the deep learning image inpainting method RFR (Recurrent Feature Reasoning) exhibits an average of 63% and 71% improvement of accuracy over Kriging and IDW under different missing rates. With regard to wind speed, there are still 36% and 50% improvements. In particular, the achieved improvement is even better for the larger missing ratio, such as under the 70% missing rate, where the accuracy of RFR is 68% and 74% higher than Kriging and IDW for temperature and also 38% and 46% higher for wind speed. In addition, the PI-RFR (Physics-Informed Recurrent Feature Reasoning) method we proposed is initialized using the spatial pattern data simulated by the numerical climate model instead of the unified average. Compared with RFR, PI-RFR has an average accuracy improvement of 10% for temperature and 9% for wind speed. When applied to reconstruct the spatial pattern based on the Antarctic site observations, where the missing rate is over 90%, the proposed method exhibits more spatial characteristics than Kriging and IDW.

54 ENVIRONMENTAL SCIENCES↗

Raindrop oscillations

A model of the change in shape of a raindrop is presented. Raindrops measured by two orthogonal cameras were classified by shape and orientation to determine the nature of the oscillation. A physical model based on potential energy was then developed to study the amplitude variation of oscillating drops. The model results show that oscillations occur about the equilibrium axis ratio, but the time average axis ratio if significantly more spherical for large amplitudes because of asymmetry in the surface potential energy. A generalization of the model to oscillations produced by turbulence yields average axis ratios that are consistent with the camera measurements. The model results for average axis ratios were applied to rainfall studies with a dual polarized radar.

Beard, K. V.↗

K and U systematics and average concentrations on the moon

The K-U, Th-U, and Th-Al systematics for lunar samples from the Apollo 11-17 missions and the Luna 16 and 20 missions are summarized. With few exceptions (granitic portion of 12013 and the Apollo 17 mare basalts) the Th/U ratio is 3.8 plus or minus 0.2. The K/U ratio for lunar samples is about 2600 and is distinctly different from the K/U ratio found in terrestrial samples and in the majority of meteorites. The K-U systematics are similar to the K-La, K-Sm, K-Ba, and K-Zr systematics, suggesting that the moon accreted approximately homogeneously. Based on a model the moon contains an average abundance of about 80 ppb U, about 200 ppm K, and about 7 times chondritic abundances for all the refractory elements such as Al, Ca, Ti, Zr, Ba, Sr, and REE. The Th-Al systematics of many lunar highland rocks are dominated by the presence of the rock-type KREEP.

Schonfeld, E.↗

Model-independent constraints on Ωm and H ( z ) from the link between geometry and growth

ABSTRACT We constrain the expansion history of the Universe and the cosmological matter density fraction in a model-independent way by exclusively making use of the relationship between background and perturbations under a minimal set of assumptions. We do so by employing a Gaussian process to model the expansion history of the Universe from present time to the recombination era. The expansion history and the cosmological matter density are then constrained using recent measurements from cosmic chronometers, Type-Ia supernovae, baryon acoustic oscillations, and redshift-space distortion data. Our results show that the evolution in the reconstructed expansion history is compatible with the Planck 2018 prediction at all redshifts. The current data considered in this study can constrain a Gaussian process on H(z) to an average $9.4 {{\ \rm per\ cent}}$ precision across redshift. We find Ωm = 0.224 ± 0.066, lower but statistically compatible with the Planck 2018 cosmology. Finally, the combination of future DESI measurements with the CMB measurement considered in this work holds the promise of $8 {{\ \rm per\ cent}}$ average constraints on a model-independent expansion history as well as a five-fold tighter Ωm constraint using the methodology developed in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine Learning Driven Contouring of High-Frequency Four-Dimensional Cardiac Ultrasound Data

Automatic boundary detection of 4D ultrasound (4DUS) cardiac data is a promising yet challenging application at the intersection of machine learning and medicine. Using recently developed murine 4DUS cardiac imaging data, we demonstrate here a set of three machine learning models that predict left ventricular wall kinematics along both the endo- and epi-cardial boundaries. Each model is fundamentally built on three key features: (1) the projection of raw US data to a lower dimensional subspace, (2) a smoothing spline basis across time, and (3) a strategic parameterization of the left ventricular boundaries. Model 1 is constructed such that boundary predictions are based on individual short-axis images, regardless of their relative position in the ventricle. Model 2 simultaneously incorporates parallel short-axis image data into their predictions. Model 3 builds on the multi-slice approach of model 2, but assists predictions with a single ground-truth position at end-diastole. To assess the performance of each model, Monte Carlo cross validation was used to assess the performance of each model on unseen data. For predicting the radial distance of the endocardium, models 1, 2, and 3 yielded average R2 values of 0.41, 0.49, and 0.71, respectively. Monte Carlo simulations of the endocardial wall showed significantly closer predictions when using model 2 versus model 1 at a rate of 48.67%, and using model 3 versus model 2 at a rate of 83.50%. These finding suggest that a machine learning approach where multi-slice data are simultaneously used as input and predictions are aided by a single user input yields the most robust performance. Subsequently, we explore the how metrics of cardiac kinematics compare between ground-truth contours and predicted boundaries. We observed negligible deviations from ground-truth when using predicted boundaries alone, except in the case of early diastolic strain rate, providing confidence for the use of such machine learning models for rapid and reliable assessments of murine cardiac function. To our knowledge, this is the first application of machine learning to murine left ventricular 4DUS data. Future work will be needed to strengthen both model performance and applicability to different cardiac disease models.

4D ultrasound↗

Validation of High-Fidelity CFD Simulations for Rocket Injector Design

Computational fluid dynamics (CFD) has the potential to improve the historical rocket injector design process by evaluating the sensitivity of performance and injector-driven thermal environments to the details of the injector geometry and key operational parameters. Methodical verification and validation efforts on a range of coaxial injector elements have shown the current production CFD capability must be improved in order to quantitatively impact the injector design process. This paper documents the status of a focused effort to compare and understand the predictive capabilities and computational requirements of a range of CFD methodologies on a set of single element injector model problems. The steady Reynolds-Average Navier-Stokes (RANS), unsteady Reynolds-Average Navier-Stokes (URANS) and three different approaches using the Large Eddy Simulation (LES) technique were used to simulate the initial model problem, a single element coaxial injector using gaseous oxygen and gaseous hydrogen propellants. While one high-fidelity LES result matches the experimental combustion chamber wall heat flux very well, there is no monotonic convergence to the data with increasing computational tool fidelity. Systematic evaluation of key flow field regions such as the flame zone, the head end recirculation zone and the downstream near wall zone has shed significant, though as of yet incomplete, light on the complex, underlying causes for the performance level of each technique. 1 Aerospace Engineer and Combustion CFD Team Leader, MS ER42, NASA MSFC, AL 35812, Senior Member, AIAA. 2 Professor and Director, Computational Combustion Laboratory, School of Aerospace Engineering, 270 Ferst Dr., Atlanta, GA 30332, Associate Fellow, AIAA. 3 Reilly Professor of Engineering, School of Mechanical Engineering, 585 Purdue Mall, West Lafayette, IN 47907, Fellow, AIAA. 4 Principal Member of Technical Staff, Combustion Research Facility, 7011 East Avenue, MS9051, Livermore, CA 94550, Associate Fellow, AIAA. 5 J. L. and G. H. McCain Endowed Chair, Mechanical Engineering, 104 Research Building East, University Park, PA 16802, Fellow, AIAA. American Institute of Aeronautics and Astronautics 1

Tucker, P. Kevin↗

Designing Digital Control Systems With Averaged Measurements

Rational criteria represent improvement over "cut-and-try" approach. Recent development in theory of control systems yields improvements in mathematical modeling and design of digital feedback controllers using time-averaged measurements. By using one of new formulations for systems with time-averaged measurements, designer takes averaging effect into account when modeling plant, eliminating need to iterate design and simulation phases.

Polites, Michael E.↗

Temperature-moisture biases in ECMWF analyses based on clear sky longwave simulations constrained by SSMI and MSU measurements and comparisons to ERBE estimates

Clear sky longwave radiation fluxes for the summer of 1988 and winter of 1989 have been simulated with a radiative transfer model that includes detailed treatment of atmospheric gas absorption. The input data to the model are humidity and temperature profiles from European Center for Medium-Range Weather Forecasts (ECMWF) analyses and surface temperature measurements from the International Satellite Cloud Climatology Project. To reduce the inherent uncertainties in humidity profiles, the ECMWF moisture fields have been adjusted, based on the climatological relationship between the moisture profile and total precipitable water (PW), so that the ECMWF PW is equal to that derived from Special Sensor Microwave Imager (SSMI) data. SSMI and ECMWF PW patterns are generally similar, but significant differences in magnitude are found in the low latitudes, particularly over the subtropical oceans in the Southern Hemisphere. At the same time the tropospheric brightness temperature (T(sub B)) at 53.74 GHz has been computed using a microwave transfer model with ECMWF temperature and humidity fields as inputs. The comparison of computed T(sub B) with the Microwave Sounding Unit (MSU) channel 2 T(sub B) revealed differences larger than 1.5 K over most of the oceans, suggesting that the ECMWF model atmosphere has a cold bias. Maximum biases greater than 3 K are found in low latitudes. After removing the bias from the ECMWF temperature field, a new set of temperature profiles yielding brightness temperature equal to MSU channel 2 T(sub B) has been obtained. Simulation results clearly demonstrate that inclusion of satellite estimates of PW and T(sub B) enhances the accuracy of the clear sky LW flux simulation, substantially reducing the differences from Earth Radiation Budget Experiment (ERBE) estimates. Global averages indicate that the model-derived top-of-atmosphere (TOA) clear sky fluxes from the adjusted moisture and temperature profiles are in very good agreement with satellite-measured ERBE values, differing by only 1.4 W/sq m. By contrast, the calculation using only ECMWF analyses without any adjustments gives a significant underestimation of TOA clear sky fluxes up to 10.7 W/sq m compared to the satellite measurements, suggesting that the constraint method can significantly improve the accuracy in the radiation budget simulation. The results also indicate that the ECMWF-based simulation errors are, at least for the two seasons studied, mainly attributed to the error in the temperature field rather than to errors in the moisture field.

Sohn, Byung-Ju↗

Model-Trained Neural Networks and Electronic Holography Demonstrated to Detect Damage in Blades

Detect Damage in Blades Electronic holography can show damaged regions in fan blades at 30 frames/sec. The electronic holograms are transformed by finite-element-model-trained artificial neural networks to visualize the damage. The trained neural networks are linked with video and graphics to visualize the bending-induced strain distribution, which is very sensitive to damage. By contrast, it is very difficult to detect damage by viewing the raw, speckled, characteristic fringe patterns. For neural-network visualization of damage, 2 frames or 2 fields are used, rather than the 12 frames normally used to compute the displacement distribution from electronic holograms. At the NASA Lewis Research Center, finite element models are used to compute displacement and strain distributions for the vibration modes of undamaged and cracked blades. A model of electronic time-averaged holography is used to transform the displacement distributions into finite-element-resolution characteristic fringe patterns. Then, a feedforward neural network is trained with the fringe-pattern/strain-pattern pairs, and the neural network, electronic holography, and video are implemented on a workstation. Now that the neural networks have been tested successfully at 30 frames/sec on undamaged and cracked cantilevers, the electronic holography and neural-network processing are being adapted for onsite damage inspection of twisted fan blades and rotormounted blades. Our conclusion is that model-trained neural nets are effective when they are trained with good models whose application is well understood. This work supports the aeromechanical testing portion of the Advanced Subsonic Technology Project.

Decker, Arthur J.↗

Annual and longitudinal variations of the Pacific North Equatorial Countercurrent

The climatological annual cycle of the Pacific North Equatorial Countercurrent (NECC) simulated by an ocean general circulation model (OGCM) was studied. The longitudinal variation of transports, degree of geostrophy, and the relationship between Ekman pumping and vertical displacement of the thermocline were emphasized. The longitudinal variation was explored using six sections along 150 deg E, 180 deg, 160 deg W, 140 deg W, 125 deg W, and 110 deg W. A primitive equation OGCM of the Pacific Ocean was run for three years and the fields used were from the third year. The fields consisted of zonal, meridional, and vertical current components and temperature and salinity averaged every three days. The model was forced with the Hellerman and Rosenstein climatological wind stress. The mean annual eastward transport (19.9 Sv) was largest at 160 deg W. The maximum-current boundaries along 160 deg W were 9.2 deg N (1.0 deg), 5.1 deg N (1.1 deg), and 187 m (90.6 m). The annual-cycle amplitude of the NECC was greatest between 160 deg W and 140 deg W. Although the NECC is geostrophic to the first order, deviations from geostrophy were found in the boreal spring and summer near the southern boundary and near the surface. Meridional local acceleration played a role between 3 deg N-5 deg N.

Lolk, Nina K.↗

Improved Regional Moment Tensor Inversion for Moderately Large Earthquakes in the Western United States Using a 3D Earth Model Based on Full Waveform Tomography

The nature of seismic sources for moderately large (moment magnitude, M w 5.0–6.5) events are commonly characterized by their moment tensor (MT) solutions and obtained by inversion of regional distance (200–1600 km) long‐period (20–50 s) waveforms. Regional MT estimates are often calculated from average plane‐layered, one‐dimensional (1D) velocity models. However, 1D model calculations can produce misfits in the arrival times and waveform shapes that introduce errors, particularly at longer distances or for shorter periods, which are necessary for analyzing lower magnitude events. Approximate Earth models (e.g., 1D) representing broad areas may be inadequate, particularly in the crust and uppermost mantle of tectonically complex regions. In this study, we show how a three‐dimensional (3D) Earth model obtained from full waveform inversion tomography can improve waveform fits and decrease phase errors. We developed a platform and workflow to perform routine 3D MT inversions and inverted MTs for 25 earthquakes in the western United States and seven nuclear explosions using an average 1D and a recent 3D Earth model, WUS256 (Rodgers et al., 2022). Using the 3D model improves waveform fits (variance reduction and phase time shifts) compared with the 1D model, and the 3D MT solutions are stable across large distances. This study shows that 3D models obtained from full waveform tomography can improve MTs and source characterization especially at far regional distances (>800 km).

Geosciences↗

A model study of the impact of emission control strategies on Los Angeles air quality

A generalized cell model is developed for the calculation of city-wide averages of photochemical smog components in Los Angeles. This model takes into account the effects of variations with time and within the city of the source strengths, the wind field, and the mixing depth. The effect of the influx of background pollution from outside the modeled volume is also included. Several control strategies for reducing automobile emissions are then introduced into the model, and their impact on predicted pollutant levels, particularly those of O3, are investigated.

Hameed, S.↗

Automated mixed traffic vehicle control and scheduling study

The operation and the expected performance of a proposed automatic guideway transit system which uses low speed automated mixed traffic vehicles (AMTVs) were analyzed. Vehicle scheduling and headway control policies were evaluated with a transit system simulation model. The effect of mixed traffic interference on the average vehicle speed was examined with a vehicle pedestrian interface model. Control parameters regulating vehicle speed were evaluated for safe stopping and passenger comfort. Some preliminary data on the cost and operation of an experimental AMTV system are included. These data were the result of a separate task conducted at JPL, and were included as background information.

Peng, T. K. C.↗

Evolution of the Earth-Moon system

The tidal evolution of the Earth-Moon system is reexamined. Several models of tidal friction are first compared in an averaged Hamiltonian formulation of the dynamics. With one of these models, full integrations of the tidally evolving Earth-Moon system are carried out in the complete, fully interacting, and chaotically evolving planetary system. Classic results on the history of the lunar orbit are confirmed by our more general model. A detailed history of the obliquity of the Earth which takes into account the evolving lunar orbit is presented.

Touma, Jihad↗

Protein solubility modeling

A thermodynamic framework (UNIQUAC model with temperature dependent parameters) is applied to model the salt-induced protein crystallization equilibrium, i.e., protein solubility. The framework introduces a term for the solubility product describing protein transfer between the liquid and solid phase and a term for the solution behavior describing deviation from ideal solution. Protein solubility is modeled as a function of salt concentration and temperature for a four-component system consisting of a protein, pseudo solvent (water and buffer), cation, and anion (salt). Two different systems, lysozyme with sodium chloride and concanavalin A with ammonium sulfate, are investigated. Comparison of the modeled and experimental protein solubility data results in an average root mean square deviation of 5.8%, demonstrating that the model closely follows the experimental behavior. Model calculations and model parameters are reviewed to examine the model and protein crystallization process. Copyright 1999 John Wiley & Sons, Inc.

Proteins/chemistry↗