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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

51 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1989

A parametric heat transfer study for cryogenic ball bearings in SSME HPOTP

A numerical modeling is to examine the effects of coolant convective heat transfer coefficient and frictional heating on the local temperature characteristics of a ball element in Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP) bearing. The present modeling uses a control-volume based, finite-difference method to solve the non-dimensionalized heat conduction equation in spherical coordinate system. The dimensionless temperature is found as a function of Biot number, heat flux ratio between the two race contacts, and location in the ball. The current results show that, for a given cooling capability, the ball temperature generally increases almost linearly with the heat input from the race-contacts. This increase is always very high at one of the two contacts. An increase in heat transfer coefficient generally reduces the ball temperature and alleviates the temperature gradient, except for the regions very close to the race contacts. For a 10-fold increase of heat transfer coefficient, temperature decrease is 35 percent for the average over entire ball, and 10 percent at the inner-race contact. The corresponding change of temperature gradient displays opposing trends between the regions immediately adjacent to the contacts and the remaining portion of the ball. The average temperature gradient in the vicinity of both contacts increases approximately 70 to 100 percent. A higher temperature gradient produces excessive thermal stress locally which may be detrimental to the material integrity. This, however, is the only unfavorable issue for an increase of heat transfer coefficient.

Chyu, Mingking K.↗

DOE OSTI2020

Thermal conductivity of two-dimensional disordered fibrous materials defined by interfiber thermal contact conductance and intrinsic conductivity of fibers

A general framework for the theoretical analysis and numerical calculations of the effective thermal conductivity of two-dimensional homogeneous and isotropic disordered fibrous materials is developed in this work based on the model of random contacts between soft-core spherocylinders. The analysis accounts for the interfiber contact conductance and intrinsic conductivity of fibers and is performed in a wide space of governing parameters that includes the fiber aspect ratio, Biot number calculated for a single thermal contact between fibers, and material density ranging from values corresponding to the percolation threshold up to those characteristic of dense fiber networks. For dense networks, exact theoretical equations for the thermal conductivity of materials composed of spherocylinders with an arbitrary aspect ratio and Biot number are derived. The effect of the intrinsic conductivity of fibers on the thermal transport in fibrous materials is found to depend on the density of contacts and can be significant in sufficiently dense fiber networks even if the Biot number for a single thermal contact is small. Semiempirical corrections to the theoretical equations are derived for small and moderate fiber densities. The power law exponent describing the approximate dependence of the conductivity on the density parameter is found to vary from values close to 1 up to values exceeding 2 when evaluated within different finite ranges of the density parameter. This finding explains the variability of scaling laws for thermal conductivity of fibrous materials suggested in the literature based on numerical simulations performed in different regions of the space of material parameters.

36 MATERIALS SCIENCE↗

NASA NTRS1995

A test of a vortex method for the computation of flap side edge noise

Upon approach to landing, a major source location of airframe noise occurs at the side edges of the part span, trailing edge flaps. In the vicinity of these flaps, a complex arrangement of spanwise flow with primary and secondary tip vortices may form. Each of these vortices is observed to become fully three-dimensional. In the present study, a numerical model is developed to investigate the noise radiated from the side edge of a flap. The inherent three-dimensionality of this flow forces us to carefully consider a numerical scheme which will be both accurate in its prediction of the flow acoustics and also computationally efficient. Vortex methods have offered a fast and efficient means of simulating many two and three-dimensional, vortex dominated flows. In vortex methods, the time development of the flow is tracked by following exclusively the vorticity containing regions. Through the Biot-Savart law, knowledge of the vorticity field enables one to obtain flow quantities at any desired location during the flow evolution. In the present study, a numerical procedure has been developed which incorporates the Lagrangian approach of vortex methods into a calculation for the noise radiated by a flow-surface interaction. In particular, the noise generated by a vortex in the presence of a flat half plane is considered. This problem serves as a basic model of flap edge flow. It also permits the direct comparison between our computed results and previous acoustic analyses performed for this problem. In our numerical simulations, the mean flow is represented by the complex potential W(z) = Aiz(exp l/2), which is obtained through conformal mapping techniques. The magnitude of the mean flow is controlled by the parameter A. This mean flow has been used in the acoustic analysis by Hardin and is considered a reasonable model of the flow field in the vicinity of the edge and away from the leading and trailing edges of the flap. To represent the primary vortex which occurs near the flap, a point vortex is introduced just below the flat half plane. Using a technique from panel methods, boundary conditions on the flap surface are satisfied by the introduction of a row of stationary point vortices along the extent of the flap. At each time step in the calculation, the strength of these vortices is chosen to eliminate the normal velocity at intermediary collocation points. The time development of the overall flow field is then tracked using standard techniques from vortex methods. Vortex trajectories obtained through this computation are in good agreement with those predicted by the analytical solution given by Hardin, thus verifying the viability of this procedure for more complex flow arrangements. For the flow acoustics, the Ffowcs Williams-Hawkings equation is numerically integrated. This equation supplies the far field acoustic pressure based upon pressures occurring along the flap surface. With our vortex method solution, surface pressures may be obtained with exceptional resolution. The Ffowcs Williams-Hawkings equation is integrated using a spatially fourth order accurate Simpson's rule. Rational function interpolation is used to obtain the surface pressures at the appropriate retarded times. Comparisons between our numerical results for the acoustic pressure and those predicted by the Hardin analysis have been made. Preliminary results indicate the need for an improved integration technique. In the future, the numerical procedure developed in this study will be applied to the case of a rectangular flap of finite thickness and ultimately modified for application to the fully three-dimensional problem.

Martin, James E.↗

NASA NTRS1982

Theoretical magnetograms based on quantitative simulation of a magnetospheric substorm

Substorm currents derived from the Rice University computer simulation of the September 19, 1976 substorm event are used to compute theoretical magnetograms as a function of universal time for various stations, integrating the Biot-Savart law over a maze of about 2700 wires and bands that carry the ring, Birkeland and horizontal ionospheric currents. A comparison of theoretical results with corresponding observations leads to a claim of general agreement, especially for stations at high and middle magnetic latitudes. Model results suggest that the ground magnetic field perturbations arise from complicated combinations of different kinds of currents, and that magnetic field disturbances due to different but related currents cancel each other out despite the inapplicability of Fukushima's (1973) theorem. It is also found that the dawn-dusk asymmetry in the horizontal magnetic field disturbance component at low latitudes is due to a net downward Birkeland current at noon, a net upward current at midnight, and, generally, antisunward-flowing electrojets.

Chen, C.-K.↗

NASA NTRS1994

Method for confining the magnetic field of the cross-tail current inside the magnetopause

A method is presented for analytically representing the magnetic field due to the cross-tail current and its closure on the magnetopause. It is an extension of a method used by Tsyganenko (1989b) to confine the dipole field inside an ellipsoidal magnetopause using a scalar potential. Given a model of the cross-tail current, the implied net magnetic field is obtained by adding to the cross-tail current field a potential field B = - del gamma, which makes all field lines divide into two disjoint groups, separated by the magnetopause (i.e., the combined field is made to have zero normal component with the magnetopause). The magnetopause is assumed to be an ellipsoid of revolution (a prolate spheroid) as an approximation to observations (Sibeck et al., 1991). This assumption permits the potential gamma to be expressed in spheroidal coordinates, expanded in spheroidal harmonics and its terms evaluated by performing inversion integrals. Finally, the field outside the magnetopause is replaced by zero, resulting in a consistent current closure along the magnetopause. This procedure can also be used to confine the modeled field of any other interior magnetic source, though the model current must always flow in closed circuits. The method is demonstrated on the T87 cross-tail current, examples illustrate the effect of changing the size and shape of the prescribed magnetopause and a comparison is made to an independent numerical scheme based on the Biot-Savart equation.

Sotirelis, T.↗

NASA NTRS1995

Transport phenomena in the micropores of plug-type phase separators

This study numerically investigates the transport phenomena within and across a porous-plug phase separator. The effect of temperature differential across a single pore and of the sidewall boundary conditions, i.e., isothermal or linear thermal gradient, are presented and discussed. The effects are quantified in terms of the evaporation mass flux across the boundary and the mean surface temperature. A two-dimensional finite element model is used to solve the continuity, momentum, and energy equations for the liquid. Temperature differentials across the pore interface of 1.0, and 1.5 K are examined and their effect on evaporation flux and mean surface temperature is shown. For isothermal side boundary conditions, the evaporation flux across the pore is directly proportional and linear with Delta T. For the case of an imposed linear thermal gradient on the side boundaries, Biot numbers of 0.0, 0.15, and 0.5 are examined. The most significant effect of Biot number is to lower the overall surface temperature and evaporation flux.

Fazah, M. M.↗

DOE OSTI2021

Simulations of biomass pyrolysis using glued-sphere CFD-DEM with 3-D intra-particle models

The shape of biomass particles influences intra-particle heat transfer, particle-particle collisions, interphase momentum, heat transfer, and pyrolysis yield. In this research, the effects of these multi-scale phenomena are captured using a glued-sphere computational fluid dynamic discrete element method (CFD-DEM) with 3-D intra-particle models. The glued-spheres resolved the shapes of particles and intra-particle temperature and species distributions. The implementation was validated using the pyrolysis data of a spherical particle and a cylindrical particle. Then, the influences of biomass shapes and sizes on pyrolysis were investigated. The results revealed the significant influence of surface areas and mass distributions on interphase heat transfer and chemical reactions. The intra-particle heat transfer is neglectable for a Biot number smaller than 0.41. Finally, the pyrolysis of irregularly shaped biomass particles in a fixed bed, fluidized bed, and spouted bed was simulated. Analysis of conversion distributions showed the largest deviations in the fixed bed and the slowest conversion rate in the spouted bed. A good balance between uniform and fast conversion is achieved in the fluidized bed. The multi-scale method proposed in this research provides a tool for the evaluation of large irregular-shaped biomass particle pyrolysis in different types of reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Interpreting T-History Data on Supercooling and Phase Change Enthalpy of a Salt-Hydrogel Complex for Building Thermal Management

Phase change materials (PCMs) are promising materials for storing and discharging latent heat in buildings to regulate the thermal environment. Measurement of salt hydrate PCM thermal performance is often challenging because of supercooling and phase segregation. Though differential scanning calorimetry (DSC) has been widely used to measure thermal properties of various phase change materials, DSC uses milligrams of sample and is unrealistic for predicting large-scale materials performance. The T-history method [1], in contrast, handles higher volumes but the analysis typically requires a lumped capacitance assumption. Here, we report and analyze T-History data on a novel salt-hydrogel complex with high viscosity, that renders low Biot number sample preparation extremely challenging. We re-examine the standard T-history method and compared against DSC data in this complex material. We report details of an experimental setup, built with necessary additions to a controlled environmental chamber. We carefully examine the validity of the lumped capacitance assumption and suggest a modified analysis using a computational approach for reliable data. Experimental parameters such as ramp rate and temperature range of the environmental chamber are considered via computational modeling to understand the effect of ramp rate on supercooling [2]. Finally, we report thermal cycling experiments on the salt-hydrogel complex. Our modified approach to measuring supercooling and enthalpy of fusion at large scales is important in better understanding the performance of phase-change materials at scale for building thermal storage. References: [1] Marin, Jose M, Belen Zalba, Luisa F Cabeza, and Harald Mehling. “Determination of Enthalpy Temperature Curves of Phase Change Materials with the Temperature-History Method: Improvement to Temperature Dependent Properties.” Measurement Science and Technology 14, no. 2 (February 1, 2003): 184–89. https://doi.org/10.1088/0957-0233/14/2/305. [2] Safari, A., R. Saidur, F.A. Sulaiman, Yan Xu, and Joe Dong. “A Review on Supercooling of Phase Change Materials in Thermal Energy Storage Systems.” Renewable and Sustainable Energy Reviews 70 (April 2017): 905–19. https://doi.org/10.1016/j.rser.2016.11.272

thermodynamics↗

NASA NTRS1993

A global analytical representation of the magnetic field produced by the region 2 Birkeland currents and the partial ring current

A quantitative model is developed of the magnetic field produced by the electric current system of region 2 Birkeland currents, closed via the partial ring current. The distribution of j-perpendicular is computed from a given axially asymmetric spatial distribution of hot isotropic magnetospheric plasma over an infinitely thin L shell in an axisymmetric purely dipolar geomagnetic field, while the field-aligned current density is found from the continuity of the net electric current. The magnetic field distribution is derived by a Biot-Savart integral over the electric current system. An assumed cosine dependence of the plasma pressure on local time makes it possible to reduce the problem of analytical representation of the B field to two dimensions. The obtained numerical fits for the partial ring current/region 2 Birkeland current magnetic field are relatively simple, continuous, and valid throughout the whole extraterrestrial space from ionospheric heights up to tens of earth radii.

Tsyganenko, N. A.↗

NASA NTRS1987

A new methodology for free wake analysis using curved vortex elements

A method using curved vortex elements was developed for helicopter rotor free wake calculations. The Basic Curve Vortex Element (BCVE) is derived from the approximate Biot-Savart integration for a parabolic arc filament. When used in conjunction with a scheme to fit the elements along a vortex filament contour, this method has a significant advantage in overall accuracy and efficiency when compared to the traditional straight-line element approach. A theoretical and numerical analysis shows that free wake flows involving close interactions between filaments should utilize curved vortex elements in order to guarantee a consistent level of accuracy. The curved element method was implemented into a forward flight free wake analysis, featuring an adaptive far wake model that utilizes free wake information to extend the vortex filaments beyond the free wake regions. The curved vortex element free wake, coupled with this far wake model, exhibited rapid convergence, even in regions where the free wake and far wake turns are interlaced. Sample calculations are presented for tip vortex motion at various advance ratios for single and multiple blade rotors. Cross-flow plots reveal that the overall downstream wake flow resembles a trailing vortex pair. A preliminary assessment shows that the rotor downwash field is insensitive to element size, even for relatively large curved elements.

Bliss, Donald B.↗

NASA NTRS1992

Mach number effects on free and wall turbulent structures in light of instability flow interactions

Our perception of fluid instabilities and turbulence is largely based on properties of incompressible vorticity distributions, omega, and the Biot-Savart Law of 'instantaneous' induction at distance. In high-speed shear layers it is the distribution of local angular momentum rho-omega that matters; modified interaction-induction between vortical elements is likely to take place only within relative Mach cones of influence, with acoustic time delay. This and the large mean density stratifications influence strongly instability vortex roll-ups, vortex mergings, and large-scale turbulent structures. Consequences of these observations on transition to turbulence and turbulence energetics are discussed in physical terms. Free shear layers and wall layers may require distinct M modeling.

Morkovin, Mark V.↗

DOE OSTI2024

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗

NASA NTRS1995

Thermocapillary flow with evaporation and condensation at low gravity. Part 2: Deformable surface

The free surface behavior of a volatile wetting liquid at low gravity is studied using scaling and numerical techniques. An open cavity model, which was applied in part 1 to investigate fluid flow and heat transfer in non-deforming pores, is used to evaluate the influence of convection on surface morphology with length scales and subcooling/superheating limits of 1 less than or equal to D less than or equal to 10(exp 2) microns and approximately 1 K, respectively. Results show that the menisci shapes of highly wetting fluids are sensitive to thermocapillary flow and to a lesser extent the recoil force associated with evaporation and condensation. With subcooling, thermocapillarity produces a suction about the pore centerline that promotes loss of mechanical equilibrium, while condensation exerts an opposing force that under some conditions offsets this destabilizing influence. With superheating, thermocapillarity and evaporation act in the same direction and mutually foster surface stability. All of these trends are magnified by high capillary and Biot numbers, and the stronger circulation intensities associated with small contact angles. These phenomena strongly depend on the thermal and interfacial equilibrium between the liquid and vapor, and have important ramifications for systems designed to maintain a pressure differential across a porous surface.

Schmidt, G. R.↗

DOE OSTI2022

Universal Anomalous Diffusion of Quantized Vortices in Ultraquantum Turbulence

In classical viscous fluids, turbulent eddies are known to be responsible for the rapid spreading of embedded particles. However, in inviscid quantum fluids where the turbulence is induced by a chaotic tangle of quantized vortices, dispersion of the particles can be achieved via a nonclassical mechanism, i.e., their binding to the evolving vortices. However, knowledge on how the vortices diffuse and spread in quantum-fluid turbulence is very limited, especially for the so-called ultraquantum turbulence (UQT) generated by a random tangle of vortices. Here we report a systematic numerical study of the apparent diffusion of vortices in UQT in superfluid helium-4 using the full Biot-Savart simulation. We reveal that the vortices in the superfluid exhibit a universal anomalous diffusion (superdiffusion) at small times, which transits to normal diffusion at large times. This behavior is found to be the result of a generic scaling property of the vortex velocity. Our simulation at finite temperatures also nicely reproduces recent experimental observations. Lastly, the knowledge obtained from this study may form the base for understanding turbulent transport and universal vortex dynamics in various quantum fluids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Monolithic Multigrid for a Reduced-Quadrature Discretization of Poroelasticity

Advanced finite-element discretizations and preconditioners for models of poroelasticity have attracted significant attention in recent years. The equations of poroelasticity offer significant challenges in both areas, due to the potentially strong coupling between unknowns in the system, saddle-point structure, and the need to account for wide ranges of parameter values, including limiting behavior such as incompressible elasticity. This paper was motivated by an attempt to develop monolithic multigrid preconditioners for the discretization developed in [C. Rodrigo et al., Comput. Methods App. Mech. Engrg, 341 (2018), pp. 467--484]; we show here why this is a difficult task and, as a result, we modify the discretization in [Rodrigo et al.] through the use of a reduced-quadrature approximation, yielding a more “solver-friendly” discretization. Local Fourier analysis is used to optimize parameters in the resulting monolithic multigrid method, allowing a fair comparison between the performance and costs of methods based on Vanka and Braess--Sarazin relaxation. Further, numerical results are presented to validate the local Fourier analysis predictions and demonstrate efficiency of the algorithms. Finally, a comparison to existing block-factorization preconditioners is also given.

97 MATHEMATICS AND COMPUTING↗