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,603 records · Page 89

Perturbations, TTVs & the (Un)reliability of Ephemerides of Kepler Planets

Fractional uncertainties quoted for orbital periods of the vast majority of planet candidates listed in recent Kepler catalogs are $<10^{-5}$, with values $\sim 10^{-6}$ (corresponding to 2 minutes per 4 years) being typical. These small uncertainties suggest ephemeris predictions for most Kepler planets are robust for decades to come. Tabulated periods list formal uncertainties in the mean times between midpoints of successive transits during the time interval in which transits were observed. However, transit timing variations (TTVs) produce errors in estimates of some planets’ orbital periods that need to be accounted for in certain dynamical investigations and ephemeris predictions. Periodic TTVs with timescales short compared to the interval of Kepler observations largely average out and do not produce significant errors in estimates of orbital periods. TTVs with timescales comparable to the interval of Kepler observations have been fit for dozens of Kepler planet candidates to estimate long-term average orbital periods (Holczer et al. 2016, ApJS 225, 9), and more detailed dynamical models have been used to estimate long-term average periods of a small number of well-studied planets. Most Kepler planets that show large TTVs are near mean-motion orbital resonances with other planets. The largest effect for planets moderately close to two-body resonances are due to rotation of the forced eccentricity vector by resonant perturbations; the timescale for this precession is typically short, so the variations tend to average out. Libration of planets locked in resonances typically occurs on timescales longer than the Kepler baseline, but few Kepler planets appear to be resonantly-locked. The more general but smaller (during the era of Kepler observations) effect is caused by secular precession of the planets' free eccentricities usually has a period much longer than the baseline of the Kepler observations, so it is not accounted for in estimates of mean periods or uncertainties. This precession causes a discrepancy between Kepler era mean orbital period and long-term mean orbital period can exist even for planets having TTVs that are too small to be observable during the epoch of Kepler observations.

Kepler↗

Nonlinear mesomechanics of composites with periodic microstructure

This work is concerned with modeling the mechanical deformation or constitutive behavior of composites comprised of a periodic microstructure under small displacement conditions at elevated temperature. A mesomechanics approach is adopted which relates the microimechanical behavior of the heterogeneous composite with its in-service macroscopic behavior. Two different methods, one based on a Fourier series approach and the other on a Green's function approach, are used in modeling the micromechanical behavior of the composite material. Although the constitutive formulations are based on a micromechanical approach, it should be stressed that the resulting equations are volume averaged to produce overall effective constitutive relations which relate the bulk, volume averaged, stress increment to the bulk, volume averaged, strain increment. As such, they are macromodels which can be used directly in nonlinear finite element programs such as MARC, ANSYS and ABAQUS or in boundary element programs such as BEST3D. In developing the volume averaged or efective macromodels from the micromechanical models, both approaches will require the evaluation of volume integrals containing the spatially varying strain distributions throughout the composite material. By assuming that the strain distributions are spatially constant within each constituent phase-or within a given subvolume within each constituent phase-of the composite material, the volume integrals can be obtained in closed form. This simplified micromodel can then be volume averaged to obtain an effective macromodel suitable for use in the MARC, ANSYS and ABAQUS nonlinear finite element programs via user constitutive subroutines such as HYPELA and CMUSER. This effective macromodel can be used in a nonlinear finite element structural analysis to obtain the strain-temperature history at those points in the structure where thermomechanical cracking and damage are expected to occur, the so called damage critical points of the structure.

Walker, Kevin P.↗

Computational Component Build-Up for the X-57 Distributed Electric Propulsion Aircraft

A computational study of the wing for the distributed electric propulsion X-57 Maxwell airplane configuration at cruise and takeoff/landing conditions was completed. Three unstructured-mesh, Navier-Stokes computational fluid dynamics methods, FUN3D, USM3D and Kestrel, were used to predict the performance buildup of components to the full X-57 configuration. The goal of the X-57 wing and distributed electric propulsion system design was to meet or exceed the required lift coefficient of 3.95 for a stall speed of 58 knots. The X-57 Maxwell airplane was designed with a small, high aspect ratio cruise wing that was designed for a high cruise lift coefficient of 0.75 at a cruise speed of 150 knots and altitude of 8,000 ft, with an angle of attack of approximately 0deg. The computational data indicates that the X-57 full aircraft drag would meet the cruise drag goal with a 25 count drag margin. The cruise configuration maximum lift coefficient is 2.07 and without including the stabilator is 1.86 at an angle of attack of 14 deg, predicted with the USM3D flow solver using the Spalart-Allmaras turbulence model. The maximum lift coefficient for the high-lift wing (with the 30deg flap deflection) without the stabilator contribution is 2.60 at an angle of attack of 13 deg. For high-lift blowing conditions with 13.7 hp/prop, the maximum lift coefficient excluding the stabilator is 4.426 at (alpha) = 13 deg. Therefore, the lift augmentation from the high-lift propellers is 1.7 and the total lift augmentation from the high-lift system (30 deg flap deflection and the high-lift blowing) is 2.38. The drag for the high-lift wing with 30 deg flap deflection is much higher than the cruise wing configuration, but the high-lift system is used only during a small portion of the flight envelope. The pitching moment is relatively constant for both blown and unblown conditions when the stabilator is excluded. Modeling the full geometry has indicated some adverse effects from the fuselage on the wing and stabilator. At high angles of attack, the solutions with the USM3D flow solver using the Spalart-Allmaras turbulence model indicates large flow separation on the wing upper surface between the two high-lift nacelles near the fuselage, and also a reduction in sectional lift on the stabilator in the first 50 percent of the stabilator semispan. However, the large flow separation near the fuselage is mostly eliminated in the solutions predicted with two codes, USM3D and Kestrel, using Hybrid Reynolds-averaged Navier Stokes/Large Eddy Simulation turbulence models.

Deere, Karen A.↗

Molecular absorption parameters in atmospheric modelling

Molecular spectroscopic parameters are compiled for a number of infrared-active molecules occurring naturally in the terrestrial atmosphere. The following molecules are included in this compilation: water vapor; carbon dioxide; ozone; nitrous oxide; carbon monoxide; methane; and oxygen. The spectral region covered extends from less than 1 micron to the far infrared, and data are presented on more than 100,000 spectral lines. The parameters included in the compilation for each line are: frequency, intensity, half-width, energy of the lower state of the transition, vibrational and rotational identifications of the upper and lower energy states, an isotopic identification, and a molecular identification. Using this data compilation, band model parameters are presented for water vapor, carbon dioxide, and ozone averaged over 20 wavenumber intervals. Using these parameters in a random model formulation, transmittance spectra are provided and compared with both degraded monochromatic calculations and laboratory data.

Mcclatchey, R. A.↗

Evapotranspiration and runoff from large land areas: Land surface hydrology for atmospheric general circulation models

A land surface hydrology parameterization for use in atmospheric GCM's is presented. The parameterization incorporates subgrid scale variability in topography, soils, soil moisture and precipitation. The framework of the model is the statistical distribution of a topography-soils index, which controls the local water balance fluxes, and is therefore taken to represent the large land area. Spatially variable water balance fluxes are integrated with respect to the topography-soils index to yield our large topography-soils distribution, and interval responses are weighted by the probability of occurrence of the interval. Grid square averaged land surface fluxes result. The model functions independently as a macroscale water balance model. Runoff ratio and evapotranspiration efficiency parameterizations are derived and are shown to depend on the spatial variability of the above mentioned properties and processes, as well as the dynamics of land surface-atmosphere interactions.

Famiglietti, J. S.↗

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management Facility (EMF) stack monitor location was qualified using a combination of scale model stacks to mitigate the risk of identifying that the sampling location does not meet the qualification criteria on the full-scale stack. The LV-S1 scale model stack was used as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling locations were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. The LV-S1, LB-S1, and LV-S2 scale model stack tests have met the criteria of the ANSI/HPS N13.1-1999 standard to demonstrate the stack sampling locations are well mixed. Verification tests of the EMF stack were performed at normal operating conditions. The minimum 1/6 DV value and the maximum 6 DV value from the scale model testing determine the range of stack flow rates for which the full-scale stack may be operated while remaining in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rate along with the average and maximum qualified stack flow rate based on the LV-S1 scale model DV values. The operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were <5°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, all of which were within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all planned fan operating configurations. This includes each combination of ACV fans with DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Approach to Improved Credibility of CFD Simulations for Rocket Injector Design

Computational fluid dynamics (CFD) has the potential to improve the historical rocket injector design process by simulating 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 an effort to understand and compare the predictive capabilities and resource requirements of a range of CFD methodologies on a set of model problem injectors. Preliminary results from a steady Reynolds-Average Navier-Stokes (RANS), an unsteady Reynolds-Average Navier Stokes (URANS) and three different Large Eddy Simulation (LES) techniques used to model a single element coaxial injector using gaseous oxygen and gaseous hydrogen propellants are presented. Initial observations are made comparing instantaneous results, corresponding time-averaged and steady-state solutions in the near -injector flow field. Significant differences in the flow fields exist, as expected, and are discussed. An important preliminary result is the identification of a fundamental mixing mechanism, accounted for by URANS and LES, but missing in the steady BANS methodology. Since propellant mixing is the core injector function, this mixing process may prove to have a profound effect on the ability to more correctly simulate injector performance and resulting thermal environments. Issues important to unifying the basis for future comparison such as solution initialization, required run time and grid resolution are addressed.

Tucker, Paul K.↗

The chemistry and mineralogy of Mars soil and dust

A single geological unit consisting of fine, apparently weathered soil material is covering large portions of the surface of Mars. This soil material has been thoroughly homogenized by global dust storms and it is plausible to assume that Mars dust is strongly correlated with it. The chemical-elemental composition of the soil was directly measured by the Viking Landers. Positive detection of Si, Al, Fe, Mg, Ca, Ti, S, Cl, and Br was achieved. Analyses of the SNC meteorites, a group of meteorites that has been suggested to be ejected Martian rocks, supply additional elemental-concentration data, broadening considerably the chemical data-base on the surface materials. A composition model for Mars soil, giving selected average elemental concentrations of major and trace elements, was recently suggested. It was constructed by combining the Viking Lander data, the SNC meteorite analyses, and other analyses. The mineralogy of the surface materials on Mars has not been directly measured yet. By use of various indirect approaches, including chemical correspondence to the surface analyses, spectral analogies, simulations of Viking Lander experiments, analyses of the SNC meteorites and various modeling efforts, the mineralogical composition was constrained to some extent. No direct analyses of soil reactivity have been done yet. Indirect evidence, mostly from the Viking biology experimental results, suggests that the soil probably has a slightly acidic reaction and is generally oxidized. Unambiguous identification of the Mars soil minerals by direct mineralogical analyses, and non-disturbed or in-situ measurements of the soil's reactivity, are of primary importance in future Mars research.

Banin, A.↗

Disturbance Dynamics in Transitional and Turbulent Flows

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. First the decay of laminar disturbances and turbulence in mean shear-free flows is studied. In laminar flows, such disturbances are linear superpositions of modes governed by the Orr-Sommerfeld equation. In turbulent flows, disturbances are described through transport equations for representative mean quantities. The link between a description based on a deterministic evolution equation and a probability based mean transport equation is established. Because an uncertainty in initial conditions exists in the laminar as well as the turbulent regime, a probability distribution must be defined even in the laminar case. Using this probability distribution, it is shown that the exponential decay of the linear modes in the laminar regime can be related to a power law decay of both the (ensemble) mean disturbance kinetic energy and the dissipation rate. The evolution of these mean disturbance quantities is then described by transport equations similar to those for the corresponding turbulent decaying flow. Second, homogeneous shear flow, where disturbances can be described by rapid distortion theory (RDT), is studied. The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined.

Grosch, Chester E.↗

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗

Modeling of turbulence for compression corner flows and internal flows

A compressible Reynolds averaged Navier-Stokes code is applied to examine the performances of basic turbulence models for unseparated and separated flows. The turbulence models considered are a zero-equation model and several two-equation models, including a new k-epsilon model with an eddy viscosity damping function depending upon the Reynolds number and the distance from the wall. The turbulence models are assessed using the following experimental flows: a channel flow, a backward facing step, and a two-dimensional compression corner at Mach 2.8. Discussions and comparisons of the experimental and numerical results are given. Mesh refinement is shown to have a significant beneficial effect on free shear layer flow results obtained with two-equation turbulence models. The results illustrate the usefulness of several turbulence models as a design tool for fluid engineering systems.

Champney, J. M.↗

Solute redistribution in dendritic solidification with diffusion in the solid

An investigation of solute redistribution during dendritic solidification with diffusion in the solid has been performed using numerical techniques. The extent of diffusion is characterized by the instantaneous and average diffusion parameters. These parameters are functions of the diffusion Fourier number, the partition ratio and the fraction solid. Numerical results are presented as an approximate model, which is used to predict the average diffusion parameter and calculate the composition of the interdendritic liquid during solidification.

Ganesan, S.↗

MIRO Computational Model

A computational model calculates the excitation of water rotational levels and emission-line spectra in a cometary coma with applications for the Micro-wave Instrument for Rosetta Orbiter (MIRO). MIRO is a millimeter-submillimeter spectrometer that will be used to study the nature of cometary nuclei, the physical processes of outgassing, and the formation of the head region of a comet (coma). The computational model is a means to interpret the data measured by MIRO. The model is based on the accelerated Monte Carlo method, which performs a random angular, spatial, and frequency sampling of the radiation field to calculate the local average intensity of the field. With the model, the water rotational level populations in the cometary coma and the line profiles for the emission from the water molecules as a function of cometary parameters (such as outgassing rate, gas temperature, and gas and electron density) and observation parameters (such as distance to the comet and beam width) are calculated.

Broderick, Daniel↗

Evaluation of Near-Tropopause Ozone Distributions in the Global Modeling Initiative Combined Stratosphere/Troposphere Model with Ozonesonde Data

The NASA Global Modeling Initiative has developed a combined stratosphere/troposphere chemistry and transport model which fully represents the processes governing atmospheric composition near the tropopause. We evaluate model ozone distributions near the tropopause, using two high vertical resolution monthly mean ozone profile climatologies constructed with ozonesonde data, one by averaging on pressure levels and the other relative to the thermal tropopause. Model ozone is high biased at the SH tropical and NH midlatitude tropopause by approx. 45% in a 4 deg. latitude x 5 deg. longitude model simulation. Increasing the resolution to 2 deg. x 2.5 deg. increases the NH tropopause high bias to approx. 60%, but decreases the tropical tropopause bias to approx. 30%, an effect of a better-resolved residual circulation. The tropopause ozone biases appear not to be due to an overly vigorous residual circulation or excessive stratosphere/troposphere exchange, but are more likely due to insufficient vertical resolution or excessive vertical diffusion near the tropopause. In the upper troposphere and lower stratosphere, model/measurement intercomparisons are strongly affected by the averaging technique. NH and tropical mean model lower stratospheric biases are less than 20%. In the upper troposphere, the 2 deg. x 2.5 deg. simulation exhibits mean high biases of approx. 20% and approx. 35% during April in the tropics and NH midlatitudes, respectively, compared to the pressure averaged climatology. However, relative-to-tropopause averaging produces upper troposphere high biases of approx. 30% and 70% in the tropics and NH midlatitudes. This is because relative-to-tropopause averaging better preserves large cross-tropopause O3 gradients, which are seen in the daily sonde data, but not in daily model profiles. The relative annual cycle of ozone near the tropopause is reproduced very well in the model Northern Hemisphere midlatitudes. In the tropics, the model amplitude of the near tropopause annual cycle is weak. This is likely due to the annual amplitude of mean vertical upwelling near the tropopause, which analysis suggests is approx. 30% weaker than in the real atmosphere.

Considine, David B.↗

Seasonal variations of the atmospheric muon neutrino spectrum measured with IceCube

This study presents an analysis of seasonal variations in the atmospheric muon neutrino flux, using 11.3 years of data from the IceCube Neutrino Observatory. By leveraging a novel spectral unfolding method, we explore the energy range from 125 GeV to 10 TeV for zenith angles from 90° to 110°, corresponding to the Antarctic atmosphere. Our findings reveal that the differential measurement of the amplitudes of the seasonal variation is consistent with an energy-dependent decrease reaching (-4.5 ± 1.2)% during Austral winter and increase to (+ 3.9 ± 1.3)% during Austral summer relative to the annual average at 10 TeV. While the unfolded flux exceeds the model predictions by up to 30%, the differential measurement of the seasonal to annual average flux remains unaffected. The measured seasonal variations of the muon neutrino spectrum are consistent with theoretical predictions using the MCEq code and the NRLMSISE-00 atmospheric model.

Astroparticle Physics and High-Energy Cosmic Pheno↗

Development of message passing-based graph convolutional networks for classifying cancer pathology reports

Abstract Background Applying graph convolutional networks (GCN) to the classification of free-form natural language texts leveraged by graph-of-words features (TextGCN) was studied and confirmed to be an effective means of describing complex natural language texts. However, the text classification models based on the TextGCN possess weaknesses in terms of memory consumption and model dissemination and distribution. In this paper, we present a fast message passing network (FastMPN), implementing a GCN with message passing architecture that provides versatility and flexibility by allowing trainable node embedding and edge weights, helping the GCN model find the better solution. We applied the FastMPN model to the task of clinical information extraction from cancer pathology reports, extracting the following six properties: main site, subsite, laterality, histology, behavior, and grade. Results We evaluated the clinical task performance of the FastMPN models in terms of micro- and macro-averaged F1 scores. A comparison was performed with the multi-task convolutional neural network (MT-CNN) model. Results show that the FastMPN model is equivalent to or better than the MT-CNN. Conclusions Our implementation revealed that our FastMPN model, which is based on the PyTorch platform, can train a large corpus (667,290 training samples) with 202,373 unique words in less than 3 minutes per epoch using one NVIDIA V100 hardware accelerator. Our experiments demonstrated that using this implementation, the clinical task performance scores of information extraction related to tumors from cancer pathology reports were highly competitive.

59 BASIC BIOLOGICAL SCIENCES↗

A theoretical investigation of noise reduction through the cylindrical fuselage of a twin-engine, propeller-driven aircraft

Interior noise in the fuselage of a twin-engine, propeller-driven aircraft with two propellers rotating in opposite directions is studied analytically. The fuselage was modeled as a stiffened cylindrical shell with simply supported ends, and the effects of stringers and frames were averaged over the shell surface. An approximate mathematical model of the propeller noise excitation was formulated which includes some of the propeller noise characteristics such as sweeping pressure waves around the sidewalls due to propeller rotation and the localized nature of the excitation with the highest levels near the propeller plane. Results are presented in the form of noise reduction, which is the difference between the levels of external and interior noise. The influence of propeller noise characteristics on the noise reduction was studied. The results indicate that the sweep velocity of the excitation around the fuselage sidewalls is critical to noise reduction.

Bhat, R. B.↗