Search NASA⌕ Search

SEARCH · Search NASA

Results for “deformation function”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

The Influence of C and Si on the Flow Behavior of NiAl Single Crystals

Alloys based on the intermetallic compound NiAl are considered potential replacements for Ni and Co-based superalloys in high temperature structural applications due to their excellent oxidation resistance, low densities, high thermal conductivities, and increased melting points. Unfortunately, NiAl exhibits low tensile ductility at room temperature and low strengths at elevated temperatures which have combined to hinder its development. Recent efforts, have revealed that NiAl in the presence of sufficient solute levels, is subject to the phenomenon of strain aging which manifests itself as: sharp yield points, abnormally low strain rate sensitivities (SRS), plateaus or peaks in yield stress and work hardening rate as a function of temperature, flow stress transients upon an upward change in strain rate, reduced tensile elongations at elevated temperatures, and serrated stress-strain curves. Though recent efforts via either alloying or the removal of interstitial impurities, have resulted in consistent room-temperature tensile elongations exceeding 5% and the elimination of serrated flow, the effects of particular substitutional and interstitial elements and the mechanisms by which they might enhance or hinder the mechanical properties remain unknown. Consequently, the purpose of the present paper is to provide a preliminary assessment of the influence of common substitutional and interstitial impurities on the deformation behavior of NiAl. To accomplish this goal a series of NiAl single crystal alloys containing various interstitial solutes were prepared and their mechanical properties were evaluated between 77 and 1100 K. Because Si is a common impurity in conventional purity single crystals grown by the Bridgman method, Si concentrations were also varied in order to determine the influence of this element.

Weaver, M. L.↗
NASA NTRS1983

Advances and trends in structural and solid mechanics; Proceedings of the Symposium, Washington, DC, October 4-7, 1982

The mechanics of materials and material characterization are considered, taking into account micromechanics, the behavior of steel structures at elevated temperatures, and an anisotropic plasticity model for inelastic multiaxial cyclic deformation. Other topics explored are related to advances and trends in finite element technology, classical analytical techniques and their computer implementation, interactive computing and computational strategies for nonlinear problems, advances and trends in numerical analysis, database management systems and CAD/CAM, space structures and vehicle crashworthiness, beams, plates and fibrous composite structures, design-oriented analysis, artificial intelligence and optimization, contact problems, random waves, and lifetime prediction. Earthquake-resistant structures and other advanced structural applications are also discussed, giving attention to cumulative damage in steel structures subjected to earthquake ground motions, and a mixed domain analysis of nuclear containment structures using impulse functions.

Noor, A. K.↗
NASA NTRS2014

Calibration and Validation of a Finite ELement Model of THor-K Anthropomorphic Test Device for Aerospace Safety Applications

The THOR anthropomorphic test device (ATD) has been developed and continuously improved by the National Highway Traffic Safety Administration to provide automotive manufacturers an advanced tool that can be used to assess the injury risk of vehicle occupants in crash tests. Recently, a series of modifications were completed to improve the biofidelity of THOR ATD [1]. The updated THOR Modification Kit (THOR-K) ATD was employed at Wright-Patterson Air Base in 22 impact tests in three configurations: vertical, lateral, and spinal [2]. Although a computational finite element (FE) model of the THOR had been previously developed [3], updates to the model were needed to incorporate the recent changes in the modification kit. The main goal of this study was to develop and validate a FE model of the THOR-K ATD. The CAD drawings of the THOR-K ATD were reviewed and FE models were developed for the updated parts. For example, the head-skin geometry was found to change significantly, so its model was re-meshed (Fig. 1a). A protocol was developed to calibrate each component identified as key to the kinematic and kinetic response of the THOR-K head/neck ATD FE model (Fig. 1b). The available ATD tests were divided in two groups: a) calibration tests where the unknown material parameters of deformable parts (e.g., head skin, pelvis foam) were optimized to match the data and b) validation tests where the model response was only compared with test data by calculating their score using CORrelation and Analysis (CORA) rating system. Finally, the whole ATD model was validated under horizontal-, vertical-, and lateral-loading conditions against data recorded in the Wright Patterson tests [2]. Overall, the final THOR-K ATD model developed in this study is shown to respond similarly to the ATD in all validation tests. This good performance indicates that the optimization performed during calibration by using the CORA score as objective function is not test specific. Therefore confidence is provided in the ATD model for uses in predicting response in test conditions not performed in this study such those observed in the spacecraft landing. Comparison studies with ATD and human models may also be performed to contribute to future changes in THOR ATD design in an effort to improve its biofidelity, which has been traditionally based on post-mortem human subject testing and designer experience.

Putnam, J. B.↗
NASA NTRS1984

High temperature ceramic interface study

Monolithic SiC and Si3N4 are susceptible to contact stress damage at static and sliding interfaces. Transformation-toughened zirconia (TTZ) was evaluated under sliding contact conditions to determine if the higher material fracture toughness would reduce the susceptibility to contact stress damage. Contact stress tests were conducted on four commercially available TTZ materials at normal loads ranging from 0.455 to 22.7 kg (1 to 50 pounds) at temperatures ranging from room temperature to 1204C (2200 F). Static and dynamic friction were measured as a function of temperature. Flexural strength measurements after these tests determined that the contact stress exposure did not reduce the strength of TTZ at contact loads of 0.455, 4.55, and 11.3 kg (1, 10, and 25 pounds). Prior testing with the lower toughness SiC and Si3N4 materials resulted in a substantial strength reduction at loads of only 4.55 and 11.3 kg (10 and 25 pounds). An increase in material toughness appears to improve ceramic material resistance to contact stress damage. Baseline material flexure strength was established and the stress rupture capability of TTZ was evaluated. Stress rupture tests determined that TTZ materials are susceptible to deformation due to creep and that aging of TTZ materials at elevated temperatures results in a reduction of material strength.

Lindberg, L. J.↗
NASA NTRS2002

A Nonlinear Viscoelastic Model for Ceramics at High Temperatures

High-temperature creep behavior of ceramics is characterized by nonlinear time-dependent responses, asymmetric behavior in tension and compression, and nucleation and coalescence of voids leading to creep rupture. Moreover, creep rupture experiments show considerable scatter or randomness in fatigue lives of nominally equal specimens. To capture the nonlinear, asymmetric time-dependent behavior, the standard linear viscoelastic solid model is modified. Nonlinearity and asymmetry are introduced in the volumetric components by using a nonlinear function similar to a hyperbolic sine function but modified to model asymmetry. The nonlinear viscoelastic model is implemented in an ABAQUS user material subroutine. To model the random formation and coalescence of voids, each element is assigned a failure strain sampled from a lognormal distribution. An element is deleted when its volumetric strain exceeds its failure strain. Element deletion has been implemented within ABAQUS. Temporal increases in strains produce a sequential loss of elements (a model for void nucleation and growth), which in turn leads to failure. Nonlinear viscoelastic model parameters are determined from uniaxial tensile and compressive creep experiments on silicon nitride. The model is then used to predict the deformation of four-point bending and ball-on-ring specimens. Simulation is used to predict statistical moments of creep rupture lives. Numerical simulation results compare well with results of experiments of four-point bending specimens. The analytical model is intended to be used to predict the creep rupture lives of ceramic parts in arbitrary stress conditions.

Powers, Lynn M.↗
NASA NTRS1999

Paramagnetic Liquid Bridge in a Gravity-Compensating Magnetic Field

Magnetic levitation was used to stabilize cylindrical columns of a paramagnetic liquid in air between two solid supports. The maximum achievable length to diameter ratio R(sub max) was approx. (3.10 +/- 0.07), very close to the Rayleigh-Plateau limit of pi. For smaller R, the stability of the column was measured as a function of the Bond number, which could be continuously varied by adjusting the strength of the magnetic field. Liquid bridges supported by two solid surfaces have been attracting scientific attention since the time of Rayleigh and Plateau. For a cylindrical bridge of length L and diameter d, it was shown theoretically that in zero gravity the maximum slenderness ratio R (identically = L/d) is pi. The stability and ultimate collapse of such bridges is of interest because of their importance in a number of industrial processes and their potential for low gravity applications. In the presence of gravity, however, the cylindrical shape of an axisymmetric bridge tends to deform, limiting its stability and decreasing the maximum achievable value of R. Theoretical studies have discussed the stability and possible shapes of axisymmetric bridges. Experiments typically are performed in either a Plateau tank, in which the bridge is surrounded by a density-matched immiscible fluid, or in a space-borne microgravity environment. It has been shown, for example, that the stability limit R can be pushed beyond pi by using flow stabilization, by acoustic radiation pressure, or by forming columns in the presence of an axial electric field. In this work, magnetic levitation was used to simulate a low gravity environment and create quasi-cylindrical liquid columns in air. Use of a magnetic field permits us to continuously vary the Bond number B identically equal to (g)(rho)d(exp 2)/4(sigma), where g is the gravitational acceleration, rho is the density of the liquid, and sigma is the surface tension of the liquid in air. The dimensionless Bond number represents the relative importance of external forces acting on the liquid column to those due to surface tension. Our central result is that in a large magnetic field gradient we could create and stabilize columns of mixtures of water and paramagnetic manganese chloride tetrahydrate (MnCl2.4H2O), achieving a length to diameter ratio very close to pi.

Mahajan, Milind P.↗
NASA NTRS1996

Thermal Excitation System for Shearography (TESS)

One of the most convenient and effective methods of stressing a part or structure for shearographic evaluation is thermal excitation. This technique involves heating the part, often convectively with a heat gun, and then monitoring with a shearography device the deformation during cooling. For a composite specimen, unbonds, delaminations, inclusions, or matrix cracking will deform during cooling differently than other more structurally sound regions and thus will appear as anomalies in the deformation field. However, one of the difficulties that cause this inspection to be dependent on the operator experience is the conventional heating process. Fanning the part with a heat gun by hand introduces a wide range of variability from person to person and from one inspection to the next. The goal of this research effort was to conduct research in the methods of thermal excitation for shearography inspection. A computerized heating system was developed for inspection of 0.61 m (24 in.) square panels. The Thermal Excitation System for Shearography (TESS) provides radiant heating with continuous digital measurement of the surface temperature profile to ensure repeatability. The TESS device functions as an accessory to any electronic shearography device.

Lansing, Matthew D.↗
NASA NTRS1972

Stresses, strains, and surface pressures in the lung caused by its weight.

In an effort to understand how the lung is deformed by its own weight, we have analyzed the distribution of regional expansion, stresses, and surface pressures in a theoretical elastic lung-shaped model using the technique of finite elements. In the upright position, the parenchyma was most expanded at the apex and least at the base. Stresses in both the vertical and lateral directions were maximal at the apex. As the lung was inflated from very low volumes to total lung capacity, parenchymal expansion and stress at the apex first decreased, then increased. This behavior can be explained by the increasing rigidity of the expanded lung which enabled it to resist distortion by its own weight. At functional residual capacity, the stress at the apex was near its minimum. The differences in intrapleural pressure down the lung were volume dependent, increasing at very low volumes. In the inverted lung, the regional differences in stress, strain, and surface pressures were less marked because of the shape of the chest.

West, J. B.↗
NASA NTRS2018

Orion Optical Navigation Progress Toward Exploration: Mission 1

Optical navigation of human spacecraft was proposed on Gemini and implemented successfully on Apollo as a means of autonomously operating the vehicle in the event of lost communication with controllers on Earth. It shares a history with the "method of lunar distances" that was used in the 18th century and gained some notoriety after its use by Captain James Cook during his 1768 Pacific voyage of the HMS Endeavor. The Orion emergency return system utilizing optical navigation has matured in design over the last several years, and is currently undergoing the final implementation and test phase in preparation for Exploration Mission 1 (EM-1) in 2019. The software development is being worked as a Government Furnished Equipment (GFE) project delivered as an application within the Core Flight Software of the Orion camera controller module. The mathematical formulation behind the initial ellipse fit in the image processing is detailed in Christian. The non-linear least squares refinement then follows the technique of Mortari as an estimation process of the planetary limb using the sigmoid function. The Orion optical navigation system uses a body fixed camera, a decision that was driven by mass and mechanism constraints. The general concept of operations involves a 2-hour pass once every 24 hours, with passes specifically placed before all maneuvers to supply accurate navigation information to guidance and targeting. The pass lengths are limited by thermal constraints on the vehicle since the OpNav attitude generally deviates from the thermally stable tail-to-sun attitude maintained during the rest of the orbit coast phase. Calibration is scheduled prior to every pass due to the unknown nature of thermal effects on the lens distortion and the mounting platform deformations between the camera and star trackers. The calibration technique is described in detail by Christian, et al. and simultaneously estimates the Brown-Conrady coefficients and the Star Tracker/Camera interlock angles. Accurate attitude information is provided by the star trackers during each pass. Figure 1 shows the various phases of lunar return navigation when the vehicle is in autonomous operation with lost ground communication. The midcourse maneuvers are placed to control the entry interface conditions to the desired corridor for safe landing. The general form of optical navigation on Orion is where still images of the Moon or Earth are processed to find the apparent angular diameter and centroid in the camera focal plane. This raw data is transformed into range and bearing angle measurements using planetary data and precise star tracker inertial attitude. The measurements are then sent to the main flight computer's Kalman filter to update the onboard state vector. The images are, of course, collected over an arc to converge the state and estimate velocity. The same basic technique was used by Apollo to satisfy loss-of-comm, but Apollo used manual crew sightings with a vehicle-integral sextant instead of autonomously processing optical imagery. The software development is past its Critical Design Review, and is progressing through test and certification for human rating. In support of this, a hardware-in-the-loop test rig was developed in the Johnson Space Center Electro-Optics Lab to exercise the OpNav system prior to integrated testing on the Orion vehicle. Figure 2 shows the rig, which the test team has dubbed OCILOT (Orion Camera In the Loop Optical Testbed). Analysis performed to date shows a delivery that satisfies an allowable entry corridor as shown in Figure 3.

Holt, Greg N.↗
NASA NTRS1976

Oceanic lithosphere and asthenosphere: The thermal and mechanical structure

A coupled thermal and mechanical solid state model of the oceanic lithosphere and asthenosphere is presented. The model includes vertical conduction of heat with a temperature dependent thermal conductivity, horizontal and vertical advection of heat, viscous dissipation or shear heating, and linear or nonlinear deformation mechanisms with temperature and pressure dependent constitutive relations between shear stress and strain rate. A constant horizontal velocity u sub 0 and temperature t sub 0 at the surface and zero horizontal velocity and constant temperature t sub infinity at great depth are required. In addition to numerical values of the thermal and mechanical properties of the medium, only the values of u sub 0, t sub 0 and t sub infinity are specified. The model determines the depth and age dependent temperature horizontal and vertical velocity, and viscosity structures of the lithosphere and asthenosphere. In particular, ocean floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of the age of the ocean floor.

Schubert, G.↗
NASA NTRS1983

Application of VLBI and satellite laser ranging to geodynamics

The NASA Crustal Dynamics Project has developed very-long baseline interferometer (VLBI) systems and satellite laser ranging (SLR) systems for geodynamics measurements. In VLBI, a radio noise signal from a distant quasar is received by two or more radio antennas and coherently recorded. These recordings are cross-correlated to determine the relative signal delays between stations which are used to derive the vector baselines between the stations. The SLR systems accurately determine the range to a retroreflector satellite as a function of time with short laser pulses. These range measurements from several stations to the same satellite are used in orbit analysis programs to determine the position of the stations and the vector baselines between the stations. Measurements with these systems have achieved precisions of a few centimeters in length for distances of several thousand km. These systems are now operating in a global network for measuring the relative motion of the N. American, Pacific, S. American, Nazca, Eurasian and Australian tectonic plates. Highly mobile VLBI and SLR systems are being operated at many sites in the active earthquake areas in western N. America in order to determine the crustal deformation and strain accumulation.

Coates, R. J.↗
NASA NTRS1992

Determining stress states using dike swarms: The Lauma Dorsa example

Initial examination of the Magellan coverage of Venus has revealed between 150 and 300 large, radially lineated landforms distributed across the planet's surface. Where the lineaments have been examined in detail, the majority fail to exhibit signatures indicative of relief at or above the resolution of the radar; however, when the sense of topographic relief may be ascertained, the lineaments commonly appear as fissures or flat-floored trenches interpreted as graben. Individual lineaments can display graben, fissure, and zero relief behavior along their length, suggesting either that these differences are a function of the resolution of the radar, or that the morphological distinctions are real but somehow genetically linked. In many instances, radial lineaments exhibiting these characteristics are directly associated with surface volcanism, including flanking and terminal flows, superimposed shield domes and pit chains, and central, calderalike topographic lows. These observable characteristics, as well as theoretical studies and comparison with similar terrestrial features, have led to the working hypothesis that many of the radial fracture systems on Venus are the surface manifestation of subsurface dikes propagating laterally from a central magma source. If this interpretation is correct, studies of terrestrial dikes suggest that the lineament directions, with localized exceptions and barring subsequent deformation, should be perpendicular to the orientation of the least compressive stress at the time of their formation. To test this hypothesis, we briefly examine a radial fracture system (63.7 degrees N, 195 degrees E) located between two deformation belts in Vinmara Planitia, and verify that the lineaments to the east behave in the expected manner. We have also chosen this feature, however, because of its proximity to Lauma Dorsa to the west. On the basis of Venera 15/16 data, both compressional and extensional origins for this deformation belt have been proposed. By examining the stratigraphy and applying our interpretation that the fracture system is linked to the presence of subsurface dikes, we present an independent evaluation of the stress state associated with Lauma Dorsa, and thus contribute to the assessment of its origin.

Grosfils, Eric B.↗
NASA NTRS1983

The dynamic response of the high-latitude thermosphere and geostrophic adjustment

The theory of geostrophic adjustment in which the adjustment process is critically dependent on the ratio of the scale size of the forcing to the Rossby radius of deformation, is applied to the qualitative explanation of certain features of the response of the high latitude thermosphere to inputs of energy and momentum. When the ratio is small, large wind speeds are efficiently generated by a momentum source, and when small, changes in the mass field associated with a heat source are more effectively produced. The Rossby radius variation with height explains the qualitatively different responses of the high latitude E and F regions during disturbed geomagnetic conditions as shown by both observations and numerical models. The present analytic solutions are in good agreement with model results concerning energy apportionment between divergent and rotational components of the wind field as a function of wavenumber.

Larsen, M. F.↗
NASA NTRS1984

Water-vapor effects on friction of magnetic tape in contact with nickel-zinc ferrite

The effects of humidity of moist nitrogen on the friction and deformation behavior of magnetic tape in contact with a nickel-zinc ferrite spherical pin were studied. The results indicate that the coefficient of friction is markedly dependent on the ambient relative humidity. Although the coefficient of friction remains low below 40-percent relative humidity, it increases rapidly with increasing relative humidity above 40 percent. The general ambient environment of the tape does not have any effect on the friction behavior if the area where the tape is in sliding contact with the ferrite pin is flooded with controlled nitrogen. The response time for the friction of the tape to humidity changes is about 10 sec. The effect of friction as a function of relative humidity on dehumidifying is very similar to that on humidifying. A surface softening of the tape due to water vapor increases the friction of the tape.

Miyoshi, K.↗
NASA NTRS1990

Elevated temperature crack growth in aluminum alloys: Tensile deformation of 2618 and FVS0812 aluminum alloys

Understanding the damage tolerance of aluminum alloys at elevated temperatures is essential for safe applications of advanced materials. The objective of this project is to investigate the time dependent subcritical cracking behavior of powder metallurgy FVS0812 and ingot metallurgy 2618 aluminum alloys at elevated temperatures. The fracture mechanics approach was applied. Sidegrooved compact tension specimens were tested at 175, 250, and 316 C under constant load. Subcritical crack growth occurred in each alloy at applied stress intensity levels (K) of between about 14 and 25 MPa/m, well below K (sub IC). Measured load, crack opening displacement and displacement rate, and crack length and growth rate (da/dt) were analyzed with several continuum fracture parameters including, the C-integral, C (sub t), and K. Elevated temperature growth rate data suggest that K is a controlling parameter during time dependent cracking. For FVS0812, da/dt is highest at 175 C when rates are expressed as a function of K. While crack growth rate is not controlled by C (sub t) at 175 C, da/dt appears to better correlate with C (sub t) at higher temperatures. Creep brittle cracking at intermediate temperatures, and perhaps related to strain aging, is augmented by time dependent transient creep plasticity at higher temperatures. The C (sub t) analysis is, however, complicated by the necessity to measure small differences in the elastic crack growth and creep contributions to the crack opening displacement rate. A microstructural study indicates that 2618 and FVS0812 are likely to be creep brittle materials, consistent with the results obtained from the fracture mechanics study. Time dependent crack growth of 2618 at 175 C is characterized by mixed transgranular and intergranular fracture. Delamination along the ribbon powder particle boundaries occurs in FVS0812 at all temperatures. The fracture mode of FVS0812 changes with temperature. At 175 C, it is characterized as dimpled rupture, and at 316 C as mixed matrix superplastic rupture and matrix-dipersoid debonding.

Leng, Yang↗
NASA NTRS1990

The effect of hydrogen and microstructure on the deformation and fracture behavior of a single crystal nickel-base superalloy

A study was conducted on the effects of internal hydrogen and microstructure on the deformation and fracture of a single crystal nickel-base superalloy. In particular, room temperature plane strain fracture toughness and tensile tests were performed on hydrogen-free and hydrogen charged samples of PWA 1480. The role of microstructure was incorporated by varying the levels of porosity and eutectic gamma/gamma prime through hot isostatic pressing and heat treatment. The room temperature behavior of PWA 1480 was unusual because precipitate shearing was not the primary deformation mechanism at all strains. At strains over 1 percent, dislocations were trapped in the gamma matrix and an attempt was made to relate this behavior to compositional differences between PWA 1480 and other superalloys. Another unique feature of the tensile behavior was cleavage of the eutectic gamma/gamma prime, which is believed to initiate the failure process. Fracture occurred on (111) planes and is likely a result of shear localization along these planes. Elimination of the eutectic gamma/gamma prime greatly improved the tensile ductility, but pososity had no effect on tensile properties. Large quantities of hydrogen (1.74 at. percent) were gas-phase charged into the material, but surprisingly this was not a function of the amount of porosity or eutectic gamma/gamma prime present. Desorption experiments suggest that the vast majority of hydrogen is at reversible lattice trapping sites. This large, uniform concentration of hydrogen dramatically reduced the tensile strain to failure, but only slightly affected the reduction in area. Available hydrogen embrittlement models were examined in light of these results and it was found that the hydrogen enhanced localized plasticity model can explain much of the tensile behavior. K(IC) fracture toughness tests were conducted, but it was necessary to also perform J(IC) tests to provide valid data.

Walston, William S.↗
NASA NTRS1995

Advanced Modeling of Micromirror Devices

The flexure-beam micromirror device (FBMD) is a phase only piston style spatial light modulator demonstrating properties which can be used for phase adaptive corrective optics. This paper presents a complete study of a square FBMD, from advanced model development through final device testing and model verification. The model relates the electrical and mechanical properties of the device by equating the electrostatic force of a parallel-plate capacitor with the counter-acting spring force of the device's support flexures. The capacitor solution is derived via the Schwartz-Christoffel transformation such that the final solution accounts for non-ideal electric fields. The complete model describes the behavior of any piston-style device, given its design geometry and material properties. It includes operational parameters such as drive frequency and temperature, as well as fringing effects, mirror surface deformations, and cross-talk from neighboring devices. The steps taken to develop this model can be applied to other micromirrors, such as the cantilever and torsion-beam designs, to produce an advanced model for any given device. The micromirror devices studied in this paper were commercially fabricated in a surface micromachining process. A microscope-based laser interferometer is used to test the device in which a beam reflected from the device modulates a fixed reference beam. The mirror displacement is determined from the relative phase which generates a continuous set of data for each selected position on the mirror surface. Plots of this data describe the localized deflection as a function of drive voltage.

Michalicek, M. Adrian↗
NASA NTRS2002

Molten-Metal Droplet Deposition on a Moving Substrate in Microgravity: Aiding the Development of Novel Technologies for Microelectronic Assembly

Driven by advancements in microelectronics manufacturing, this research investigates the oblique (non-axisymmetric) impact of liquid-metal droplets on flat substrates. The problem of interest is relevant to the development of the novel technology of on-demand dispension (printing) of microscopic solder deposits for the surface mounting of microelectronic devices. The technology, known as solder jetting, features on-demand deposition of miniature solder droplets (30 to 120 microns in diameter) in very fine, very accurate patterns using techniques analogous to those developed for the ink-jet printing industry. Despite its promise, severe limitations exist currently with regards to the throughput rates of the technology; some of these limitations are largely due to the lack of the capability for reliable prediction of solder bump positioning and shapes, especially under ballistic deposition conditions where the droplet impact phenomena are inherently three-dimensional. The study consists of a theoretical and an experimental component. The theoretical work uses a finite element formulation to simulate numerically the non-axisymmetric (3-D) fluid mechanics and heat transfer phenomena of a liquid solder droplet impacting at an angle alpha on a flat substrate. The work focuses on the pre-solidification regime. The modeling of the most challenging fluid mechanics part of the process has been completed successfully. It is based upon the full laminar Navier-Stokes equations employing a Lagrangian frame of reference. Due to the large droplet deformation, the surface (skin) as well as the volumetric mesh have to be regenerated during the calculations in order to maintain the high accuracy of the numerical scheme. The pressure and velocity fields are then interpolated on the newly created mesh. The numerical predictions are being tested against experiments, for cases where wetting phenomena are not important. For the impact parameters used in the example shown (We = 2.38, Fr = 16300, Re = 157), the droplet rolls along the substrate, but its shape remains practically axisymmetric for all impact angles within the range from 0 to 60 deg. Interestingly, the substrate/droplet contact area during the recoiling phase of the impact is not a monotonically decreasing function of time. The experimental component of the research tests the numerical predictions and provides necessary input data (contact angles) for the theoretical model. The experiments are performed in microgravity (2.2s drop tower of the NASA GRC) in order to allow for the use of mm-size solder droplets, which make feasible the performance of accurate measurements, while maintaining similitude of the relevant fluid dynamic groups (Re, Fr, We, Ste). Preliminary oblique impact experiments have been performed using water droplets in normal gravity.

Megaridis, C. M.↗