Search NASA⌕ Search

SEARCH · Search NASA

Results for “convention equation”

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 37 records · Page 2

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

Optical antenna gain. I - Transmitting antennas

The gain of centrally obscured optical transmitting antennas is analyzed in detail. The calculations, resulting in near- and far-field antenna gain patterns, assume a circular antenna illuminated by a laser operating in the TEM-00 mode. A simple polynomial equation is derived for matching the incident source distribution to a general antenna configuration for maximum on-axis gain. An interpretation of the resultant gain curves allows a number of auxiliary design curves to be drawn that display the losses in antenna gain due to pointing errors and the cone angle of the beam in the far field as a function of antenna aperture size and its central obscuration. The results are presented in a series of graphs that allow the rapid and accurate evaluation of the antenna gain which may then be substituted into the conventional range equation.

Klein, B. J.↗

Non-linear duct acoustics and its application to fan noise

Quite often the scalar Helmholtz equation is assumed to be the fundamental equation of duct acoustics and various relations are derived from it. An investigation is, therefore, conducted regarding the underlying assumptions leading to the Helmholtz equation. It is found that serious errors, even well below 186 dB, are possible if the assumption of linearity were to be made. In duct acoustics these errors are highest near the cutoff points and in the case of ducts with transonic flow. It is pointed out that serious errors could arise if a nonlinear problem is analyzed by linear methods. At times, an entire phenomenon remains unpredictable by linear methods. Attention is given to the characteristics of nonlinear systems, the phenomenon of multiple pure tones in the case of noise from the aircraft engine fans, the limitations of the conventional acoustics equation, duct shape and nonlinearity, and the numerical simulation of a fan with 8 blades.

Vaidya, P. G.↗

Simple phenomenological modeling of transition-region capacitance of forward-biased p-n junction diodes and transistor diodes

The derivation of a simple expression for the capacitance C(V) associated with the transition region of a p-n junction under a forward bias is derived by phenomenological reasoning. The treatment of C(V) is based on the conventional Shockley equations, and simpler expressions for C(V) result that are in general accord with the previous analytical and numerical results. C(V) consists of two components resulting from changes in majority carrier concentration and from free hole and electron accumulation in the space-charge region. The space-charge region is conceived as the intrinsic region of an n-i-p structure for a space-charge region markedly wider than the extrinsic Debye lengths at its edges. This region is excited in the sense that the forward bias creates hole and electron densities orders of magnitude larger than those in equilibrium. The recent Shirts-Gordon (1979) modeling of the space-charge region using a dielectric response function is contrasted with the more conventional Schottky-Shockley modeling.

Lindholm, F. A.↗

GPU acceleration of rank-reduced coupled-cluster singles and doubles

Here, we have developed a graphical processing unit (GPU) accelerated implementation of our recently introduced rank-reduced coupled-cluster singles and doubles (RR-CCSD) method. RR-CCSD introduces a low-rank approximation of the doubles amplitudes. This is combined with a low-rank approximation of the electron repulsion integrals via Cholesky decomposition. The result of these two low-rank approximations is the replacement of the usual fourth-order CCSD tensors with products of second- and third-order tensors. In our implementation, only a single fourth-order tensor must be constructed as an intermediate during the solution of the amplitude equations. Owing in large part to the compression of the doubles amplitudes, the GPU-accelerated implementation shows excellent parallel efficiency (95% on eight GPUs). Our implementation can solve the RR-CCSD equations for up to 400 electrons and 1550 basis functions—roughly 50% larger than the largest canonical CCSD computations that have been performed on any hardware. In addition to increased scalability, the RR-CCSD computations are faster than the corresponding CCSD computations for all but the smallest molecules. We test the accuracy of RR-CCSD for a variety of chemical systems including up to 1000 basis functions and determine that accuracy to better than 0.1% error in the correlation energy can be achieved with roughly 95% compression of the ov space for the largest systems considered. We also demonstrate that conformational energies can be predicted to be within 0.1 kcal mol -1 with efficient compression applied to the wavefunction. Finally, we find that low-rank approximations of the CCSD doubles amplitudes used in the similarity transformation of the Hamiltonian prior to a conventional equation-of-motion CCSD computation will not introduce significant errors (on the order of a few hundredths of an electronvolt) into the resulting excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗

Deciphering coulombic loss in lithium-ion batteries and beyond

Lithium-ion batteries are pivotal for modern energy storage, yet accurately predicting their lifespan remains a critical challenge. While descriptors like coulombic efficiency are widely used to assess battery longevity, the unclear physical origins of coulombic losses cause semi-quantitative correlation with capacity, complicating battery development. Here, we combine high-precision leakage current and open-circuit-voltage measurements with charge conservation principles to explore microscopic charge consumptions at electrode-electrolyte interfaces across diverse chemistries. We demonstrate that coulombic loss arises from a synergy between local charge neutrality and global charge compensation, reconciling its quantitative correlation to capacity. Contrary to conventional assumptions equating coulombic loss with irreversible capacity loss, this framework resolves systematic overestimations and paradoxical phenomena in existing chemistries. Our findings establish physics-informed criteria for accelerated lifespan evaluation and guide rational design of long-life lithium-ion batteries and beyond.

Batteries↗

Liquid solution centrifugation for safe, scalable, and efficient isotope separation

A general method of separating isotopes by centrifuging dissolved chemical compounds in a liquid is introduced. This technique can be applied to almost all elements and leads to large separation factors. The method has been demonstrated in several isotopic systems including Ca, Mo, O, and Li with single-stage selectivities of 1.046 to 1.067 per neutron mass difference (e.g., 1.43 in 40 Ca/ 48 Ca), which are beyond the capabilities of various conventional methods. Equations are derived to model the process, and the results agree with those of the experiments. The scalability of the technique has been demonstrated by a three-stage enrichment of 48 Ca with a total 40 Ca/ 48 Ca selectivity of 2.43, and the scalability is more broadly supported through analogies to gas centrifuge, whereby countercurrent centrifugation can further multiply the separation factor by 5 to 10 times per stage in a continuous process. Optimal centrifuge conditions and solutions can achieve both high-throughput and highly efficient isotope separation.

07 ISOTOPE AND RADIATION SOURCES↗

Transmitter and receiver antenna gain analysis for laser radar and communication systems

A comprehensive and fairly self-contained study of centrally obscured optical transmitting and receiving antennas is presented and is intended for use by the laser radar and communication systems designer. The material is presented in a format which allows the rapid and accurate evaluation of antenna gain. The Fresnel approximation to scalar wave theory is reviewed and the antenna analysis proceeds in terms of the power gain. Conventional range equations may then be used to calculate the power budget. The transmitter calculations, resulting in near and far field antenna gain patterns, assumes the antenna is illuminated by a laser operating in the fundamental cavity mode. A simple equation is derived for matching the incident source distribution to a general antenna configuration for maximum on-axis gain. An interpretation of the resultant gain curves allows a number of auxiliary design curves to be drawn which display the losses in antenna gain due to pointing errors and the cone angle of the outgoing beam as a function of antenna size and central obscuration. The use of telescope defocusing as an approach to spreading the beam for target acquisition is compared to some alternate methods.

Klein, B. J.↗

Photoconductor-metal contact at higher densities.

The formation of space charge in the barrier region is discussed. This region extends to a distance, x-sub o, equal to 200 A plus twice the mean free path of majority carriers. The conventional transport equation can be used only for x greater than x-sub o, and the carrier density at x-sub o represents a boundary condition for the bulk. Its change as a function of applied voltage, temperature, and light intensity in photoconducting CdS is discussed. The time dependence of the space-charge formation in the region between 0 and x-sub o is analyzed. It is shown under which conditions the individuality of the metal contact is observable.

Dussel, G. A.↗

Buckling of periodic structures

Equations are developed for the buckling of a general lattice structure that has repetitive geometry. Equilibrium at a typical node is expressed using finite element techniques, and the only assumption is that the response is periodic. By basing the stiffness matrix on the exact solution of the beam column equation, accurate results are obtained for complex buckling behavior that would require a very large system of equations using conventional techniques. The present method requires the eigenvalues of only a 6x6 determinant. The results are used to study the buckling of isogrid cylinders, three-element truss columns and polygonal rings. Details of the analysis including expressions for all terms in the governing stability determinant are given.

Anderson, M. S.↗

On the theory of cosmic-ray-mediated shocks with variable compression ratio

Cosmic-ray-mediated shocks may accelerate enough cosmic rays to high enough energies that they escape the shock, carrying an appreciable amount of energy before being convected to downstream infinity. Under such conditions, it is noted, the overall compression ratio cannot be determined from the conservation equations as in conventional hydrodynamic treatments, and the standard equations for shock acceleration admit arbitrarily high compression ratios. A procedure is outlined for obtaining the structure of high Mach number, cosmic-ray-mediated shocks, including their overall compresion ratio, around a low Mach number viscous subshock. Analytic solutions are obtained by quardrature for an energy-dependent diffusion coefficient in the limit of extreme sensitivity to energy, which, unlike previous solutions, include the finite thermal pressure of the preshock gas.

Eichler, D.↗

The finite-difference time-domain (FD-TD) method for electromagnetic scattering and interaction problems

The formulation and recent applications of the finite-difference time-domain (FD-TD) method for the numerical modeling of electromagnetic scattering and interaction problems are considered. It is shown that improvements in FD-TD modeling concepts and software implementation often make it a preferable choice for structures which cannot be easily treated by conventional integral equations and asymptotic approaches. Recent FD-TD modeling validations in research areas including coupling to wires and wire bundles in free space and cavities, scattering from surfaces in relativistic motion, inverse scattering, and radiation condition theory, are reviewed. Finally, the advantages and disadvantages of FD-TD, and guidelines concerning when FD-TD should and should not be used in high-frequency electromagnetic modeling problems, are summarized.

Taflove, A.↗

Effect of general relativity on a near-Earth satellite in the geocentric and barycentric reference frames

Whether one uses a solar-system barycentric frame or a geocentric frame when including the general theory of relativity in orbit determinations for near-Earth satellites, the results should be equivalent to some limiting accuracy. The purpose of this paper is to clarify the effects of relativity in each frame and to demonstrate their equivalence through the analysis of real laser-tracking data. A correction to the conventional barycentric equations of motion is shown to be required.

Ries, J. C.↗

Efficient load measurements using singular value decomposition

Various basic research was performed on efficient load measurement estimation techniques for aircraft structure analysis. An overview is presented of the load measurement problem. Two basic equivalent approaches to load measurement evaluations were considered. Under approach 1, the load values are modeled as depending linearly on the measured values. Under approach 2, the measured values depend linearly on the load values. By using the modern Singular Value Decomposition method, it was shown that under all conditions of the number of loads and number of gages, approach 1 is equivalent to approach 2. By using the conventional normal equation (linear regression) approach, approach 1 is only valid when the number of loads is equal to or greater than the number of gages, while approach 2 is the reverse. Furthermore, except for the case of the number of loads equals the number of gages, the load prediction formulas under the two approaches are not equivalent.

Yao, Kung↗

Determination of joint drives for stable end-effector motion in flexible robotic systems

The prescribed tasks in high speed robotic systems are severely deteriorated because of their manipulator dynamic deflections. On the other hand conventional dynamic modeling techniques fail to reveal appropriate control forces in flexible systems. The conventional dynamic equations of motion for systems subject to kinematical constraints are modified by a new concept of control force representation. The directions of the control forces are selected such that they correspond to the joint degrees of freedom. Then the joint control forces and torques that yield unperturbed prescribed motions are solved simultaneously with the system motion. A flexible manipulator is presented to illustrate the methods proposed.

Ider, Sitki Kemal↗

The effects of general relativity on near-earth satellites

Whether one uses a solar system barycentric frame or a geocentric frame when including the general theory of relativity in orbit determination for near-earth satellites, the results should be equivalent to some limiting accuracy. The purpose of this paper is to clarify the effects of relativity in each frame and to demonstrate their equivalence through the analysis of three years of laser tracking data taken on the Lageos satellite. It is demonstrated that the simpler formulation in the geocentric frame is adequate for the purpose of near-earth satellite orbit determination. A correction to the conventional barycentric equations of motion is shown to be required.

Ries, J. C.↗

On two-dimensional relativistic stellar winds

Stellar wind is of great interest, because many of the astrophysical systems possess wind like structures. For a long time only nonrelativistic winds have been studied. However, recently it was proposed by Kennel et al. that the relativistic pulsar wind with the plasma, consisting of electrons and positrons, can be responsible for the observed features of the Crab nebula. The study of Kennel et al. revealed the inconsistency of the assumption of the wind zero temperature with the observational data. It has been shown that only for high relativistic temperatures can high Mach numbers be reached, which allows the possibility of a shock formation. The analysis of Kennel et al. was extended into the 2-D case by using conventional magnetohydrodynamic equations for a relativistic plasma with an isotropic relativistic temperature. The state equation is assumed polytropic.

Gedalin, M. E.↗