SEARCH · Engineering Papers

Results for “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 757 records · Page 42

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Modal Testing of a Flexible Wing on a Dynamically Active Test Fixture Using Fixed Base Correction Method - IFASD 2019

In modal testing and finite element model correlation, analysts desire modal results using free-free or rigid boundary conditions to ease comparisons of test versus analytical data. It is often expensive both in cost and schedule to build and test with boundary conditions that replicate the free-free or rigid boundaries. Static test fixtures for load testing are often large, heavy, and unyielding, but do not provide adequate boundaries for modal tests because they are dynamically too flexible and often contain natural frequencies within the test article frequency range of interest. The dynamic coupling between the test article and test fixture complicates the model updating process because significant effort needs to be spent on modeling the test fixture and boundary conditions in addition to the test article. If the modal results could be corrected for fixture coupling, then setups used for other structural testing could be adequate for modal testing and would allow significant schedule and cost savings by eliminating a unique setup for only modal testing. To simplify future modal tests, this paper reports on a fixed base correction method that was investigated during modal testing of a full-scale, half-span, flexible wing cantilevered from a static test fixture. The results of this fixed base correction approach look very promising. The method aided in producing similar wing modal characteristics for two different physical boundary configurations of a dynamically active test fixture.

modal testing↗

NASA NTRS1979

A fast Cauchy-Riemann solver

The inhomogeneous Cauchy-Riemann equations in a rectangle are discretized by a finite difference approximation. Several different boundary conditions are treated explicitly, leading to algorithms which have overall second-order accuracy. All boundary conditions with either u or v prescribed along a side of the rectangle can be treated by similar methods. The algorithms presented here have nearly minimal time and storage requirements and seem suitable for development into a general-purpose direct Cauchy-Riemann solver for arbitrary boundary conditions.

Ghil, M.↗

DOE OSTI2023

A Medium-/Low-Voltage Joint State Estimator Through Linear Uncertainty Propagation

Traditionally, distribution system state estimations (DSSE) are challenged by the lack of measurements at both primary and secondary sides of the system. The widely available cable television (CATV) voltage sensors installed in low-voltage (LV) networks bring opportunities to achieve higher quality DSSE covering a broader area of the distribution network. This study proposes a medium-/low-voltage (MV/LV) joint distribution system state estimation approach using the untapped CATV measurements. It aims at addressing the need for system situational awareness at the grid edge while improving the estimation accuracy at both the primary and secondary sides compared to its disjointed counterpart. Linearized measurement functions and boundary condition uncertainty propagation rules are derived to ensure the computational efficiency and accuracy of the joint state estimator. Numerical experiments are conducted on an IEEE test feeder to demonstrate the efficacy of the proposed method and the value of CATV measurements.

joint state estimation↗

NASA NTRS1993

A novel pulse technique for thermal diffusivity determination of high-temperature levitated materials - Inverse analysis and experimental set-up

In order to determine the thermal diffusivity of materials at high temperatures, a two-step extended flash technique, which is applicable to levitated spherical samples, is proposed. The containerless flash method is modeled as an axisymmetric transient conduction heat transfer problem within the sphere subjected to a radiative boundary condition on its surface. The problem of nonlinearity which arises from the radiative heat transfer boundary condition is solved by replacing it with the measured time-dependent surface temperature data, thus giving rise to an 'Inverse Analysis'. Upon obtaining the analytic solution for the temperature field, the determination of the thermal diffusivity turns into a minimization problem, whereby the difference between the actual surface boundary condition and the analytic solution is minimized. In performing the proposed experiments, there is a need to undertake a cool-down experiment which should be governed by the lumped parameter analysis.

Shen, F.↗

NASA NTRS2020

Passive Aeroelastic Tailored Wing Modal Test Using the Fixed Base Correction Method

In modal testing and finite element model correlation, analysts desire modal results using free-free or rigid boundary conditions to ease comparisons of test versus analytical data. It is often expensive both in cost and schedule to build and test with boundary conditions that replicate the free-free or rigid boundaries. Static test fixtures for load testing are often large, heavy, and unyielding, but do not provide adequate boundaries for modal tests because they are dynamically too flexible and often contain natural frequencies within the frequency range of interest of the test article. Dynamic coupling between the test article and test fixture complicates the model updating process because significant effort is required to model the test fixture and boundary conditions in addition to the test article. If there were a way to correct the modal results for fixture coupling, then setups used for other structural testing could be adequate for modal testing. In the case described in this paper, a partial static loads testing setup was used, which allowed significant schedule and cost savings by eliminating a unique setup for a modal test. A fixed base correction technique was investigated during modal testing of a flexible wing cantilevered from part of a static test fixture. The technique was successfully used to measure the wing modes de-coupled from the dynamically active test fixture. The technique is promising for future aircraft applications, but more research is needed.

structural dynamics↗

NASA NTRS1989

Further comparison of two synoptic surface wind and pressure analysis methods

The geostrophic stream and potential functions of Sangster and the derived flat pressure field of Pielke and Cram are further compared. A simple numerical experiment with an idealized mountain-atmosphere system compares the two methods and their sensitivity to lateral boundary conditions. The flat pressure field of Pielke and Cram is essentially equivalent to the stream-function of Sangster, although the two methods have different boundary condition formulations and resulting different sensitivities to these boundary conditions.

Cram, Jennifer M.↗

NASA NTRS1998

The Origin of Antarctic Precipitation: A Modeling Approach

Isotope concentrations in polar ice cores have long been used to estimate paleotemperatures. Underlying the use of this "isotope paleothermometer" is the assumption that the relationship between surface temperature and isotope concentration over time at a single geographical point is the same as that observed over space during the present-day climate. The validity of this assumption may in fact be compromised by several factors related to climate change. The specific factor studied in this paper involves the evaporative sources for polar precipitation. Climatic changes in the relative strengths of these sources would imply a need for a recalibration of the paleothermometer. To quantify such changes, we performed two GCM simulations, one of present-day climate and the other of the climate during the Last Glacial Maximum (LGM), roughly 18000 years ago. Evaporative sources of Antarctic precipitation were established using special tracer diagnostics. Results suggest that polar precipitation during the LGM does indeed consist of (relatively) more water from tropical oceans, a direct reflection of the LGM's increased equator-to-pole temperature gradient and its increased sea ice extent, which reduces high latitude evaporation. This result implies that an uncalibrated ice core paleothermometer would produce LGM temperatures that are biased slightly low. Because LGM boundary conditions are still under debate, we performed a third GCM simulation using a modified set of LGM boundary conditions. Using this simulation gives some qualitatively similar results, though the tropical contribution is not quite as high. Uncertainties in the LGM boundary conditions does hamper success in calibrating the paleothermometer.

Delaygue, Gilles↗

NASA NTRS2014

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗

NASA NTRS1990

A critical evaluation of two-equation models for near wall turbulence

A variety of two-equation turbulence models,including several versions of the K-epsilon model as well as the K-omega model, are analyzed critically for near wall turbulent flows from a theoretical and computational standpoint. It is shown that the K-epsilon model has two major problems associated with it: the lack of natural boundary conditions for the dissipation rate and the appearance of higher-order correlations in the balance of terms for the dissipation rate at the wall. In so far as the former problem is concerned, either physically inconsistent boundary conditions have been used or the boundary conditions for the dissipation rate have been tied to higher-order derivatives of the turbulent kinetic energy which leads to numerical stiffness. The K-omega model can alleviate these problems since the asymptotic behavior of omega is known in more detail and since its near wall balance involves only exact viscous terms. However, the modeled form of the omega equation that is used in the literature is incomplete-an exact viscous term is missing which causes the model to behave in an asymptotically inconsistent manner. By including this viscous term and by introducing new wall damping functions with improved asymptotic behavior, a new K-tau model (where tau is identical with 1/omega is turbulent time scale) is developed. It is demonstrated that this new model is computationally robust and yields improved predictions for turbulent boundary layers.

Speziale, Charles G.↗

NASA NTRS1994

A mathematical model for simulating noise suppression of lined ejectors

A mathematical model containing the essential features embodied in the noise suppression of lined ejectors is presented. Although some simplification of the physics is necessary to render the model mathematically tractable, the current model is the most versatile and technologically advanced at the current time. A system of linearized equations and the boundary conditions governing the sound field are derived starting from the equations of fluid dynamics. A nonreflecting boundary condition is developed. In view of the complex nature of the equations, a parametric study requires the use of numerical techniques and modern computers. A finite element algorithm that solves the differential equations coupled with the boundary condition is then introduced. The numerical method results in a matrix equation with several hundred thousand degrees of freedom that is solved efficiently on a supercomputer. The model is validated by comparing results either with exact solutions or with approximate solutions from other works. In each case, excellent correlations are obtained. The usefulness of the model as an optimization tool and the importance of variable impedance liners as a mechanism for achieving broadband suppression within a lined ejector are demonstrated.

Watson, Willie R.↗

NASA NTRS1977

Computer program for calculating aerodynamic characteristics of upper-surface-blowing and over-wing-blowing configurations

The program is based on the inviscid wing-jet interaction theory of Lan and Campbell, and the jet entrainment theory of Lan. In the interaction theory, the flow perturbations are computed both inside and outside the jet, separately, and then matched on the jet surface to satisfy the jet boundary conditions. The jet Mach number is allowed to be different from the free stream value (Mach number nonuniformity). These jet boundary conditions require that the static pressure be continuous across the jet surface which must always remain as a stream surface. These conditions, as well as the wing-surface tangency condition, are satisified only in the linearized sense. The detailed formulation of these boundary conditions is based on the quasi-vortex-lattice method of Lan.

Lan, C. E.↗

NASA NTRS2018

Simulation of Sound Absorption by Scattering Bodies Treated with Acoustic Liners Using a Time-Domain Boundary Element Method

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction, and are characterized by a complex valued frequency-dependent impedance, Z(w). Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering of geometric bodies treated with acoustic liners. This work uses an admittance boundary condition where the admittance, Y(w), is defined to be the inverse of impedance, i.e., Y(w) = 1/Z(w). An admittance boundary condition will be derived and coupled with a time domain boundary integral equation. The solution will be obtained iteratively using spatial and temporal basis functions and will allow for acoustic scattering problems to be modeled with geometries consisting of both unlined and soft surfaces. Stability will be demonstrated through eigenvalue analysis.

Pizzo, Michelle E.↗

NASA NTRS1992

Numerical analysis of hypersonic turbulent film cooling flows

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

Chen, Y. S.↗

DOE OSTI2024

Proton deflectometry analysis in magnetized plasmas: Magnetic field reconstruction in one dimension

Proton deflectometry is used in magnetized high-energy-density plasmas to observe electromagnetic fields. Here we describe a reconstruction algorithm to recover the electromagnetic fields from proton fluence data in 1-D. The algorithm is verified against analytic solutions and applied to example data. Next, we study the role of source fluence uncertainty for 1-D reconstructions. We show that reconstruction boundary conditions can be used to constrain the source fluence profile and use this to develop a reconstruction using a specified pair of boundary conditions on the magnetic field. From these considerations, we experimentally demonstrate a hybrid mesh-fluence reconstruction technique where fields are reconstructed from fluence data in an interior region with boundary conditions supplied by direct mesh measurements at the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2020

Hall2De Simulations of a Magnetic Nozzle

The two-dimensional, fluid-based Hall thruster code Hall2De is adapted to simulate a low temperature propulsive magnetic nozzle. Electron cyclotron resonance (ECR) thrusters developed at the Office National d’Etudes et de Recherches Aérospatiales (ONERA) in France and the University of Michigan are investigated. For each device, the simulated domain is defined from the nozzle throat to a location 10 thruster radii downstream. A numerical mesh is defined based on the measured magnetic field, and the governing fluid equations for singlycharge ions, electrons, and neutrals are solved. The boundary conditions for the throat of the magnetic nozzle are based on measured values of electron temperature and ion velocity. Multiple boundary conditions for the electron temperature in the downstream regions of the domain are considered. The simulated ion velocity and potential agrees well with experiment, while the simulated divergence efficiencies and thrusts are 25% lower and 45-55% higher than experiment. When a Neumann boundary condition for electron temperature is applied to the far-plume boundaries, the nozzle was seen to be purely isothermal. This suggests that the thermal conductivity along field lines in real devices is much lower than classical simulations predict, and anomalous resistivity may be needed to reproduce the temperature distribution observed in experiment. The larger beam divergence observed in our simulation may also be related to the lack of anomalous resistivity.

Jorns, Benjamin A.↗

NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗

NASA NTRS1991

Unsteady-flow-field predictions for oscillating cascades

The unsteady flow field around an oscillating cascade of flat plates with zero stagger was studied by using a time marching Euler code. This case had an exact solution based on linear theory and served as a model problem for studying pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step size was shown for a moderate reduced frequency. Results show that an approximate nonreflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer to the airfoils than when reflective boundaries are used. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates captures the unsteady flow field better than when uniform grids are used as long as the 'Courant Friedrichs Levy' (CFL) number is less than 1 for a sufficient portion of the grid. Finally, a solution based on an optimization of grid, CFL number, and boundary conditions shows good agreement with linear theory.

Huff, Dennis L.↗

NASA NTRS1993

Matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation for aeroassisted plane-change maneuvers

In this paper we develop a general procedure for constructing a matched asymptotic expansion of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is inappropriate since it is not uniformly valid over a narrow range of the independent variable. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished - one where the left boundary condition coincides with or lies to the right of the singular region and one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure, and its potential application to aeroassisted plane change is described.

Calise, Anthony J.↗