SEARCH · Engineering Papers

Results for “numerical solution”

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 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Notch Sensitivity of Woven Ceramic Matrix Composites Under Tensile Loading: An Experimental, Analytical, and Finite Element Study

The stress concentrations associated with circular notches and subjected to uniform tensile loading in woven ceramic matrix composites (CMCs) have been investigated for high-efficient turbine engine applications. The CMC's were composed of Nicalon silicon carbide woven fabric in SiNC matrix manufactured through polymer impregnation process (PIP). Several combinations of hole diameter/plate width ratios and ply orientations were considered in this study. In the first part, the stress concentrations were calculated measuring strain distributions surrounding the hole using strain gages at different locations of the specimens during the initial portion of the stress-strain curve before any microdamage developed. The stress concentration was also calculated analytically using Lekhnitskii's solution for orthotropic plates. A finite-width correction factor for anisotropic and orthotropic composite plate was considered. The stress distributions surrounding the circular hole of a CMC's plate were further studied using finite element analysis. Both solid and shell elements were considered. The experimental results were compared with both the analytical and finite element solutions. Extensive optical and scanning electron microscopic examinations were carried out for identifying the fracture behavior and failure mechanisms of both the notched and notched specimens. The stress concentration factors (SCF) determined by analytical method overpredicted the experimental results. But the numerical solution underpredicted the experimental SCF. Stress concentration factors are shown to increase with enlarged hole size and the effects of ply orientations on stress concentration factors are observed to be negligible. In all the cases, the crack initiated at the notch edge and propagated along the width towards the edge of the specimens.

Haque, A.↗

DOE OSTI2021

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗

NASA NTRS1972

Analytical treatment of gas flows through multilayer insulation, project 1

A theoretical investigation of gas flow inside a multilayer insulation system was made for the case of the broadside pumping process. A set of simultaneous first-order differential equations for the temperature and pressure of the gas molecules through the perforations on the insulation layers. A modified Runge-Kutta method was used for numerical experiment. The numerical stability problem was also investigated. It was shown that when the relaxation time is less than the time period over which the gas properties change appreciably, the set of differential equations can be replaced by a set of algebraic equations for solution. Numerical examples were given and comparisons with experimental data were made.

Lin, J. T.↗

NASA NTRS1989

Interaction between Tollmien-Schlichting waves and free-stream disturbances in boundary-layer flows

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small freestream disturbance and a small localized variation of the surface geometry has been carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T-S waves generated by the interaction between the freestream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the Tollmien-Schlichting waves.

Bodonyi, R. J.↗

NASA NTRS1995

Flowfield and acoustic characteristics of telescope cavity in SOFIA platform

Unsteady three-dimensional flowfields are calculated for the Stratospheric Observatory For Infrared Astronomy (SOFIA) at both free-flight cruise and wind tunnel conditions with a view to help in the design process of an acoustically quiet telescope cavity and to understand the flow physics of a three dimensional cavity. The calculation method is based on the numerical solution of thin layer Navier-Stokes equations on a Chimera overset grid system. The Boeing 747-200 aircraft is examined as one option for the SOFIA platform. The flowfield domain is composed of 45 grids consisting of over 4.1 million points. Numerical simulations are performed for both wind tunnel and free-flight cruise conditions at one freestream condition of M(infinity) = 0.85, alpha = 2.5 deg. Comparison of results from wind tunnel simulation show good agreement with experimental data for time-averaged surface pressures, drag for the empennage, and sound pressure levels and power spectra at various locations within the cavity and on the telescope. The presence of the open cavity induces an incremental drag increase, an increased acoustic radiation, and an increase in unsteady pressure loads on the telescope. Its impact on the effectiveness of aircraft control surfaces appears minimal.

Srinivasan, G. R.↗

NASA NTRS1999

Active Control of Flow Separation Over an Airfoil

Designing an aircraft without conventional control surfaces is of interest to aerospace community. In this direction, smart actuator devices such as synthetic jets have been proposed to provide aircraft maneuverability instead of control surfaces. In this article, a numerical study is performed to investigate the effects of unsteady suction and blowing on airfoils. The unsteady suction and blowing is introduced at the leading edge of the airfoil in the form of tangential jet. Numerical solutions are obtained using Reynolds-Averaged viscous compressible Navier-Stokes equations. Unsteady suction and blowing is investigated as a means of separation control to obtain lift on airfoils. The effect of blowing coefficients on lift and drag is investigated. The numerical simulations are compared with experiments from the Tel-Aviv University (TAU). These results indicate that unsteady suction and blowing can be used as a means of separation control to generate lift on airfoils.

Ravindran, S. S.↗

NASA NTRSDate not supplied

Numerical Methods for Jet Noise Predictions Using the Generalized Acoustic Analogy

Numerical methods used in several implementations of the Generalized Acoustic Analogy for the prediction of turbulent jet noise are discussed. The equation governing sound propagation in the analogy formulation is presented along with analytical simplifications to reduce the computational burden of obtaining solutions. Numerical methods for the resulting boundary-value problem for two classes of jet flows, axisymmetric and non-axisymmetric mean flows, are discussed in detail. Example results from problems of the two classed of jets are presented.

jet noise↗

NASA NTRSDate not supplied

Phased Microphone Array on Aircraft Fuselage

This paper studies some challenging features for phased microphone arrays on aircraft fuselage, including the effects of sound scattering such as refraction by nonuniform boundary layer flows, reflection by curved surfaces, and diffraction by sharp edges and smooth geometry. The objective is to show if and how source locations and source amplitudes can be accurately captured by the array data analysis, in the presence of these features. The study starts with canonical problems to separately examine the individual features, making use of their respective analytical solutions, which provide exact results for all quantities needed for the array data processing. To simulate realistic flight environments, a generic aircraft geometry is then utilized, involving multiple scattering mechanisms and multiple noise sources to model aircraft engine noise. For this application, numerical solutions are computed to provide noise propagation and scattering from the sources to the microphones, as well as the quantities needed for array data analysis. The study reveals the effects of various features on fuselage-mounted arrays and provides guidelines to account for these effects in array data analysis.

Yueping Guo↗

DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Toward real-time optimization through model reduction and model discrepancy sensitivities

Optimization problems arise in a range of scenarios, from optimal control to model parameter estimation. In many applications, such as the development of digital twins, it is essential to solve these optimization problems within wall-clock-time limitations. However, this is often unattainable for complex systems, such as those modeled by nonlinear partial differential equations. One strategy for mitigating this issue is to construct a reduced-order model (ROM) that enables more rapid optimization. In particular, the use of nonintrusive ROMs—those that do not require access to the full-order model at evaluation time—is popular because they facilitate the computation of optimization solutions within the wall-clock time requirements. However, the optimization solution will be unreliable if the iterates move outside the ROM training data. This article proposes the use of hyper-differential sensitivity analysis with respect to model discrepancy (HDSA-MD) as a computationally efficient tool to augment ROM-constrained optimization and improve its reliability. The proposed approach consists of two phases: (i) an offline phase where several full-order model evaluations are computed to train the ROM, and (ii) an online phase where a ROM-constrained optimization problem is solved, a limited number of full-order model evaluations are computed, and HDSA-MD is used to enhance the optimization solution. Numerical results are demonstrated for two examples, atmospheric contaminant control and wildfire ignition location estimation, in which a ROM is trained offline using inaccurate atmospheric data. In conclusion, the HDSA-MD update yields a significant improvement in the ROM-constrained optimization solution using only one full-order model evaluation online with corrected atmospheric data.

PDE-constrained optimization↗

NASA NTRS1973

Analysis of gas flow through a multilayer insulation system.

A theoretical investigation of gas flow inside a multilayer insulation system has been made for the case of the broadside pumping process. A set of simultaneous first-order differential equations for the temperature and pressure of the gas mixture was obtained by considering the diffusion mechanism of the gas molecules through the perforations on the insulation layers. A modified Runge-Kutta method was used for numerical experiment. The numerical stability problem was investigated. It has been shown that when the relaxation time is small compared with the time period over which the gas properties change appreciably, the set of differential equations can be replaced by a set of algebraic equations for solution. Numerical examples were given and comparison with experimental data was made.

Lin, J. T.↗

NASA NTRS1973

Thermal engineering research

Shih, C. C.↗

NASA NTRS1980

Radiative transfer theory for passive microwave remote sensing of a two-layer random medium with cylindrical structures

A model of the vegetation layer as a two-layer random medium with a small correlation length l sub rho in the horizontal direction, and a large correlation length l sub z in the vertical direction, is presented for fields with cylindrical structures. As l sub z approaches infinity, closed form solutions are derived for the brightness temperatures; the kernels in the scattering terms of the radiative transfer equations result in delta functions showing that forward scattering is dominant over all other directions. The results are compared with the Gaussian quadrature method data for numerical solution of the radiative transfer equations.

Chuang, S. L.↗

NASA NTRS1990

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

NASA NTRS1991

The free compressible viscous vortex

The present study investigates the effects of compressibility on free (unsteady) viscous heat-conducting vortices. Analytical solutions are found in the limit of large but finite Reynolds number and small but finite Mach number. It is shown that the spreading of the vortex causes a radial flow. This flow is given by the solution of an ordinary differential equation, which gives the dependence of the radial velocity on the tangential velocity, density, and temperature profiles of the vortex. Estimates of the radial velocity found by solving this equation are found to be in good agreement with numerical solutions of the full equations. The equations for the viscous evolution are expanded in powers of Mach number to obtain detailed analytical solutions. It is shown that swirling axisymmetric compressible flows generate negative radial velocities far from the vortex core owing to viscous effects, regardless of the initial distributions of vorticity, density, and entropy.

Colonius, Tim↗

NASA NTRS1992

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

Bodonyi, R. J.↗

NASA NTRS1994

A note on the accuracy of spectral method applied to nonlinear conservation laws

Fourier spectral method can achieve exponential accuracy both on the approximation level and for solving partial differential equations if the solutions are analytic. For a linear partial differential equation with a discontinuous solution, Fourier spectral method produces poor point-wise accuracy without post-processing, but still maintains exponential accuracy for all moments against analytic functions. In this note we assess the accuracy of Fourier spectral method applied to nonlinear conservation laws through a numerical case study. We find that the moments with respect to analytic functions are no longer very accurate. However the numerical solution does contain accurate information which can be extracted by a post-processing based on Gegenbauer polynomials.

Shu, Chi-Wang↗

NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗