Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable geometry”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Probabilistic Micromechanics and Macromechanics for Ceramic Matrix Composites

The properties of ceramic matrix composites (CMC's) are known to display a considerable amount of scatter due to variations in fiber/matrix properties, interphase properties, interphase bonding, amount of matrix voids, and many geometry- or fabrication-related parameters, such as ply thickness and ply orientation. This paper summarizes preliminary studies in which formal probabilistic descriptions of the material-behavior- and fabrication-related parameters were incorporated into micromechanics and macromechanics for CMC'S. In this process two existing methodologies, namely CMC micromechanics and macromechanics analysis and a fast probability integration (FPI) technique are synergistically coupled to obtain the probabilistic composite behavior or response. Preliminary results in the form of cumulative probability distributions and information on the probability sensitivities of the response to primitive variables for a unidirectional silicon carbide/reaction-bonded silicon nitride (SiC/RBSN) CMC are presented. The cumulative distribution functions are computed for composite moduli, thermal expansion coefficients, thermal conductivities, and longitudinal tensile strength at room temperature. The variations in the constituent properties that directly affect these composite properties are accounted for via assumed probabilistic distributions. Collectively, the results show that the present technique provides valuable information about the composite properties and sensitivity factors, which is useful to design or test engineers. Furthermore, the present methodology is computationally more efficient than a standard Monte-Carlo simulation technique; and the agreement between the two solutions is excellent, as shown via select examples.

Murthy, Pappu L. N.↗
NASA NTRS2001

Automated Approach to Very High-Order Aeroacoustic Computations

Computational aeroacoustics requires efficient, high-resolution simulation tools. For smooth problems, this is best accomplished with very high-order in space and time methods on small stencils. However, the complexity of highly accurate numerical methods can inhibit their practical application, especially in irregular geometries. This complexity is reduced by using a special form of Hermite divided-difference spatial interpolation on Cartesian grids, and a Cauchy-Kowalewski recursion procedure for time advancement. In addition, a stencil constraint tree reduces the complexity of interpolating grid points that am located near wall boundaries. These procedures are used to develop automatically and to implement very high-order methods (> 15) for solving the linearized Euler equations that can achieve less than one grid point per wavelength resolution away from boundaries by including spatial derivatives of the primitive variables at each grid point. The accuracy of stable surface treatments is currently limited to 11th order for grid aligned boundaries and to 2nd order for irregular boundaries.

Dyson, Rodger W.↗
NASA NTRS2012

Polarimetric Retrievals of Surface and Cirrus Clouds Properties in the Region Affected by the Deepwater Horizon Oil Spill

In 2010, the Goddard Institute for Space Studies (GISS) Research Scanning Polarimeter (RSP) performed several aerial surveys over the region affected by the oil spill caused by the explosion of the Deepwater Horizon offshore platform. The instrument was deployed on the NASA Langley B200 aircraft together with the High Spectral Resolution Lidar (HSRL), which provides information on the distribution of the aerosol layers beneath the aircraft, including an accurate estimate of aerosol optical depth. This work illustrates the merits of polarization measurements in detecting variations of ocean surface properties linked to the presence of an oil slick. In particular, we make use of the degree of linear polarization in the glint region, which is severely affected by variations in the refractive index but insensitive to the waviness of the water surface. Alterations in the surface optical properties are therefore expected to directly affect the polarization response of the RSP channel at 2264 nm, where both molecular and aerosol scattering are negligible and virtually all of the observed signal is generated via Fresnel reflection at the surface. The glint profile at this wavelength is fitted with a model which can optimally estimate refractive index, wind speed and direction, together with aircraft attitude variations affecting the viewing geometry. The retrieved refractive index markedly increases over oil-contaminated waters, while the apparent wind speed is significantly lower than in adjacent uncontaminated areas, suggesting that the slick dampens high-frequency components of the ocean wave spectrum. The constraint on surface reflectance provided by the short-wave infrared channels is a cornerstone of established procedures to retrieve atmospheric aerosol microphysical parameters based on the inversion of the RSP multispectral measurements. This retrieval, which benefits from the ancillary information provided by the HSRL, was in this specific case hampered by prohibitive variability in atmospheric conditions (very inhomogeneous aerosol distribution and cloud cover). Although the results presented for the surface are essentially unaffected, we discuss the results obtained by typing algorithms in sorting the complex mix of aerosol types, and show evidence of oriented ice in cirrus clouds present in the area. In this context, polarization measurements at 1880 nm were used to infer ice habit and cirrus optical depth, which was found in the subvisual/threshold-visible regime, confirming the utility of the aforementioned RSP channel for the remote sensing of even thin cold clouds.

Polarization↗
NASA NTRS2018

Parametric Analysis of Aircraft WingWeight Using Low-Order Physics-Based Analysis

In the conceptual aircraft design phase, prediction of the empty weight typically relies on empirically-based regression equations which execute quickly and require little detailed information about the internal structural layout. Since they are based on existing aircraft, however, empirical methods can lose their validity for newer technologies and unconventional configurations. Designers can transition to higher-order, physics-based analysis methods to improve the accuracy of the weight prediction, but at the cost of complex model setup and increased computational time. This paper describes a methodology for low-order aero-structural analysis of conceptual aircraft configurations that increases the use of physics-based analysis in conceptual design, but is less complex and time-consuming than higher-order methods such as finite-element analysis. The methodology uses Vehicle Sketch Pad (OpenVSP) to model the aircraft geometry, and ASWING to perform the aero-structural analysis. The internal forces and moments from the ASWING analysis are post-processed to calculate the resulting direct and shear stresses in the structure, and the thickness distributions of the aircraft components are varied to match the maximum von Mises stress at each cross section to the material allowable. To offset the increased computational time relative to empirical weight equations, a process is studied which uses parametric variation to develop a regression equation relating the weight of the aircraft wing to major design variables. This new weight equation is similar to existing empirical equations, but is built using the more physics-based methodology; the new equation could be used to augment or replace portions of the empirical database to improve the validity of the wing weight prediction for unconventional configurations and advanced technologies.

Olson, Erik D.↗
NASA NTRS2006

The Character and Variability of the Eta Carinae Wind Lines

The binarity of Eta Carinae has been debated for a long time. We have searched for more evidence for a companion star in a spectroscopic investigation of the Eta Carinae stellar wind lines, using moderate spectral and high angular resolution HST/STIS data. Over Eta Carinae's 5.54 year spectroscopic period many of the observable wind lines in the NUV/Optical spectral region exhibit peculiar line profiles with unusual velocity shifts relative to the system velocity. Some of the lines are exclusively blue-shifted over the entire cycle. Their ionization/excitation imply formation not in the stellar wind but rather in the interface between the two massive stars. We have analyzed velocity and intensity variations over the spectroscopic period and interpreted what the variations tell us about the geometry of the nebular structure close to Eta Carinae.

Nielsen, K. E.↗
NASA NTRS1974

Geodesic synchrotron radiation in the Kerr geometry by the method of asymptotically factorized Green's functions

The scalar, electromagnetic, and gravitational geodesic-synchrotron-radiation (GSR) spectra are determined for the case of a test particle moving on a highly relativistic circular orbit about a rotating (Kerr) black hole. It is found that the spectral shape depends only weakly on the value of the angular-momentum parameter (a/M) of the black hole, but the total radiated power drops unexpectedly for a value of at least 0.95 and vanishes as the value approaches unity. A spin-dependent factor (involving the inner product of the polarization of a radiated quantum with the source) is isolated to explain the dependence of the spectral shape on the spin of the radiated field. Although the scalar wave equation is solved by separation of variables, this procedure is avoided for the vector and tensor cases by postulating a sum-over-states expansion for the Green's function similar to that found to hold in the scalar case. The terms in this sum, significant for GSR, can then be evaluated in the geometric-optics approximation without requiring the use of vector or tensor spherical harmonics.

Chrzanowski, P. L.↗
NASA NTRS2005

Tracing Material Flow Paths in Friction Stir Welds

Heat and mechanical work are coupled in the friction stir welding process. The process variables are RPM, translational weld speed, and downward plunge force. The strain-temperature history of a metal element at each point on the cross-section of the weld is determined by the process variables plus the individual flow path taken by the particular filament of metal flowing around the tool and ending on flat point. The strain-temperature history determines the properties of a metal element on the weld cross-section. The strain-temperature history is carefully controlled in metal processes where direct control is feasible. Indirect estimates of the flow paths and the strain-temperature histories of filaments comprising friction stir welds can be made from a model, if the model is good enough. This paper describes marker studies of flow path geometries for various process parameters. Observed geometries are compared with geometries estimated from models.

Sanders, Johnny↗
NASA NTRS2013

Aerodynamic Investigation of Incidence Angle Effects in a Large Scale Transonic Turbine Cascade

Aerodynamic measurements showing the effects of large incidence angle variations on an HPT turbine blade set are presented. Measurements were made in NASA's Transonic Turbine Blade Cascade Facility which has been used in previous studies to acquire detailed aerodynamic and heat transfer measurements for CFD code validation. The current study supports the development of variable-speed power turbine (VSPT) speed-change technology for the NASA Large Civil Tilt Rotor (LCTR) vehicle. In order to maintain acceptable main rotor propulsive efficiency, the VSPT operates over a nearly 50 percent speed range from takeoff to altitude cruise. This results in 50deg or more variations in VSPT blade incidence angles. The cascade facility has the ability to operate over a wide range of Reynolds numbers and Mach numbers, but had to be modified in order to accommodate the negative incidence angle variation required by the LCTR VSPT operation. Using existing blade geometry with previously acquired aerodynamic data, the tunnel was re-baselined and the new incidence angle range was exercised. Midspan exit total pressure and flow angle measurements were obtained at seven inlet flow angles. For each inlet angle, data were obtained at five flow conditions with inlet Reynolds numbers varying from 6.83×10(exp 5) to 0.85×10(exp 5) and two isentropic exit Mach numbers of 0.74 and 0.34. The midspan flowfield measurements were acquired using a three-hole pneumatic probe located in a survey plane 8.6 percent axial chord downstream of the blade trailing edge plane and covering three blade passages. Blade and endwall static pressure distributions were also acquired for each flow condition.

McVetta, Ashlie B.↗
NASA NTRS2012

Aerodynamic Investigation of Incidence Angle Effects in a Large Scale Transonic Turbine Cascade

Aerodynamic measurements showing the effects of large incidence angle variations on an HPT turbine blade set are presented. Measurements were made in NASA's Transonic Turbine Blade Cascade Facility which has been used in previous studies to acquire detailed aerodynamic and heat transfer measurements for CFD code validation. The current study supports the development of variable-speed power turbine (VSPT) speed-change technology for the NASA Large Civil Tilt Rotor (LCTR) vehicle. In order to maintain acceptable main rotor propulsive efficiency, the VSPT operates over a nearly 50% speed range from takeoff to altitude cruise. This results in 50 degrees or more variations in VSPT blade incidence angles. The cascade facility has the ability to operate over a wide range of Reynolds numbers and Mach numbers, but had to be modified in order to accommodate the negative incidence angle variation required by the LCTR VSPT operation. Using existing blade geometry with previously acquired aerodynamic data, the tunnel was re-baselined and the new incidence angle range was exercised. Midspan exit total pressure and flow angle measurements were obtained at seven inlet flow angles. For each inlet angle, data were obtained at five flow conditions with inlet Reynolds numbers varying from 6.83 × 10(exp 5) to 0.85 ×10(exp 5) and two isentropic exit Mach numbers of 0.74 and 0.34. The midspan flowfield measurements were acquired using a three-hole pneumatic probe located in a survey plane 8.6% axial chord downstream of the blade trailing edge plane and covering three blade passages. Blade and endwall static pressure distributions were also acquired for each flow condition

McVetta, Ashlie B.↗
NASA NTRS2014

Aerodynamic Investigation of Incidence Angle Effects in a Large Scale Transonic Turbine Cascade

Aerodynamic measurements showing the effects of large incidence angle variations on an HPT turbine blade set are presented. Measurements were made in NASA's Transonic Turbine Blade Cascade Facility which has been used in previous studies to acquire detailed aerodynamic and heat transfer measurements for CFD code validation. The current study supports the development of variable-speed power turbine (VSPT) speed-change technology for the NASA Large Civil Tilt Rotor (LCTR) vehicle. In order to maintain acceptable main rotor propulsive efficiency, the VSPT operates over a nearly 50 percent speed range from takeoff to altitude cruise. This results in 50 deg or more variations in VSPT blade incidence angles. The cascade facility has the ability to operate over a wide range of Reynolds numbers and Mach numbers, but had to be modified in order to accommodate the negative incidence angle variation required by the LCTR VSPT operation. Using existing blade geometry with previously acquired aerodynamic data, the tunnel was re-baselined and the new incidence angle range was exercised. Midspan exit total pressure and flow angle measurements were obtained at seven inlet flow angles. For each inlet angle, data were obtained at five flow conditions with inlet Reynolds numbers varying from 6.83×10 (exp 5) to 0.85×10(exp 5) and two isentropic exit Mach numbers of 0.74 and 0.34. The midspan flowfield measurements were acquired using a three-hole pneumatic probe located in a survey plane 8.6 percent axial chord downstream of the blade trailing edge plane and covering three blade passages. Blade and endwall static pressure distributions were also acquired for each flow condition.

McVetta, Ashlie B.↗
NASA NTRS1975

Design, fabrication and acoustic tests of a 36 inch (0.914 meter) statorless turbotip fan

The LF336/E is a 36 inch (0.914 meter) diameter fan designed to operate in a rotor-alone configuration. Design features required for modification of the existing LF336/A rotor-stator fan into the LF336/E statorless fan configuration are discussed. Tests of the statorless fan identified an aerodynamic performance deficiency due to inaccurate accounting of the fan exit swirl during the aerodynamic design. This performance deficiency, related to fan exit static pressure levels, produced about a 20 percent thrust loss. A study was then conducted for further evaluation of the fan exit flow fields typical of statorless fan systems. This study showed that through proper selection of fan design variables such as pressure ratio, radius ratio, and swirl distributions, performance of a statorless fan configuration could be improved with levels of thrust approaching the conventional rotor-stator fan system. Acoustic measurements were taken for the statorless fan system at both GE and NASA, and when compared to other lift fan systems, showed noise levels comparable to the quietest lift fan configuration which included rotor-stator spacing and acoustic treatment. The statorless fan system was also used to determine effects of rotor leading edge serrations on noise generations. A cascade test program identified the serration geometry based on minimum pressure losses, wake turbulence levels and noise generations.

Smith, E. G.↗
NASA NTRS1991

CCD photometry of 2060 Chiron, 1991 January

Observations of 2060 Chiron was performed on 7 to 8 Jan. 1991 with the Mt. Palomar 1.52 m telescope in the Gunn-R passband. On-chip field stars were used to perform differential reductions. The repeatability of the 5.9 hour light curve was excellent, both within a night and from night to night. No evidence for short-term secular variations similar to those seen last year by both Luu and Jewitt (1990) and Buratti and Dunbar (1991) is seen in the new light curve. Chiron's rotational light curve appears strikingly similar to that obtained a year earlier by Luu and Jewitt (1990), both in amplitude and shape. Both light curves show strongly correlated changes over a timescale of perhaps 15 minutes. These same features were marginally visible in the 1986 light curve. Such behavior is believed to be evidence that Chiron may be more aspherical than the 4 percent intensity variation might otherwise indicate, and favors a viewing geometry where the subearth latitude is rather low. Chiron was much fainter in 1985, when a partial light curve was obtained by Marcialis. Due to the lower sampling rate of these early data, no conclusions can be made regarding the high-frequency light curve structure back then. All three of these light curves differ significantly from that obtained by Buratti and Dunbar (1991), one week before the observations of Luu and Jewitt. The Chiron field was calibrated using Landolt standards on Ut 15 Mar. 1991. A mean R-magnitude of 15.6 + or - 0.1 was found. Variability of 2060 Chiron was demonstrated over timescales of minutes, hours, and years. An intense campaign was urged to monitor the photometric behavior of Chiron throughout the 1990s.

Buratti, B. J.↗
NASA NTRS2022

Simulated Lunar Surface Hydration Measurements using Multispectral Lidar at 3 µm

Accurately measuring the variability of spectroscopic signatures of hydration (H2O + OH) on the illuminated lunar surface at 3 µm as a function of latitude, lunar time of day, and composition is crucial to determining the generation and destruction mechanisms of OH species and understanding the global water cycle. A prime complication in analysis of the spectroscopic feature is the accurate removal of thermal emission, which can modify or even eliminate the hydration feature depending on the data processing methods used and assumptions made. An orbital multispectral lidar, with laser illumination at key diagnostic wavelengths, would provide uniform, zero-phase geometry, complete latitude and time of day coverage from a circular polar orbit, and is agnostic to the thermal state of the surface. We have performed measurement simulations of a four-wavelength multispectral lidar using spectral mixtures of hydrated mid-ocean ridge basalt (MORB) glasses and lunar regolith endmembers to assess the lidar performance in measuring hydration signatures on the lunar surface. Our results show a feasible system with wavelengths at 1.5 µm, 2.65 µm, 2.8 µm, and 3.1 µm can measure lunar hydration with a precision of 52 ppm (1σ) or better. These results, combined with the uniform measurement capabilities of multispectral lidar make it a valuable spectroscopic technique for elucidating mechanisms of OH/H2O generation, migration, and destruction.

D. R. Cremons↗
NASA NTRS1992

Adaptive multigrid domain decomposition solutions for viscous interacting flows

Several viscous incompressible flows with strong pressure interaction and/or axial flow reversal are considered with an adaptive multigrid domain decomposition procedure. Specific examples include the triple deck structure surrounding the trailing edge of a flat plate, the flow recirculation in a trough geometry, and the flow in a rearward facing step channel. For the latter case, there are multiple recirculation zones, of different character, for laminar and turbulent flow conditions. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable, non-staggered grid computation.

Rubin, Stanley G.↗
NASA NTRS2016

Terrestrial Testing of the CapiBRIC, a Microgravity Optimized Brine Processor

Utilizing geometry based static phase separation exhibited in the radial vaned capillary drying tray, a system was conceived to recover water from brine. This technology has been named the Capillary BRIC; abbreviated CapiBRIC. The CapiBRIC utilizes a capillary drying tray within a drying chamber. Water is recovered from clean water vapor evaporating from the free surface leaving waste brine solids behind. A novel approach of optimizing the containment geometry to support passive capillary flow and static phase separation provides the opportunity for a low power system that is not as susceptible to fouling as membranes or other technologies employing physical barriers across the free brine surface to achieve phase separation in microgravity. Having been optimized for operation in microgravity, full-scale testing of the CapiBRIC as designed cannot be performed on the ground as the force of gravity would dominate over the capillary forces. However, subscale units relevant to full-scale design were used to characterize fill rates, containment stability, and interaction with a variable volume reservoir in the PSU Dryden Drop Tower (DDT) facility. PSU also using tested units scaled such that capillary forces dominated in a 1-g environment to characterize evaporation from a free-surface in 1-g upward, sideways and downward orientations. In order to augment the subscale testing performed by PSU, a full scale 1-g analogue of the CapiBRIC drying unit was initiated to help validate performance predictions regarding expected water recovery ratio, estimated processing time, and interface definitions for inlets, outlets, and internal processes, including vent gas composition. This paper describes the design, development and test of the terrestrial CapiBRIC prototypes.

Sargusingh, Miriam J.↗
NASA NTRS2018

Reduction in 317-780 nm Radiance Reflected from the Sunlit Earth During the Eclipse of 21 August 2017

Ten wavelength channels of calibrated radiance image data from the sunlit Earth are obtained every 65 min during Northern Hemisphere summer from the EPIC (Earth Polychromatic Imaging Camera) instrument on the DSCOVR (Deep Space Climate Observatory) satellite located near the Earth-Sun Lagrange 1 point (L1), about 1.5 million km from the Earth. The L1 location permitted seven observations of the Moon's shadow on the Earth for about 3 h during the 21 August 2017 eclipse. Two of the observations were timed to coincide with totality over Casper, Wyoming, and Columbia, Missouri. Since the solar irradiances within five channels (λi = 388, 443, 551, 680, and 780 nm) are not strongly absorbed in the atmosphere, they can be used for characterizing the eclipse reduction in reflected radiances for the Earth's sunlit face containing the eclipse shadow. Five channels (λi = 317:5, 325, 340, 688, and 764 nm) that are partially absorbed in the atmosphere give consistent reductions compared to the non-absorbed channels. This indicates that cloud reflectivities dominate the 317.5-780 nm radiances reflected back to space from the sunlit Earth's disk with a significant contribution from Rayleigh scattering for the shorter wavelengths. An estimated reduction of 10% was obtained for spectrally integrated radiance (387 to 781 nm) reflected from the sunlit Earth towards L1 for two sets of observations on 21 August 2017, while the shadow was in the vicinity of Casper, Wyoming (42.8666° N, 106.3131° W; centered on 17:44:50 UTC), and Columbia, Missouri (38.9517° N, 92.3341° W; centered on 18:14:50 UTC). In contrast, when non-eclipse days (20 and 23 August) are compared for each wavelength channel, the change in reflected light is much smaller (less than 1% for 443 nm compared to 9% (Casper) and 8% (Columbia) during the eclipse). Also measured was the ratio R(sub EN)(λi) of reflected radiance on adjacent non-eclipse days divided by radiances centered in the eclipse totality region with the same geometry for all 10 wavelength channels. The measured R(sub EN)(443 nm) was smaller for Columbia (169) than for Casper (935), because Columbia had more cloud cover than Casper. R(sub EN)(λi) forms a useful test of a 3-D radiative transfer models for an eclipse in the presence of optically thin clouds. Specific values measured at Casper with thin clouds are R(sub EN)(340 nm)=475, R(sub EN)(388 nm)=3500, REN(443 nm)=935, REN(551 nm)=5455, REN(680 nm)=220, and R(sub EN)(780 nm)=395. Some of the variability is caused by changing cloud amounts within the moving region of totality during the 2.7 min needed to measure all 10 wavelength channels.

DSCOVR↗
NASA NTRSDate not supplied

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models.The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar convexwall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University. The favorable pressure gradient is induced bythe streamwise surface curvature. The DNS results will be compared with experimental measurements and RANS computations for the same flow conditions and wall geometry. The DNS datasets will also be analyzed to provide an improved physical understanding of hypersonic boundary-layer turbulence subject to mechanical nonequilibrium, including the validity of Morkovin’s hypothesis and the various available compressibility transformations, as well as to obtain turbulence statistics that are relevant to RANS modeling, including the distribution ofeddy viscosity, the budgets of turbulent kinetic energy and Reynolds stresses, and the velocity-temperature correlations. The RANS predictions will be obtained by using both the commonly used one- and two-equation models that assume a constant turbulent Prandtl number and an algebraic, variable Prandtl number model for the turbulent energy flux.

turbulent boundary layers↗
NASA NTRS2012

Journal and Wave Bearing Impedance Calculation Software

The wave bearing software suite is a MALTA application that computes bearing properties for user-specified wave bearing conditions, as well as plain journal bearings. Wave bearings are fluid film journal bearings with multi-lobed wave patterns around the circumference of the bearing surface. In this software suite, the dynamic coefficients are outputted in a way for easy implementation in a finite element model used in rotor dynamics analysis. The software has a graphical user interface (GUI) for inputting bearing geometry parameters, and uses MATLAB s structure interface for ease of interpreting data. This innovation was developed to provide the stiffness and damping components of wave bearing impedances. The computational method for computing bearing coefficients was originally designed for plain journal bearings and tilting pad bearings. Modifications to include a wave bearing profile consisted of changing the film thickness profile given by an equation, and writing an algorithm to locate the integration limits for each fluid region. Careful consideration was needed to implement the correct integration limits while computing the dynamic coefficients, depending on the form of the input/output variables specified in the algorithm.

Hanford, Amanda↗