SEARCH · Engineering Papers

Results for “weak boundary conditions”

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.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2004

On the Final Mass of Giant Planets

In the core accretion model of giant planet formation, when the core reaches critical mass, hydrostatic equilibrium is no longer possible and gas accretion ensues. If the envelope is radiative, the critical core mass is nearly independent of the boundary conditions and is roughly M(sub crit) ~ 10Mass of the Earth (with weak dependence on the rate of planetesimal accretion M(sub core) and the disk opacity k). Given that such a core may form at the present location of Jupiter in a time comparable to its Type I migration time (10(exp 5) - 10(exp 6) years) provided that the nebula was significantly enhanced in solids with respect to the MMSN and stall at this location in a weakly turbulent (alpha approximately less than 10(exp -4) disk, it may be appropriate to assume that such objects inevitably form and drive the evolution of late-phase T Tauri star disks. Here we investigate the final masses of giant planets in disks with one or more than one such cores. Although the presence of several planets would lead to Type II migration (due to the effective viscosity resulting from the planetary tidal torques), we ignore this complication for now and simply assume that each core has stalled at its location in the disk. Once a core has achieved critical mass, its gaseous accretion is governed by the given Kelvin-Helmholtz timescale.

Estrada, P. R.↗

NASA NTRS2001

A Hybrid Numerical Analysis Method for Structural Health Monitoring

A new hybrid surface-integral-finite-element numerical scheme has been developed to model a three-dimensional crack propagating through a thin, multi-layered coating. The finite element method was used to model the physical state of the coating (far field), and the surface integral method was used to model the fatigue crack growth. The two formulations are coupled through the need to satisfy boundary conditions on the crack surface and the external boundary. The coupling is sufficiently weak that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume can be set up independently. Thus when modeling crack growth, the finite element mesh can remain fixed for the duration of the simulation as the crack mesh is advanced. This method was implemented to evaluate the feasibility of fabricating a structural health monitoring system for real-time detection of surface cracks propagating in engine components. In this work, the authors formulate the hybrid surface-integral-finite-element method and discuss the mechanical issues of implementing a structural health monitoring system in an aircraft engine environment.

Forth, Scott C.↗

NASA NTRS1977

Effect of flight on jet noise from supersonic underexpanded flows

Experiments on underexpanded cold jet flows from a convergent nozzle under simulated flight conditions have shown that a large periodic spinning motion of the jet can occur with greatly enhanced broadband noise production. Shadowgraph pictures indicate that this oscillatory jet motion accompanies the generation of random weak shock waves at the source. These waves appear to be generated at the point downstream of the nozzle exit where the shock cells in the jet begin to disappear. The weak shock waves propagate upstream and have been identified to be the cause of enhanced broadband jet noise production in flight. In addition, the results show that the boundary layer flow conditions over the outside of the primary nozzle (simulating engine cowl flow in flight) have a key role in the production of these random weak shock waves.

Sarohia, V.↗

DOE OSTI2020

Attribution of NAO Predictive Skill Beyond 2 Weeks in Boreal Winter

Weeks 3–6 averaged winter North Atlantic Oscillation (NAO) predictive skill in a state-of-the-art coupled climate prediction system is attributed to two principle sources: upper and lower boundary conditions linked to the stratosphere and El Niño-Southern Oscillation (ENSO), respectively. A 20-member ensemble of 45-day reforecasts over 1999–2015 is utilized, together with uninitialized simulations with the atmospheric component of the prediction system forced with observed radiative forcing and lower boundary conditions. NAO forecast skill for lead times out to 6 weeks is higher following extreme stratospheric polar vortex conditions (weak and strong vortex events) compared to neutral states. Enhanced weeks 3-6 NAO predictive skill for weak vortex events results primarily from stratospheric downward coupling to the troposphere, while enhanced skill for strong vortex events can be partly attributed to lower boundary forcing related to the ENSO phenomenon. Implications for forecast system development and improvement are discussed.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1988

Oscillations of drops in zero gravity with weak viscous effects

Nonlinear oscillations and other motions of large axially symmetric liquid drops in zero gravity are studied numerically by a boundary-integral method. The effect of small viscosity is included in the computations by retaining first-order viscous terms in the normal stress boundary condition. This is accomplished by making use of a partial solution of the boundary-layer equations which describe the weak vortical surface layer. Small viscosity is found to have a relatively large effect on resonant mode coupling phenomena.

Lundgren, T. S.↗

NASA NTRS1987

Solitary waves in the resonant phenomenon between a surface gravity wave packet and an internal gravity wave

A two-layer inviscid incompressible fluid system of intermediate depth is considered. A multiple-scales perturbation technique is applied to the basic equations and boundary conditions for a two-layer fluid system to derive a system of weakly nonlinear partial integrodifferential equations governing the resonant interaction between a surface gravity wave packet and an internal gravity wave at an intermediate depth, providing a bridge between the existing shallow and deep fluid theories. The convolution integral term in these equations accounts for the dispersion in the lower-layer fluid. An iterative fast Fourier transform scheme is developed to find solitary wave solutions to this system of equations. The overtaking collision of two pairs of solitary waves, simulated using a spectral method, is found to be inelastic. It is found that the amplitude of the solitary waves changes slightly after the collision. The phase shifts these solitary waves undergo was calculated numerically.

Sepulveda, Nicasio↗

DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗

NASA NTRS2005

Sensitivity of Stratospheric Retrievals from Radio Occultations on Upper Boundary Conditions

The main uncertainty in the stratospheric retrievals from GPS radio occultation (RO) measurements comes from the lack of reliable measurements in the upper stratosphere and above where the bending due to the neutral atmosphere is weak and residual ionospheric effects are strong. In this work, we quantify the bias and uncertainty of the refractivity and temperature retrievals due to different upper boundary strategies using a simulation study. We use lidar refractivity and temperature profiles as the input states in generating the synthetic occultations. Random noise levels commensurate with the CHAMP RO measurements are then added to the simulated data. Through this study, the sensitivity of stratospheric retrievals to upper boundary methods and parameters are examined. Such error characterizations are important prerequisites towards the effective use of GPS RO data in climate monitoring.

GPS↗

NASA NTRS1991

In situ backscatter measurements over Pacific Ocean using two coherent focused CO2 lidars

Preliminary results of CS lidar measurements over the Pacific Ocean using two CO2 Doppler lidars operating at 9.1 and 10.6 microns are presented. SNR calculations during real-time and postprocessing were performed using a multistep signal search and signal power estimation algorithm optimized for weak signal conditions. Typical aerosol structure observed in the troposphere was found to consist of a boundary layer, a cloud-pumped layer, and a clean background layer; backscatter from cloud-free air was found to vary over three to four orders of magnitude or more.

Rothermel, Jeffry↗

NASA NTRS1986

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations--potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomeclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗

NASA NTRS1989

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations: potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomenclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗

DOE OSTI2023

Atmospheric Turbulent Intermittency Over the Arctic Sea‐Ice Surface During the MOSAiC Expedition

Turbulent motions in the Arctic stable boundary layer are characterized by intermittency, but they are rarely investigated due to limited observations, in particular over the sea‐ice surface. In the present study, we explore the characteristics of turbulent intermittency over the Arctic sea‐ice surface using data collected during the Multidisciplinary drifting Observation for the Study of Arctic Climate expedition from October 2019 to September 2020. We first develop a new algorithm, which performs well in identifying the spectral gap over the Arctic sea‐ice surface. Then the characteristics of intermittency are investigated. It is found that the strength of intermittency increases under the conditions of light surface wind speed, small surface wind speed gradient, and strong surface air temperature gradient. The momentum flux, sensible heat flux, and latent heat flux calculated by raw eddy‐covariance fluctuations are overestimated by 3%, 10%, and 24%, respectively, because submesoscale motions are included. Furthermore, the characteristics of the atmospheric boundary layer structure under various intermittency conditions reveal that strong low‐level jets are favorable to surface turbulent motions that result in weak intermittency, while strong temperature inversions above the surface layer suppress surface turbulent motions and lead to strong intermittency.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1986

Analyses of broadband noise mechanisms of rotors

The various source mechanisms which generate broadband noise on a range of rotors are reviewed. Analyses of these mechanisms are presented and compared to existing experimental data. The sources considered are load fluctuations due to inflow turbulence, due to turbulent blade boundary layers passing the trailing edge, and due to tip vortex formation turbulence. Vortex shedding noise due to laminar boundary layers and blunt trailing edges is not considered in detail as it can be avoided in most cases. Present analyses are adequate to predict the spectra from a wide variety of experiments on fans, helicopter rotors, and wind turbines to within about 5 to 10 dB. Better knowledge of the inflow turbulence improves the accuracy of the predictions. Inflow turbulence noise depends strongly on ambient conditions and dominates at low frequencies. Trailing edge and tip vortex noise are important at higher frequencies if inflow turbulence is weak. Boundary layer trailing edge noise increases slowly with angle of attack but not as rapidly as tip vortex formation noise. Tip noise can be important at high angles of attack for wide chord, square edge tips.

George, A. R.↗

DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

NASA NTRS1982

On the instability of shock waves attached to wedges and cones

A numerical study is made of the stability of weak and strong shocks attached to wedges and cones. The stability is shown to depend on the boundary conditions imposed, in good agreement with the experimental observations and the principle of minimum entropy. The behavior of the solutions is explained by quasi-static considerations.

Salas, M. D.↗

NASA NTRS1974

Unsteady temperature distribution in volume reflectors

The development of unsteady temperature distributions in semitransparent scattering materials under the influence of arbitrary incident radiative flux with other than specified temperature boundary conditions is considered. Specifically, such a problem is considered for an approximate radiative transfer model which allows analytic solutions. The model represents the unsteady, one-dimensional, radiative heating of a weakly absorbing, semitransparent medium with specified conductive heat flux at one boundary. The radiation field is modeled using an approximation to the radiative heat flux obtained by solution of the Kubelka-Munk differential equations.

Weston, K. C.↗

DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗