Search NASA⌕ Search

SEARCH · Search NASA

Results for “Biot 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Microwave heating and joining of ceramic cylinders: A mathematical model

A thin cylindrical ceramic sample is placed in a single mode microwave applicator in such a way that the electric field strength is allowed to vary along its axis. The sample can either be a single rod or two rods butted together. We present a simple mathematical model which describes the microwave heating process. It is built on the assumption that the Biot number of the material is small, and that the electric field is known and uniform throughout the cylinder's cross-section. The model takes the form of a nonlinear parabolic equation of reaction-diffusion type, with a spatially varying reaction term that corresponds to the spatial variation of the electromagnetic field strength in the waveguide. The equation is analyzed and a solution is found which develops a hot spot near the center of the cylindrical sample and which then propagates outwards until it stabilizes. The propagation and stabilization phenomenon concentrates the microwave energy in a localized region about the center where elevated temperatures may be desirable.

Booty, Michael R.↗
DOE OSTI2026

Integration of the Biot–Gassmann Fluid Substitution Method and Machine Learning-Based Velocity–Stress Relationship for Estimating In Situ Stresses

Recent advancements have shown that in situ stresses can be reliably estimated through an integrated machine/deep learning (ML/DL)-based framework, which relies on models trained and validated using true triaxial ultrasonic velocity (TUV) experimental data that involve measurements of ultrasonic velocity in saturated rocks under varying stress configurations. However, when the goal is to interpret lower frequency measurements, it may be more appropriate to run experiments on dry rocks and then obtain Biot–Gassmann-derived equivalent saturated velocities (low-frequency approximation) and employ these quantities for training ML/DL models to predict in situ stress. Whether the dispersion effect of frequency on the velocity–stress relationship substantially impacts in situ stress prediction is an important and unresolved question. This work presents an enhancement of ML/DL-based workflow by training and implementing ML/DL models using equivalent saturated acoustic velocities (low-frequency) obtained by applying Biot–Gassmann fluid substitution on the ultrasonic velocities of dry cores. The models were trained on TUV data sets derived from three subsurface cores extracted from the geothermal well 16B(78)-32 at the Utah FORGE site. Each core was subjected to 75 unique stress configurations for velocity measurement in the dry state. The ML/DL trained on the TUV data set with equivalent saturated velocities demonstrated promising performance to predict in situ stress in subsurface geological rocks using velocity–stress relationships with R 2 of 0.86, 0.971, and 0.975 and root mean squared error (RMSE) of 2.59, 1.92, and 1.80 for validation/testing phases of vertical, minimum horizontal, and maximum horizontal stress models, respectively. Additionally, interpretation and explanation by Shapley additive explanations (SHAP) analysis further improved scientific validation and model reliability for estimating in situ stresses.

colloids↗
NASA NTRS1999

The Influence of Convection on Magnetotail Variability

This study investigates the evolution of the magnetotail's magnetic field with the aid of a self-consistent two-dimensional model. In this model the plasma mantle continuously supplies particles to the magnetotail, the ion current periodically updates the magnetic field using the Biot-Savart law. The simulated magnetotail evolves into a quasi-steady state, characterized by the periodic motion of the model's near-Earth X-line. This variability results from the nonadiabatic acceleration of ions in the current sheet and their rapid loss from the tail. The characteristic time scale of variability in the magnetotail is on the order of 4 - 5 minutes. We also investigate how the magnetotail's topology responds to increased convection electric fields, and show examples of observations of variability in the magnetotail.

Peroomian, Vahe↗
NASA NTRS2014

The Large-Scale Current System During Auroral Substorms

We present an empirical model of the equivalent current system in the ionosphere during the peak of a classical bulge-type auroral sub storm. This model is derived from measurements made by ~110 ground magnetometer stations during 116 substorms. The data are temporally and spatially organized using global auroral images obtained by the Polar Visible Imaging System Earth Camera. The empirical equivalent current system displays three key features: a poleward shift of the westward electrojet connecting the postmidnight and premidnight components; a polar cap swirl; and significantly different magnitudes of the postmidnight and premidnight westward electrojets. This leads us to propose a two-wedge current system linking the ionosphere to the magnetosphere. The bulge current wedge is located in the premidnight region just equatorward of the open-closed field line boundary while another three-dimensional current system is located in the postmidnight region well within the auroral oval. We use Biot and Savart calculations and Tsyganenko mapping and show that this new model is a likely solution for the large-scale current system.

magnetic storms↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2022

Current distribution monitoring enables quench and damage detection in superconducting fusion magnets

Abstract Fusion magnets made from high temperature superconducting ReBCO CORC ® cables are typically protected with quench detection systems that use voltage or temperature measurements to trigger current extraction processes. Although small coils with low inductances have been demonstrated, magnet protection remains a challenge and magnets are typically operated with little knowledge of the intrinsic performance parameters. We propose a protection framework based on current distribution monitoring in fusion cables with limited inter-cable current sharing. By employing inverse Biot-Savart techniques to distributed Hall probe arrays around CORC ® Cable-In-Conduit-Conductor (CICC) terminations, individual cable currents are recreated and used to extract the parameters of a predictive model. These parameters are shown to be of value for detecting conductor damage and defining safe magnet operating limits. The trained model is then used to predict cable current distributions in real-time, and departures between predictions and inverse Biot-Savart recreated current distributions are used to generate quench triggers. The methodology shows promise for quality control, operational planning and real-time quench detection in bundled CORC ® cables for compact fusion reactors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗
DOE OSTI2023

Time-lapse VSP integration and calibration of subsurface stress field utilizing machine learning approaches: A case study of the morrow B formation, FWU

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse vertical seismic profiling (VSP) integration. The selected study site is a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

58 GEOSCIENCES↗
NASA NTRS1995

The effects of vortex modeling on blade-vortex interaction noise prediction

The use of a blade vortex interaction noise prediction scheme, based on CAMRAD/JA, FPR and RAPP, quantifies the effects of errors and assumptions in the modeling of the helicopter's shed vortex on the acoustic predictions. CAMRAD/JA computes the wake geometry and inflow angles that are used in FPR to solve for the aerodynamic surface pressures. RAPP uses these surface pressures to predict the acoustic pressure. Both CAMRAD/JA and FPR utilize the Biot-Savart Law to determine the influence of the vortical velocities on the blade loading and both codes use an algebraic vortex model for the solid body rotation of the vortex core. Large changes in the specification of the vortex core size do not change the inplane wake geometry calculated by CAMRAD/JA and only slightly affect the out-of-plane wake geometry. However, the aerodynamic surface pressure calculated by FPR changes in both magnitude and character with small changes to the core size used by the FPR calculations. This in turn affects the acoustic predictions. Shifting the CAMRAD/JA wake geometry away from the rotor plane by 1/4 chord produces drastic changes in the acoustic predictions indicating that the prediction of acoustic pressure is extremely sensitive to the miss distance between the vortex and the blade and that this distance must be calculated as accurately as possible for acceptable noise predictions. The inclusion or exclusion of a vortex in the FPR-RAPP calculation allows for the determination of the relative importance of that vortex as a BVI noise source.

Gallman, Judith M.↗
NASA NTRS1999

Intrinsic Variability in the Quiet-Time Magnetotail

This study investigates the evolution of the magnetotail's magnetic field with the aid of a self-consistent two-dimensional model in which the ion current periodically updates the magnetic field. The plasma mantle supplies particles continuously to the magnetotail, and the perturbation magnetic field is calculated from the ion current using the Biot-Savart law. The simulated magnetotail evolves into a quasi-steady state, characterized by the periodic motion of the near-Earth X-line in the model. This variability is caused by the nonadiabatic acceleration of ions in the current sheet and their rapid loss from the tail. Particularly noteworthy is the value found for the characteristic time scale of variability in the magnetotail. on the order of 4 - 5 minutes.

Peroomian, Vahe↗
DOE OSTI2026

Flash Heating of a Hollow Rod to Determine Thermal Properties

In a nuclear reactor, prolonged exposure to high neutron flux levels can change material properties of reactor components. Due to the danger posed to personnel by highly radioactive components such as spent fuel, noncontact methods must be used to determine these properties. The objective of the current research is to test the viability of using a flash-heating method to determine the thermal diffusivity of a silicon carbide (SiC) tube, which could be used as cladding for a nuclear fuel rod. The flash-heating experiment was conducted at Oak Ridge National Laboratory and temperature results were analyzed as a two-Dimensional cross section which assumed axial symmetry. The model was fitted iteratively to the temperature measurements using nonlinear regression which required values for thermal diffusivity, Biot number and a term to determine the magnitude of the heat absorbed by the flash. The method could be expanded in the future to estimate the thermal diffusivity of spent nuclear fuel rods in an effort to increase the certainty of new reactor designs using precise thermal properties throughout the core life. A statistical analysis of the results of this work is provided as part of the analysis.

McMasters, Robert L [Virginia Military Institute, ↗
NASA NTRS2013

Development of Polyimide Foam for Aircraft Sidewall Applications

In this paper, the use of polyimide foam as a lining in double panel applications is considered. It is being investigated here as a replacement for aircraft grade glass fiber and has a number of attractive functional attributes, not the least of which is its high fire resistance. The test configuration studied here consisted of two 1mm (0.04 in.) thick, flat aluminum panels separated by 12.7 cm (5.0 in.) with a 7.6 cm (3.0 in.) thick layer of foam centered in that space. Random incidence transmission loss measurements were conducted on this buildup, and conventional poro-elastic models were used to predict the performance of the lining material. Results from two densities of foam are considered. The Biot parameters of the foam were determined by a combination of direct measurement (for density, flow resistivity and Young s modulus) and inverse characterization procedures (for porosity, tortuosity, viscous and thermal characteristic length, Poisson s ratio and loss factor). The inverse characterization procedure involved matching normal incidence standing wave tube measurements of absorption coefficient and transmission loss of the isolated foam with finite element predictions. When the foam parameters determined in this way were used to predict the performance of the complete double panel system, reasonable agreement was obtained between the measured transmission loss and predictions made using a commercial statistical energy analysis code.

Silcox, Richard↗
DOE OSTI2021

Impacts of Anisotropic Porosity on Heat Transfer and Off-Gassing during Biomass Pyrolysis

The pore structure of biogenic materials imbues the ability to deliver water and nutrients through a plant from root to leaf. This anisotropic pore granularity can also play a significant role in processes such as biomass pyrolysis that are used to convert these materials into useful products like heat, fuel, and chemicals. Evolutions in modeling of biomass pyrolysis as well as imaging of pore structures allow for further insights into the concerted physics of phase change-induced off-gassing, heat transfer, and chemical reactions. In this work, we report a biomass single particle model which incorporates these physics to explore the impact of implementing anisotropic permeability and diffusivity on the conversion time and yields predicted for pyrolysis of oak and pine particles. Simulation results showed that anisotropic permeability impacts predicted conversion time more than 2 times when the Biot number is above 0.1 and pyrolysis numbers (Py 1 , Py 2 ) are less than 20. Pore structure significantly impacts predicted pyrolytic conversion time (>8 times) when the Biot number is above 1 and the pyrolysis number is below 1, i.e., the “conduction controlled” regime. Therefore, these nondimensional numbers reflect that when internal heat conduction limits pyrolysis performance, internal pyrolysis off-gassing further retards effective heat transfer rates as a closely coupled phenomenon. Overall, this study highlights physically meaningful opportunities to improve particle-scale pyrolysis modeling and experimental validation relevant to a variety of feedstock identities and preparations, guiding the future design of pyrolyzers for efficient biomass conversion.

09 BIOMASS FUELS↗
DOE OSTI2021

Biomass pyrolysis devolatilization kinetics of herbaceous and woody feedstocks

Devolatilization kinetics were determined using a modified micropyrolyzer reactor for several biomass feedstocks: switchgrass, corn stover, red oak, and pine. Here, the micropyrolyzer was directly coupled to a flame ionization detector (FID) to track the release of volatiles from the biomass. Time series data from these experiments was analyzed to determine apparent devolatilization rates. Care was taken to assure the experiments were isothermal and kinetically limited calculating a Biot number less than 0.1 and pyrolysis numbers greater than 10, which simplifies the derivation of devolatilization rates. A single, first order reaction was able to model devolatilization rates at temperatures up to 500 °C. No correlation was found between the inorganic content of the biomass and its rate of devolatilization. Apparent activation energies were in the range of 54.9–88.4 kJ mol –1 . The rate coefficient at 500 °C was calculated as 1.90–5.14 s –1 for the four feedstocks.

09 BIOMASS FUELS↗
NASA NTRS1979

Three-dimensional simulation of the free shear layer using the vortex-in-cell method

We present numerical simulations of the evolution of a mixing layer from an initial state of uniform vorticity with simple two- and three-dimensional small perturbations. A new method for tracing a large number of three-dimensional vortex filaments is used in the simulations. Vortex tracing by Biot-Savart interaction originally implied ideal (non-viscous) flow, but we use a 3-d mesh, Fourier transforms and filtering for vortex tracing, which implies 'modeling' of subgrid scale motion and hence some viscosity. Streamwise perturbations lead to the usual roll-up of vortex patterns with spanwise uniformity maintained. Remarkably, spanwise perturbations generate streamwise distortions of the vortex filaments and the combination of both perturbations leads to patterns with interesting features discernable in the movies and in the records of enstrophy and energy for the three components of the flow.

Couet, B.↗
NASA NTRS1994

Effects of leading and trailing edge flaps on the aerodynamics of airfoil/vortex interactions

A numerical procedure has been developed for predicting the two-dimensional parallel interaction between a free convecting vortex and a NACA 0012 airfoil having leading and trailing edge integral-type flaps. Special emphasis is placed on the unsteady flap motion effects which result in alleviating the interaction at subcritical and supercritical onset flows. The numerical procedure described here is based on the implicit finite-difference solutions to the unsteady two-dimensional full potential equation. Vortex-induced effects are computed using the Biot-Savart Law with allowance for a finite core radius. The vortex-induced velocities at the surface of the airfoil are incorporated into the potential flow model via the use of the velocity transpiration approach. Flap motion effects are also modeled using the transpiration approach. For subcritical interactions, our results indicate that trailing edge flaps can be used to alleviate the impulsive loads experienced by the airfoil. For supercritical interactions, our results demonstrate the necessity of using a leading edge flap, rather than a trailing edge flap, to alleviate the interaction. Results for various time-dependent flap motions and their effect on the predicted temporal sectional loads, differential pressures, and the free vortex trajectories are presented

Hassan, Ahmed A.↗
NASA NTRS1987

Biot theory and acoustical properties of high porosity fibrous materials and plastic foams

Experimental values of acoustic wave propagation constant and characteristic impedance in fibrous materials, and normal absorption for two plastic foams, were compared to theoretical predictions obtained with Biot's theory. The best agreement was observed for fibrous materials between Biot's theory and Delany and Bazley experiments for a nearly zero mass coupling parameter. For foams, the lambda/4 structure resonance effect on absorption was calculated by using four-pole modelling of the medium. A significant mass coupling parameter is then necessary for obtaining agreement between the behavior of the measured absorption coefficients and the theoretical predictions. It is shown how the formalism used for predicting foams absorption coefficients may be used for studying the acoustic behavior of multi-layered media.

Allard, J.↗
NASA NTRS1983

Thermal Modeling of Bridgman Crystal Growth

Heat Flow modeled for moving or stationary rod shaped sample inside directional-solidification furnace. Program effectively models one-dimensional heat flow in translating or motionless rod-shaped sample inside of directionalsolidification furnace in which adiabatic zone separates hot zone and cold zone. Applicable to systems for which Biot numbers in hot and cold zones are less than unity.

Cothran, E.↗