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,153 records · Page 64

An adaptive knowledge-based data-driven approach for turbulence modeling using ensemble learning technique under complex flow configuration: 3D PWR sub-channel with DNS data

This work describes a new approach to increase the accuracy of Reynolds-averaged Navier–Stokes (RANS) in modeling turbulence flow leveraging the machine learning technique. Traditionally, different turbulence models for Reynolds stress are developed for different flow patterns based on human knowledge. Each turbulence model has a certain application domain and prediction uncertainty. In recent years, with the rapid improvements of machine learning techniques, researchers start to develop an approach to compensate for the prediction discrepancy of traditional turbulence models with statistical models and data. However, the approach has deficiencies in several aspects. For example, the amount of human knowledge introduced to the statistical model couldn’t be controlled, which makes the statistical model learn from a very naïve stage and limits its application. In this work, a new approach is developed to address those deficiencies. Here, the new approach uses the “ensemble learning” technique to control the amount of human knowledge introduced into the statistical model. Therefore, the new approach could be adaptive to the multiple application domains. In conclusion, according to the results of case study, the new approach shows higher accuracy than both traditional turbulence models and the previous machine learning approach.

42 ENGINEERING↗

X-Ray Polarization of the Magnetar 1E 1841−045

We report on IXPE and NuSTAR observations beginning 40 days after the 2024 outburst onset of magnetar 1E 1841−045, marking the first IXPE observation of a magnetar in an enhanced state. Our spectropolarimetric analysis indicates that both a blackbody (BB) plus double power-law (PL) and a double blackbody plus power-law spectral model fit the phase-averaged intensity data well, with a hard PL tail (Γ = 1.19 and 1.35, respectively) dominating above ≈5 keV. For the former model, we find the soft PL (the dominant component at soft energies) exhibits a polarization degree (PD) of ≈30% while the hard PL displays a PD of ≈40%. Similarly, the cool BB of the 2BB+PL model possesses a PD of ≈15% and a hard PL PD of ≈57%. For both models, each component has a polarization angle (PA) compatible with celestial north. Model-independent polarization analysis supports these results, wherein the PD increases from ≈15% to ≈70% in the 2–3 keV and 6–8 keV ranges, respectively, while the PA remains nearly constant. We find marginal evidence for phase-dependent variability of the polarization properties, namely a higher PD at phases coinciding with the hard X-ray pulse peak. We compare the hard X-ray PL to the expectation from resonant inverse Compton scattering (RICS) and secondary pair cascade synchrotron radiation from primary high-energy RICS photons; both present reasonable spectropolarimetric agreement with the data, albeit the latter does so more naturally. We suggest that the soft PL X-ray component may originate from a Comptonized corona in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbulence Modeling Workshop

Advances in turbulence modeling are needed in order to calculate high Reynolds number flows near the onset of separation and beyond. To this end, the participants in this workshop made the following recommendations. (1) A national/international database and standards for turbulence modeling assessment should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be taken of other efforts already under-way, such as that of the European Research Community on Flow, Turbulence, and Combustion (ERCOFTAC) consortium. Carefully selected "unit" experiments will be needed, as well as advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least understood and is a key component of two-equation and higher models. Second priority should be given to the development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes (RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in turbulence modeling development, validation, and implementation should include these approaches as well.

Rubinstein, R.↗

Are existing magnetospheric models excessively stretched?

The distribution of the average north-south magnetic field component Bz in the vicinity of the neutral sheet has been investigated. This component is crucial for mappings between the nightside polar ionosphere and the equatorial magnetosphere. Data sets consisting of about 0.5 R(E) averages of magnetic field observations by the IMP/HEOS and ISEE spacecraft have been compared to the field predicted by the Tsyganenko (1987, 1989) models T87 and T89. It was found that both T87 and T89 underestimate Bz in the near tail region by as much as a factor of 2. Modified versions of the T87 model, incorporating plasma sheet warping, were obtained by fitting the model parameters via nonlinear least squares to the ISEE data set and yielded Bz values in agreement with the ISEE data. The study reveals an enormous scatter among the observed baseline values of Bz (on a time scale of 10-20 min), as well as intrinsic biases imposed by the mathematical structure of tail models, and these two factors (especially the first one) greatly limit the accuracy of model predictions of tail Bz.

Peredo, Mauricio↗

Intermediate Fidelity Solid State Recorder Modeling for NISAR

At 26Tb/day return average, NISAR (National Aero- nautics and Space Administration (NASA) and the In- dian Space Research Organization (ISRO) Synthetic Aperture Radar) will stress ground station facilities with a variable cadence of data generation observing nearly all of earth's land and ice to produce global time- series maps of cryospheric, solid-Earth and ecosystem phenomena. Between the instrument data generation source and downlink sink lies the primary science data bu er, the SSR (Solid State Recorder). The SSR has deferred, discrete le deletion events that occur only after a le is played back to all required ground segments (sometimes both NASA and ISRO). A forward dispatch, advancing frontier repropagation playback scheduler is proposed as a model of the SSR's greedy playback scheduling behavior. Schedules based on the existing low delity model and new intermediate delity model are compared to assess the impact of deferred, discrete deletion on planning. The low delity model is found to produce unfeasible schedules that exceed SSR capacity, causing 3.87% of all observations in a 12 day simulation to be removed by repair actions.

Doubleday, Joshua R.↗

Using Climate Models to Learn: Investigating the Phenomenon of Increasing Surface Air Temperatures Using a Global Climate Modeling Approach

Since 1880, Earth’s average global temperature has increased by approximately 1.4 °F (0.8 °C). Approximately two-thirds of the warming has occurred since 1975 (Wuebbles, Easterling, Hayhoe, et al. 2017). Despite consensus about global climate change (GCC) within the scientific community, students struggle to comprehend - why temperatures are increasing globally - at what rate this change is occurring, and - how climate scientists predict Earth’s future climate. These challenges are due, in part, to the abstraction inherent to understanding GCC, a phenomenon that occurs over long periods of time (past and future), and over large regions (regional and global). With the adoption of the Next Generation Science Standards (NGSS Lead States 2013), GCC is emphasized at the secondary school level, with a focus on analyzing and understanding the evidence for GCC (HS-ESS3-1), climate modeling (HS-ESS2-6, HS-ESS3-5), and the impact of humans on the environment (HS-LS2-7, HS-ESS3-6).

climate models↗

A Tropospheric Emission Spectrometer HDO/H2O Retrieval Simulator for Climate Models

Retrievals of the isotopic composition of water vapor from the Aura Tropospheric Emission Spectrometer (TES) have unique value in constraining moist processes in climate models. Accurate comparison between simulated and retrieved values requires that model profiles that would be poorly retrieved are excluded, and that an instrument operator be applied to the remaining profiles. Typically, this is done by sampling model output at satellite measurement points and using the quality flags and averaging kernels from individual retrievals at specific places and times. This approach is not reliable when the model meteorological conditions influencing retrieval sensitivity are different from those observed by the instrument at short time scales, which will be the case for free-running climate simulations. In this study, we describe an alternative, categorical approach to applying the instrument operator, implemented within the NASA GISS ModelE general circulation model. Retrieval quality and averaging kernel structure are predicted empirically from model conditions, rather than obtained from collocated satellite observations. This approach can be used for arbitrary model configurations, and requires no agreement between satellite-retrieved and model meteorology at short time scales. To test this approach, nudged simUlations were conducted using both the retrieval-based and categorical operators. Cloud cover, surface temperature and free-tropospheric moisture content were the most important predictors of retrieval quality and averaging kernel structure. There was good agreement between the D fields after applying the retrieval-based and more detailed categorical operators, with increases of up to 30 over the ocean and decreases of up to 40 over land relative to the raw model fields. The categorical operator performed better over the ocean than over land, and requires further refinement for use outside of the tropics. After applying the TES operator, ModelE had D biases of 8 over ocean and 34 over land compared to TES D, which were less than the biases using raw model D fields.

Field, R. D.↗

Microbial seasonality promotes soil respiratory carbon emission in natural ecosystems: A modeling study

Seasonality is a key feature of the biosphere and the seasonal dynamics of soil carbon (C) emissions represent a fundamental mechanism regulating the terrestrial–climate interaction. We applied a microbial explicit model—CLM-Microbe—to evaluate the impacts of microbial seasonality on soil C cycling in terrestrial ecosystems. Additionally, the CLM-Microbe model was validated in simulating belowground respiratory fluxes, that is, microbial respiration, root respiration, and soil respiration at the site level. On average, the CLM-Microbe model explained 72% (n = 19, p < 0.0001), 65% (n = 19, p < 0.0001), and 71% (n = 18, p < 0.0001) of the variation in microbial respiration, root respiration, and soil respiration, respectively. We then compared the model simulations of soil respiratory fluxes and soil organic C content in top 1 m between the CLM-Microbe model with (CLM-Microbe) and without (CLM-Microbe_wos) seasonal dynamics of soil microbial biomass in natural biomes. Removing soil microbial seasonality reduced model performance in simulating microbial respiration and soil respiration, but led to slight differences in simulating root respiration. Compared with the CLM-Microbe, the CLM-Microbe_wos underestimated the annual flux of microbial respiration by 0.6%–32% and annual flux of soil respiration by 0.4%–29% in natural biomes. Correspondingly, the CLM-Microbe_wos estimated higher soil organic C content in top 1 m (0.2%–7%) except for the sites in Arctic and boreal regions. Our findings suggest that soil microbial seasonality enhances soil respiratory C emissions, leading to a decline in SOC storage. An explicit representation of soil microbial seasonality represents a critical improvement for projecting soil C decomposition and reducing the uncertainties in global C cycle projection under the changing climate.

59 BASIC BIOLOGICAL SCIENCES↗

A rapid and accurate method for calculation of stratospheric photolysis rates with molecular scattering

A method for calculating the photodissociation rates needed for photochemical modeling of the stratosphere, which includes the effects of molecular scattering, is described. The procedure is based on Sokolov's method of averaging functional correction. The radiation model and approximations used to calculate the radiation field are examined. The approximated diffuse fields and photolysis rates are compared with exact data. It is observed that the approximate solutions differ from the exact result by 10 percent or less at altitudes above 15 km; the photolysis rates differ from the exact rates by less than 5 percent for altitudes above 10 km and all zenith angles, and by less than 1 percent for altitudes above 15 km.

Boughner, Robert E.↗

X-59 Sonic Boom Test Results from the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel

A wind tunnel test was conducted to investigate near-field sonic boom pressure signatures from a model of the X-59 Low-Boom Flight Demonstrator aircraft. A 1.62%-scale model of the aircraft in the C612A configuration was fabricated for the wind tunnel test, which took place in the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel in September and October 2021. The model had provisions for two different mounting options: a swept blade strut that attached at the top of model ahead of the inlet, and rear-entry sting that was made as one piece with a dummy nacelle, and which had a 2”-long cylindrical segment aft of the nozzle exit before tapering up in size. The blade strut allowed for a clean aft end of the model for evaluation of the shocks from that region, while the sting avoided the significant distortions of the flow and shocks from the blade strut along the top of the model. Both the sting and the strut had adapters that attached to a force balance. The model had alternate parts for ±0.5° deflections of the flaps, ailerons, and stabilator, and ±1° deflections of the T-tail horizontal surface. Off-body static pressure measurements of the flow field below the model were made by use of a pressure rail which had 420 orifices along its tip. The model was positioned at various heights from the rail by vertical movement of the wind tunnel strut, and at various longitudinal stations relative to the rail by means of a linear actuator mounted between the tunnel strut and the balance. The longitudinal positioning allowed multiple pressure signatures to be obtained along different portions of the rail. These signatures were aligned by accounting for the model longitudinal movement and then averaged to take out the effects of tunnel flow distortions and the interference of the rail flow field and shocks on the model pressure signatures. The test was run at approximate Mach numbers of 1.36, 1.4, and 1.47, and the model was set at various angles of attack and roll relative to the rail. Plots of the model signatures for all the variations of Mach number, model angles, control deflections, and height relative to the rail are provided throughout the report. Repeatability was generally very good and gave confidence in the quality of the measurements. The signatures measured at various heights from the rail provided insight into the aging of the model shocks as they propagated from 1.2 to 3 body lengths from the model. Off-track signatures up to 45° from centerline obtained by rolling the model gave indications of the shock flow fields across the width of the sonic boom carpet. The deflections of the various control surfaces allowed assessment of the boom sensitivity to the control surface movements.

Sonic boom↗

Clouds and Convective Self-Aggregation in a Multi-Model Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective56equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climatecsensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique amongst intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than un-aggregated simulations. Plain Language Summary This study investigates tropical clouds and climate using results from more than 30 different numerical models set up in a simplified framework. The dataset of model simulations is unique in that it includes a wide range of model types configured in a consistent manner. We address some of the biggest open questions in climate science, including how cloud properties change with warming and the role that the tendency of clouds to form clusters plays in determining the average climate and how climate changes. While there are large differences in how the different models simulate average temperature, humidity, and cloudiness, in a majority of models, the amount of high clouds decreases as climate warms. Nearly all models simulate a tendency for clouds to cluster together. There is agreement that when the clouds are clustered, the atmosphere is drier with fewer clouds overall. We don’t find a conclusive result for how cloud clustering changes as the climate warms.

Allison A. Wing↗

Turbomachinery Flows Modeled

Last year, researchers at the NASA Lewis Research Center used the average passage code APNASA to complete the largest three-dimensional simulation of a multistage axial flow compressor to date. Consisting of 29 blade rows, the configuration is typical of those found in aeroengines today. The simulation, which was executed on the High Performance Computing and Communications (HPCC) Program IBM SP2 parallel computer located at the NASA Ames Research Center, took nearly 90 hr to complete. Since the completion of this activity, a fine-grain, parallel version of APNASA has been written by a team of researchers from General Electric, NASA Lewis, and NYMA. Timing studies performed on the SP2 have shown that, with eight processors assigned to each blade row, the simulation time is reduced by a factor of six. For this configuration, the simulation time would be 15 hr. The reduction in computing time indicates that an overnight turnaround of a multistage configuration simulation is feasible. In addition, average passage forms of two-equation turbulence models were formulated. These models are currently being incorporated into APNASA.

Adamczyk, John J.↗

Analysis of Approximations in Modeling of BWR Bundle Void Distributions

In boiling water reactors, complex heterogeneous bundle designs, control blades adjacent to the corner of bundles, and the presence of boiling can lead to complex internal void distributions. A few approximations exist to model these void distributions. They could be modeled using a 1D axial solver in which each axial node is assumed to be at an average void, or each pin cell could be modeled with its own void concentration. In the latter case, the void could be discretized in pin-centered or coolant-centered channels. The goal of this project was to quantify the effect of using the different approximations for modeling internal void distributions on neutronics calculations. Using 3D void distributions calculated with CTF, Monte Carlo Neutral Particle (MCNP) transport code models were created for GE-9 and GE-14 lattices. For each model, the internal void distribution from CTF at a given axial node was selected, and a lattice calculation was carried out with MCNP. Comparisons between models using a lattice-averaged void, or using a void distribution in coolant-centered channels, showed large differences in reactivity which in some cases were well above 1,000 pcm, and it also showed differences in normalized fission rates greater than 20%. It was also found that using a lattice average void can lead to a significant difference in the estimation of the worth of a control blade. The differences found when comparing results from models using pin-centered and coolant-centered channels were up to 200 pcm in reactivity and up to 1.4% in the normalized fission rates. In addition to these two sets of comparisons, MCNP models were set up so that each subchannel had a saturated liquid component around the fuel pins and a saturated vapor component in the center to approximate annular flow. In comparison to the models using coolant-centered subchannels, up to 1–3% differences in normalized fission rates could be found.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Martian general circulation experiment with large topography

A three-layer general circulation model of the Martian atmosphere is described, and the assumptions governing the model are discussed. The simulated, zonally averaged circulation is found to have only limited sensitivity to differences between this model and an earlier general circulation model; this circulation compares reasonably well with observations. It is also found that the meridional mass flow produced by the seasonal condensation of CO2 in the winter polar region has a major influence on the circulation; owing to the weak influence of atmospheric heat transport, however, the mass flow is governed almost entirely by radiation. Quasi-barotropic stationary waves, which are forced kinematically by the topography and which resemble topographically forced terrestrial planetary waves, are generated by the model in the winter hemisphere region of strong eastward flow, while baroclinic stationary waves are thermally forced by topography in the tropics and summer subtropics. It is also concluded that transient baroclinically unstable waves, of somewhat lower dominant wavenumber than those found on the earth, are generated in winter midlatitudes; their amplitudes, wavenumbers, and phase speeds closely agree with what has been deduced from the Viking lander observations.

Pollack, J. B.↗

On representation of fracture profiles by fractional integrals of a Wiener process

An enhancement of the crack diffusion model of brittle fracture is proposed. In the model, certain probabilities are evaluated as averages over the set of all virtual crack paths, and the latter are modeled as realization of a Wiener process. However, commonly observed fracture profiles are smoother. In fact, one can show that their fractal dimension is smaller than the value 3/2 for a Wiener process. Accordingly, it is proposed to represent fracture profiles by graphs of fractional integrals of a Wiener process. This allows one to employ additional fractographic information obtained from fracture profiles, e.g. their fractal dimension, in a predictive formalism, evaluating probabilities of failure.

Kunin, B.↗

Nonlinear System Identification for Aeroelastic Systems with Application to Experimental Data

Representation and identification of a nonlinear aeroelastic pitch-plunge system as a model of the Nonlinear AutoRegressive, Moving Average eXogenous (NARMAX) class is considered. A nonlinear difference equation describing this aircraft model is derived theoretically and shown to be of the NARMAX form. Identification methods for NARMAX models are applied to aeroelastic dynamics and its properties demonstrated via continuous-time simulations of experimental conditions. Simulation results show that (1) the outputs of the NARMAX model closely match those generated using continuous-time methods, and (2) NARMAX identification methods applied to aeroelastic dynamics provide accurate discrete-time parameter estimates. Application of NARMAX identification to experimental pitch-plunge dynamics data gives a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

Electron Flux Models for Different Energies at Geostationary Orbit

Forecast models were derived for energetic electrons at all energy ranges sampled by the third-generation Geostationary Operational Environmental Satellites (GOES). These models were based on Multi-Input Single-Output Nonlinear Autoregressive Moving Average with Exogenous inputs methodologies. The model inputs include the solar wind velocity, density and pressure, the fraction of time that the interplanetary magnetic field (IMF) was southward, the IMF contribution of a solar wind-magnetosphere coupling function proposed by Boynton et al. (2011b), and the Dst index. As such, this study has deduced five new 1 h resolution models for the low-energy electrons measured by GOES (30-50 keV, 50-100 keV, 100-200 keV, 200-350 keV, and 350-600 keV) and extended the existing >800 keV and >2 MeV Geostationary Earth Orbit electron fluxes models to forecast at a 1 h resolution. All of these models were shown to provide accurate forecasts, with prediction efficiencies ranging between 66.9% and 82.3%.

Boynton, R. J.↗

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗