Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Contact enforcement”

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.

35 records · Page 2

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗

A mortar thermomechanical contact computational framework for nuclear fuel performance simulation

Nuclear fuel performance simulations involve the modeling of complex physical phenomena, ranging from fission gas release to fuel swelling and other temperature-induced effects. For light-water reactors (LWRs), swelling of the fuel and the pressure it imposes on the clad when they come into contact causes permanent clad deformation. Accurately characterizing the fuel-cladding interaction, which involves multiple physics, is essential to accurately simulate the fuel/cladding system. Thermomechanical modeling of this problem using a variationally consistent enforcement (e.g., a mortar approach) has been shown to improve the quality of results and facilitate convergence. Here, we present a general multiphysics computational framework for solving nuclear fuel problems using a mortar approach in BISON, a nuclear fuel performance code. In this study analyses show that using the mortar approach, which enables variationally consistent constraint enforcement, improves the quality of results as compared to the more commonly used node-on-face enforcement for representative LWR nuclear fuel simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Frictionless Contact of Multilayered Composite Half Planes Containing Layers With Complex Eigenvalues

A previously developed local-global stiffness matrix methodology for the response of a composite half plane, arbitrarily layered with isotropic, orthotropic or monoclinic plies, to indentation by a rigid parabolic punch is further extended to accommodate the presence of layers with complex eigenvalues (e.g., honeycomb or piezoelectric layers). First, a generalized plane deformation solution for the displacement field in an orthotropic layer or half plane characterized by complex eigenvalues is obtained using Fourier transforms. A local stiffness matrix in the transform domain is subsequently constructed for this class of layers and half planes, which is then assembled into a global stiffness matrix for the entire multilayered half plane by enforcing continuity conditions along the interfaces. Application of the mixed boundary condition on the top surface of the half plane indented by a rigid punch results in an integral equation for the unknown pressure in the contact region. The integral possesses a divergent kernel which is decomposed into Cauchy-type and regular parts using the asymptotic properties of the local stiffness matrix and a relationship between Fourier and finite Hilbert transform of the contact pressure. The solution of the resulting singular integral equation is obtained using a collocation technique based on the properties of orthogonal polynomials developed by Erdogan and Gupta. Examples are presented that illustrate the important influence of low transverse properties of layers with complex eigenvalues, such as those exhibited by honeycomb, on the load versus contact length response and contact pressure distributions for half planes containing typical composite materials.

Zhang, Wang↗

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗

Maintenance cost study of rotary wing aircraft, phase 2

The Navy's maintenance and materials management data base was used in a study to determine the feasibility of predicting unscheduled maintenance costs for the dynamic systems of military rotary wing aircraft. The major operational and design variables were identified and the direct maintenance man hours per flight hour were obtained by step-wise multiple regression analysis. Five nonmilitary helicopter users were contacted to supply data on which variables were important factors in civil applications. These uses included offshore oil exploration and support, police and fire department rescue and enforcement, logging and heavy equipment movement, and U.S. Army military operations. The equations developed were highly effective in predicting unscheduled direct maintenance man hours per flying hours for military aircraft, but less effective for commercial or public service helicopters, probably because of the longer mission durations and the much higher utilization of civil users.

Source record↗

Automated and Adaptive Mission Planning for Orbital Express

The Orbital Express space mission was a Defense Advanced Research Projects Agency (DARPA) lead demonstration of on-orbit satellite servicing scenarios, autonomous rendezvous, fluid transfers of hydrazine propellant, and robotic arm transfers of Orbital Replacement Unit (ORU) components. Boeing's Autonomous Space Transport Robotic Operations (ASTRO) vehicle provided the servicing to the Ball Aerospace's Next Generation Serviceable Satellite (NextSat) client. For communication opportunities, operations used the high-bandwidth ground-based Air Force Satellite Control Network (AFSCN) along with the relatively low-bandwidth GEO-Synchronous space-borne Tracking and Data Relay Satellite System (TDRSS) network. Mission operations were conducted out of the RDT&E Support Complex (RSC) at the Kirtland Air Force Base in New Mexico. All mission objectives were met successfully: The first of several autonomous rendezvous was demonstrated on May 5, 2007; autonomous free-flyer capture was demonstrated on June 22, 2007; the fluid and ORU transfers throughout the mission were successful. Planning operations for the mission were conducted by a team of personnel including Flight Directors, who were responsible for verifying the steps and contacts within the procedures, the Rendezvous Planners who would compute the locations and visibilities of the spacecraft, the Scenario Resource Planners (SRPs), who were concerned with assignment of communications windows, monitoring of resources, and sending commands to the ASTRO spacecraft, and the Mission planners who would interface with the real-time operations environment, process planning products and coordinate activities with the SRP. The SRP position was staffed by JPL personnel who used the Automated Scheduling and Planning ENvironment (ASPEN) to model and enforce mission and satellite constraints. The lifecycle of a plan began three weeks outside its execution on-board. During the planning timeframe, many aspects could change the plan, causing the need for re-planning. These variable factors, ranging from shifting contact times to ground-station closures and required maintenance times, are discussed along with the flexibility of the ASPEN tool to accommodate changes to procedures and the daily or long-range plan, which contributed to the success of the mission. This paper will present an introduction to ASPEN, a more in-depth discussion on its use on the Orbital Express mission, and other relative work. A description of ground operations after the SRP deliveries were made is included, and we briefly discuss lessons learned from the planning perspective and future work.

scheduling↗

A Finite Element Approach for Simplified 2D Nonlinear Dynamic Contact/Impact Analysis

In this paper, a simplified numerical approach for finite element dynamic analysis of an inelastic solid structure subjected to solid object impact is presented. The approach approximates the impacting solid as the selected multiple nodes, for which mass of the impactor is distributed. The node-to-segment contact formulation with the penalty constraint technique incorporated is employed to impose contact conditions between the nodes and the surface of the receiver structure. The node-to-segment algorithm is integrated into Newton–Raphson time integration scheme and the Lagrange multiplier technique is applied to enforce the identical displacements for the selected nodes throughout the analysis process. The approach is verified using two-dimensional plane strain models considering elastic-perfectly-plastic material behavior. The results obtained using the proposed approach are in a good agreement with those simulated using a commercial finite element code, ABAQUS dynamic/implicit, in terms of displacements and stress distribution fields. The proposed approach is shown to be computationally superior to general finite element method-based contact/impact analysis without significantly sacrificing the accuracy.

Seungwook Seok↗

Embedding space approach to Lorentzian CFT amplitudes and causal spherical functions

Conformal field theory in a Minkowski setting is discussed in an embedding space approach, paying special attention to causality constraints for four-point amplitudes. The physics of dilatation and Lorentz boost is emphasized in specifying the noncompact maximal Abelian subgroup of S O ( d , 2 ) . Reduction of a conformal field theory four-point amplitudes as functions of cross ratios is shown to be equivalent to enforcing H bi-invariance, i.e., F ( h g h ′ ) = F ( g ) , with g ∈ S O ( d , 2 ) and H an appropriate subgroup. Causality is imposed by introducing appropriate semigroups. Causal zonal spherical functions are constructed, making contact with Minkowski conformal blocks introduced previously. Published by the American Physical Society 2024

Agarwal, Pulkit (ORCID:0000000346581691)↗

A Flexible Multibody Approach to Space Launch System Liftoff Pad Separation with Umbilical Disconnect

A flexible multibody dynamic framework is applied to the Space Launch System (SLS) liftoff Coupled Loads Analysis (CLA), enabling computational efficiencies and systematic inclusion of components and various nonlinearities. The objective of this work is to simulate the SLS liftoff transient pad separation event inclusive of umbilical disconnects. The vehicle stacking preloads and cryogenically induced preloads, due to the fueling of the Core Stage (CS), are integrated via specialized deformed coupling procedure to ensure full capture of the liftoff transient event. The Vehicle Stabilizer System (VSS) hydraulic struts are nonlinear and provide amplitude dependent damping, significantly contributing to the liftoff transient event. Ground Extensible Columns (EC) provide additional stiffness and preloads on the Mobile Launcher (ML), also significantly contributing to the liftoff transient. An advanced Henkel-Mar algorithm, with separation and re-contact, thus tracking the potential post separation re-contact between the booster aft skirt and ML, the umbilicals and vehicle, and ground EC separation. For this system level nonlinear dynamic simulation, the flexible multibody framework provides a systematic approach for the addition of component specific algorithms to enforce nonlinearities. This multibody framework is similar to an object-oriented programming language integrating classes of associated subroutines. As such, the complexities associated with the nonlinear CLA are greatly reduced, facilitating simplified bookkeeping as well as accelerating the execution of the nonlinear time-domain simulations.

Joel W Sills Jr.↗

Improved Fuel Particle and Compact Matrix Fission Product Release Prediction with Improved Mass Transfer Model

This study was motivated by the need to incorporate more physics-based mass transfer models in BISON in order to improve fuel particle and compact matrix fission product release prediction. This document details the integration of new modeling capabilities in BISON, including (1) development of sorption mass transfer models; improvement of thermal contact model by adding neon to the gas gap inventory, (3) implementation of PCEA and IG-110 graphite thermal models, (4) validation using Advanced Gas Reactor (AGR)-3/4 compact and capsule rings data. These new capabilities have been shown to enforce the desired physics with satisfactory accuracy. BISON’s predictions of fission product release of AGR-3/4 compacts compare favorably with PARFUME and the concentration profiles across capsule rings show reasonable trend to the experiment measurement.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Developing nanometer scale, atomically precise metallo-catalysts with Molecular Lego

We propose to create atomically precise, highly robust, nanometer scale macromolecules that organize biomimetic metal binding groups in designed three-dimensional pockets to accelerate specific polymerization reactions. These atomically precise catalysts will create atomically precise polymers with improved mechanical and environmental properties from renewable building blocks. We seek to create and study nanometer scale catalysts that mimic metallo-enzymes in their activity and selectivity because these larger catalysts will be able to make more extensive non-covalent contacts with the transition state of the growing polymer. These catalysts will have extended lifetimes because they will better protect the metals and higher activity because additional reactive groups can be brought in close contact with the active site to enhance reactivity. They will be far more tolerant to extremes of temperature and will not denature in non-aqueous solvents because they are held together by multiple, strong, covalent bonds. Our approach to catalysts could be used to develop catalysts to replace noble metal-based catalysts when targeting other reactions, because the geometry of the active site enforced by the Molecular Lego scaffolding will access new reactivity, mimicking how nature uses earth abundant metals for much of its catalysis. We will initially target Lewis acid-based catalysts that assemble aliphatic polyesters with improved stereoselectivity, tacticity and complex alternation of monomers. Aliphatic polyesters have received growing attention as attractive, environmentally benign, and sustainable alternatives to polymers developed from petroleum feedstocks. Aliphatic polyesters are environmentally friendly because they undergo facile hydrolytic degradation to benign products, and they have excellent properties and high biocompatibility. We will assemble these catalysts using unique “Spiroligomers” (aka Molecular Lego) chemistry developed in the Schafmeister group combined with the organometallic chemistry expertise of the Dobereiner group. Polymers are an attractive synthetic target of nanoscale atomically precise catalysts because their resulting structure, stereochemistry and tacticity has profound impact on their properties (crystalline vs amorphous) and provides a readout of the catalytic mechanism. Polymer synthesis also serves as an example of “atomically precise manufacturing” where atomically precise nanoscale catalysts construct polymers with atomic precision at increasing length scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Object impedance control for cooperative manipulation - Theory and experimental results

The dynamic control module being developed in the Dynamic and Strategic Control of Cooperative Manipulators (DASCCOM) project at the Stanford University Aerospace Robotics Laboratory is described. First, the cooperative manipulation problem is analyzed from a systems perspective, and the desirable features of a control system for cooperative manipulation are discussed. Next, a control policy is developed that enforces a controlled impedance not of the individual arm endpoints, but of the manipulated object itself. A parallel implementation for a multiprocessor system is presented. The controller fully compensates for the system dynamics and directly controls the object internal forces. Most importantly, it presents a simple, powerful, intuitive interface to the strategic controller. Experimental results for a dual two-link arm robotic system are presented to verify the controllers performance, for both free-motion slews and environmental contact.

Schneider, Stanley A.↗

Object impedance control for cooperative manipulation - Theory and experimental results

This paper presents the dynamic control module of the Dynamic and Strategic Control of Cooperating Manipulators (DASCCOM) project at Stanford University's Aerospace Robotics Laboratory. First, the cooperative manipulation problem is analyzed from a systems perspective, and the desirable features of a control system for cooperative manipulation are discussed. Next, a control policy is developed that enforces a controlled impedance not of the individual arm endpoints, but of the manipulated object itself. A parallel implementation for a multiprocessor system is presented. The controller fully compensates for the system dynamics and directly controls the object internal forces. Most importantly, it presents a simple, powerful, intuitive interface to higher level strategic control modules. Experimental results from a dual two-link-arm robotic system are used to compare the object impedance controller with other strategies, both for free-motion slews and environmental contact.

Schneider, Stanley A.↗

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗

"Electrostructural Phase Changes" In Charged Particulate Clouds: Planetary and Astrophysical Implications

There is empirical evidence that freely-suspended triboelectrostatically charged particulate clouds of dielectric materials undergo rapid conversion from (nominally) monodispersed "aerosols" to a system of well-defined grain aggregates after grain motion or fluid turbulence ceases within the cloud. In United States Microgravity Laboratory Space Shuttle experiments USML-1 and USML-2, it was found that ballistically-energized grain dispersions would rapidly convert into populations of filamentary aggregates after natural fluid (air) damping of grain motion. Unless continuously disrupted mechanically, it was impossible to maintain a non-aggregated state for the grain clouds of sand-size materials. Similarly, ground- based experiments with very fine dust-size material produced the same results: rapid, impulsive "collapse" of the dispersed grains into well-defined filamentary structures. In both ground-based and microgravity experiments, the chains or filaments were created by long-range dipole electrostatic forces and dipole-induced dielectric interactions, not by monopole interactions. Maintenance of the structures was assisted by short-range static boundary adhesion forces and van der Waals interactions. When the aggregate containers in the USML experiments were disturbed after aggregate formation, the quiescently disposed filaments would rearrange themselves into fractal bundles and tighter clusters as a result of enforced encounters with one another. The long-range dipole interactions that bring the grains together into aggregates are a product of randomly-distributed monopole charges on the grain surfaces. In computer simulations, it has been shown that when the force vectors of all the random charges (of both sign) on a grain are resolved mathematically by assuming Coulombic interaction between them, the net result is a dipole moment on individual grains, even though the grains are electrically neutral insofar as there is no predominance, on their surface, of one charge sign over another. The random charges of both sign derive from natural grain-to-grain interactions that produce triboelectrification via charge exchange every time grain surfaces make contact with one another. The conversion from a random distribution of grains (upon which there are randomly distributed charges) into an organization of electrostatically-ordered aggregates, can be regarded (within the framework of granular-material science) as an "electrical or Coulombic phase change" of the particulate cloud. It is not totally dissimilar from the more normal phase-change concept in which, for example, a gas with long free-path-molecules suddenly becomes a solid as a result of structural ordering of the molecules (notably, also the result of electronic forces, albeit at a different scale). In both the gas-to-solid case, and the aerosol-to-aggregate case, the same materials and charges are present before and after the phase change, but their arrangement now has a higher degree of order and a lower-energy configuration. An input of energy into the system is required to reverse the situation. The aggregates in the USML experiments were observed to undergo at least two phase changes as noted above. The point about phase changes, and by implication, the "electrostructural" reorganizations in particulate clouds, is the following: (a) they can occur very rapidly, almost spontaneously, above a critical cloud density, (b) in going from a higher energy state to a lower energy state, they convert to a denser system, (c) energy must be required to reverse the situation, implying that energy is released during the high-to-low energy phase change. In applying this information to natural particulate clouds, some inferences can be made (it is stressed that reference is still to dielectric materials attracted by dipole forces). There are several natural settings to which the USML observations apply, and to which the phase-change implications likewise apply. Dense clouds of triboelectrically-charged, kinetically-energized grains are to be found in volcanic eruptions (particularly on earth), aeolian dust storms (particularly on Mars), meteorite impact ejecta curtains (on all planets), in "immature" debris rings around planets (e.g., that from which our own Moon may have condensed), and in gravitationally collapsing protoplanetary dust/planetesimal debris disks where dielectric granules are being increasingly brought into collisional relationships with one another (increasing both charge exchange and physical proximity). It is noted that in many of these cases, the degree of electrical charging on the grains is likely to be much higher than that in the USML experiments where charging was not enhanced above the "normal", naturally encountered level for the particular materials (quartz, glass, and various silicate minerals). Application of the phase-change concept suggests that volcanic, aeolian, and impact debris clouds may, under certain circumstances, undergo rapid, impulsive, or even catastrophic collapse into a denser state that will lead to rapid precipitation or fall-out of suspended particulates. Although this idea has been suggested previously by the author , the phase change concept possibly permits some new insights into cloud-system behavior. For example, in a protoplanetary debris disk, the work of gravity may suddenly be enhanced by electrostatically driven collapse of the system when materials reach a critical intergranular spacing or grain number density. This might reduce the rate of planet formation by orders of magnitude (considering the ratio of g-forces to electrostatic ones for very small grains in close proximity), and indeed, it might drive the collapse into a system configuration that would perhaps not be created by gravity alone. Additional information is contained in the original.

Marshall, J. R.↗

The CE/SE Method: a CFD Framework for the Challenges of the New Millennium

The space-time conservation element and solution element (CE/SE) method, which was originated and is continuously being developed at NASA Glenn Research Center, is a high-resolution, genuinely multidimensional and unstructured-mesh compatible numerical method for solving conservation laws. Since its inception in 1991, the CE/SE method has been used to obtain highly accurate numerical solutions for 1D, 2D and 3D flow problems involving shocks, contact discontinuities, acoustic waves, vortices, shock/acoustic waves/vortices interactions, shock/boundary layers interactions and chemical reactions. Without the aid of preconditioning or other special techniques, it has been applied to both steady and unsteady flows with speeds ranging from Mach number = 0.00288 to 10. In addition, the method has unique features that allow for (i) the use of very simple non-reflecting boundary conditions, and (ii) a unified wall boundary treatment for viscous and inviscid flows. The CE/SE method was developed with the conviction that, with a solid foundation in physics, a robust, coherent and accurate numerical framework can be built without involving overly complex mathematics. As a result, the method was constructed using a set of design principles that facilitate simplicity, robustness and accuracy. The most important among them are: (i) enforcing both local and global flux conservation in space and time, with flux evaluation at an interface being an integral part of the solution procedure and requiring no interpolation or extrapolation; (ii) unifying space and time and treating them as a single entity; and (iii) requiring that a numerical scheme be built from a nondissipative core scheme such that the numerical dissipation can be effectively controlled and, as a result, will not overwhelm the physical dissipation. Part I of the workshop will be devoted to a discussion of these principles along with a description of how the ID, 2D and 3D CE/SE schemes are constructed. In Part II, various applications of the CE/SE method, particularly those involving chemical reactions and acoustics, will be presented. The workshop will be concluded with a sketch of the future research directions.

Chang, Sin-Chung↗

nBn Infrared Detector Containing Graded Absorption Layer

It has been proposed to modify the basic structure of an nBn infrared photodetector so that a plain electron-donor- type (n-type) semiconductor contact layer would be replaced by a graded n-type III V alloy semiconductor layer (i.e., ternary or quarternary) with appropriate doping gradient. The abbreviation nBn refers to one aspect of the unmodified basic device structure: There is an electron-barrier ("B" ) layer between two n-type ("n" ) layers, as shown in the upper part of the figure. One of the n-type layers is the aforementioned photon-absorption layer; the other n-type layer, denoted the contact layer, collects the photocurrent. The basic unmodified device structure utilizes minority-charge-carrier conduction, such that, for reasons too complex to explain within the space available for this article, the dark current at a given temperature can be orders of magnitude lower (and, consequently, signal-to-noise ratios can be greater) than in infrared detectors of other types. Thus, to obtain a given level of performance, less cooling (and, consequently, less cooling equipment and less cooling power) is needed. [In principle, one could obtain the same advantages by means of a structure that would be called pBp because it would include a barrier layer between two electron-acceptor- type (p-type) layers.] The proposed modifications could make it practical to utilize nBn photodetectors in conjunction with readily available, compact thermoelectric coolers in diverse infrared- imaging applications that could include planetary exploration, industrial quality control, monitoring pollution, firefighting, law enforcement, and medical diagnosis.

Gunapala, Sarath D.↗