Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Numerical evaluation of the Rayleigh integral for planar radiators using the FFT

Rayleigh's integral formula is evaluated numerically for planar radiators of any shape, with any specified velocity in the source plane using the fast Fourier transfrom algorithm. The major advantage of this technique is its speed of computation - over 400 times faster than a straightforward two-dimensional numerical integration. The technique is developed for computation of the radiated pressure in the nearfield of the source and can be easily extended to provide, with little computation time, the vector intensity in the nearfield. Computations with the FFT of the nearfield pressure of baffled rectangular plates with clamped and free boundaries are compared with the 'exact' solution to illuminate any errors. The bias errors, introduced by the FFT, are investigated and a technique is developed to significantly reduce them.

Williams, E. G.↗
NASA NTRSDate not supplied

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗
NASA NTRS1996

Artificial Boundary Conditions Based on the Difference Potentials Method

While numerically solving a problem initially formulated on an unbounded domain, one typically truncates this domain, which necessitates setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The issue of setting the ABC's appears to be most significant in many areas of scientific computing, for example, in problems originating from acoustics, electrodynamics, solid mechanics, and fluid dynamics. In particular, in computational fluid dynamics (where external problems present a wide class of practically important formulations) the proper treatment of external boundaries may have a profound impact on the overall quality and performance of numerical algorithms. Most of the currently used techniques for setting the ABC's can basically be classified into two groups. The methods from the first group (global ABC's) usually provide high accuracy and robustness of the numerical procedure but often appear to be fairly cumbersome and (computationally) expensive. The methods from the second group (local ABC's) are, as a rule, algorithmically simple, numerically cheap, and geometrically universal; however, they usually lack accuracy of computations. In this paper we first present a survey and provide a comparative assessment of different existing methods for constructing the ABC's. Then, we describe a relatively new ABC's technique of ours and review the corresponding results. This new technique, in our opinion, is currently one of the most promising in the field. It enables one to construct such ABC's that combine the advantages relevant to the two aforementioned classes of existing methods. Our approach is based on application of the difference potentials method attributable to V. S. Ryaben'kii. This approach allows us to obtain highly accurate ABC's in the form of certain (nonlocal) boundary operator equations. The operators involved are analogous to the pseudodifferential boundary projections first introduced by A. P. Calderon and then also studied by R. T. Seeley. The apparatus of the boundary pseudodifferential equations, which has formerly been used mostly in the qualitative theory of integral equations and PDE'S, is now effectively employed for developing numerical methods in the different fields of scientific computing.

Tsynkov, Semyon V.↗
DOE OSTI2021

Phase-weighted slant stacking for surface wave dispersion measurement

SUMMARY Surface wave retrieval from ambient noise records using seismic interferometry techniques has been widely used for multiscale shear wave velocity (Vs) imaging. One key step during Vs imaging is the generation of dispersion spectra and the extraction of a reliable dispersion curve from the retrieved surface waves. However, the sparse array geometry usually affects the ability for high-frequency (>1 Hz) seismic signals’ acquisition. Dispersion measurements are degraded by array response due to sparse sampling and often present smeared dispersion spectra with sidelobe artefacts. Previous studies usually focus on interferograms’ domain (e.g. cross-correlation function) and attempt to enhance coherent signals before dispersion measurement. We propose an alternative technique to explicitly deblur dispersion spectra through use of a phase-weighted slant-stacking algorithm. Numerical examples demonstrate the strength of the proposed technique to attenuate array responses as well as incoherent noise. Three different field examples prove the flexibility and superiority of the proposed technique: the first data set consists of ambient noise records acquired using a nodal seismometer array; the second data set utilizes distributed acoustic sensing (DAS) and a marine fibre-optic cable to acquire a similar ambient noise data set; the last data set is a vibrator-based active-source surface wave data. The enhanced dispersion measurements provide cleaner and higher-resolution spectra without distortions which will assist both human interpreters as well as ML algorithms in efficiently picking curves for subsequent Vs inversion.

58 GEOSCIENCES↗
NASA NTRS1987

A finite volume Euler calculation of the aerodynamics of transonic airfoil-vortex interaction

Unsteady inviscid transonic airfoil-vortex interaction is numerically analyzed by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables is used to form the numerical dissipation terms to stabilize the scheme. A velocity decomposition technique is applied to alleviate the problem of vortex diffusion by the numerical dissipation terms and to treat the interaction of a Rankine vortex with an airfoil accurately. Results obtained are compared with available numerical data.

Damodaran, Murali↗
NASA NTRS1990

Microwave holographic metrology for antenna diagnosis

Advances in antenna diagnostic methodologies have been very significant in recent years. In particular, microwave holographic diagnostic techniques have been applied very successfully in improving the performance of reflector and array antennas. These techniques use the knowledge of the measured amplitude and phase of the antenna radiated fields and then take advantage of the existing Fourier transform relationships between the radiated fields and the effective aperture or current distribution to eventually determine the reflector surface or array excitation coefficients anomalies. In this paper an overview of the recent developments in applying microwave holography is presented. The theoretical, numerical and measurement aspects of this technique is detailed by providing representative results.

Rahmat-Samii, Y.↗
NASA NTRS1992

Homogeneous quantum electrodynamic turbulence

The electromagnetic field equations and Dirac equations for oppositely charged wave functions are numerically time-integrated using a spatial Fourier method. The numerical approach used, a spectral transform technique, is based on a continuum representation of physical space. The coupled classical field equations contain a dimensionless parameter which sets the strength of the nonlinear interaction (as the parameter increases, interaction volume decreases). For a parameter value of unity, highly nonlinear behavior in the time-evolution of an individual wave function, analogous to ideal fluid turbulence, is observed. In the truncated Fourier representation which is numerically implemented here, the quantum turbulence is homogeneous but anisotropic and manifests itself in the nonlinear evolution of equilibrium modal spatial spectra for the probability density of each particle and also for the electromagnetic energy density. The results show that nonlinearly interacting fermionic wave functions quickly approach a multi-mode, dynamic equilibrium state, and that this state can be determined by numerical means.

Shebalin, John V.↗
DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Phase-space Energization of Ions in Oblique Shocks

Abstract Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique’s utility for simulation, laboratory, and spacecraft analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

Phase-Space Energization of Ions in Oblique Shocks

Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique's utility for simulation, laboratory, and spacecraft analysis.

Plasma astrophysics↗
NASA NTRS1986

Optimal control of quaternion propagation errors in spacecraft navigation

Optimal control techniques are used to drive the numerical error (truncation, roundoff, commutation) in computing the quaternion vector to zero. The normalization of the quaternion is carried out by appropriate choice of a performance index, which can be optimized. The error equations are derived from Friedland's (1978) theoretical development, and a matrix Riccati equation results for the computation of the gain matrix. Simulation results show that a high precision of the order of 10 to the -12th can be obtained using this technique in meeting the q(T)q=1 constraint. The performance of the estimator in the presence of the feedback control that maintains the normalization, is studied.

Vathsal, S.↗
DOE OSTI2020

Many-Body Effects in FeN4 Center Embedded in Graphene

We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Advances and trends in computational structural mechanics

The development status and applicational range of techniques in computational structural mechanics (CSM) are evaluated with a view to advances in computational models for material behavior, discrete-element technology, quality assessment, the control of numerical simulations of structural response, hybrid analysis techniques, techniques for large-scale optimization, and the impact of new computing systems on CSM. Primary pacers of CSM development encompass prediction and analysis of novel materials for structural components, computational strategies for large-scale structural calculations, and the assessment of response prediction reliability together with its adaptive improvement.

Noor, Ahmed K.↗
NASA NTRS1990

Upwind techniques for flows with multiple translational temperatures

The numerical simulation of high-speed flows out of chemical and thermodynamic equilibrium is considered. A general model for the species internal energy is utilized, which allows for multiple translational temperatures. A study of the characteristic wave speeds for the case of two translational temperatures, one associated with the average thermal speed of the heavy particles and the other associated with the translational temperature of free electrons, yields some very interesting results, with acoustic wave speeds of the symmetric form u + or - a no longer appearing. An asymptotic analysis in terms of small electron mass fractions yields approximate wave speed which recover the symmetric form, with the addition of a new wave speed associated with the electron energy equation.

Cinnella, Pasquale↗
NASA NTRS1993

An evaluation of excitation techniques for time domain based flutter data processing

Numerically simulated, wind tunnel and flight test experimental data were used to evaluate various types of system excitation sources for subcritical flutter testing. Emphasis was placed upon the determination of modal frequencies and damping from relatively short time series records. Parameter identification based upon digital time series models was used to predict modal characteristics from subcritical test data. Various types of excitation were considered. These were: 1) impulsive loading and the resulting transient free vibration response, 2) random forced response, in which the 'mechanical' forcing function could be directly measured, 3) sine sweep forced response, in which the 'mechanical' forcing function could be directly measured, and 4) forced response from unsteady aerodynamic excitation or turbulence, in which the excitation could only be indirectly inferred from other measurements.

Pinkelman, J. K.↗
DOE OSTI2020

Comparison of electro-thermal performance of advanced cooling techniques for electric vehicle motors

Permanent magnet synchronous motors are commonly used in the powertrain of electric vehicles because of their high power and torque densities. Attempts to increase power and torque densities beyond the state-of-the-art often suffer from thermal limitations of the adopted winding-wire insulation class. In this paper, overall effectiveness of four different cooling technologies, namely conventional stator jacket cooling, embedded circular and rectangular cooling channels within stator core, and direct winding heat exchangers have been studied numerically. Here, by realizing temperature dependent magnetic and material properties, winding and core losses, a two-way coupling algorithm has been utilized to study the detailed electro-magnetic and thermal performance of BMW i3 motor with aforementioned cooling techniques. The numerical results illustrated that the direct winding heat exchangers approach provides better electro-thermal performance in comparison to the other cooling techniques analyzed. On the contrary, embedded rectangular and circular channel in stator core adversely affect the electro-magnetic performance and overall efficiency of the motor, despite their superior thermal performance compared to the jacket cooling

42 ENGINEERING↗
NASA NTRS2000

Magnetic Susceptibility Effects and Lorentz Damping in Diamagnetic Fluids

A great number of crystals (semi-conductor and protein) grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity and g-jitter. Both static and dynamic (rotating or travelling wave) magnetic fields can be used to reduce the effects of convection in materials processing. In semi-conductor melts, due to their relatively high electrical conductivity, the induced Lorentz force can be effectively used to curtail convective effects. In melts/solutions with reduced electrical conductivity, such as aqueous solutions used in solution crystal growth, protein crystal growth and/or model fluid experiments for simulating melt growth, however, the variation of the magnetic susceptibility with temperature and/or concentration can be utilized to better damp fluid convection than the Lorentz force method. This paper presents a comprehensive, comparative numerical study of the relative damping effects using static magnetic fields and gradients in a simple geometry subjected to a thermal gradient. The governing equations are formulated in general terms and then simplified for the numerical calculations. Operational regimes, based on the best damping technique for different melts/solutions are identified based on fluid properties. Comparisons are provided between the numerical results and available results from experiments in surveyed literature.

Ramachandran, Narayanan↗
NASA NTRS1982

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗