Search NASA⌕ Search

SEARCH · Search NASA

Results for “code patterns”

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 505 records · Page 28

Concerning the Perturbations Generated by Flux Transfer Events

The interaction of the solar wind with the Earth's magnetosphere is often highly unsteady. Bursts of magnetic reconnection at multiple locations on the dayside equatorial magnetopause generate flux transfer events or FTEs: twisted ropes of interconnected magnetosheath and magnetospheric magnetic field lines. Once formed, the events move antisunward, displacing and perturbing the ambient media. This talk explores the perturbations predicted by both global numerical simulations and analytical models. Results from global hybrid code simulation confirm the predictions of analytical models indicating that the events generate standing forward slow mode waves as their speeds relative to the magnetosheath flow approach the Alfven velocity. Geometric considerations lead to the conclusion that events generated by component reconnection on the dayside magnetopause move poleward and exhibit strong signatures during intervals of southward IMF orientation, but move towards the flanks and exhibit weak signatures during intervals of northward IMF orientation. Changing event orientations and magnetosheath/magnetospheric magnetic field orientation can enhance the amplitudes of events reaching the flanks. Although the orientations of events on the flanks inferred from multispacecraft timing techniques are consistent with the predictions of the component reconnection model, occurrence patterns versus latitude and IMF orientation require an explanation in terms of both the component and antiparallel reconnection models.

Sibeck, David G.↗

Modification of a Macromechanical Finite-Element Based Model for Impact Analysis of Triaxially-Braided Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. For the current analytical approach, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The constitutive model requires stiffness and strength properties of an equivalent unidirectional composite. Simplified micromechanics methods are used to determine the equivalent stiffness properties, and results from coupon level tests on the braided composite are utilized to back out the required strength properties. Simulations of quasi-static coupon tests of several representative braided composites are conducted to demonstrate the correlation of the model. Impact simulations of a represented braided composites are conducted to demonstrate the capability of the model to predict the penetration velocity and damage patterns obtained experimentally.

Goldberg, Robert K.↗

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗

A new Mars radiation environment model with visualization

A new model for the radiation environment to be found on the planet Mars due to Galactic Cosmic Rays (OCR) has been developed at the NASA Langley Research Center. Solar modulated primary particles rescaled for Mars conditions are transported through the Martian atmosphere, with temporal properties modeled with variable timescales, down to the surface, with altitude and backscattering patterns taken into account. The Martian atmosphere has been modeled by using the Mars Global Reference Atmospheric Model--version 2001 (Mars-GRAM 2001). The altitude to compute the atmospheric thickness profile has been determined by using a model for the topography based on the data provided by the Mars Orbiter Laser Altimeter (MOLA) instrument on board the Mars Global Surveyor (MGS) spacecraft. The Mars surface composition has been modeled based on averages over the measurements obtained from orbiting spacecraft and at various landing sites, taking into account the possible volatile inventory (e.g., CO2 ice, H2O ice) along with its time variation throughout the Martian year. Particle transport has been performed with the HZETRN heavy ion code. The Mars Radiation Environment Model has been made available worldwide through the Space Ionizing Radiation Effects and Shielding Tools (SIREST) website, a project of NASA Langley Research Center. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

NASA Center LaRC↗

Generational Analysis Reveals that TGF-Beta1 Inhibits the Rate of Angiogenesis in Vivo by Selective Decrease in the Number of New Vessels

Quantitative analysis of vascular generational branching demonstrated that transforming growth factor-beta1 (TGF-beta1), a multifunctional cytokine and angiogenic regulator, strongly inhibited angiogenesis in the arterial tree of the developing quail chorioallantoic membrane (CAM) by inhibition of the normal increase in the number of new, small vessels. The cytokine was applied uniformly in solution at embryonic day 7 (E7) to the CAMs of quail embryos cultured in petri dishes. After 24 h the rate of arterial growth was inhibited by as much as 105% as a function of increasing TGF-beta1 concentration. Inhibition of the rate of angiogenesis in the arterial tree by TGF-beta1 relative to controls was measured in digital images by three well-correlated, computerized methods. The first computerized method, direct measurement by the computer code VESGEN of vascular morphological parameters according to branching generations G(sub 1) through G(sub greater than or equal to 5), revealed that TGF-beta1 selectively inhibited the increase in the number density of small vessels, N(sub v greater than or equal to 5), (382 plus or minus 85 per square centimeter) for specimens treated with 1 microgram TGF-beta1/CAM for 24 h, compared to 583 plus or minus 99 per square centimeter for controls), but did not significantly affect other parameters such as average vessel length or vessel diameter. The second and third methods, the fractal dimension (D(sub f)) and grid intersection (rho (sub v)), are statistical descriptors of spatial pattern and density. According to D(sub f) and rho(sub v), arterial density increased in control specimens from 1.382 plus or minus 0.007 and 662 plus or minus 52 per square centimeters at E7 (0 h) to 1.439 plus or minus 0.013 and 884 plus or minus 55 per square centimeters at E8 (24 h), compared to 1.379 plus or minus 0.039 and 650 plus or minus 111 per square centimeter for specimens treated with 1 microgram TGF-beta1/CAM for 24 h. TGF-beta1 therefore regulates vascular pattern and the rate of angiogenesis in a unique fingerprint manner, as do other major angiogenic regulators that include VEGF, FGF-2 (bFGF), and angiostatin. TGF-beta1 did not stimulate angiogenesis significantly at low cytokine concentrations, which suggests that this quail CAM model of angiogenesis is not associated with an inflammatory response.

Parsons-Wingerter, Patricia↗

Visual scanning with or without spatial uncertainty and time-sharing performance

An experiment is reported that examines the pattern of task interference between visual scanning as a sequential and selective attention process and other concurrent spatial or verbal processing tasks. A distinction is proposed between visual scanning with or without spatial uncertainty regarding the possible differential effects of these two types of scanning on interference with other concurrent processes. The experiment required the subject to perform a simulated primary tracking task, which was time-shared with a secondary spatial or verbal decision task. The relevant information that was needed to perform the decision tasks were displayed with or without spatial uncertainty. The experiment employed a 2 x 2 x 2 design with types of scanning (with or without spatial uncertainty), expected scanning distance (low/high), and codes of concurrent processing (spatial/verbal) as the three experimental factors. The results provide strong evidence that visual scanning as a spatial exploratory activity produces greater task interference with concurrent spatial tasks than with concurrent verbal tasks. Furthermore, spatial uncertainty in visual scanning is identified to be the crucial factor in producing this differential effect.

Liu, Yili↗

Advanced electromagnetic methods for aerospace vehicles

The Advanced Helicopter Electromagnetics is centered on issues that advance technology related to helicopter electromagnetics. Progress was made on three major topics: composite materials; precipitation static corona discharge; and antenna technology. In composite materials, the research has focused on the measurements of their electrical properties, and the modeling of material discontinuities and their effect on the radiation pattern of antennas mounted on or near material surfaces. The electrical properties were used to model antenna performance when mounted on composite materials. Since helicopter platforms include several antenna systems at VHF and UHF bands, measuring techniques are being explored that can be used to measure the properties at these bands. The effort on corona discharge and precipitation static was directed toward the development of a new two dimensional Voltage Finite Difference Time Domain computer program. Results indicate the feasibility of using potentials for simulating electromagnetic problems in the cases where potentials become primary sources. In antenna technology the focus was on Polarization Diverse Conformal Microstrip Antennas, Cavity Backed Slot Antennas, and Varactor Tuned Circular Patch Antennas. Numerical codes were developed for the analysis of two probe fed rectangular and circular microstrip patch antennas fed by resistive and reactive power divider networks.

Balanis, Constantine A.↗

Pilot Injector Redesign to Reduce N+3 Cycle Emissions for a Gas-Turbine Combustor

An overview is given of an effort for the use of CFD analysis to complement design and configuration definition of third generation Lean-Direct Injection combustion concepts (LDI-3) for NASA's N plus 3 program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for a three-cup, nineteen-element flametube array with redesigned pilot injectors to improve spray and emissions characteristics when compared to a previous LDI-3 design. All computations were performed with a consistent approach to mesh-generation, spray modeling, ignition and kinetics modeling for a "medium-power" cycle condition. Computational predictions of the aerodynamics of a new pre-filming pilot injector were used to arrive at an optimized aerothermal design that meets effective area and fuel-air mixing criteria. The newly designed pilot injectors were shown to provide considerable improvements in aerodynamic stability, flame-tube pattern factor and NOx emissions, when compared to the original design.

Nitrogen Oxides↗

Interaction of two glancing, crossing shock waves with a turbulent boundary-layer at various Mach numbers

A preliminary experimental investigation was conducted to study two crossing, glancing shock waves of equal strengths, interacting with the boundary-layer developed on a supersonic wind tunnel wall. This study was performed at several Mach numbers between 2.5 and 4.0. The shock waves were created by fins (shock generators), spanning the tunnel test section, that were set at angles varying from 4 to 12 degrees. The data acquired are wall static pressure measurements, and qualitative information in the form of oil flow and schlieren visualizations. The principle aim is two-fold. First, a fundamental understanding of the physics underlying this flow phenomena is desired. Also, a comprehensive data set is needed for computational fluid dynamic code validation. Results indicate that for small shock generator angles, the boundary-layer remains attached throughout the flow field. However, with increasing shock strengths (increasing generator angles), boundary layer separation does occur and becomes progressively more severe as the generator angles are increased further. The location of the separation, which starts well downstream of the shock crossing point, moves upstream as shock strengths are increased. At the highest generator angles, the separation appears to begin coincident with the generator leading edges and engulfs most of the area between the generators. This phenomena occurs very near the 'unstart' limit for the generators. The wall pressures at the lower generator angles are nominally consistent with the flow geometries (i.e. shock patterns) although significantly affected by the boundary-layer upstream influence. As separation occurs, the wall pressures exhibit a gradient that is mainly axial in direction in the vicinity of the separation. At the limiting conditions the wall pressure gradients are primarily in the axial direction throughout.

Hingst, Warren R.↗

Pilot Injector Redesign to Reduce N+3 Cycle Emissions for a Gas-Turbine Combustor

An overview is given of an effort for the use of CFD analysis to complement design and configuration definition of third generation Lean-Direct Injection combustion concepts (LDI-3) for NASA’s N+3 program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for a three-cup, nineteen-element flame tube array with redesigned pilot injectors to improve spray and emissions characteristics when compared to a previous LDI-3 design. All computations were performed with a consistent approach to mesh-generation, spray modeling, ignition and kinetics modeling for a ‘medium-power’ cycle condition. Computational predictions of the aerodynamics of a new pre-filming pilot injector were used to arrive at an optimized aerothermal design that meets effective area and fuel-air mixing criteria. The newly designed pilot injectors were shown to provide considerable improvements in aerodynamic stability, flame-tube pattern factor and NOx emissions, when compared to the original design.

nitrogen oxides injectors↗

Software Tool for Tracking & Mapping the NASA Orion AA-2 Test Flight Ejectable Data Recorders in Real Time

On 2 July 2019, the NASA Ascent Abort 2 flight took place off the Florida coast to test the emergency systems to separate the Orion Crew Module (CM) from the future Space Launch System rocket in the event of a malfunction. During this high-altitude test, instrumentation data was recorded on twelve customized buoyant Ejectable Data Recorders (EDRs) and subsequently jettisoned from the CM in mid-air. Upon release, the EDRs activated their GPS-Iridium beacon systems and began transmitting Short Burst Data (SBD) messages via the Iridium satellite network to relay their individual location and system health information. To locate, track and retrieve each EDR from the ocean surface in real-time, multiple open-source programming tools (Python and Linux shells) were developed for parsing the incoming Iridium binary SBD messages. For this, a Linux laptop was used to receive the Iridium-generated emails containing the SBD messages and autonomously execute the parsing tools. The received SBD data contained location, timestamp and health status information that was translated, saved, and subsequently used for simultaneously generating a continuously updated color-coded tabular display summary and unique KML files used with Google Earth to track their locations. Once their locations were known, dedicated recovery vessels retrieved all EDRs from the ocean. An additional tool was also developed in order to generate 5- and 10-minute geolocation predictions for each EDR by deriving the displacement distance, elapsed time, displacement heading and velocity based on the latest known information available. The recovery vessels were also tracked with the use of a separate commercial GPS beacon system. After jettison, 67% of the EDRs transmitted valid data by the time they were retrieved from the ocean. However, the real-time information presented by the plotting tool allowed for the ready depiction of EDR dispersal patterns and reference drift trajectories, which contributed to the recovery of all twelve EDRs and the AA-2 flight data. Lastly, the available data showed that the distance between the software’s reported drift/predicted locations and the recovery locations did not exceed 38 meters, therefore demonstrating the advantages of this software tool for supporting real-time tracking and recovery efforts of beacon devices.

Moxey, Lucas↗

Boundary Layer Development on a Turbine Blade in a Linear Cascade

Several different boundary-layer development patterns for flow over the suction surface of a turbine airfoil in a linear cascade were studied and documented using a sliding surface hot-film sensor. The state of the boundary layer, whether laminar, transitional or turbulent, was determined at numerous locations along the airfoil suction surface from leading to trailing edge. Boundary-layer transition from laminar to turbulent flow through laminar separation and turbulent reattachment, or through a combination of bypass transition and strong and weak separation and turbulent reattachment, or through solely bypass transition without separation, was observed and benchmark data were recorded. Surface flow visualization and numerical boundary-layer analysis results are consistent with the hot-film data. Flow and geometry information necessary for nmerical code operation is available.

Halstead, Dave↗

A Sensitivity Study on the Effects of Particle Chemistry, Asphericity and Size on the Mass Extinction Efficiency of Mineral Dust in the Earth's Atmosphere: From the Near to Thermal IR

To determine a plausible range of mass extinction efficiencies (MEE) of terrestrial atmospheric dust from the near to thermal IR, sensitivity analyses are performed over an extended range of dust microphysical and chemistry perturbations. The IR values are subsequently compared to those in the near-IR, to evaluate spectral relationships in their optical properties. Synthesized size distributions consistent with measurements, model particle size, while composition is defined by the refractive indices of minerals routinely observed in dust, including the widely used OPAC/Hess parameterization. Single-scattering properties of representative dust particle shapes are calculated using the T-matrix, Discrete Dipole Approximation and Lorenz-Mie light-scattering codes. For the parameterizations examined, MEE ranges from nearly zero to 1.2 square meters per gram, with the higher values associated with non-spheres composed of quartz and gypsum. At near-IR wavelengths, MEE for non-spheres generally exceeds those for spheres, while in the thermal IR, shape-induced changes in MEE strongly depend on volume median diameter (VMD) and wavelength, particularly for MEE evaluated at the mineral resonant frequencies. MEE spectral distributions appear to follow particle geometry and are evidence for shape dependency in the optical properties. It is also shown that non-spheres best reproduce the positions of prominent absorption peaks found in silicates. Generally, angular particles exhibit wider and more symmetric MEE spectral distribution patterns from 8-10 micrometers than those with smooth surfaces, likely due to their edge-effects. Lastly, MEE ratios allow for inferring dust optical properties across the visible-IR spectrum. We conclude the MEE of dust aerosol are significant for the parameter space investigated, and are a key component for remote sensing applications and the study of direct aerosol radiative effects.

Hansell, R. A., Jr.↗

Nacelle Aerodynamic and Inertial Loads (NAIL) project

A flight test survey of pressures measured on wing, pylon, and nacelle surfaces and of the operating loads on Boeing 747/Pratt & Whitney JT9D-7A nacelles was made to provide information on airflow patterns surrounding the propulsion system installations and to clarify processes responsible for inservice deterioration of fuel economy. Airloads at takeoff rotation were found to be larger than at any other normal service condition because of the combined effects of high angle of attack and high engine airflow. Inertial loads were smaller than previous estimates indicated. A procedure is given for estimating inlet airloads at low speeds and high angles of attack for any underwing high bypass ratio turbofan installation approximately resembling the one tested. Flight procedure modifications are suggested that may result in better fuel economy retention in service. Pressures were recorded on the core cowls and pylons of both engine installations and on adjacent wing surfaces for use in development of computer codes for analysis of installed propulsion system aerodynamic drag interference effects.

Source record↗

Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics

A computational flow field and predicted jet noise source analysis is presented for asymmetrical fan chevrons on a modern separate flow nozzle at take off conditions. The propulsion airframe aeroacoustic asymmetric fan nozzle is designed with an azimuthally varying chevron pattern with longer chevrons close to the pylon. A baseline round nozzle without chevrons and a reference nozzle with azimuthally uniform chevrons are also studied. The intent of the asymmetric fan chevron nozzle was to improve the noise reduction potential by creating a favorable propulsion airframe aeroacoustic interaction effect between the pylon and chevron nozzle. This favorable interaction and improved noise reduction was observed in model scale tests and flight test data and has been reported in other studies. The goal of this study was to identify the fundamental flow and noise source mechanisms. The flow simulation uses the asymptotically steady, compressible Reynolds averaged Navier-Stokes equations on a structured grid. Flow computations are performed using the parallel, multi-block, structured grid code PAB3D. Local noise sources were mapped and integrated computationally using the Jet3D code based upon the Lighthill Acoustic Analogy with anisotropic Reynolds stress modeling. In this study, trends of noise reduction were correctly predicted. Jet3D was also utilized to produce noise source maps that were then correlated to local flow features. The flow studies show that asymmetry of the longer fan chevrons near the pylon work to reduce the strength of the secondary flow induced by the pylon itself, such that the asymmetric merging of the fan and core shear layers is significantly delayed. The effect is to reduce the peak turbulence kinetic energy and shift it downstream, reducing overall noise production. This combined flow and noise prediction approach has yielded considerable understanding of the physics of a fan chevron nozzle designed to include propulsion airframe aeroacoustic interaction effects.

Massey, Steven J.↗

The Effect of Wake Passing on Turbine Blade Film Cooling

The effect of upstream blade row wake passing on the showerhead film cooling performance of a downstream turbine blade has been investigated through a combination of experimental and computational studies. The experiments were performed in a steady-flow annular turbine cascade facility equipped with an upstream rotating row of cylindrical rods to produce a periodic wake field similar to that found in an actual turbine. Spanwise, chordwise, and temporal resolution of the blade surface temperature were achieved through the use of an array of nickel thin-film surface gauges covering one unit cell of showerhead film hole pattern. Film effectiveness and Nusselt number values were determined for a test matrix of various injectants, injectant blowing ratios, and wake Strouhal numbers. Results indicated a demonstratable reduction in film effectiveness with increasing Strouhal number, as well as the expected increase in film effectiveness with blowing ratio. An equation was developed to correlate the span-average film effectiveness data. The primary effect of wake unsteadiness was found to be correlated well by a chordwise-constant decrement of 0.094-St. Measurable spanwise film effectiveness variations were found near the showerhead region, but meaningful unsteady variations and downstream spanwise variations were not found. Nusselt numbers were less sensitive to wake and injection changes. Computations were performed using a three-dimensional turbulent Navier-Stokes code which was modified to model wake passing and film cooling. Unsteady computations were found to agree well with steady computations provided the proper time-average blowing ratio and pressure/suction surface flow split are matched. The remaining differences were isolated to be due to the enhanced mixing in the unsteady solution caused by the wake sweeping normally on the pressure surface. Steady computations were found to be in excellent agreement with experimental Nusselt numbers, but to overpredict experimental film effectiveness values. This is likely due to the inability to match actual hole exit velocity profiles and the absence of a credible turbulence model for film cooling.

Heidmann, James David↗

Rapid Prototyping Integrated With Nondestructive Evaluation and Finite Element Analysis

Most reverse engineering approaches involve imaging or digitizing an object then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. Rapid prototyping (RP) refers to the practical ability to build high-quality physical prototypes directly from computer aided design (CAD) files. Using rapid prototyping, full-scale models or patterns can be built using a variety of materials in a fraction of the time required by more traditional prototyping techniques (refs. 1 and 2). Many software packages have been developed and are being designed to tackle the reverse engineering and rapid prototyping issues just mentioned. For example, image processing and three-dimensional reconstruction visualization software such as Velocity2 (ref. 3) are being used to carry out the construction process of three-dimensional volume models and the subsequent generation of a stereolithography file that is suitable for CAD applications. Producing three-dimensional models of objects from computed tomography (CT) scans is becoming a valuable nondestructive evaluation methodology (ref. 4). Real components can be rendered and subjected to temperature and stress tests using structural engineering software codes. For this to be achieved, accurate high-resolution images have to be obtained via CT scans and then processed, converted into a traditional file format, and translated into finite element models. Prototyping a three-dimensional volume of a composite structure by reading in a series of two-dimensional images generated via CT and by using and integrating commercial software (e.g. Velocity2, MSC/PATRAN (ref. 5), and Hypermesh (ref. 6)) is being applied successfully at the NASA Glenn Research Center. The building process from structural modeling to the analysis level is outlined in reference 7. Subsequently, a stress analysis of a composite cooling panel under combined thermomechanical loading conditions was performed to validate this process.

Abdul-Aziz, Ali↗

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid (s/1) interface shape and convection patterns during solidification of pure Al and Al-4 wt% Cu alloy. The objective of the numerical task was to predict the influence of convective velocity on an insoluble particle near a s/l interface. These predictions would then be used to define the minimum gravity level (g) required to investigate the fundamental physics of interaction between a particle and a s/I interface. To satisfy this objective, steady state calculations were performed for different gravity levels and orientations with the gravity vector. ne furnace configuration used in this analysis is the proposed International Space Station Furnace, Quench Module Insert (QMI) 1. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. General model of binary alloy solidification was based on the finite element code FIDAP. It was found that for the worst case orientation of 90 degrees with the gravity vector and a g level of 10(exp -4)g(sub o) (g(sub o) = 9.8 m/s(exp 2)) the dominant forces acting on the particle would be the fundamental drag and interfacial forces.

Bune, Andris V.↗