Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate 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.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Propagation and Breaking at High Altitudes of Gravity Waves Excited by Tropospheric Forcing

An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

Prusa, Joseph M.↗
NASA NTRS1995

Aerothermodynamic Measurement and Prediction for Modified Orbiter at Mach 6 and 10

Detailed heat-transfer rate distributions measured laterally over the windward surface of an orbiter-like configuration using thin-film resistance heat-transfer gauges and globally using the newly developed relative intensity, two-color thermographic phosphor technique are presented for Mach 6 and 10 in air. The angle of attack was varied from 0 to 40 deg, and the freestream Reynolds number based on the model length was varied from 4 x 10(exp 5) to 6 x 10(exp 6) at Mach 6, corresponding to laminar, transitional, and turbulent boundary layers; the Reynolds number at Mach 10 was 4 x 10(exp 5), corresponding to laminar flow. The primary objective of the present study was to provide detailed benchmark heat-transfer data for the calibration of computational fluid-dynamics codes. Predictions from a Navier-Stokes solver referred to as the Langley aerothermodynamic upwind relaxation algorithm and an approximate boundary-layer solving method known as the axisymmetric analog three-dimensional boundary layer code are compared with measurement. In general, predicted laminar heat-transfer rates are in good agreement with measurements.

Micol, John R.↗
NASA NTRS1994

Coupled electromechanical response of composite beams with embedded piezoelectric sensors and actuators

Unified mechanics are developed with the capability to model both sensory and active composite laminates with embedded piezoelectric layers. A discrete-layer formulation enables analysis of both global and local electromechanical response. The mechanics include the contributions from elastic, piezoelectric, and dielectric components. The incorporation of electric potential into the state variables permits representation of general electromechanical boundary conditions. Approximate finite element solutions for the static and free-vibration analysis of beams are presented. Applications on composite beams demonstrate the capability to represent either sensory or active structures and to model the complicated stress-strain fields, the interactions between passive/active layers, interfacial phenomena between sensors and composite plies, and critical damage modes in the material. The capability to predict the dynamic characteristics under various electrical boundary conditions is also demonstrated.

Saravanos, D. A.↗
NASA NTRS1994

Self-consistent steady state model of the low-latitude boundary layer

Results are presented from a steady state numerical model of the low-latitude boundary layer (LLBL) on closed field lines and its coupling to the dayside auroral ionosphere. In the model the boundary layer approximation is used, the result being that inertial forces are taken into account only in the main flow direction (-x) where they are balanced by pressure forces, j x B forces, and viscous forces. Motion in the transverse directions (y and z) is treated kinematically, the force balances in these two directions being purely static. Computationally, the model is two dimensional, describing the motion of plasma and frozen-in magnetic field in the equatorial (xy) plane but allowing for lowest-order polynomial variation of some quantities with the coordinate (z) perpendicular to that plane. The plasma expands and compresses isentropically; the magnetic field is calculated self-consistently, which leads to approximately parabolic field line shape in planes parallel to the magnetopause (the xz plane), with maximum field curvature near the magnetopause edge of the LLBL. Coupling to the ionosphere via region 1 field-aligned currents is included. The effects of the ionosphere are represented by two parallel resistive plates at fixed height above and below the equatorial plane. The model can be used to investigate the influence of various physical parameters, for example, viscous and magnetic Reynolds numbers, and of boundary conditions at the magnetopause and in the magnetosphere on the LLBL development in the -x direction. Special attention is given to viscous effects which, under suitable circumstances, lead to a region 1 current that first increases and then reduces with increasing longitude away from local noon. Asymptotic matching of the antisunward motion of the cool LLBL plasma to sunward convection of hot plasma in the magnetosphere is illustrated along with the entrainment of magnetospheric plasma by the antisunward LLBL flow.

Drakou, E.↗
NASA NTRS1979

Perturbation solutions of combustion instability problems

A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

Googerdy, A.↗
NASA NTRS1988

Propagation of quasiplane waves along an impedance boundary

The parabolic approximation for the acoustic equations of motion is applied to the study of the sound field generated by a time harmonic plane wave at grazing incidence to a finite impedance boundary. The resulting equations possess a solution which may be expressed in terms of the complementary error function. Asymptotic expansion of this solution for field points near the boundary provides results compatible with those for a point source on the boundary for both the soft boundary (finite impedance) and hard boundary (the limit in which the impedance becomes infinite) cases. The presence of a surface wave in the solution is also established.

Mcaninch, G. L.↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High-throughput characterization of Lu-doped zirconia

Lu-doped zirconia (LuDZ) was reported in the 1960s to be a poor oxygen ion conductor. However, we identified this system as worthy of reexamination based on the results of a machine learning-based screen for promising oxygen ion conductors, further supported by a physical argument: because Lu is the heaviest and smallest of the lanthanides, it should offer the possibility of high conductivity due to its high polarizability and size compatibility with the zirconia structure. Here the relationship between dopant concentration and conductivity in the LuDZ system was determined using a high throughput experimental methodology in which spatially resolved conductivity was measured across a compositionally graded film, spanning from Zr 0.93 Lu 0.07 O 1.965 to Zr 0.82 Lu 0.18 O 1.91 . The conductivity broadly peaks with composition at 15 cation % Lu, coinciding with the approximate phase boundary of the fully stabilized cubic phase. The composition-dependent conductivity measurements reveal a surprising increase in activation energy for ion transport with increasing Lu content with no change in this trend at the composition of the conductivity maximum. The pre-exponential factor increases sharply with Lu content, but displays a distinct plateau once the conductivity maximum composition is attained. At 500 °C, the conductivity of Zr 0.85 Lu 0.15 O 1.925 is 2 × 10 –3 S cm –1 , approaching that of optimally doped Sc-stabilized zirconia, which has the highest reported conductivity among zirconia compounds. In comparison to other oxide ion conductors, Lu-stabilized zirconia offers the advantage over ceria of negligible electronic conductivity and is easier to process than (La,Sr)(Ga,Mg)O 3-δ .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Self-supervised physics-informed generative networks for phase retrieval from a single X-ray hologram

X-ray phase contrast imaging significantly improves the visualization of structures with weak or uniform absorption, broadening its applications across a wide range of scientific disciplines. Propagation-based phase contrast is particularly suitable for time- or dose-critical in vivo/in situ/operando (tomography) experiments because it requires only a single intensity measurement. However, the phase information of the wave field is lost during the measurement and must be recovered. Conventional algebraic and iterative methods often rely on specific approximations or boundary conditions that may not be met by many samples or experimental setups. In addition, they require manual tuning of reconstruction parameters by experts, making them less adaptable for complex or variable conditions. Here we present a self-learning approach for solving the inverse problem of phase retrieval in the near-field regime of Fresnel theory using a single intensity measurement (hologram). A physics-informed generative adversarial network is employed to reconstruct both the phase and absorbance of the unpropagated wave field in the sample plane from a single hologram. Unlike most state-of-the-art deep learning approaches for phase retrieval, our approach does not require paired, unpaired, or simulated training data. This significantly broadens the applicability of our approach, as acquiring or generating suitable training data remains a major challenge due to the wide variability in sample types and experimental configurations. The algorithm demonstrates robust and consistent performance across diverse imaging conditions and sample types, delivering quantitative, high-quality reconstructions for both simulated data and experimental datasets acquired at beamline P05 at PETRA III (DESY, Hamburg), operated by Helmholtz-Zentrum Hereon. Furthermore, it enables the simultaneous retrieval of both phase and absorption information.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Prediction of far flow field in trailing vortices.

A finite-difference machine code is brought to bear on the wake-vortex problem in the quasi-cylindrical boundary-layer approximation. A turbulent-energy model containing new features is developed. Parameters of the model are evaluated by comparison of calculated velocities and turbulent intensities with measurements in an axisymmetric wake. Comparisons are made with a previous calculation of the decay of an isolated vortex and with wind tunnel and flight measurements in trailing vortices. A self-similar solution develops at large axial distance that decays with the square root of distance. A slower decay occurs in the preceding transition region.

Baldwin, B. S.↗
NASA NTRS1973

General atmospheric circulation driven by polar and diurnal surface temperature variations.

Described is a global circulation model for the Venus atmosphere that includes the effects of both polar cooling and diurnal temperature variation. It is based on a linearized Boussinesq approximation and boundary conditions derived from theoretical and empirical considerations. The time-dependent, three-dimensional flow field is deduced without any a priori assumptions about its configuration. Results show that the mean atmospheric motions are essentially zonal in a narrow belt near the equator and change to become meridional over most of the globe. The circulation pattern is not symmetrical and rotates about the polar axis of the planet with the period of the solar day.

Bohachevsky, I. O.↗
NASA NTRS1974

Decay of far flow field in trailing vortices

A finite difference machine code is used in the wake vortex problems in the quasi-cylindrical boundary layer approximation. A turbulent energy model containing new features is developed that accounts for the major effects disclosed by more advanced models in which the parameters are not yet established. Several puzzles that arose in previous theoretical investigations of wake vortices are resolved.

Baldwin, B. S.↗
NASA NTRS1975

Investigation B: Laminar confined coaxial entrance flow with heat generation

The results of a parametric study on the entrance flow region in a gas core nuclear reactor are presented. The physical system is modeled as laminar confined, coaxial flow with heat generation in the inner fluid. The governing equations include the boundary layer approximations and the assumptions of only radial radiative transport of energy represented as an energy diffusion term. The Von Mises transformation and a zeta transformation are used to transform the equations into nonlinear nonhomogeneous convective-diffusion equations. A unique combination of forward and backward difference equations which yields accurate results at moderate computational times, is used in the numerical method. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Bobba, G. K. M.↗