Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partitioned methods”

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 343 records · Page 19

The Pilot Lunar Geologic Mapping Project: Summary Results and Recommendations from the Copernicus Quadrangle

The first systematic lunar geologic maps were completed at 1:1M scale for the lunar near side during the 1960s using telescopic and Lunar Orbiter (LO) photographs [1-3]. The program under which these maps were completed established precedents for map base, scale, projection, and boundaries in order to avoid widely discrepant products. A variety of geologic maps were subsequently produced for various purposes, including 1:5M scale global maps [4-9] and large scale maps of high scientific interest (including the Apollo landing sites) [10]. Since that time, lunar science has benefitted from an abundance of surface information, including high resolution images and diverse compositional data sets, which have yielded a host of topical planetary investigations. The existing suite of lunar geologic maps and topical studies provide exceptional context in which to unravel the geologic history of the Moon. However, there has been no systematic approach to lunar geologic mapping since the flight of post-Apollo scientific orbiters. Geologic maps provide a spatial and temporal framework wherein observations can be reliably benchmarked and compared. As such, a lack of a systematic mapping program means that modern (post- Apollo) data sets, their scientific ramifications, and the lunar scientists who investigate these data, are all marginalized in regard to geologic mapping. Marginalization weakens the overall understanding of the geologic evolution of the Moon and unnecessarily partitions lunar research. To bridge these deficiencies, we began a pilot geologic mapping project in 2005 as a means to assess the interest, relevance, and technical methods required for a renewed lunar geologic mapping program [11]. Herein, we provide a summary of the pilot geologic mapping project, which focused on the geologic materials and stratigraphic relationships within the Copernicus quadrangle (0-30degN, 0-45degW).

Skinner, J. A., Jr.↗

The Germanium Dichotomy in Martian Meteorites

Germanium is a moderately volatile and siderophile element that follows silicon in its compatibility during partial melting of planetary mantles. Despite its obvious usefulness in planetary geochemistry germanium is not analyzed routinely, with there being only three prior studies reporting germanium abundances in Martian meteorites. The broad range (1-3 ppm) observed in Martian igneous rocks is in stark contrast to the narrow range of germanium observed in terrestrial basalts (1.5 plus or minus 0.1 ppm). The germanium data from these studies indicates that nakhlites contain 2-3 ppm germanium, while shergottites contain approximately 1 ppm germanium, a dichotomy with important implications for core formation models. There have been no reliable germanium abundances on chassignites. The ancient meteoritic breccia, NWA 7533 (and paired meteorites) contains numerous clasts, some pristine and some impact melt rocks, that are being studied individually. Because germanium is depleted in the Martian crust relative to chondritic impactors, it has proven useful as an indicator of meteoritic contamination of impact melt clasts in NWA 7533. The germanium/silicon ratio can be applied to minerals that might not partition nickel and iridium, like feldspars. We report germanium in minerals from the 3 known chassignites, 2 nakhlites and 5 shergottites by LAICP- MS using a method optimized for precise germanium analysis.

Humayun, M.↗

Efficient computation of volume in flow predictions

An efficient method for calculating cell volumes for time-dependent three-dimensional flow predictions by finite volume calculations is presented. Eight arbitrary corner points are considered and the shape face is divided into two planar triangles. The volume is then dependent on the orientation of the partitioning. In the case of a hexahedron, it is noted that any open surface with a boundary that is a closed curve possesses a surface vector independent of the surface shape. Expressions are defined for the surface vector, which is independent of the partitioning surface diagonal used to quantify the volume. Using a decomposition of the cell volume involving two corners, with each the vertex of three diagonals and six corners which are vertices of one diagonal, gives portions which are tetrahedra. The resultant mesh is can be used for time-dependent finite volume calculations one requires less computer time than previous methods.

Vinokur, M.↗

Correlations in cosmic density fields

A method is proposed to place constraints on the functional form of the high-order correlation functions zeta(sub n) that arise in cosmic density fields at large scales. This technique is based on a mass-in-cell statistic and a difference of mass in partitions of a cell. The relationship between these measures is sensitive to the formal structure of the zeta(sub n) as well as their amplitudes. This relationship is quantified in several theoretical models of structure, based on the hierarchical clustering paradigm. The results lead to a test for specific types of hierarchical clustering that is sensitive to correlations of all orders. The method is applied to examples of simulated large-scaled structure dominated by cold dark matter. In the preliminary study, the hierarchical paradigm appears to be a realistic approximation over a broad range of the scales. Furthermore, there is evidence that graphs of low-order vertices are dominant. On the basis of simulated data a phenomological model is specified that gives a good representation of clustering from linear scales to the strongly clustered regime (zeta(sub 2) approximately 500).

Bromley, B. C.↗

A computational procedure for large rotational motions in multibody dynamics

A computational procedure suitable for the solution of equations of motion for multibody systems is presented. The present procedure adopts a differential partitioning of the translational motions and the rotational motions. The translational equations of motion are then treated by either a conventional explicit or an implicit direct integration method. A principle feature of this procedure is a nonlinearly implicit algorithm for updating rotations via the Euler four-parameter representation. This procedure is applied to the rolling of a sphere through a specific trajectory, which shows that it yields robust solutions.

Park, K. C.↗

QSPIN: A High Level Java API for Quantum Computing Experimentation

QSPIN is a high level Java language API for experimentation in QC models used in the calculation of Ising spin glass ground states and related quadratic unconstrained binary optimization (QUBO) problems. The Java API is intended to facilitate research in advanced QC algorithms such as hybrid quantum-classical solvers, automatic selection of constraint and optimization parameters, and techniques for the correction and mitigation of model and solution errors. QSPIN includes high level solver objects tailored to the D-Wave quantum annealing architecture that implement hybrid quantum-classical algorithms [Booth et al.] for solving large problems on small quantum devices, elimination of variables via roof duality, and classical computing optimization methods such as GPU accelerated simulated annealing and tabu search for comparison. A test suite of documented NP-complete applications ranging from graph coloring, covering, and partitioning to integer programming and scheduling are provided to demonstrate current capabilities.

Quantu↗

Laser induced photoluminescence spectroscopy of cometary radicals

Flash photolysis together with laser excitation of the product fragments was used in laboratory studies of cometary radicals. The LIPS method has been applied to the CN radical to determine: (1) Radiative lifetimes of individual rotational levels in the zeroth vibrational level of the B state; (2) energy partitioning during photodissociation of C2N2; and (3) vibrational and rotational excitation during formation of CN radicals in the photodissociation of dicyanoacetylene.

Jackson, W. M.↗

Hybrid state vector methods for structural dynamic and aeroelastic boundary value problems

A computational technique is developed that is suitable for performing preliminary design aeroelastic and structural dynamic analyses of large aspect ratio lifting surfaces. The method proves to be quite general and can be adapted to solving various two point boundary value problems. The solution method, which is applicable to both fixed and rotating wing configurations, is based upon a formulation of the structural equilibrium equations in terms of a hybrid state vector containing generalized force and displacement variables. A mixed variational formulation is presented that conveniently yields a useful form for these state vector differential equations. Solutions to these equations are obtained by employing an integrating matrix method. The application of an integrating matrix provides a discretization of the differential equations that only requires solutions of standard linear matrix systems. It is demonstrated that matrix partitioning can be used to reduce the order of the required solutions. Results are presented for several example problems in structural dynamics and aeroelasticity to verify the technique and to demonstrate its use. These problems examine various types of loading and boundary conditions and include aeroelastic analyses of lifting surfaces constructed from anisotropic composite materials.

Lehman, L. L.↗

Flux vector splitting and approximate Newton methods

In the present investigation, the basic approach is employed to view an iterative scheme as Newton's method or as a modified Newton's method. Attention is given to various modified Newton methods which can arise from differencing schemes for the Euler equations. Flux vector splitting is considered as the basic spatial differencing technique. This technique is based on the partition of a flux vector into groups which have certain properties. The Euler equations fluxes can be split into two groups, the first group having a flux Jacobian with all positive eigenvalues, and the second group having a flux Jacobian with all negative eigenvalues. Flux vector splitting based on a velocity-sound speed split is considered along with the use of numerical techniques to analyze nonlinear systems, and the steady Euler equations for quasi-one-dimensional flow in a nozzle. Results are given for steady flows with shocks.

Jespersen, D. C.↗

Microgreens Root Zone/Shoot Zone Partitioned Planting Box

To enable sustainable food production in future human exploration missions, plant growth is being studied by the Space Crop Production Team at KSC. Microgreens are good candidates for food supplements and contain specific nutrients that are lacking in the prepackaged diet, including vitamin C and vitamin K. Because they are densely sown, typical growth methods do not allow the ability to distinguish between the levels of evaporation from the rooting substrate and transpiration from the leaves. With larger plants, the root and shoot zone can be separated to distinguish these fluxes and accurate transpiration measurements of plant canopies are feasible. Furthermore, separation of the root and shoot zone may also be beneficial when harvesting microgreens in microgravity as it may reduce microbial contamination of the edible biomass by the roots, which have high microbial loads. Using a root and shoot separator box when harvesting may help with microbial contamination, but these tiny plants are challenging to handle in microgravity, so harvest management remains an open question. The innovation proposed here is a microgreen root/shoot partitioned planting box, which offers a solution to these challenges with accurate gas exchange measurements and a safe microgreen harvest in low gravity environments. Being able to measure transpiration of a microgreen canopy will be important for modeling plant growth in reduced gravity environments, so the first objective was to develop a planting unit with a seal between the root zone and the shoot zone. With an unsuitable harvesting technique, freshly harvested microgreens may add debris to the cabin, so the second objective was to test different harvesting techniques and management approaches associated with this innovation. These two objectives were pursued in parallel since many goals were the same: develop a planting unit that 1) separates the shoots from the roots, 2) allows acceptable germination rate, and 3) allows for seedlings to emerge and develop. What differed was the need to have a seal, which was only applicable for our gas exchange goal, and the need to have an embedded harvesting mechanism and bagging method which only applied to our harvesting goal. Testing of the various harvesting mechanisms and bagging methods was performed during a series of parabolic flights. All parabolic flight procedures took place inside a rented secondary containment chamber (e.g., glovebox) that was developed by the University of Louisville specifically for experiments involving fluids and other materials that may become airborne during reduced gravity flight. Three different harvesting methods and two different bagging collection methods were tested for microgreens. A third bagging method was initially tested but found to be unsuitable. Human factors were also taken into consideration, to identify which harvest and bagging collection methods would be easiest to use with favorable results in microgravity. Three parabolic flight tests were performed in total, one in November 2021 and two in December 2021.

Gioia Massa↗

Studies in tilt rotor VTOL aircraft aeroelasticity, volume 2

Two methods for natural mode vibration analysis are discussed. The first consists of a direct approach based on a finite element representation of the complete structure as an entity. The mass and stiffness matrices for the complete structure are assembled by properly combining the mass and stiffness matrices of the individual elements into which the structure has been divided. The second approach is that of component mode synthesis. This method is based on the concept of synthesizing the natural modes of the complete structure from modes of conveniently difined substructures, or components, into which the structure has been partitioned. In this way the expedient of reducing the system degrees of freedom, and thus the size of the eigenvalue problem, can be introduced by partial modal synthesis.

Kvaternik, R. G.↗

Progressive content-based retrieval of image and video with adaptive and iterative refinement

A method and apparatus for minimizing the time required to obtain results for a content based query in a data base. More specifically, with this invention, the data base is partitioned into a plurality of groups. Then, a schedule or sequence of groups is assigned to each of the operations of the query, where the schedule represents the order in which an operation of the query will be applied to the groups in the schedule. Each schedule is arranged so that each application of the operation operates on the group which will yield intermediate results that are closest to final results.

Li, Chung-Sheng↗

Electron Energization and Thermal to Non-Thermal Energy Partition During Earth's Magnetotail Reconnection

Electrons in earth's magnetotail are energized significantly both in the form of heating and in the form of acceleration to non-thermal energies. While magnetic reconnection is considered to play an important role in this energization, it still remains unclear how electrons are energized and how energy is partitioned between thermal and non-thermal components. Here, we show, based on in situ observations by NASA's magnetospheric multiscale mission combined with multi-component spectral fitting methods, that the average electron energy ε¯(or equivalently temperature) is substantially higher when the locally averaged electric field magnitude |𝐸| is also higher. While this result is consistent with the classification of “plasma-sheet” and “tail-lobe” reconnection during which reconnection is considered to occur on closed and open magnetic field lines, respectively, it further suggests that a stochastic Fermi acceleration in 3D, reconnection-driven turbulence is essential for the production and confinement of energetic electrons in the reconnection region. The puzzle is that the non-thermal power-law component can be quite small even when the electric field is large and the bulk population is significantly heated. The fraction of non-thermal electron energies varies from sample to sample between ∼20% and ∼60%, regardless of the electric field magnitude. Interestingly, these values of non-thermal fractions are similar to those obtained for the above-the-looptop hard x-ray coronal sources for solar flares.

M Oka↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Energy loss partitioning during ballistic impact of polymer composites

The objective of this study is to determine the energy dissipation processes in polymer-matrix composites during impact of ballistic projectiles. These processes include heat, fiber deformation and breakage, matrix deformation and fracture, and interfacial delamination. In this study, experimental measurements were made, using specialized specimen designs and test methods, to isolate the energy consumed by each of these processes during impact in the ballistic range. Using these experiments, relationships between material parameters and energy dissipation were examined. Composites with the same matrix but reinforced with Kevlar, PE, and graphite fabric were included in this study. These fibers were selected based on the differences in their intrinsic properties. Matrix cracking was found to be one of the most important energy absorption mechanisms during impact, especially in ductile samples such as Spectra-900 PE and Kevlar-49 reinforced polymer. On the contrary, delamination dominated the energy dissipation in brittle composites such as graphite reinforced materials. The contribution from frictional forces was also investigated and the energy partitioning among the different processes evaluated.

Zee, Ralph H.↗

An alternative to Guyan reduction of finite-element models

Structural modeling is a key part of structural system identification for large space structures. Finite-element structural models are commonly used in practice because of their general applicability and availability. The initial models generated by using a standard computer program such as NASTRAN, ANSYS, SUPERB, STARDYNE, STRUDL, etc., generally contain tens of thousands of degrees of freedom. The models must be reduced for purposes of identification. Not only does the magnitude of the identification effort grow exponentially as a function of the number of degrees of freedom, but numerical procedures may also break down because of accumulated round-off errors. Guyan reduction is usually applied after a static condensation. Misapplication of Guyan reduction can lead to serious modeling errors. It is quite unfortunate and disappointing, since the accuracy of the original detailed finite-element model one tries very hard to achieve is lost by the reduction. First, why and how Guyan reduction always causes loss of accuracy is examined. An alternative approach is then introduced. The alternative can be thought of as an improvement of Guyan reduction, the Rayleigh-Ritz method, and in particular the recent algorithm of Wilson, Yuan, and Dickens. Unlike Guyan reduction, the use of the alternative does not need any special insight, experience, or skill for partitioning the structural degrees of freedom. In addition to model condensation, this alternative approach can also be used for predicting analytically, quickly, and economically, what are those structural modes that are excitable by a force actuator at a given trial location. That is, in the excitation of the structural modes for identification, it can be used for guiding the placement of the force actuators.

Lin, Jiguan Gene↗

Flight Guidance System Requirements Specification

This report describes a requirements specification written in the RSML-e language for the mode logic of a Flight Guidance System of a typical regional jet aircraft. This model was created as one of the first steps in a five-year project sponsored by the NASA Langley Research Center, Rockwell Collins Inc., and the Critical Systems Research Group of the University of Minnesota to develop new methods and tools to improve the safety of avionics designs. This model will be used to demonstrate the application of a variety of methods and techniques, including safety analysis of system and subsystem requirements, verification of key properties using theorem provers and model checkers, identification of potential sources mode confusion in system designs, partitioning of applications based on the criticality of system hazards, and autogeneration of avionics quality code. While this model is representative of the mode logic of a typical regional jet aircraft, it does not describe an actual or planned product. Several aspects of a full Flight Guidance System, such as recovery from failed sensors, have been omitted, and no claims are made regarding the accuracy or completeness of this specification.

Miller, Steven P.↗

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗