Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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 919 records · Page 51

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗

The fluid dynamic and shear environment in the NASA/JSC rotating-wall perfused-vessel bioreactor

The rotating-wall perfused-vessel (RWPV) bioreactor, used for both microgravity and Earth-based cell science experiments, is characterized in terms of the fluid dynamic and fluid shear stress environment. A numerical model of the flow field is developed and verified with laser Doppler velocimeter measurements. The effects of changes in operating conditions, including rotation rates and fluid perfusion rates, are investigated with the numerical model. The operating conditions typically used for ground-based experiments (equal rotation of the inner and outer cylinders) leads to flow patterns with relatively poor mass distribution characteristics. Approximately 50% of the inlet-perfused fluid bypasses the bulk of the fluid volume and flows to the perfusion exit. For operating conditions typical in microgravity, small differential rotation rates between the inner and outer cylinders lead to greatly improved flow distribution patterns and very low fluid shear stress levels over a large percentage of the fluid volume. Differences in flow patterns for the different operating conditions are explored. Large differences in the hydrodynamic environments for operating conditions typical of true microgravity and ground-based "microgravity simulations" are demonstrated.

Bioreactors↗

Computer Modeling of Non-Isothermal Crystallization

A realistic computer model for simulating isothermal and non-isothermal phase transformations proceeding by homogeneous and heterogeneous nucleation and interface-limited growth is presented. A new treatment for particle size effects on the crystallization kinetics is developed and is incorporated into the numerical model. Time-dependent nucleation rates, size-dependent growth rates, and surface crystallization are also included. Model predictions are compared with experimental measurements of DSC/DTA peak parameters for the crystallization of lithium disilicate glass as a function of particle size, Pt doping levels, and water content. The quantitative agreement that is demonstrated indicates that the numerical model can be used to extract key kinetic data from easily obtained calorimetric data. The model can also be used to probe nucleation and growth behavior in regimes that are otherwise inaccessible. Based on a fit to data, an earlier prediction that the time-dependent nucleation rate in a DSC/DTA scan can rise above the steady-state value at a temperature higher than the peak in the steady-state rate is demonstrated.

Kelton, K. F.↗

Modelling the pole tide and its effect on the earth's rotation

The pole tide is the response of the ocean to incremental centrifugal forces associated with the Chandler wobble. The tide has a potentially important effect on the period and damping of the wobble, but it is at present not well constrained by observations. Here, both analytical and numerical models for the pole tide are constructed. The analytical models consider the tide first in a global ocean and then in an enclosed basin on a beta-plane. The results are found to approach equilibrium linearly with decreasing frequency and inversely with increasing basin depth. The numerical models solve Laplace's tidal equations over the world's oceans using realistic continental boundaries and bottom topography. The results indicate that the effects of non-equilibrium portion of the deep ocean tide on the Chandler wobble period and damping are negligible.

Carton, J. A.↗

Revision of Faraday rotation measure constraints on the primordial magnetic field using the IllustrisTNG simulation

ABSTRACT Previously derived Faraday rotation constraints on the volume-filling intergalactic magnetic field (IGMF) have used analytical models that made a range of simplifying assumptions about magnetic field evolution in the intergalactic medium and did not consider the effect of baryonic feedback on large-scale structures. In this work, we revise existing Faraday rotation constraints on the IGMF using a numerical model of the intergalactic medium from the IllustrisTNG cosmological simulation that includes a sophisticated model of the baryonic feedback. We use the IllustrisTNG model to calculate the rotation measure and compare the resulting mean and median of the absolute value of the rotation measure with data from the NRAO VLA Sky Survey (NVSS). The numerical model of the intergalactic medium includes a full magnetohydrodynamic model of the compressed primordial magnetic field as well as a model of the regions where the magnetic field is not primordial, but is rather produced by the process of baryonic feedback. Separating these two types of regions, we are able to assess the influence of the primordial magnetic field on the Faraday rotation signal. We find that by correcting for regions of compressed primordial field and accounting for the fact that part of the intergalactic medium is occupied by magnetic fields spread by baryonic feedback processes rather than by the primordial field relaxes the Faraday rotation bound by a factor of ≃3. This results in B0 < 1.8 × 10−9 G for large correlation length IGMFs.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗

Rotorcraft Transmission Noise Path Model, Including Distributed Fluid Film Bearing Impedance Modeling

A computational approach for simulating the effects of rolling element and journal bearings on the vibration and sound transmission through gearboxes has been demonstrated. The approach, using ARL/Penn State s CHAMP methodology, uses Component Mode Synthesis of housing and shafting modes computed using Finite Element (FE) models to allow for rapid adjustment of bearing impedances in gearbox models. The approach has been demonstrated on NASA GRC s test gearbox with three different bearing configurations: in the first condition, traditional rolling element (ball and roller) bearings were installed, and in the second and third conditions, the traditional bearings were replaced with journal and wave bearings (wave bearings are journal bearings with a multi-lobed wave pattern on the bearing surface). A methodology for computing the stiffnesses and damping in journal and wave bearings has been presented, and demonstrated for the journal and wave bearings used in the NASA GRC test gearbox. The FE model of the gearbox, along with the rolling element bearing coupling impedances, was analyzed to compute dynamic transfer functions between forces applied to the meshing gears and accelerations on the gearbox housing, including several locations near the bearings. A Boundary Element (BE) acoustic model was used to compute the sound radiated by the gearbox. Measurements of the Gear Mesh Frequency (GMF) tones were made by NASA GRC at several operational speeds for the rolling element and journal bearing gearbox configurations. Both the measurements and the CHAMP numerical model indicate that the journal bearings reduce vibration and noise for the second harmonic of the gear meshing tones, but show no clear benefit to using journal bearings to reduce the amplitudes of the fundamental gear meshing tones. Also, the numerical model shows that the gearbox vibrations and radiated sound are similar for journal and wave bearing configurations.

Hambric, Stephen A.↗

The Future of Land-Based Wind Turbine Rotor Technology: The Perspective from NREL

The Big Adaptive Rotor (BAR) Project looks at potential innovative pathways for future land-based wind turbine technology. BAR is led by the National Renewable Energy Laboratory and is sponsored by the Wind Energy Technology Office of the US Department of Energy. This talk will describe innovations and challenges faced by the teams working on BAR. The innovations include numerical models predicting the stability and performance of highly flexible blades, downwind rotors to maximize wind farm power, and controlled bending of 100-meter-long blades during rail transport. The challenges of developing predictive numerical models and formulating successful value propositions supporting new technologies such as distributed aero control devices will also be discussed.

BAR↗

Integrated Model of the Eye/Optic Nerve Head Biomechanical Environment

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a concern for long-duration space flight. Previously, it has been suggested that ocular changes observed in VIIP syndrome are related to the cephalad fluid shift that results in altered fluid pressures [1]. We are investigating the impact of changes in intracranial pressure (ICP) using a combination of numerical models, which simulate the effects of various environment conditions, including finite element (FE) models of the posterior eye. The specific interest is to understand how altered pressures due to gravitational changes affect the biomechanical environment of tissues of the posterior eye and optic nerve sheath. METHODS: Additional description of the numerical modeling is provided in the IWS abstract by Nelson et al. In brief, to simulate the effects of a cephalad fluid shift on the cardiovascular and ocular systems, we utilized a lumped-parameter compartment model of these systems. The outputs of this lumped-parameter model then inform boundary conditions (pressures) for a finite element model of the optic nerve head (Figure 1). As an example, we show here a simulation of postural change from supine to 15 degree head-down tilt (HDT), with primary outcomes being the predicted change in strains at the optic nerve head (ONH) region, specifically in the lamina cribrosa (LC), retrolaminar optic nerve, and prelaminar neural tissue (PLNT). The strain field can be decomposed into three orthogonal components, denoted as the first, second and third principal strains. We compare the peak tensile (first principal) and compressive (third principal) strains, since elevated strain alters cell phenotype and induces tissue remodeling. RESULTS AND CONCLUSIONS: Our lumped-parameter model predicted an IOP increase of c. 7 mmHg after 21 minutes of 15 degree HDT, which agreed with previous reports of IOP in HDT [1]. The corresponding FEM simulations predicted a relative increase in the magnitudes of the peak tensile and compressive strains in the lamina cribrosa of 42 and 43, respectively (Fig. 2). The corresponding changes in the optic nerve strains were 17 and 39, while in the PLNT they were 47 and 43. These magnitudes of relative elevations in peak strains may induce a phenotypic response in resident mechano-responsive resident cells [2]. This approach may be expanded to investigate other environmental changes (e.g. parabolic flight). Through our VIIP SCHOLAR project, we will validate and improve these integrated models by measuring patient-specific changes in optic nerve sheath geometry in patients with idiopathic intracranial hypertension before and after lumbar puncture and CSF removal.

biomechanics↗

Towards Resolving the Risk of Turbine Collision on Fish: FY21-22 Seedling Report Phase I

The most important concern for permitting tidal and river turbines is the collision risk of marine animals with the turbine blades. Our understanding of the risk to individual fish from colliding with turbine blades is poor; if these collisions were to occur, it is unknown whether fish will sustain recoverable injuries or be killed. Equally unknown is the impact these collisions might have on populations, particularly for threatened, endangered, or commercially managed fish species. In addition to observations of interactions of fish with turbines, numerical models need to be developed or expanded to predict impacts on fish populations. These models could replace expensive and technically challenging monitoring programs in high energy, often murky, tidal or river waters. In this project, a multi-pronged approach was taken to understand the state of knowledge for collision risk of fish with turbines and build a pathway forward. This approach included a workshop with experts, a literature review of modeling and empirical studies, the development of a research framework, and the identification of case studies to address through future work. This project is a first step towards the future development of effective and robust numerical models for assessing the collision risk of fish around turbines. This will achieve a balanced and improved estimate of the severity of collision risk to fish at the population scale.

Garavelli, Lysel↗

Solute transport and the prediction of breakaway oxidation in gamma + beta Ni-Cr-Al alloys

The Al transport and the condition leading to breakaway oxidation during the cyclic oxidation of gamma + beta NiCrAl alloys have been studied. The Al concentration/distance profiles were measured after various cyclic oxidation exposures at 1200 C. It was observed that cyclic oxidation results in a decreasing Al concentration at the oxide/metal interface, maintaining a constant flux of Al to the Al2O3 scale. It was also observed that breakaway oxidation occurs when the Al concentration at the oxide/metal interface approaches zero. A numerical model was developed to simulate the diffusional transport of Al and to predict breakaway oxidation in gamma + beta NiCrAl alloys undergoing cyclic oxidation. In a comparison of two alloys with similar oxide spalling characteristics, the numerical model was shown to predict correctly the onset of breakaway oxidation in the higher Al-content alloy.

Nesbitt, J. A.↗

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) eV and 1 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass less than about 10(exp -2) eV, axions produced deep inside the neutron star simply stream out; in a previous paper this case has been addressed. Remarkably, for an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become 'trapped' and radiated from an axion sphere. In this paper the trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a leakage approximation scheme for axion-energy transport. The axion opacity is computed due to inverse nucleon-nucleon, axion bremsstrahlung, and numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande 2 (K2) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The earlier estimate is confirmed and refined of the axion mass above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration--the most sensitive barometer of axion cooling--it is concluded that for an axion mass of greater than about 0.3 eV, axion emission would not have had a significant effect on the neutrino bursts detected by K2 and IMB. The present work, together with the previous work, strongly suggests that an axion with mass in the interval 10(exp -3) eV to 0.3 eV is excluded by SN 1987A.

Burrows, Adam↗

On-Orbit Propellant Motion Resulting from an Impulsive Acceleration

In-space docking and separation maneuvers of spacecraft that have large fluid mass fractions may cause undesirable spacecraft motion in response to the impulsive-acceleration-induced fluid motion. An example of this potential low-gravity fluid management problem arose during the development of the shuttle/Centaur vehicle. Experimentally verified numerical modeling techniques were developed to establish the propellant dynamics, and subsequent vehicle motion, associated with the separation of the Centaur vehicle from the shuttle cargo bay. Although the shuttle/Centaur development activity has been suspended, the numerical modeling techniques are available to predict on-orbit liquid motion resulting from impulsive accelerations for other missions and spacecraft.

Hochstein, John I.↗

Hexagonal Distributed Embedded Energy Converters (HexDEECs)

The HexDEEC is a small, characteristic length approximating a centimeter, energy transducer that converts the dynamic deformations of its elastomer housing into electricity through a variable capacitance charging-discharging cycle. This device is a type of Distributed Embedded Energy Converter Technology (DEEC-Tec), a new domain for marine renewable energy research that utilizes a conglomeration of small distributed embedded energy converters (DEECs) that, in aggregate, form larger metamaterial frameworks. These resulting DEEC-Tec metamaterials can then, in turn, be used to construct flexible ocean wave energy converters called flexWECs, which can utilize a broad band of ocean wave frequencies and lack highly loaded rigid bodies. These systems also provide new avenues of wave energy harvesting such as actively transforming topologies (e.g., shape and form) and morphologies (e.g., stiffness and damping throughout its entire structure) in real time. Presented, is one specific type of DEEC: the HexDEEC, which is currently being developed by the United States National Renewable Energy Laboratory. This transducer shows promise in aiding the adoption and further development of the DEEC-Tec domain. The following presentation focuses on the promise of this technology and current work being done to analyze the performance of an individual HexDEEC design. The HexDEEC is composed of a hyperelastic hexagonal housing, nominally silicon rubber, with six electrodes on its inner faces. The upper three electrodes share the same charge while the lower three electrodes oppose the upper electrode charges. Externally, the HexDEEC has two arms extending away from the middle vertices of the hexagon. Via principles governing the relationship between electrical capacitance and electrical potential (voltage and charge), electricity is generated when the HexDEEC's arms are dynamically pulled or released under tensile loading, as doing so causes the distance between the upper and lower sets of electrodes to change - varying the energy converter's overall capacitance. Analytical and numerical modeling is being used to evaluate the mechanics and electrical energy generated by the HexDEEC. Equations to describe the capacitance and electrostatic forces acting on this unique system have been developed and implemented into the numerical modeling software STAR-CCM+, along with models to describe its hyperelastic material, such as the Mooney-Rivlin 3-parameter model. So far, an initial design has been analyzed and we plan to further optimize it to increase power production. Individual HexDEECs have been fabricated by drawing uncured liquid silicon rubber into molds via vacuum pressure. To simplify manufacturing, HexDEEC sub-components - e.g., electrodes, wires - can be placed within those molds such that they are directly embedded into the hexagonal housing during the curing process. Furthermore, DEEC-Tec metamaterials can be created by interweaving or sequentially layering multiple HexDEEC strands together. The HexDEEC based metamaterial could then generate electricity through its gross deformations. Ultimately, HexDEECs represent a specific type of energy transducer that can be leveraged, by the DEEC-Tec domain, to create metamaterials used to construct novel flexWECs.

DEEC-Tec↗

The NASA Ames 16-Inch Shock Tunnel Nozzle Simulations and Experimental Comparison

The 16-Inch Shock Tunnel at NASA Ames Research Center is a unique test facility used for hypersonic propulsion testing. To provide information necessary to understand the hypersonic testing of the combustor model, computational simulations of the facility nozzle were performed and results are compared with available experimental data, namely static pressure along the nozzle walls and pitot pressure at the exit of the nozzle section. Both quasi-one-dimensional and axisymmetric approaches were used to study the numerous modeling issues involved. The facility nozzle flow was examined for three hypersonic test conditions, and the computational results are presented in detail. The effects of variations in reservoir conditions, boundary layer growth, and parameters of numerical modeling are explored.

TokarcikPolsky, S.↗

Methods for computing comet core temperatures

The temperature profile within the comet nucleus provides the key to an understanding of the history of the volatiles within a comet. Certain difficulties arise in connection with current cometary temperature models. It is shown that the constraint of zero net heat flow can be used to derive general analytical expressions which will allow for the determination of comet core temperature for a spherically symmetric comet, taking into account information about the surface temperature and the thermal conductivity. The obtained results are compared with the expression for comet core temperatures considered by Klinger (1981). Attention is given to analytical results, an example case, and numerical models. The formalization developed makes it possible to determine the core temperature on the basis of the numerical models of the surface temperature.

Mckay, C. P.↗

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) and 10 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become trapped and radiated from an axion sphere. The trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a diffusion approximation for axion-energy transport. The axion opacity due to inverse nucleon-nucleon, axion bremsstrahlung is computed; and then the numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande II (KII) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The estimate of the axion mass is confirmed above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration - the most sensitive barometer of axion cooling - it is concluded that for an axion mass greater than about 3 eV axion emission would not have had a significant effect on the neutrino bursts detected by KII and IMB. It is strongly suggested that an axion with mass in the interval 10(exp -3) to 3 eV is excluded by the observation of neutrinos from SN 1987A.

Burrows, Adam↗