Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tension tracking”

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 37 records · Page 2

Compact Telescoping Surface Array for Mars Solar Power

This paper describes a lightweight, large-area solar array concept for Mars surface power called the Compact Telescoping Surface Array (CTSA). The design is derived from the Compact Telescoping Array (CTA) proposed in 2015 for high-power spacecraft. The CTSA deploys horizontally from Mars landers, provides 1000 sq.m of solar cell area, and generates about 50-80 kW daytime and 10 kW nighttime power (from energy storage) near the equator with clear skies. The total mass is about 1500 kg, and the stowed volume is about 10 cu.m, equivalent to 200 W/kg and 30 kW/cu.m at 1 astronomical unit (AU) from the sun. These challenging mass and volume goals require compact, lightweight structures to deploy and tension thin solar cell blankets in the Mars 0.38 g gravity field. Other features include single-axis rotation for daily or seasonal sun tracking, dust shedding, and feathering in high winds and sufficient strength for sustained 40 m/s winds in dust storms and brief 100 m/s winds in dust devils (dynamic pressure of 120 Pa). Based on preliminary design and analysis reported in the paper, CTSA appears to be a feasible approach for 1000 sq.m-class Mars solar arrays that autonomously deploy from landers prior to human arrival.

Pappa, Richard↗

Development of a tether deployment monitoring system

The feasibility of monitoring the position of a tether during orbital deployment is discussed. Attention is given to the techniques for attaching reflectors to the tether at periodic intervals during the deployment process, which could then be tracked visually or with a radar unit that is pointed towards the targets by the control system. Consideration is also given to a method for tether attitude and tension detection, and a launching device to attach small instrumentation platforms onto the tether during deployment. The configuration of a tether deployment monitoring system (TEDEMS) and the equipment developed to perform the above tasks are described, along with an analysis of requirements for the construction of a family of space-borne TEDEMS suitable for use on the Space Shuttle, Expendable Launch Vehicles, or the Space Station. The potential uses of the TEDEMS for launching small tether-borne experimental packages are discussed. Design and block diagrams are included.

Ibanez, P.↗

Tension and compression fatigue response of unnotched 3D braided composites

The unnotched compression and tension fatigue response of a 3-D braided composite was measured. Both gross compressive stress and tensile stress were plotted against cycles to failure to evaluate the fatigue life of these materials. Damage initiation and growth was monitored visually and by tracking compliance change during cycle loading. The intent was to establish by what means the strength of a 3-D architecture will start to degrade, at what point will it degrade beyond an acceptable level, and how this material will typically fail.

Portanova, M. A.↗

Remobilizing the Interface of Thermocapillary Driven Bubbles Retarded By the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow that propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillary to direct the movement of bubbles in space is the fact that surfactant impurities, which are unavoidably present in the continuous phase, can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. (This can also be the case if kinetic desorption of surfactant at the back end of the bubble is much slower than convection.) For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like and bubbles translate as solid particles). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature-induced tension gradient, and can decrease to near zero the bubble's thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have be retarded by the adsorption onto the bubble surface of a surfactant impurity. Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble surface much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition, the rapid kinetic exchange and high bulk concentration maintain a saturated surface with uniform surface concentrations. This prevents retarding surface tension gradients and keeps the thermocapillary velocity high. In our reports over the first 2 years, we presented numerical simulations of the bubble motion and surfactant transport which verified theoretically the concept of remobilization, and the development of an apparatus to track and measure the velocity of rising bubbles in a glycerol/water surfactant solution. This year, we detail experimental observations of remobilization. Two polyethylene oxide surfactants were studied, C12E6 (CH3(CH2)11(OCH2)6OH) and C10E8 (CH3(CH2)4(OCH2CH2)8OH). Measurements of the kinetic exchange for these surfactants show that the one with the longer hydrophobe chain C12E6 has a lower rate of kinetic exchange. In addition, this surfactant is much less soluble in the glycerol/water mixture because of the shorter ethoxylate chain. As a result, we found that C12E6 had only a very limited ability to remobilize rising bubbles because of the limited kinetic exchange and reduced solubility. However, C10E8, with its higher solubility and more rapid exchange was found to dramatically remobilize rising bubbles. We also compared our theoretical calculations to the experimental measurements of velocity for both the non-remobilizing and remobilizing surfactants and found excellent agreement. We further observed that for C10E8 at high concentrations, which exceeded the critical micelle concentrations, additional remobilization was measured. In this case the rapid exchange of monomer between micelle and surfactant provides an additional mechanism for maintaining a uniform surface concentrations.

Palaparthi, Ravi↗

Several Developments in Space Tethers

Five reports address different aspects of development of tethers to be deployed from spacecraft in orbit around the Earth. The first report discusses proposed optoelectronic tracking of retroreflective objects located at intervals or of retroreflective coats along the entire length of a tether to measure lateral motions. The second report describes digitally controlled spooling machinery that retracts or extends a tape tether at controlled speed and tension in the spool isolated from uncontrolled tension on the outside. The third report discusses part of this machinery that pivots to accommodate misalignments between the deployed and spooled portions of the tether and contains rollers used to exert tension and speed control. The fourth report discusses aspects of designs of proposed electrodynamic tethers, which would be electrically conductive and would interact with the magnetic field of the Earth to exert forces to modify orbits of deploying spacecraft. The fifth report discusses electrical aspects of designs of electrodynamic tape tethers, including the use of solar cells or motional electromagnetic force to generate currents in tethers and the use of electron emitters and electron and ion collectors at opposite ends of tethers to make electrical contact with the thin plasma in surrounding space.

Santangelo, Andrew↗

Design, Development and Test Challenges: Separation Mechanisms for the Orion Pad Abort-1 Flight Test

On May 6, 2010, NASA launched the first successful integrated flight test, Pad Abort-1, of the Orion Project from the White Sands Missile Range in Las Cruces, New Mexico. This test demonstrated the ability to perform an emergency pad abort of a full-scale 4.8 m diameter, 8200 kg crew capsule. During development of the critical separation mechanisms for this flight test, various challenges were overcome related to environments definition, installation complications, separation joint retraction speed, thruster ordnance development issues, load path validation and significant design loads increases. The Launch Abort System retention and release (LAS R&R) mechanism consisted of 6 discrete structural connections between the LAS and the crew module (CM) simulator, each of which had a preloaded tension tie, Superbolt torque-nut and frangible nut. During the flight test, the frangible nuts were pyrotechnically split, permitting the CM to separate from the LAS. The LAS separation event was the driving case in the shock environment for many co-located hardware items. During development testing, it was necessary to measure the source shock during the separation event so the predicted shock environment could be validated and used for certification testing of multiple hardware items. The Lockheed Martin test team measured the source separation shock due to the LAS R&R function, which dramatically decreased the predicted environment by 90% at 100 Hz. During development testing a hydraulic tensioner was used to preload the joint; however, the joint relaxation with the tensioner proved unsatisfactory so the design was modified to include a Superbolt torque-nut. The observed preload creep during lab testing was 4% after 30 days, with 2.5% occurring in the first 24 hours. The conversion of strain energy (preload) to kinetic energy (retraction) was measured to be 50-75%. Design features and careful monitoring of multiple strain gauges on each tension tie allowed a pure tensile load to be applied after stacking at the launch pad. Following installation, preload in each joint was monitored for 24 hours. Due to unforeseen complications and the influence of temperature on the portable data acquisition system, the team encountered difficulty in tracking the joint relaxation. In some cases, bond-line failure of the strain gauges occurred.

Dinsel, Alison↗

Calibration of the 11x11-Foot Transonic Wind Tunnel at NASA Ames Research Center

A static pipe calibration characterizing the streamwise static pressure distribution was conducted at the 11-By 11-Foot Transonic Wind Tunnel at NASA Ames Research Center. This data is used to determine the local Mach number within the test section and evaluate buoyancy corrections to axial force measurements. The 60 foot long, 6 inch diameter pipe contained 444 static pressure taps spanning the test section and nozzle regions of the tunnel. The forward end of the pipe extends into the settling chamber and is held by four cables mounting to the tunnel shell, and the aft end is fixed on the institutional model support system. A hydraulic cylinder at the aft end of the pipe provides tension on the system to reduce vibration and to keep the pipe as level as possible throughout the test section. The previous calibration was improved upon by using pressure scanners with greater accuracy, ensuring a uniform pressure tube length for each tap to control pneumatic lag, optically tracking any streamwise movement of the pipe, and more tightly controlling the tunnel condition set points. Typically this calibration is conducted with the pipe on tunnel centerline and 33 inches below centerline for sting-mounted models and semi-span (i.e. floor-mounted) models respectively, however schedule demands permitted only the centerline calibration. The semi-span calibration is planned to be completed in the summer of 2020. Immediately following the static pipe calibration, a shorter, 9 foot static pipe used as the calibration check standard was installed to obtain its first post-calibration pressure dataset. This short static pipe consists of 148 static pressure taps distributed along the pipe section and one total pressure tap at the end of an ogive nose. Performing the calibration test and its check standard back-to-back allows this dataset to establish a reliable baseline for future calibration check standard testing. Over time the use of a calibration check standard offers the ability to assess the stability of the calibration through statistical process control in an efficient and costeffective manner, thereby potentially increasing the time required between full tunnel calibrations.

White, Bethany↗

Modeling and Simulation of the Angel Upper Limb Offload Device: Branching Into New Methods

BACKGROUND: The Active Response Gravity Offload System (ARGOS) provides an analog environment for extravehicular activity (EVA) testing and training. Discomfort has been observed during longer suited test sessions. While the subject’s core is offloaded during surface EVA evaluations, his/her arms experience full Earth gravity and can become overly fatigued, especially during suited tests which involve reaching and prolonged arm extensions. A device (ARGOS Negation of Gravitational Effects on the Limbs: ANGEL) to offload the weight of the arms and suit sleeves is being developed by JSC’s Flight Systems Branch of the Software, Robotics, and Simulation Division. Previously we have shared preliminary modeling of that device and kinematics based on motion capture data. Here we present an alternative approach to determine device kinematics by calculating ANGEL component angles with an OpenSim plugin. We compare calculated angles to inverse kinematics (IK) derived ones with the goal of validating the model. This new method can be further informative for device design and analytically testing different configurations to achieve desired reduced gravity conditions (e.g., lunar gravity (Lg) or Martian gravity (Mg)). We have compared calculated angles with IK-derived angles in tests with a shirt-sleeve subject positioned in a test stand with a Mark-III Hard Upper Torso (HUT) and Portable Life Support System (PLSS) mockup and arm weights to emulate the weight of the suit sleeve as well as a suited subject in ARGOS with a Mark-III suit. A variety of upper body tasks were completed in the former and full-body tasks in the latter. METHODS AND RESULTS: To model the offload device, we augment the OpenSim human model topology with the offload mechanism components and joints, using CAD models to represent the mechanism graphically. The joint angles of the device are calculated in the OpenSim plugin by modeling how the components configure themselves under the offloading spring tension given a particular IK-derived arm position. There are four ANGEL components with a total of 5 degrees of freedom (DOFs), each component has a single DOF except for the cuff which is modeled as 2 DOFs. The sickle/yaw bracket and cuff rotation angles are determined statically based on the assumptions that the sickle will track the attachment point of the cuff and that the cuff will rotate such that the attachment point is at its highest point. The cuff tilt, linker and V-bracket angles are then determined by optimizing their positions to approach a mechanical equilibrium. The calculated linker angle is compared to three different methods of determining the linker line-of-force kinematically (from V-bracket to center-cuff, cuff highest point or marker-derived position). Given the joint angles of the device, the spring force and resulting force on the arm is computed by the plugin and applied as an external load in inverse dynamics (ID) to enable study of overall shoulder joint torques as well as offload achieved. We verify the calculated joint angles by using the inverse kinematic data. The average difference in angles is the smallest for the V-bracket and linker, around 1 to 5 degrees for most trials. The resulting offload and shoulder torque are comparable between calculated and IK-derived angles. In summary, we have developed a method to calculate the joint angles of an exoskeleton-like upper limb offloading device currently in development. We have also developed a custom plugin which will be a valuable tool to optimize device configurations for a desired gravitational environment, probe the offload achieved for motions recorded outside of our test suite, and inform future design improvements.

L B Nilsson↗

Modeling and Simulation of The Angel Upper Limb Offload Device: Branching into New Methods

BACKGROUND: The Active Response Gravity Offload System (ARGOS) provides an analog environment for extravehicular activity (EVA) testing and training. Discomfort has been observed during longer suited test sessions. While the subject’s core is offloaded during surface EVA evaluations, his/her arms experience full Earth gravity and can become overly fatigued, especially during suited tests which involve reaching and prolonged arm extensions. A device (ARGOS Negation of Gravitational Effects on the Limbs: ANGEL) to offload the weight of the arms and suit sleeves is being developed by JSC’s Flight Systems Branch of the Software, Robotics, and Simulation Division. Previously we have shared preliminary modeling of that device and kinematics based on motion capture data. Here we present an alternative approach to determine device kinematics by calculating ANGEL component angles with an OpenSim plugin. We compare calculated angles to inverse kinematics (IK) derived ones with the goal of validating the model. This new method can be further informative for device design and analytically testing different configurations to achieve desired reduced gravity conditions (e.g., lunar gravity (Lg) or Martian gravity (Mg)). We have compared calculated angles with IK-derived angles in tests with a shirt-sleeve subject positioned in a test stand with a Mark-III Hard Upper Torso (HUT) and Portable Life Support System (PLSS) mockup and arm weights to emulate the weight of the suit sleeve as well as a suited subject in ARGOS with a Mark-III suit. A variety of upper body tasks were completed in the former and full-body tasks in the latter. METHODS AND RESULTS: To model the offload device, we augment the OpenSim human model topology with the offload mechanism components and joints, using CAD models to represent the mechanism graphically. The joint angles of the device are calculated in the OpenSim plugin by modeling how the components configure themselves under the offloading spring tension given a particular IK-derived arm position. There are four ANGEL components with a total of 5 degrees of freedom (DOFs), each component has a single DOF except for the cuff which is modeled as 2 DOFs. The sickle/yaw bracket and cuff rotation angles are determined statically based on the assumptions that the sickle will track the attachment point of the cuff and that the cuff will rotate such that the attachment point is at its highest point. The cuff tilt, linker and V-bracket angles are then determined by optimizing their positions to approach a mechanical equilibrium. The calculated linker angle is compared to three different methods of determining the linker line-of-force kinematically (from V-bracket to center-cuff, cuff highest point or marker-derived position). Given the joint angles of the device, the spring force and resulting force on the arm is computed by the plugin and applied as an external load in inverse dynamics (ID) to enable study of overall shoulder joint torques as well as offload achieved. We verify the calculated joint angles by using the inverse kinematic data. The average difference in angles is the smallest for the V-bracket and linker, around 1 to 5 degrees for most trials. The resulting offload and shoulder torque are comparable between calculated and IK-derived angles. In summary, we have developed a method to calculate the joint angles of an exoskeleton-like upper limb offloading device currently in development. We have also developed a custom plugin which will be a valuable tool to optimize device configurations for a desired gravitational environment, probe the offload achieved for motions recorded outside of our test suite, and inform future design improvements.

L B Nilsson↗

SEDS1 mission software verification using a signal simulator

The first flight of the Small Expendable Deployer System (SEDS1) is schedule to fly as the secondary payload of a Delta 2 in March, 1993. The objective of the SEDS1 mission is to collect data to validate the concept of tethered satellite systems and to verify computer simulations used to predict their behavior. SEDS1 will deploy a 50 lb. instrumented satellite as an end mass using a 20 km tether. Langley Research Center is providing the end mass instrumentation, while the Marshall Space Flight Center is designing and building the deployer. The objective of the experiment is to test the SEDS design concept by demonstrating that the system will satisfactorily deploy the full 20 km tether without stopping prematurely, come to a smooth stop on the application of a brake, and cut the tether at the proper time after it swings to the local vertical. Also, SEDS1 will collect data which will be used to test the accuracy of tether dynamics models used to stimulate this type of deployment. The experiment will last about 1.5 hours and complete approximately 1.5 orbits. Radar tracking of the Delta II and end mass is planned. In addition, the SEDS1 on-board computer will continuously record, store, and transmit mission data over the Delta II S-band telemetry system. The Data System will count tether windings as the tether unwinds, log the times of each turn and other mission events, monitor tether tension, and record the temperature of system components. A summary of the measurements taken during the SEDS1 are shown. The Data System will also control the tether brake and cutter mechanisms. Preliminary versions of two major sections of the flight software, the data telemetry modules and the data collection modules, were developed and tested under the 1990 NASA/ASEE Summer Faculty Fellowship Program. To facilitate the debugging of these software modules, a prototype SEDS Data System was programmed to simulate turn count signals. During the 1991 summer program, the concept of simulating signals produced by the SEDS electronics systems and circuits was expanded and more precisely defined. During the 1992 summer program, the SEDS signal simulator was programmed to test the requirements of the SEDS Mission software, and this simulator will be used in the formal verification of the SEDS Mission Software. The formal test procedures specification was written which incorporates the use of the signal simulator to test the SEDS Mission Software and which incorporates procedures for testing the other major component of the SEDS software, the Monitor Software.

Pierson, William E.↗

Dendritic growth in the presence of convection

The motion of the freezing front between a dendritic crystal and a supercooled liquid is studied using an interface evolution equation derived from a boundary integral transformation of the transient convective-diffusion equation. A new steady-state theory is introduced that incorporates the effects of convection in dendritic growth. It is shown that in the absence of capillary effects the shape of the crystal-melt interface is a paraboloid of revolution, similar to that found in situations where diffusion is the sole heat transfer mechanism. A relation between the supercooling, the product of the tip velocity and tip radius, and the strength of the flow is derived which reduces to the well-known Ivantsov theory in the absence of convection. A non-linear interface-tracking algorithm is developed and used to study the temporal and spatial evolution of the dendritic interface. The important role of capillarity and convection on the interface dynamics is established and the response of the interface to finite amplitude disturbances is examined for the first time. Tip splitting is identified as the dominant destabilization mechanism in the limit of zero surface tension. Finite surface tension leads to interface stabilization, irrespective of the magnitude and structure of the external perturbations. Finally, convection significantly decreases the magnitude of the freezing velocity.

Beaghton, Pantelis John↗

Tethered Vehicle Control and Tracking System

A kite system includes a kite and a ground station. The ground station includes a sensor that can be utilized to determine an angular position and velocity of the kite relative to the ground station. A controller utilizes a fuzzy logic control system to autonomously fly the kite. The system may include a ground station having powered winding units that generate power as the lines to the kite are unreeled. The control system may be configured to fly the kite in a crosswind trajectory to increase line tension for power generation. The sensors for determining the position of the kite are preferably ground-based.

North, David D.↗

Tethered Vehicle Control and Tracking System

A kite system includes a kite and a ground station. The ground station includes a sensor that can be utilized to determine an angular position and velocity of the kite relative to the ground station. A controller utilizes a fuzzy logic control system to autonomously fly the kite. The system may include a ground station having powered winding units that generate power as the lines to the kite are unreeled. The control system may be configured to fly the kite in a crosswind trajectory to increase line tension for power generation. The sensors for determining the position of the kite are preferably ground-based.

North, David D.↗

Stereo Imaging Velocimetry

Stereo imaging velocimetry (SIV) will permit the collection of quantitative, three-dimensional flow data from any optically transparent fluid that can be seeded with tracer particles. This includes such diverse experiments as the study of multiphase flow, bubble nucleation and migration, pool combustion, and crystal growth. This technique will be useful to the microgravity science community as our investigations of fluid behavior in reduced-gravity environments enhance our knowledge of heat transfer, surface tension, concentration-gradient-driven anomalies, and residual effects from g-jitter. In its proposed configuration, the NASA Lewis Research Center's Stereo Imaging Velocimeter will consist of at least two charged coupled device (CCD) cameras, oriented at some relative angle with respect to each other. The cameras will observe a fluid experiment that has been seeded with tracer particles that are neutrally buoyant to permit accurate flow tracking. Except for the tracer particles, this measurement technique will be nonintrusive. Velocity accuracies will be on the order of 1 to 5 percent of full field. Each camera will make a two-dimensional record of the motion of the seed particles in the observation volume. Three-dimensional data will be obtained by computationally combining the two-dimensional information.

Source record↗

Ground-based PIV and numerical flow visualization results from the surface tension driven convection experiment

The Surface Tension Driven Convection Experiment (STDCE) is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the United States Microgravity Laboratory-1 (USML-1) Spacelab mission planned for June, 1992. One of the components of data collected during the experiment is a video record of the flow field. This qualitative data is then quantified using an all electric, two dimensional Particle Image Velocimetry (PIV) technique called Particle Displacement Tracking (PDT), which uses a simple space domain particle tracking algorithm. Results using the ground based STDCE hardware, with a radiant flux heating mode, and the PDT system are compared to numerical solutions obtained by solving the axisymmetric Navier Stokes equations with a deformable free surface. The PDT technique is successful in producing a velocity vector field and corresponding stream function from the raw video data which satisfactorily represents the physical flow. A numerical program is used to compute the velocity field and corresponding stream function under identical conditions. Both the PDT system and numerical results were compared to a streak photograph, used as a benchmark, with good correlation.

Pline, Alexander D.↗

Ground-based PIV and numerical flow visualization results from the Surface Tension Driven Convection Experiment

The Surface Tension Driven Convection Experiment (STDCE) is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the United States Microgravity Laboratory-1 (USML-1) Spacelab mission planned for June, 1992. One of the components of data collected during the experiment is a video record of the flow field. This qualitative data is then quantified using an all electric, two dimensional Particle Image Velocimetry (PIV) technique called Particle Displacement Tracking (PDT), which uses a simple space domain particle tracking algorithm. Results using the ground based STDCE hardware, with a radiant flux heating mode, and the PDT system are compared to numerical solutions obtained by solving the axisymmetric Navier Stokes equations with a deformable free surface. The PDT technique is successful in producing a velocity vector field and corresponding stream function from the raw video data which satisfactorily represents the physical flow. A numerical program is used to compute the velocity field and corresponding stream function under identical conditions. Both the PDT system and numerical results were compared to a streak photograph, used as a benchmark, with good correlation.

Pline, Alexander D.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

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

Chen, Yen-Sen↗

Representing Matrix Cracks Through Decomposition of the Deformation Gradient Tensor in Continuum Damage Mechanics Methods

A method is presented to represent the large-deformation kinematics of intraply matrix cracks and delaminations in continuum damage mechanics (CDM) constitutive material models. The method involves the additive decomposition of the deformation gradient tensor into 'crack' and 'bulk material' components. The response of the intact bulk material is represented by a reduced deformation gradient tensor, and the opening of an embedded cohesive interface is represented by a normalized cohesive displacement-jump vector. The rotation of the embedded interface is tracked as the material deforms and as the crack opens. The distribution of the total local deformation between the bulk material and the cohesive interface components is determined by minimizing the difference between the cohesive stress and the bulk material stress projected onto the cohesive interface. The improvements to the accuracy of CDM models that incorporate the presented method over existing approaches are demonstrated for a single element subjected to simple shear deformation and for a finite element model of a unidirectional open-hole tension specimen. The material model is implemented as a VUMAT user subroutine for the Abaqus/Explicit finite element software. The presented deformation gradient decomposition method reduces the artificial load transfer across matrix cracks subjected to large shearing deformations, and avoids the spurious secondary failure modes that often occur in analyses based on conventional progressive damage models.

Leone, Frank A., Jr.↗