Search NASA⌕ Search

SEARCH · Search NASA

Results for “Immersed Moving Boundary”

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.

29 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
DOE OSTI2022

A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles

This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.

42 ENGINEERING↗
DOE OSTI2025

On the motion of compact objects in relativistic viscous fluids

We present a world-line effective field theory of compact objects moving relativistically through a viscous fluid. The theory is valid when velocity gradients are small compared to the inverse size of the object. Working within the EFT eliminates the need to solve a boundary value problem by turning all interactions between the fluid and the object into a source term in the action. We use the EFT to derive the relativistic equations of motion for a compact object immersed in a viscous fluid in a curved background, when the relative velocity of the object and the fluid is small compared to the speed of light.

astrophysical black holes↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1959

An Investigation of the Effect of a Highly Favorable Pressure Gradient on Boundary-Layer Transition as Caused by Various Types of Roughnesses on a 10-foot-Diameter Hemisphere at Subsonic Speeds

Tests were made on a 10-foot-diameter hemispherical nose at Reynolds numbers up to 10 x 10(exp 6) and at a maximum Mach number of about 0.1 to determine the effects of a highly favorable pressure gradient on boundary-layer transition caused by roughness. Both two-dimensional and three-dimensional roughness particles were used, and the transition of the boundary layer was determined by hot-wire anemometers. The roughness Reynolds number for transition R(sub k,t) caused by three-dimensional particles such as Carborundum grains, spherical particles, and rimmed craters was found. The results show that for particles immersed in the boundary layer, R(sub k,t) is independent of the particle size or position on the hemispherical nose and depends mainly on the height-to-width ratio of the particle. The values of R(sub k,t) found on the hemispherical nose compare closely with those previously found on a flat plate and on airfoils with roughness. For two-dimensional roughness, the ratio of roughness height to boundary-layer displacement thickness necessary to cause transition was found to increase appreciably as the roughness was moved forward on the nose. Also included in the investigation were studies of the spread of turbulence behind a single particle of roughness and the effect of holes such as pressure orifices.

Peterson, John B., Jr.↗
NASA NTRS2002

A New Time Domain Formulation for Broadband Noise Predictions

A new analytic result in acoustics called "Formulation 1B," proposed by Farassat, is used to compute the loading noise from an unsteady surface pressure distribution on a thin airfoil in the time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the loading source term. The formulation contains a far field surface integral that depends on the time derivative and the surface gradient of the pressure on the airfoil, as well as a contour integral on the boundary of the airfoil surface. As a first test case, the new formulation is used to compute the noise radiated from a flat plate, moving through a sinusoidal gust of constant frequency. The unsteady surface pressure for this test case is analytically specified from a result based on linear airfoil theory. This test case is used to examine the velocity scaling properties of Formulation 1B and to demonstrate its equivalence to Formulation 1A of Farassat. The new acoustic formulation, again with an analytic surface pressure, is then used to predict broadband noise radiated from an airfoil immersed in homogeneous, isotropic turbulence. The results are compared with experimental data previously reported by Paterson and Amiet. Good agreement between predictions and measurements is obtained. Finally, an alternative form of Formulation 1B is described for statistical analysis of broadband noise.

Casper, J.↗
NASA NTRS2002

A New Time Domain Formulation for Broadband Noise Predictions

A new analytic result in acoustics called "Formulation 1B," proposed by Farassat, is used to compute the loading noise from an unsteady surface pressure distribution on a thin airfoil in the time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the loading source term. The formulation contains a far field surface integral that depends on the time derivative and the surface gradient of the pressure on the airfoil, as well as a contour integral on the boundary of the airfoil surface. As a first test case, the new formulation is used to compute the noise radiated from a flat plate, moving through a sinusoidal gust of constant frequency. The unsteady surface pressure for this test case is analytically specied from a result based on linear airfoil theory. This test case is used to examine the velocity scaling properties of Formulation 1B and to demonstrate its equivalence to Formulation 1A of Farassat. The new acoustic formulation, again with an analytic surface pressure, is then used to predict broadband noise radiated from an airfoil immersed in homogeneous, isotropic turbulence. The results are compared with experimental data previously reported by Paterson and Amiet. Good agreement between predictions and measurements is obtained. Finally, an alternative form of Formulation 1B is described for statistical analysis of broadband noise.

Casper, Jay H.↗
NASA NTRS2002

Broadband Noise Predictions Based on a New Aeroacoustic Formulation

A new analytic result in acoustics called 'Formulation 1B,' proposed by Farassat, is used to compute the loading noise from an unsteady surface pressure distribution on a thin airfoil in the time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the loading source term. The formulation contains a far-field surface integral that depends on the time derivative and the surface gradient of the pressure on the airfoil, as well as a contour integral on the boundary of the airfoil surface. As a first test case, the new formulation is used to compute the noise radiated from a flat plate, moving through a sinusoidal gust of constant frequency. The unsteady surface pressure for this test case is specified analytically from a result that is based on linear airfoil theory. This test case is used to examine the velocity scaling properties of Formulation 1B, and to demonstrate its equivalence to Formulation 1A, of Farassat. The new acoustic formulation, again with an analytic surface pressure, is then used to predict broadband noise radiated from an airfoil immersed in homogeneous turbulence. The results are compared with experimental data previously reported by Paterson and Amiet. Good agreement between predictions and measurements is obtained. The predicted results also agree very well with those of Paterson and Amiet, who used a frequency-domain approach. Finally, an alternative form of Formulation 1B is described for statistical analysis of broadband noise.

Casper, J.↗
NASA NTRS2018

Understanding Space Weather Part III: the Sun's Domain

The Sun exports a continuous outflow of plasma into interplanetary space: the solar wind. The solar wind primarily comprises two components: high- and slow-speed flows. These move with velocities ranging from 200 to 800 km/s depending on the source of the particular flow. As well as its speed, the density, temperature, and even the composition of the solar wind change. Adding to its intrinsic variability, there are embedded transients resulting from flares and coronal mass ejections that further complicate its dynamics and space weather impacts. The solar wind interacts differently with each of the solar system objects it encounters based on their magnetic and atmospheric properties. Even more complex processes occur as the solar wind encounters the interstellar medium, at the outer boundaries of the Sun's domain. The solar wind stretches to beyond 100 au (where 1 au defined to be equal to 149 597 870 700 m) from the Sun, which means that Earth is essentially immersed in the very hot solar atmosphere, and that leads to many space weather impacts on life and society. The specific space weather impacts on Earth will be discussed in detail in the next two papers in this series.

Strong, Keith↗
DOE OSTI2025

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Enhanced Boiling on Micro-Configured Composite Surfaces Under Microgravity Conditions

In order to accommodate the growing thermal management needs of future space platforms, several two-phase active thermal control systems (ATCSs) have evolved and were included in the designs of space stations. Compared to the pumped single-phase liquid loops used in the conventional Space Transportation System and Spacelab, ATCSs offer significant benefits that may be realized by adopting a two-phase fluid-loop system. Alternately, dynamic power systems (DPSs), based on the Rankine cycle, seem inevitably to be required to supply the electrical power requirements of expanding space activities. Boiling heat transfer is one of the key technologies for both ATCSs and DPSs. Nucleate boiling near critical heat flux (CHF) can transport very large thermal loads with much smaller device size and much lower pumping power. However, boiling performance deteriorates in a reduced gravity environment and operation in the CHF regime is precarious because any slight overload will cause the heat transfer to suddenly move to the film boiling regime, which in turn, will result in burnout of the heat transfer surfaces. New materials, such as micro-configured metal-graphite composites, can provide a solution for boiling enhancement. It has been shown experimentally that this type of material manifests outstanding boiling heat transfer performance and their CHF is also extended to higher values. Due to the high thermal conductivity of graphite fiber (up to 1,200 W/m-K in the fiber direction), the composite surfaces are non-isothermal during the boiling process. The composite surfaces are believed to have a much wider safe operating region (a more uniform boiling curve in the CHF regime) because non-isothermal surfaces have been found to be less sensitive to variations of wall superheat in the CHF regime. The thermocapillary forces formed by the temperature difference between the fiber tips and the metal matrix play a more important role than the buoyancy in the bubble detachment, for the bubble detachment manifests itself by a necking process which should not be weakened by reduced gravity. In addition, the composite surfaces introduce no extra pressure drop, no fouling and do not impose significant primary or maintenance costs. All of these suggest that this type of composite is an ideal material for the challenge of accounting for both reliability and economy of the relevant components applied in the ATCSs, the DPSs and other devices in future space missions. The aim of the proposed work is to experimentally investigate high nucleate pool boiling performance on a micro-configured metal-graphite composite surface and to determine the mechanisms of the nucleate boiling heat transfer both experimentally and theoretically. Freon-113 and water will be used as the test liquids to investigate wettability effects on boiling characteristics. The Cu-Gr and Al-Gr composites with various volume fractions of graphite fibers will be tested to obtain the heat transfer characteristic data in the nucleate boiling region and in the CHF regime. In the experiments, the bubble emission and coalescence processes will be recorded by a video camera with a magnifying borescope probe immersed in the working fluid. The temperature profile in the thermal boundary layer on the composite surfaces will be measured by a group of micro thermocouples consisting of four ultra fine micro thermocouples. This instrument was developed and successfully used to measure the temperature profile of evaporating liquid thin layers by the proposers in a study performed at the NASA/Lewis Research Center. A two tier model to explain the nucleate boiling process and the performance enhancement on the composite surfaces has been suggested by the authors. According to the model, the thicknesses of the microlayer and the macrolayer underneath the bubbles and mushrooms, can be estimated by the geometry of the composite surface. The experimental results will be compared to the predictions from the model, and in turn, to revise and improve it.

Zhang, Nengli↗