PROPOSED STUDIES IN THE FIELD OF MULTIVARIATE APPROXIMATIONS
Numerical analysis of functions with independent variables - guidance systems
MATHEMATICS↗
SEARCH · Engineering Papers
Results for “Numerical approximation & analysis”
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.
Chemical control of rate and onset temperature of nadimide polymerization
The chemistry of norbornenyl capped imide compounds (nadimides) is briefly reviewed with emphasis on the contribution of Diels-Alder reversion in controlling the rate and onset of the thermal polymerization reaction. Control of onset temperature of the cure exotherm by adjusting the concentration of maleimide is demonstrated using selected model compounds. The effects of nitrophenyl compounds as free radical retarders on nadimide reactivity are discussed. A simple copolymerization model is proposed for the overall nadimide cure reaction. An approximate numerical analysis is carried out to demonstrate the ability of the model to simulate the trends observed for both maleimide and nitrophenyl additions.
Lauver, R. W.↗
Synchrotron emissivity near the electron cyclotron and upper hybrid frequencies
The spontaneous synchrotron emissivity from a high density, weakly relativistic plasma in thermal equilibrium is discussed. Thermal effects on the dielectric polarization of the plasma are included in the analysis, and the result is compared with computations of the emissivity based upon the cold plasma approximation. A numerical analysis is performed for frequencies in the vicinity of both the electron cyclotron frequency and the upper hybrid resonance. Significant modifications to the synchrotron emissivity in the cold plasma approximation are found in both of these frequency regimes.
Freund, H. P.↗
Analytical and numerical modeling of optical second harmonic generation in anisotropic crystals using #SHAARP package
Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. #SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.
36 MATERIALS SCIENCE↗
Ramsey envelope modulation in NV diamond magnetometry
Nitrogen-vacancy (NV) spin ensembles in diamond provide an advanced magnetic sensing platform, with applications in both the physical and life sciences. The development of isotopically engineered 15 NV diamond offers advantages over naturally occurring 14 NV for magnetometry, due to its simpler hyperfine structure. However, for sensing modalities requiring a bias magnetic field not aligned with the sensing NV axis, the absence of a quadrupole moment in the 15 N nuclear spin leads to pronounced envelope modulation effects in time-dependent measurements of 15 NV spin evolution. While such behavior in spin echo experiments are well studied, analogous effects in Ramsey measurements and the implications for magnetometry remain underexplored. Here, we derive the modulated 15 NV Ramsey response to a misaligned bias field, using a simple vector description of the effective magnetic field on the nuclear spin. The predicted modulation properties are then compared to experimental results, revealing significant magnetic sensitivity loss if unaddressed. We demonstrate that double-quantum coherences of the NV S=1 electronic spin states dramatically suppress these envelope modulations, while additionally proving resilient to other parasitic effects such as strain heterogeneity and temperature shifts.
36 MATERIALS SCIENCE↗
Volume extrapolation via eigenvector continuation
Here, we develop an extension of eigenvector continuation (EC) that makes it possible to extrapolate simulations of quantum systems in finite periodic boxes across large ranges of box sizes. The formal justification for this approach, which we call finite-volume eigenvector continuation (FVEC), is provided by matching periodic functions at different box sizes. As concrete FVEC implementation we use a discrete variable representation based on plane-wave states and present several applications calculated within this framework.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Universal relations for neutron star 𝑓-mode and 𝑔-mode oscillations
Among the various oscillation modes of neutron stars, 𝑓- and 𝑔- modes are the most likely to be ultimately observed in binary neutron star mergers due to their relatively large coupling and shared frequencies with tidal excitations. The 𝑓-mode frequency and damping time are known to correlate in normal neutron stars with their compactness, and previous fits to hadronic stars are extended and shown to be valid for an extremely broad sampling of equations of state using a piecewise polytropic parametrization scheme for hadrons and a constant sound-speed parametrization for quark matter. Separate fits applicable to quark (self-bound) stars are improved. Much more significant correlations exist with tidal deformability, and therefore with moment of inertia and quadrupole moment. It is conclusively demonstrated that these correlations are the same for all types of stars, whether hadronic, hybrid, or pure quark, and its accuracy is quantified. A novel 1-node branch of the 𝑓-mode that occurs in low-mass hybrid stars in a narrow mass range just beyond the critical mass necessary for a phase transition to appear is identified. This 1-node branch shows the largest, but still small, deviations from the universal correlation for any configuration. It is characterized by a nonmonotonic relation between neutron star mass and 𝑓-mode frequency, in contrast to the behavior otherwise observed in normal, quark and hybrid stars. The 𝑔-mode only exists in matter with a nonbarotropic equation of state involving temperature, chemical potential or composition (such as being out of beta equilibrium), or a phase transition in barotropic matter. Here, the 𝑔-mode therefore could serve as a probe for studying phase transitions in hybrid stars. In contrast with the 𝑓-mode, 𝑔-mode frequencies do not correlate well with tidal deformability, but depend strongly on properties of the transition (the density and the magnitude of the discontinuity) at the transition. Imposing causality and maximum mass constraints, a fit involving neutron star and phase transition properties is found and the 𝑔-mode frequency is determined to have an upper bound of about 1.25 kHz. However, if the sound speed 𝑐𝑠 in the inner core at densities above the phase transition density is restricted to 𝑐$^{2}_{𝑠}$ ≤1/3, 𝑔-mode frequencies can only reach about 0.8 kHz, which are significantly lower than 𝑓-mode frequencies (1.3–2.8 kHz). 𝑔-mode gravitational wave damping times are found to be extremely long, >10 4 s (102 s) in the inner core with 𝑐$^{2}_{𝑠}$ ≤1/3 (1), in comparison with 𝑓-mode damping times (0.1–1 s).
Composition of astronomical objects↗
Spectral density reconstruction with Chebyshev polynomials
Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semiexclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. Here, in the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to obtain the local density of states for graphene in a magnetic field as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
CSM research: Methods and application studies
Computational mechanics is that discipline of applied science and engineering devoted to the study of physical phenomena by means of computational methods based on mathematical modeling and simulation, utilizing digital computers. The discipline combines theoretical and applied mechanics, approximation theory, numerical analysis, and computer science. Computational mechanics has had a major impact on engineering analysis and design. When applied to structural mechanics, the discipline is referred to herein as computational structural mechanics. Complex structures being considered by NASA for the 1990's include composite primary aircraft structures and the space station. These structures will be much more difficult to analyze than today's structures and necessitate a major upgrade in computerized structural analysis technology. NASA has initiated a research activity in structural analysis called Computational Structural Mechanics (CSM). The broad objective of the CSM activity is to develop advanced structural analysis technology that will exploit modern and emerging computers, such as those with vector and/or parallel processing capabilities. Here, the current research directions for the Methods and Application Studies Team of the Langley CSM activity are described.
Knight, Norman F., Jr.↗
Computers boost structural technology
Derived from matrix methods of structural analysis and finite element methods developed over the last three decades, computational structures technology (CST) blends computer science, numerical analysis, and approximation theory into structural analysis and synthesis. Recent significant advances in CST include stochastic-based modeling, strategies for performing large-scale structural calculations on new computing systems, and the integration of CST with other disciplinary modules for multidisciplinary analysis and design. New methodologies have been developed at NASA for integrated fluid-thermal structural analysis and integrated aerodynamic-structure-control design. The need for multiple views of data for different modules also led to the development of a number of sophisticated data-base management systems. For CST to play a role in the future development of structures technology and in the multidisciplinary design of future flight vehicles, major advances and computational tools are needed in a number of key areas.
Noor, Ahmed K.↗
Computational structures technology and UVA Center for CST
Rapid advances in computer hardware have had a profound effect on various engineering and mechanics disciplines, including the materials, structures, and dynamics disciplines. A new technology, computational structures technology (CST), has recently emerged as an insightful blend between material modeling, structural and dynamic analysis and synthesis on the one hand, and other disciplines such as computer science, numerical analysis, and approximation theory, on the other hand. CST is an outgrowth of finite element methods developed over the last three decades. The focus of this presentation is on some aspects of CST which can impact future airframes and propulsion systems, as well as on the newly established University of Virginia (UVA) Center for CST. The background and goals for CST are described along with the motivations for developing CST, and a brief discussion is made on computational material modeling. We look at the future in terms of technical needs, computing environment, and research directions. The newly established UVA Center for CST is described. One of the research projects of the Center is described, and a brief summary of the presentation is given.
Noor, Ahmed K.↗
Computational structures technology
Computational structures technology (CST), which has emerged from FEM developments, is a fusion of materials modeling, structural and dynamic analysis and synthesis methods, on the one hand, and numerical analysis and approximation theory, on the other. In addition to computational materials modeling, CST encompasses computational methods for predicting the response, performance, failure, and service life of structures and their components, as well as automated methods for structural synthesis and optimization.
Noor, Ahmed K.↗
A scaling law for energy transfer by inelastic electron-molecule collisions in mixtures
The equation governing the electron energy distribution in the presence of a spatially uniform electric field in a weakly ionized gas was reformulated into an integral equation for the logarithmic slope of the distribution function. For gas mixtures in which the dominant electron energy loss mechanism is by vibrational excitation of the molecules, this equation is suitable for approximate analysis and exact numerical solution by iteration. Superelastic collisions are easily included in this formulation, and do not seriously effect the convergence of the numerical scheme. The approximate analytical results are only qualitatively correct, but suggest appropriate parameters which correlate the exact numerical results very well. The distribution function as well as certain gross properties such as net energy transfer into vibration, mean energy, and drift velocity depend primarily on a single nondimensional parameter involving only E/N and the cross sections.
Bienkowski, G. K.↗
Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations
Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.
Beam, R. M.↗
Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations
Previously announced in STAR as N81-33870
Beam, R. M.↗
Derivation and evaluation of an approximate analysis for three-dimensional viscous subsonic flow with large secondary velocities
An approximate analysis is presented for calculating three-dimensional, low Mach number, laminar viscous flows in curved passages with large secondary flows and corner boundary layers. The analysis is based on the decomposition of the overall velocity field into inviscid and viscous components with the overall velocity being determined from superposition. An incompressible vorticity transport equation is used to estimate inviscid secondary flow velocities to be used as corrections to the potential flow velocity field. A parabolized streamwise momentum equation coupled to an adiabatic energy equation and global continuity equation is used to obtain an approximate viscous correction to the pressure and longitudinal velocity fields. A collateral flow assumption is invoked to estimate the viscous correction to the transverse velocity fields. The approximate analysis is solved numerically using an implicit ADI solution for the viscous pressure and velocity fields. An iterative ADI procedure is used to solve for the inviscid secondary vorticity and velocity fields. This method was applied to computing the flow within a turbine vane passage with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, and for a constant radius curved rectangular duct with R/D = 12 and 14 and with inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls.
Anderson, O. L.↗
Turbulent flow field predictions in sharply curved turn around ducts
In this investigation, two-dimensional turbulent flow of incompressible Newtonian fluids in sharply curved 180 deg turn around ducts is studied. Results of an approximate numerical flow field analysis utilizing an orthogonal, body-fitted, curvilinear coordinate system are compared to results based on a traditional cylindrical reference frame. Qualitative indication of general streamfield characteristics as well as quantitative benchmarks for the planning of future experimentation are provided. In addition, preliminary results of an augmented kappa-epsilon turbulence model analysis, which explicitly accounts for the effects of streamline curvature and pressure strain in internal turbulent flows, are presented. Specific model difficulties are discussed and comparisons with standard kappa-esilon model predictions are included.
Santi, L. M.↗
Turbulent flow field predictions in sharply curved turn-around ducts
Two-dimensional turbulent flow of incompressible Newtonian fluids in sharply curved 180 deg turn around ducts is studied. Results of an approximate numerical flow field analysis utilizing an orthogonal, body-fitted, curvilinear coordinate system are compared to results based on a traditional cylindrical reference frame. Qualitative induction of general streamfield characteristics as well as quantitative benchmarks for the planning of future experimentation are provided. In addition, preliminary results of an augmented kappa - epsilon turbulence model analysis, which explicitly accounts for the effects of streamline curvature and pressure strain in internal turbulent flows, are presented. Specific model difficulties are discussed and comparisons with standard kappa - epsilon model predictions are included.
Santi, L. M.↗