Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geometric discontinuities”

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 73 records · Page 4

Path Profiles of Cn2 Derived from Radiometer Temperature Measurements and Geometrical Ray Tracing

Atmospheric turbulence has significant impairments on the operation of Free-Space Optical (FSO) communication systems, in particular temporal and spatial intensity fluctuations at the receiving aperture resulting in power surges and fades, changes in angle of arrival, spatial coherence degradation, etc. The refractive index structure parameter Cn2 is a statistical measure of the strength of turbulence in the atmosphere and is highly dependent upon vertical height. Therefore to understand atmospheric turbulence effects on vertical FSO communication links such as space-to-ground links, it is necessary to specify Cn2 profiles along the atmospheric propagation path. To avoid the limitations on the applicability of classical approaches, propagation simulation through geometrical ray tracing is applied. This is achieved by considering the atmosphere along the optical propagation path as a spatial distribution of spherical bubbles with varying relative refractive index deviations representing turbulent eddies. The relative deviations of the refractive index are statistically determined from altitude-dependent and time-varying temperature fluctuations, as measured by a microwave profiling radiometer. For each representative atmosphere ray paths are analyzed using geometrical optics, which is particularly advantageous in situations of strong turbulence where there is severe wavefront distortion and discontinuity. The refractive index structure parameter is then determined as a function of height and time.

Atmospheric turbulence↗

Constraints on sequential discontinuities from the geometry of on-shell spaces

We present several classes of constraints on the discontinuities of Feynman integrals that go beyond the Steinmann relations. These constraints follow from a geometric formulation of the Landau equations that was advocated by Pham, in which the singularities of Feynman integrals correspond to critical points of maps between on-shell spaces. To establish our results, we review elements of Picard-Lefschetz theory, which connect the homotopy properties of the space of complexified external momenta to the homology of the combined space of on-shell internal and external momenta. An important concept that emerges from this analysis is the question of whether or not a pair of Landau singularities is compatible — namely, whether or not the Landau equations for the two singularities can be satisfied simultaneously. Under conditions we describe, sequential discontinuities with respect to non-compatible Landau singularities must vanish. Although we only rigorously prove results for Feynman integrals with generic masses in this paper, we expect the geometric and algebraic insights that we gain will also assist in the analysis of more general Feynman integrals.

97 MATHEMATICS AND COMPUTING↗

Free Surface Phenomena Under Low- and Zero-Gravity Conditions

In a low- or zero-gravity environment the free surface of a liquid can behave in striking, unexpected ways. For example, in a partially filled container, a free surface that is well behaved under terrestrial conditions can rise to an arbitrarily large height or even fail to exist when gravity is absent. An answer is sought to the central mathematical questions: Under what conditions can free surfaces exist and what are their properties? -- and experimental questions of what means can be devised to observe and to measure the surfaces quantitatively. Current activity focuses on material selection and the design of optical diagnostic methods for in-space experiments. These experiments will test mathematical predictions of discontinuous transition from existence to nonexistence of capillary free surfaces in certain geometric configurations under zero gravity conditions. Liquids and container materials suitable for achieving the required contact-angle range and optical properties are being investigated, as well as the effects of contaminants and fluid motion on contact angle. The optical diagnostic technique investigation concerns laser-induced fluorescence, including emphasis on data on acquisition, sources and detectors, reliability, data management, and computer control.

Concus, P.↗

XMM-Newton and TESS observations of the highly variable polar V496 UMa

We studied the temporal and spectral behavior of V496 UMa from the optical to the X-ray regimes. Methods. We used archival XMM-Newton and TESS observations obtained in 2017 and 2019, respectively, to perform a spectral and timing analysis of the highly variable magnetic CV. Results. The light curves of the TESS and XMM-Newton satellites reveal a double-humped pattern modulated with the periodicity of 91.058467 ± 0.00001 min. V496 UMa displays a two-pole accretion geometry in the high accretion state. X-ray spectra from these regions are composed of thermal plasma radiation and soft blackbody components with almost identical temperatures and a total accretion rate of $\dot{M}$ = 1.4(8) × 10 -11 M ⊙ yr -1 . The X-ray centers of the humps show longitudinal shifts of -18° and 4°, and shifts around photometric phase zero of -172° and -186°, for the main hump and second hump, respectively. The long-term ZTF light curves reveal high and low accretions states. Low-state ZTF and SDSS photometric data are consistent with an 0.8 M ⊙ white dwarf at 10 000 K and a main-sequence donor star with a spectral type of M5.0 at a Gaia determined distance of 758 pc. Conclusions. V496 UMa is a very bright polar in X-rays when it is in the high state. Due to its unusual geometric structure, mass accretion onto the second accretion pole is interrupted occasionally. This discontinuous behavior does not follow a certain pattern in time and has been observed so far only in the high state. The X-ray light curves display clear evidence of an accretion stream at the photometric phase of φ = 0.81, which does not show up in optical light curves. An accurate period was derived using the combined TESS and XMM-Newton data, which differs by 3.8 σ from published results.

79 ASTRONOMY AND ASTROPHYSICS↗

I-BIEM, an iterative boundary integral equation method for computer solutions of current distribution problems with complex boundaries: A new algorithm. I - Theoretical

A new algorithm for an iterative computation of solutions of Laplace's or Poisson's equations in two dimensions, using Green's second identity, is presented. This algorithm converges strongly and geometrically and can be applied to curved, irregular, or moving boundaries with nonlinear and/or discontinuous boundary conditions. It has been implemented in Pascal on a number of micro- and minicomputers and applied to several geometries. Cases with known analytic solutions have been tested. Convergence to within 0.1 percent to 0.01 percent of the theoretical values are obtained in a few minutes on a microcomputer.

Cahan, B. D.↗

Distributed optical proximity sensor system - HexEYE

A novel optical proximity sensor system, capable of measuring the distance to the orientation and the discontinuity at a local area of an object surface, was designed and evaluated for robotic applications. The sensor system gets its name, Hexagonal Eye, (HexEYE) from its shape where seven identical hexagonal sensing units were configured hexagonally into a compact geometric structure. The seven sensing units were grouped into multiple combinations of three sensing units to measure the surface orientation as well as surface discontinuity. The distinctive features of HexEYE are an order of magnitude increased in distance sensitivity by optimizing the curvature of a conical mirror, the enhanced measurement accuracy based on multiple levels of sensor fusion, and the compactness in size due to a sensing mechanism based on the Gaussian lens law. A prototype of single sensing unit has been built and was evaluated experimentally.

Lee, Sukhan↗

Alfven wave refraction by interplanetary inhomogeneities

Pioneer 6 magnetic data reveals that the propagation direction of Alfven waves in the interplanetary medium is strongly oriented along the ambient field. Magnetic fluctuations of frequencies up to 1/30 sec in the spacecraft frame are shown to satisfy a necessary condition for Alfven wave normal. It appears from this analysis that geometrical hydromagnetics may satisfactorily describe deviation of the wave normal from the background field. The rotational discontinuity is likely also to propagate along the field lines.

Daily, W. D.↗

Alfven wave refraction by interplanetary inhomogeneities.

Pioneer 6 magnetic data reveal that the propagation direction of Alfven waves in the interplanetary medium is strongly oriented along the ambient field. Magnetic fluctuations of frequencies up to 1/30 sec in the spacecraft frame are shown to satisfy a necessary condition for Alfven waves, and a variance matrix analysis is used to determine the Alfven wave normal. It appears from this analysis that geometrical hydromagnetics may satisfactorily describe deviation of the wave normal from the background field. The rotational discontinuity is likely also to propagate along the field lines.

Daily, W. D.↗

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗

Computer modeling of electromagnetic problems using the geometrical theory of diffraction

Some applications of the geometrical theory of diffraction (GTD), a high frequency ray optical solution to electromagnetic problems, are presented. GTD extends geometric optics, which does not take into account the diffractions occurring at edges, vertices, and various other discontinuities. Diffraction solutions, analysis of basic structures, construction of more complex structures, and coupling using GTD are discussed.

Burnside, W. D.↗

An efficient high-frequency analysis of modal reflection and transmission coefficients for a class of waveguide discontinuities

A high-frequency analysis of electromagnetic modal reflection and transmission coefficients is presented for waveguide discontinuities formed by joining different waveguide sections. The analysis uses an extended version of the concept of geometrical theory of diffraction based equivalent edge currents in conjunction with the reciprocity theorem to describe interior scattering effects. If the waveguide modes and their associated modal rays can be found explicitly, general two- and three-dimensional waveguide geometries can be analyzed. Expressions are developed for two-dimensional reflection and transmission coefficients. Numerical results are given for a flanged, semi-infinite parallel plate waveguide and for the junction between two linearly tapered waveguides.

Pathak, P. H.↗

A bubblelike coronal mass ejection flux rope in the solar wind

A resolution to the question of whether coronal mass ejections are loops or bubbles is proposed and applied to the geometrical analysis of a solar wind event detected at 1 AU by ISEE 1 and 3. The discontinuity orientations, the size determined by time of passage, and the magnetic cloud signature are fit into the topology of a flux rope loop distorted by expansion into a thick rope with comparable dimensions in both the ecliptic and meridional planes. The looped rope fills a bubblelike cavity, thus preserving both types of proposed coronal mass ejection geometries. Other interesting features of the data include an apparent separation by the rope core of bidirectionally streaming protons in the leading section from electrons in the trailing section, possible vortical flow within the magnetic cloud, and a well-defined filamentary structure behind the shock.

Crooker, N. U.↗

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗

Nonlinear dynamics, bifurcations, and multi-stability in a vibro-impact system with geometric and multi-segmented freeplay nonlinearities

Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. Here, the effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.

42 ENGINEERING↗

Estimation of subpixel vegetation cover using red-infrared scattergrams

The bulk properties of discontinuous vegetation canopies are estimated at subpixel scales by applying the method of moments to a linear stochastic geometric model of canopy-soil reflectance and one set of multispectral observations without ground truth. The procedure involves the formulation of conditional moments for subsets of pixels that possess similar properties and can be identified through their common orientation in red-infrared scattergrams. The analysis is facilitated by assuming geometric similarity among the canopy elements and by formulating a sampling scale ratio in terms of the bulk geometric scales of the canopy and the pixel. Three versions of the method are demonstrated using two simulated scenes and an actual forested watershed for which aerial radiometric data and corresponding ground truth were obtained.

Jasinski, Michael F.↗

Evaluation of a strained-coordinate perturbation procedure - Nonlinear subsonic and transonic flows

An evaluation is made of a perturbation method devised to obtain highly accurate approximations to families of strongly nonlinear solutions which are either continuous or discontinuous, and which represent variations in some arbitrary parameter. The method first defines a unit perturbation by using two nonlinear solutions which differ from one another by a nominal change in some geometric or flow parameter, then employs this unit perturbation to predict a family of related nonlinear solutions over a range of parameter variation. Coordinate straining is incorporated into this perturbation method for determining the unit perturbation to account for the movement of discontinuities and maxima of high-gradient regions due to the perturbation. Attention is given to transonic and subsonic flows. Comparisons of the perturbation results with the corresponding 'exact' nonlinear solutions show a remarkable accuracy and range of validity of the perturbation method across the spectrum of examples considered.

Stahara, S. S.↗

Unification of color postprocessing techniques for 3-dimensional computational mechanics

To facilitate the understanding of complex three-dimensional numerical models, advanced interactive color postprocessing techniques are introduced. These techniques are sufficiently flexible so that postprocessing difficulties arising from model size, geometric complexity, response variation, and analysis type can be adequately overcome. Finite element, finite difference, and boundary element models may be evaluated with the prototype postprocessor. Elements may be removed from parent models to be studied as independent subobjects. Discontinuous responses may be contoured including responses which become singular, and nonlinear color scales may be input by the user for the enhancement of the contouring operation. Hit testing can be performed to extract precise geometric, response, mesh, or material information from the database. In addition, stress intensity factors may be contoured along the crack front of a fracture model. Stepwise analyses can be studied, and the user can recontour responses repeatedly, as if he were paging through the response sets. As a system, these tools allow effective interpretation of complex analysis results.

Bailey, Bruce Charles↗

Georgia Tech Accelerated, Compressed, and Regularized Compute of Kinetic-based PDEs (Final Report)

This report summarizes the collaborative effort between Lawrence Livermore National Laboratory and Georgia Tech to enhance the BoBa library for tensor train computation in PDE solvers, with a target on kinetic equations and their continuum limits. We aimed to reduce computational cost and memory usage by replacing traditional array-based computations with tensor trains. We examined the compressibility of time-evolving solutions to the Euler equations with discontinuities. We also explored using the first invsicid and linear regularization of the compressible flow equations via the information geometric regularization (IGR). We explored this in a tensor train formulation. To identify that inverse terms in the IGR equations pose problems for tensor train formulations and investigate efficient methods for batched inversion of tensor trains.

97 MATHEMATICS AND COMPUTING↗