Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element contact”

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 397 records · Page 22

Higher order finite element analysis of thick composite laminates

A higher order, sub-parametric, laminated, 3D solid finite element was used for the analysis of very thick laminated composite plates. The geometry of this element is defined by four nodes in the X-Y plane which define a prism of material through the thickness of the laminate. There are twenty-four degrees of freedom at each node; translations at the upper and lower surfaces of the laminate in each of the three coordinate directions, and the derivatives of these translations with respect to each coordinate. This choice of degrees of freedom leads to displacement and strain compatibility at the corners. Stacking sequence effects are accounted for by explicitly integrating the strain energy density through the thickness of the element. The laminated solid element was combined with a gap-contact element to analyze thick laminated composite lugs loaded through flexible pins. The resulting model accounts for pin bending effects that produce non-uniform bearing stresses through the thickness of the lug. A thick composite lug experimental test program was performed, and provided data that was used to validate the analytical model. Two lug geometries and three stacking sequences were tested.

Goering, J.↗

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗

Electrostatic Switching in Vertically Oriented Nanotubes for Nonvolatile Memory Applications

We have demonstrated electrostatic switching in vertically oriented nanotubes or nanofibers, where a nanoprobe was used as the actuating electrode inside an SEM. When the nanoprobe was manipulated to be in close proximity to a single tube, switching voltages between 10 V - 40 V were observed, depending on the geometrical parameters. The turn-on transitions appeared to be much sharper than the turn-off transitions which were limited by the tube-to-probe contact resistances. In many cases, stiction forces at these dimensions were dominant, since the tube appeared stuck to the probe even after the voltage returned to 0 V, suggesting that such structures are promising for nonvolatile memory applications. The stiction effects, to some extent, can be adjusted by engineering the switch geometry appropriately. Nanoscale mechanical measurements were also conducted on the tubes using a custom-built anoindentor inside an SEM, from which preliminary material parameters, such as the elastic modulus, were extracted. The mechanical measurements also revealed that the tubes appear to be well adhered to the substrate. The material parameters gathered from the mechanical measurements were then used in developing an electrostatic model of the switch using a commercially available finite-element simulator. The calculated pull-in voltages appeared to be in agreement to the experimentally obtained switching voltages to first order.

plasma-enhanced chemical vapor deposition (PECVD).↗

Computational mechanics needs study

In order to assess the needs in computational mechanics over the next decade, we formulated a questionnaire and contacted computational mechanics researchers and users in industry, government, and academia. As expected, we found a wide variety of computational mechanics usage and research. This report outlines the activity discussed with those contacts, as well as that in our own organizations. It should be noted that most of the contacts were made before the recent decline of the defense industry. Therefore, areas which are strongly defense-oriented may decrease in relative importance. In order to facilitate updating of this study, names of a few key researchers in each area are included as starting points for future literature surveys. These lists of names are not intended to represent those persons doing the best research in that area, nor are they intended to be comprehensive. They are, as previously stated, offered as starting points for future literature searches. Overall, there is currently a broad activity in computational mechanics in this country, with the breadth and depth increasing as more sophisticated software and faster computers become more available. The needs and desires of the workers in this field are as diverse as their background and organizational products. There seems to be some degree of software development in any organization (although the level of activity is highly variable from one organization to another) which has any research component in its mission. It seems, however, that there is considerable use of commercial software in almost all organizations. In most industrial research organizations, it appears that very little actual software development is contracted out, but that most is done in-house, using a mixture of funding sources. Government agencies vary widely in the ratio of in-house to out-house ratio. There is a considerable amount of experimental verification in most, but not all, organizations. Generally, the amount of experimental verification is more than we expected. Of all the survey contacts, one or two believe that the resources they are allocated are sufficient, but most do not. Some believe they have only half the resources they need. Some see their resource deficits as short-term, while others see it as a trend which will continue or perhaps worsen. The pessimism is stronger in the defense and aerospace industry. When considering only the nonlinear development efforts, there appears to be an even mix of geometric and material nonlinearity. There is not much particular emphasis in linear analysis unless it is in extension of current analysis capabilities to larger problems. The primary exception is concern about modeling of composites, where proven methodologies have trailed element and computer hardware development. Most of the people we spoke to use finite element techniques, but there is some finite difference and boundary element work ongoing. There is also some interest in multiple methods. Coupling of finite elements and boundary elements appears to be of high interest, since the two analysis types are complementary.

Griffin, O. Hayden, Jr.↗

Compatibility of Segments of Thermoelectric Generators

A method of calculating (usually for the purpose of maximizing) the power-conversion efficiency of a segmented thermoelectric generator is based on equations derived from the fundamental equations of thermoelectricity. Because it is directly traceable to first principles, the method provides physical explanations in addition to predictions of phenomena involved in segmentation. In comparison with the finite-element method used heretofore to predict (without being able to explain) the behavior of a segmented thermoelectric generator, this method is much simpler to implement in practice: in particular, the efficiency of a segmented thermoelectric generator can be estimated by evaluating equations using only hand-held calculator with this method. In addition, the method provides for determination of cascading ratios. The concept of cascading is illustrated in the figure and the definition of the cascading ratio is defined in the figure caption. An important aspect of the method is its approach to the issue of compatibility among segments, in combination with introduction of the concept of compatibility within a segment. Prior approaches involved the use of only averaged material properties. Two materials in direct contact could be examined for compatibility with each other, but there was no general framework for analysis of compatibility. The present method establishes such a framework. The mathematical derivation of the method begins with the definition of reduced efficiency of a thermoelectric generator as the ratio between (1) its thermal-to-electric power-conversion efficiency and (2) its Carnot efficiency (the maximum efficiency theoretically attainable, given its hot- and cold-side temperatures). The derivation involves calculation of the reduced efficiency of a model thermoelectric generator for which the hot-side temperature is only infinitesimally greater than the cold-side temperature. The derivation includes consideration of the ratio (u) between the electric current and heat-conduction power and leads to the concept of compatibility factor (s) for a given thermoelectric material, defined as the value of u that maximizes the reduced efficiency of the aforementioned model thermoelectric generator.

Snyder, G. Jeffrey↗

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗

Passive heat removal in horizontally oriented micro-HTGRs

There are novel, horizontally oriented high-temperature gas-cooled micro-reactor (HTGR) designs under consideration because of their portability, simplicity, and reliable operation. These features of the micro-HTGR design can meet the fission battery attributes, such as economic, standardize, installed, unattended and reliable. Although there have been HTGR related safety studies, micro-HTGRs when horizontally oriented are expected to exhibit significantly different thermal physics. During the unavailability of helium circulation, the internal reactor core is designed to cool by block-to-block conduction and radiation, and the reactor vessel surface is cooled by the ambient air. This scenario is anticipated during the transport of the micro-HTGR in a shipping container. The conduction and radiation between the prismatic micro-HTGR blocks in the core can be influenced by variances in the thermal contacts. This work investigated the conduction within a simulated horizontal HTGR core. An experimental setup was used to validate a numerical model. The experimental setup consisted of a hexagonal assembly with scaled prismatic blocks placed within a high-temperature vacuum environment. The gaps between the blocks were well controlled and monitored. The experimental setup was designed in such a way that the temperature variation in the axial direction was minimal, such that the experiment could be observed as a 2D (r, ) heat transfer problem. The pressurized tube was equipped with Infra-Red (IR) transparent windows, and an IR camera was used to capture the spatio-temporal temperature evolution in the hexagonal block assembly during the cooldown experiments. The experimental scenario was computationally modeled with a finite element analysis (FEA) program. Once validated, the computational model was used to investigate the impact of gap conductance on overall decay heat removal. Using a conservative estimate for gap conductance value (100 W/m 2 -K) between the prismatic blocks, there is a negligible increase in temperature observed during decay heat generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Extraction of low-energy constants of single- and double-β decays from lattice QCD: A sensitivity analysis

Lattice quantum chromodynamics (LQCD) has the promise of constraining low-energy constants (LECs) of nuclear effective field theories (EFTs) from first-principles calculations that incorporate the dynamics of quarks and gluons. Given the Euclidean and finite-volume nature of LQCD outputs, complex mappings are developed in recent years to obtain the Minkowski and infinite-volume counterparts of LQCD observables. In particular, as LQCD is moving toward computing a set of important few-nucleon matrix elements at the physical values of the quark masses, it is important to investigate whether the anticipated precision of LQCD spectra and matrix elements will be sufficient to guarantee tighter constraints on the relevant LECs than those already obtained from phenomenology, considering the nontrivial mappings involved. With a focus on the leading-order LECs of the pionless EFT, L 1,A and $g$$^{NN}_{v}$, which parametrize, respectively, the strength of the isovector axial two-body current in a single-β decay (and other related processes such pp fusion), and of the isotensor contact two-body operator in the neutrinoless double-β decay within the light neutrino exchange scenario, the expected uncertainty on future extractions of L 1,A and $g$$^{NN}_{v}$ are examined using synthetic data at the physical values of the quark masses. It is observed that achieving small uncertainties in L 1,A will be challenging, and (sub)percent-level precision in the two-nucleon spectra and matrix elements is essential in reducing the uncertainty on this LEC compared to the existing constraints. On the other hand, the short-distance coupling of the neutrinoless double-β decay, $g$$^{NN}_{v}$, is shown to be less sensitive to uncertainties on both LQCD energies and the matrix element, and can likely be constrained with percent-level precision in the upcoming LQCD calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of Creep Properties of Ultra-High-Temperature Materials by a Novel Non-Contact Technique

A non-contact technique for measuring the creep properties of materials has been developed and validated as part of a collaboration among the University of Massachusetts, NASA Marshall Space Flight Center Electrostatic Levitation Facility (ESL), and the University of Tennessee. This novel method has several advantages over conventional creep testing. The sample is deformed by the centripetal acceleration from the rapid rotation, and the deformed shapes are analyzed to determine the strain. Since there is no contact with grips, there is no theoretical maximum temperature and no concern about chemical compatibility. Materials may be tested at the service temperature even for extreme environments such as rocket nozzles, or above the service temperature for accelerated testing of materials for applications such as jet engines or turbopumps for liquid-fueled engines. The creep measurements have been demonstrated to 2400 C with niobium, while the test facility, the NASA MSFC ESL, has processed materials up to 3400 C. Furthermore, the ESL creep method employs a distribution of stress to determine the stress exponent from a single test, versus the many tests required by conventional methods. Determination of the stress exponent from the ESL creep tests requires very precise measurement of the surface shape of the deformed sample for comparison to deformations predicted by finite element models for different stress exponents. An error analysis shows that the stress exponent can be determined to about 1% accuracy with the current methods and apparatus. The creep properties of single-crystal niobium at 1985 C showed excellent agreement with conventional tests performed according to ASTM Standard E-139. Tests on other metals, ceramics, and composites relevant to rocket propulsion and turbine engines are underway.

Hyers, Robert W.↗

A Thermostructural Analysis of a Diboride Composite Leading Edge

In an effort to support the design of zirconium diboride composite leading edges for hypersonic vehicles, a finite element model (FEM) of a prototype leading edge was created and finite element analysis (FEA) was employed to assess its thermal and structural response to aerothermal boundary conditions. Unidirectional material properties for the structural components of the leading edge, a continuous fiber reinforced diboride composite, were computed with COSTAR. These properties agree well with those experimentally measured. To verify the analytical approach taken with COSMOS/M, an independent FEA of one of the leading edge assembly components was also done with COSTAR. Good agreement was obtained between the two codes. Both showed that a unidirectional lay-up had the best margin of safety for a simple loading case. Both located the maximum stress in the same region and ply. The magnitudes agreed within 4 percent. Trajectory based aerothermal heating was then applied to the leading edge assembly FEM created with COSMOS/M to determine steady state temperature response, displacement, stresses, and contact forces due to thermal expansion and thermal strains. Results show that the leading edge stagnation line temperature reached 4700 F. The maximum computed failure index for the laminated composite components peaks at 4.2, and is located at the bolt flange in layer 2 of the side bracket. The temperature gradient in the tip causes a compressive stress of 279 ksi along its width and substantial tensile stresses within its depth.

Kowalski, Tom↗

Simulation System for Training in Laparoscopic Surgery

A computer-based simulation system creates a visual and haptic virtual environment for training a medical practitioner in laparoscopic surgery. Heretofore, it has been common practice to perform training in partial laparoscopic surgical procedures by use of a laparoscopic training box that encloses a pair of laparoscopic tools, objects to be manipulated by the tools, and an endoscopic video camera. However, the surgical procedures simulated by use of a training box are usually poor imitations of the actual ones. The present computer-based system improves training by presenting a more realistic simulated environment to the trainee. The system includes a computer monitor that displays a real-time image of the affected interior region of the patient, showing laparoscopic instruments interacting with organs and tissues, as would be viewed by use of an endoscopic video camera and displayed to a surgeon during a laparoscopic operation. The system also includes laparoscopic tools that the trainee manipulates while observing the image on the computer monitor (see figure). The instrumentation on the tools consists of (1) position and orientation sensors that provide input data for the simulation and (2) actuators that provide force feedback to simulate the contact forces between the tools and tissues. The simulation software includes components that model the geometries of surgical tools, components that model the geometries and physical behaviors of soft tissues, and components that detect collisions between them. Using the measured positions and orientations of the tools, the software detects whether they are in contact with tissues. In the event of contact, the deformations of the tissues and contact forces are computed by use of the geometric and physical models. The image on the computer screen shows tissues deformed accordingly, while the actuators apply the corresponding forces to the distal ends of the tools. For the purpose of demonstration, the system has been set up to simulate the insertion of a flexible catheter in a bile duct. [As thus configured, the system can also be used to simulate other endoscopic procedures (e.g., bronchoscopy and colonoscopy) that include the insertion of flexible tubes into flexible ducts.] A hybrid approach has been followed in developing the software for real-time simulation of the visual and haptic interactions (1) between forceps and the catheter, (2) between the forceps and the duct, and (3) between the catheter and the duct. The deformations of the duct are simulated by finite-element and modalanalysis procedures, using only the most significant vibration modes of the duct for computing deformations and interaction forces. The catheter is modeled as a set of virtual particles uniformly distributed along the center line of the catheter and connected to each other via linear and torsional springs and damping elements. The interactions between the forceps and the duct as well as the catheter are simulated by use of a ray-based haptic-interaction- simulating technique in which the forceps are modeled as connected line segments.

Basdogan, Cagatay↗

Structural, Criticality, and Radiation Dose Calculations to Support SNF Loading into a DOE Standard Canister

The DOE Standard Canister Demonstration Project includes the development of an internal support structure (ISS) for the 4.6-m-long and 45.7-cm-diameter canister. This study presents structural evaluations, criticality safety assessments, and dose rate calculations conducted within the scope of the ISS design process to support smooth canister loading operations and to ensure safe storage, transportation, and disposal of Peach Bottom 1 Core II (PB2) and Fort St. Vrain (FSV) SNF, currently stored at INL’s CPP-603 facility. The ISS includes a 316L stainless-steel basket that holds 12 PB2 SNF rods. The PB2 basket rests on top of six individual, 78.7-cm-long, 316L stainless-steel columns that are equally spaced and welded to the inner canister wall at the lower end of the shell. These columns represent the FSV basket and can hold one FSV SNF element. An A92014 T6 aluminum spacer disc is bolted to the bottom plate of the PB2 basket to vertically restrain the FSV SNF element after the PB2 basket is placed inside the canister above the FSV basket. The structural evaluations of the ISS followed applicable ASME BPVC.III.3 guidelines and included finite element (FE) analyses of the PB2 basket structure; analyses of welds and bolds; buckling analyses of selected components, and acceptability assessments of the expected basket deformations under loading operations. The criticality safety assessments used the Monte Carlo N-Particle (MCNP) software architecture Version 6.2, including ENDF/B-V continuous-energy cross-section libraries, considering intact SNF in a single storage overpack or two multistorage overpack configurations and intact or failed SNF configured for disposal. The dose rate computations are based on source terms taken from the DOE Spent Fuel Database. The isotopic composition was decay corrected for the year 2022 using the ORIGEN module in the SCALE suite. A 19-group photon spectrum and a 27-group neutron-source spectrum were generated and used in MCNP to calculate estimated dose-equivalent rates, both on DOE Standard Canister contact and at a radial distance of 1 m from the canister surface. The results of this study indicate a structurally sound system that can uphold its criticality safety functions throughout its intended operational phases. Furthermore, this study provides confidence that sufficient radiological protection is technically achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Structural, Criticality, and Radiation Dose Calculations to support SNF Loading into a DOE Standard Canister

The DOE Standard Canister Demonstration Project includes the development of an internal support structure (ISS) for the 4.6-m long and 45.7-cm diameter canister. This study presents structural evaluations, criticality safety assessments, and dose rate calculations conducted within the scope of the ISS design process to support smooth canister loading operations and to ensure safe storage, transportation, and disposal of Peach Bottom 1 Core II (PB2) and Fort St. Vrain (FSV) SNF, currently stored at the CPP-603 facility of the INL. The ISS includes a 316L stainless steel basket that holds twelve PB2 SNF rods. Further, six individual, 78.7-cm long, 316L stainless steel columns are equally spaced and welded to the inner canister wall at the lower end of the shell. These columns represent the FSV basket and can hold one FSV SNF element. An A92014 T6 aluminum spacer disc is bolted to the bottom plate of the PB2 basket to vertically restrain the FSV SNF element after the PB2 basket is placed inside the canister on top of the FSV basket. The structural evaluations of the ISS followed applicable ASME BPVC.III.3 guidelines and included finite element (FE) analyses of the PB2 basket structure; analyses of welds and bolds; buckling analyses of selected components, and acceptability assessments of the expected basket deformations under loading operations. The criticality safety assessments used the Monte Carlo N-Particle (MCNP) software architecture version 6.2 including ENDF/B-V continuous energy cross-section libraries, considering intact SNF in a single storage overpack or two multi-storage overpack configurations, and intact or failed SNF configured for disposal. The dose rate computations are based on source terms taken from the DOE Spent Fuel Database. The isotopic composition was decay corrected for the year 2022 using the ORIGEN module in the SCALE suite. A 19-group photon spectrum and 27-group neutron-source spectra were generated and used in MCNP to calculate estimated dose-equivalent rates, both on DOE Standard Canister contact and at a radial distance of from the canister surface. The results of this study indicate a structurally sound system that can uphold its criticality safety functions throughout its intended operational phases. Further, they increase confidence that sufficient radiological protection is technically achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measurement of Photovoltaic Module Deformation Dynamics During Hail Impact Using Digital Image Correlation

Stereo high-speed video of photovoltaic modules undergoing laboratory hail tests was processed using digital image correlation to determine module surface deformation during and immediately following impact. The purpose of this work was to demonstrate a methodology for characterizing module impact response differences as a function of construction and incident hail parameters. Video capture and digital image analysis were able to capture out-of-plane module deformation to a resolution of ±0.1 mm at 11 kHz on an in-plane grid of 10 × 10 mm over the area of a 1 × 2 m commercial photovoltaic module. With lighting and optical adjustments, the technique was adaptable to arbitrary module designs, including size, backsheet color, and cell interconnection. Furthermore, impacts were observed to produce an initially localized dimple in the glass surface, with peak deflection proportional to the square root of incident energy. Subsequent deformation propagation and dissipation were also captured, along with behavior for instances when the module glass fractured. Natural frequencies of the module were identifiable by analyzing module oscillations postimpact. Limitations of the measurement technique were that the impacting ice ball obscured the data field immediately surrounding the point of contact, and both ice and glass fracture events occurred within 100 μs, which was not resolvable at the chosen frame rate. Increasing the frame rate and visualizing the back surface of the impact could be applied to avoid these issues. Applications for these data include validating computational models for hail impacts, identifying the natural frequencies of a module, and identifying damage initiation mechanisms.

14 SOLAR ENERGY↗

Interface Problem Formulation Improvements with Application to Nuclear Fuel Performance Analysis

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation Program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This report documents new developments and robustness improvements in mechanical and thermal (gap heat transfer) contact formulations. The improvements range from the migration of industrial level (“assessment”) nuclear fuel model setups to the usage of mortar formulations, the addition of frictional contact to one-dimensional layered representations of fuel and cladding components, and the addition of the Petrov-Galerkin approach to dual mortar, which improves performance on curved, relatively coarse meshes. In addition, the Lagrange-multiplier enforcement of mechanical mortar contact constraints has been extended to two additional types of enforcement: penalty and augmented Lagrange-Uzawa. We show that the latter approach yields the same interface results as dual mortar in the Multiphysics Object-Oriented Simulation Environment, with the advantage of not worsening the condition number of the system matrix—thereby enabling the use of some general implementations of iterative preconditioners, at the expense of additional system evaluations (i.e., augmentations).

42 ENGINEERING↗

Estimating Hardness from the USDC Tool-Bit Temperature Rise

A method of real-time quantification of the hardness of a rock or similar material involves measurement of the temperature, as a function of time, of the tool bit of an ultrasonic/sonic drill corer (USDC) that is being used to drill into the material. The method is based on the idea that, other things being about equal, the rate of rise of temperature and the maximum temperature reached during drilling increase with the hardness of the drilled material. In this method, the temperature is measured by means of a thermocouple embedded in the USDC tool bit near the drilling tip. The hardness of the drilled material can then be determined through correlation of the temperature-rise-versus-time data with time-dependent temperature rises determined in finite-element simulations of, and/or experiments on, drilling at various known rates of advance or known power levels through materials of known hardness. The figure presents an example of empirical temperature-versus-time data for a particular 3.6-mm USDC bit, driven at an average power somewhat below 40 W, drilling through materials of various hardness levels. The temperature readings from within a USDC tool bit can also be used for purposes other than estimating the hardness of the drilled material. For example, they can be especially useful as feedback to control the driving power to prevent thermal damage to the drilled material, the drill bit, or both. In the case of drilling through ice, the temperature readings could be used as a guide to maintaining sufficient drive power to prevent jamming of the drill by preventing refreezing of melted ice in contact with the drill.

Bar-Cohen, Yoseph↗

Design Considerations for Lightweight Space Radiators Based on Fabrication and Test Experience With a Carbon-Carbon Composite Prototype Heat Pipe

This report discusses the design implications for spacecraft radiators made possible by the successful fabrication and proof-of-concept testing of a graphite-fiber-carbon-matrix composite (i.e., carbon-carbon (C-C)) heat pipe. The prototype heat pipe, or space radiator element, consists of a C-C composite shell with integrally woven fins. It has a thin-walled furnace-brazed metallic (Nb-1%Zr) liner with end caps for containment of the potassium working fluid. A short extension of this liner, at increased wall thickness beyond the C-C shell, forms the heat pipe evaporator section which is in thermal contact with the radiator fluid that needs to be cooled. From geometric and thermal transport properties of the C-C composite heat pipe tested, a specific radiator mass of 1.45 kg/sq m can be derived. This is less than one-fourth the specific mass of present day satellite radiators. The report also discusses the advantage of segmented space radiator designs utilizing heat pipe elements, or segments, in their survivability to micrometeoroid damage. This survivability is further raised by the use of condenser sections with attached fins, which also improve the radiation heat transfer rate. Since the problem of heat radiation from a fin does not lend itself to a closed analytical solution, a derivation of the governing differential equation and boundary conditions is given in appendix A, along with solutions for rectangular and parabolic fin profile geometries obtained by use of a finite difference computer code written by the author.

Juhasz, Albert J.↗

Design Considerations for Lightweight Space Radiators Based on Fabrication and Test Experience with a Carbon-Carbon Composite Prototype Heat Pipe

This report discusses the design implications for spacecraft radiators made possible by the successful fabrication and Proof-of-concept testing of a graphite-fiber-carbon-matrix composite (i.e., carbon-carbon (C-C)) heat pipe. The proto-type heat pipe, or space radiator element, consists of a C-C composite shell with integrally woven fins. It has a thin-walled furnace-brazed metallic (Nb-1%Zr) liner with end caps for containment of the potassium working fluid. A short extension of this liner, at increased wall thickness beyond the C-C shell, forms the heat pipe evaporator section which is in thermal contact with the radiator fluid that needs to be cooled. From geometric and thermal transport properties of the C-C composite heat pipe tested, a specific radiator mass of 1.45 kg/m2 can be derived. This is less than one-fourth the specific mass of present day satellite radiators. The report also discusses the advantage of segmented space radiator designs utilizing heat pipe elements, or segments, in their survivability to micro-meteoroid damage. This survivability is further raised by the use of condenser sections with attached fins, which also improve the radiation heat transfer rate. Since the problem of heat radiation from a fin does not lend itself to a closed analytical solution, a derivation of the governing differential equation and boundary conditions is given in appendix A, along with solutions for rectangular and parabolic fin profile geometries obtained by use of a finite difference computer code written by the author.

Juhasz, Albert J.↗