Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “missing values”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Scheduling for Urban Air Mobility using Safe Learning

This work considers the scheduling problem for Urban Air Mobility (UAM) vehicles travelling between origin-destination pairs with both hard and soft trip deadlines. Each route is described by a discrete probability distribution over trip completion times (or delay) and over interarrival times of requests (or demand) for the route along with a fixed hard or soft deadline. Soft deadlines carry a cost that is incurred when the deadline is missed. An online, safe scheduler is developed that ensures that hard deadlines are never missed and that average cost of missing soft deadlines is minimized. The system is modelled as a Markov Decision Process (MDP) and safe model based learning is used to find the probabilistic distributions over route delays and demand. Monte Carlo Tree Search (MCTS) Earliest Deadline First (EDF) is used to safely explore the learned models in an online fashion and develop a near-optimal non-preemptive scheduling policy. These results are compared with Value Iteration (VI) and MCTS (Random) scheduling solutions.

Urban Air Mobility↗
DOE OSTI2026

Uniting Theory and Experiment to Deliver Flexible MOFS for Superior Methane (NG) Storage

The objective of the project was to use previous insights developed through synthesis and quantitative modeling of rigid metal–organic frameworks (MOFs) in an established synergistic theoretical/experimental team to create, modify, and evaluate flexible MOFs (FlexMOFs) for natural gas (NG) storage and release at practically useful pressures and transform the NG storage economy. The goal was reduced pressure absorbed natural gas (ANG) FlexMOF storage at operating pressures less than 100 bar with physisorption exploiting the favorable thermodynamics and kinetics of flexible porous material opening in response to adsorption. Hydrogen behavior in this context was also considered. The specific aim of the project was to design and develop a standard computational modeling methodology for, first, detailed atomistic retrodiction of FlexMOF gating behavior and ultimately prediction of the effects of functionalization and/or substitution on structural transitioning. The project was also geared towards the establishment of the interaction of methane with the framework and binding sites using modeling and the FlexMOF will be both internally and external validated as a SMART metric. This has value to the scientific community both from the obvious standpoint of providing a better understanding of the promising test case systems (CdIF-13 and the MIL-53(Al) series of MOFs), but also in providing an avenue of obtaining insight into these FlexMOF systems in general, which is of particular interest given the tendency of structure-function correlation to lag behind synthesis methodology (making the latter a hit or miss proposition for applications). This gap remains significant for FlexMOF systems whose gate opening behaviors complicate computational examination. The computational methodologies are relatively inexpensive in terms of both money and computational resources enhancing the general viability of these methodologies. The ultimate gain to the public will be in the application of these techniques to design systems for natural gas storage and use to cut down on green-house emissions.

03 NATURAL GAS↗
NASA NTRS1977

Mariner Venus Mercury 1973 S/X-band experiment

The S/X-band experiment on the Mariner Venus/Mercury 1973 spacecraft constituted a unique opportunity to demonstrate the capability of an X-band downlink coherent with the normal S-band downlink. This was both a technological and scientific experiment, and the results indicated that it was successful in both cases. Analysis of the tracking data shows that the new S/X data type was capable of reducing the miss distance at the planet Mercury by 80% (post-processed data). The use of S/X electron content was demonstrated by comparison with Faraday rotation data. An X-band turnaround telemetry experiment showed the feasibility of a planetary X-band link. In the science area, the model atmospheric environment of Venus was refined. The ionosphere of the planet was measured to a higher accuracy than before, and the value of the dual-frequency link for measuring the scale size of turbulence was demonstrated. The estimate of the scale size was increased from 100 m to above 5 km.

Levy, G. S.↗
NASA NTRS2014

Statistics Clinic

Do you have elevated p‐values? Is the data analysis process getting you down? Do you experience anxiety when you need to respond to criticism of statistical methods in your manuscript? You may be suffering from Insufficient Statistical Support Syndrome (ISSS). For symptomatic relief of ISSS, come for a free consultation with JSC biostatisticians at our help desk during the poster sessions at the HRP Investigators Workshop. Get answers to common questions about sample size, missing data, multiple testing, when to trust the results of your analyses and more. Side effects may include sudden loss of statistics anxiety, improved interpretation of your data, and increased confidence in your results.

Feiveson, Alan H.↗
NASA NTRS2014

Data Preservation, Information Preservation, and Lifecyle of Information Management at NASA GES DISC

Data lifecycle management awareness is common today; planners are more likely to consider lifecycle issues at mission start. NASA remote sensing missions are typically subject to life cycle management plans of the Distributed Active Archive Center (DAAC), and NASA invests in these national centers for the long-term safeguarding and benefit of future generations. As stewards of older missions, it is incumbent upon us to ensure that a comprehensive enough set of information is being preserved to prevent the risk for information loss. This risk is greater when the original data experts have moved on or are no longer available. Preservation of items like documentation related to processing algorithms, pre-flight calibration data, or input-output configuration parameters used in product generation, are examples of digital artifacts that are sometimes not fully preserved. This is the grey area of information preservation; the importance of these items is not always clear and requires careful consideration. Missing important metadata about intermediate steps used to derive a product could lead to serious challenges in the reproducibility of results or conclusions. Organizations are rapidly recognizing that the focus of life-cycle preservation needs to be enlarged from the strict raw data to the more encompassing arena of information lifecycle management. By understanding what constitutes information, and the complexities involved, we are better equipped to deliver longer lasting value about the original data and derived knowledge (information) from them. The NASA Earth Science Data Preservation Content Specification is an attempt to define the content necessary for long-term preservation. It requires new lifecycle infrastructure approach along with content repositories to accommodate artifacts other than just raw data. The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) setup an open-source Preservation System capable of long-term archive of digital content to augment its raw data holding. This repository is being used for such missions as HIRDLS, UARS, TOMS, OMI, among others. We will provide a status of this implementation; report on challenges, lessons learned, and detail our plans for future evolution to include other missions and services.

data management↗
NASA NTRS2016

Tropical Cyclone Activity in the North Atlantic Basin During the Weather Satellite Era, 1960-2014

This Technical Publication (TP) represents an extension of previous work concerning the tropical cyclone activity in the North Atlantic basin during the weather satellite era, 1960-2014, in particular, that of an article published in The Journal of the Alabama Academy of Science. With the launch of the TIROS-1 polar-orbiting satellite in April 1960, a new era of global weather observation and monitoring began. Prior to this, the conditions of the North Atlantic basin were determined only from ship reports, island reports, and long-range aircraft reconnaissance. Consequently, storms that formed far from land, away from shipping lanes, and beyond the reach of aircraft possibly could be missed altogether, thereby leading to an underestimate of the true number of tropical cyclones forming in the basin. Additionally, new analysis techniques have come into use which sometimes has led to the inclusion of one or more storms at the end of a nominal hurricane season that otherwise would not have been included. In this TP, examined are the yearly (or seasonal) and 10-year moving average (10-year moving average) values of the (1) first storm day (FSD), last storm day (LSD), and length of season (LOS); (2) frequencies of tropical cyclones (by class); (3) average peak 1-minute sustained wind speed ( ) and average lowest pressure ( ); (4) average genesis location in terms of north latitudinal ( ) and west longitudinal ( ) positions; (5) sum and average power dissipation index ( ); (6) sum and average accumulated cyclone energy ( ); (7) sum and average number of storm days ( ); (8) sum of the number of hurricane days (NHD) and number of major hurricane days (NMHD); (9) net tropical cyclone activity index (NTCA); (10) largest individual storm (LIS) PWS, LP, PDI, ACE, NSD, NHD, NMHD; and (11) number of category 4 and 5 hurricanes (N4/5). Also examined are the December-May (D-M) and June-November (J-N) averages and 10-year moving average values of several climatic factors, including the (1) oceanic Nino index ( ); (2) Atlantic multi-decadal oscillation ( ) index; (3) Atlantic meridional mode ( ) index; (4) global land-ocean temperature index ( ); and (5) quasi-biennial oscillation ( ) index. Lastly, the associational aspects (using both linear and nonparametric statistical tests) between selected tropical cyclone parameters and the climatic factors are examined based on their 10-year moving average trend values.

satellite↗
NASA NTRS1993

Multiple steady states in atmospheric chemistry

The equations describing the distributions and concentrations of trace species are nonlinear and may thus possess more than one solution. This paper develops methods for searching for multiple physical solutions to chemical continuity equations and applies these to subsets of equations describing tropospheric chemistry. The calculations are carried out with a box model and use two basic strategies. The first strategy is a 'search' method. This involves fixing model parameters at specified values, choosing a wide range of initial guesses at a solution, and using a Newton-Raphson technique to determine if different initial points converge to different solutions. The second strategy involves a set of techniques known as homotopy methods. These do not require an initial guess, are globally convergent, and are guaranteed, in principle, to find all solutions of the continuity equations. The first method is efficient but essentially 'hit or miss' in the sense that it cannot guarantee that all solutions which may exist will be found. The second method is computationally burdensome but can, in principle, determine all the solutions of a photochemical system. Multiple solutions have been found for models that contain a basic complement of photochemical reactions involving O(x), HO(x), NO(x), and CH4. In the present calculations, transitions occur between stable branches of a multiple solution set as a control parameter is varied. These transitions are manifestations of hysteresis phenomena in the photochemical system and may be triggered by increasing the NO flux or decreasing the CH4 flux from current mean tropospheric levels.

Stewart, Richard W.↗
NASA NTRS1997

Evidence from Quasi-Periodic Oscillations for a Millisecond Pulsar in the Low Mass X-Ray Binary 4U 0614+091

We have detected quasi-periodic oscillations (QPOs) near 1 kHz from the low mass X-ray binary 4U 0614+091 in observations with RXTE. The observations span several months and sample the source over a large range of X-ray luminosity. In every interval QPOs are present above 400 Hz with fractional RMS amplitudes from 3 to 12% over the full PCA band. At high count rates, two high frequency QPOs are detected simultaneously. The difference of their frequency centroids is consistent with a constant value of 323 Hz in all observations. During one interval a third signal is detected at 328 +/- 2 Hz. This suggests the system has a stable 'clock' which is most likely the neutron star with spin period 3.1 msec. Thus, our observations of 4U 0614+091 and those of 4U 1728-34 provide the first evidence for millisecond pulsars within low-mass X-ray binary systems and reveal the 'missing-link' between millisecond radiopulsars and the late stages of binary evolution in low mass X-ray binaries. The constant difference of the high frequency QPOs sug,,ests a beat-frequency interpretation. In this model, the high frequency QPO is associated with the Keplerian frequency of the inner accretion disk and the lower frequency QPO is a 'beat' between the differential rotation frequency of the inner disk and the spinning neutron star. Assuming the high frequency QPO is a Keplerian orbital frequency for the accretion disk, we find a maximum mass of 1.9 solar mass and a maximum radius of 17 km for the neutron star.

Ford, E.↗
NASA NTRS2005

Real-Time Simulation for Verification and Validation of Diagnostic and Prognostic Algorithms

To verify that a health management system (HMS) performs as expected, a virtual system simulation capability, including interaction with the associated platform or vehicle, very likely will need to be developed. The rationale for developing this capability is discussed and includes the limited capability to seed faults into the actual target system due to the risk of potential damage to high value hardware. The capability envisioned would accurately reproduce the propagation of a fault or failure as observed by sensors located at strategic locations on and around the target system and would also accurately reproduce the control system and vehicle response. In this way, HMS operation can be exercised over a broad range of conditions to verify that it meets requirements for accurate, timely response to actual faults with adequate margin against false and missed detections. An overview is also presented of a real-time rocket propulsion health management system laboratory which is available for future rocket engine programs. The health management elements and approaches of this lab are directly applicable for future space systems. In this paper the various components are discussed and the general fault detection, diagnosis, isolation and the response (FDIR) concept is presented. Additionally, the complexities of V&V (Verification and Validation) for advanced algorithms and the simulation capabilities required to meet the changing state-of-the-art in HMS are discussed.

Aguilar, Robet↗
NASA NTRS2010

MISSE 2 PEACE Polymers Experiment Atomic Oxygen Erosion Yield Error Analysis

Atomic oxygen erosion of polymers in low Earth orbit (LEO) poses a serious threat to spacecraft performance and durability. To address this, 40 different polymer samples and a sample of pyrolytic graphite, collectively called the PEACE (Polymer Erosion and Contamination Experiment) Polymers, were exposed to the LEO space environment on the exterior of the International Space Station (ISS) for nearly 4 years as part of the Materials International Space Station Experiment 1 & 2 (MISSE 1 & 2). The purpose of the PEACE Polymers experiment was to obtain accurate mass loss measurements in space to combine with ground measurements in order to accurately calculate the atomic oxygen erosion yields of a wide variety of polymeric materials exposed to the LEO space environment for a long period of time. Error calculations were performed in order to determine the accuracy of the mass measurements and therefore of the erosion yield values. The standard deviation, or error, of each factor was incorporated into the fractional uncertainty of the erosion yield for each of three different situations, depending on the post-flight weighing procedure. The resulting error calculations showed the erosion yield values to be very accurate, with an average error of 3.30 percent.

McCarthy, Catherine E.↗
NASA NTRS2013

Refinement of the Robert-Bonamy Formalism: Considering Effects from the Line Coupling

Since it was developed in 1979, the Robert-Bonamy (RB) formalism has been widely used in calculating pressure broadened half-widths and induced shifts for many molecular systems. However, this formalism contains several approximations whose applicability has not been thoroughly justified. One of them is that lines of interest are well isolated. When these authors developed the formalism, they have relied on this assumption twice. First, in calculating the spectral density F(ω), they have only considered the diagonal matrix elements of the relaxation operator. Due to this simplification, effects from the line mixing are ignored. Second, when they applied the linked cluster theorem to remove the cutoff, they have assumed the matrix elements of the operator exp(-iS(sub 1) - S(sub 2)) can be replaced by the exponential of the matrix elements of -iS(sub 1) - S(sub 2). With this replacement, effects from the line coupling are also ignored. Although both these two simplifications relied on the same approximation, their validity criteria are completely different and the latter is more stringent than the former. As a result, in many cases where the line mixing becomes negligible, significant effects from the line coupling have been completely missed. In the present study, we have developed a new method to evaluate the matrix elements of exp(-iS(sub 1) - S(sub 2)) and have refined the RB formalism such that line coupling can be taken into account. Our numerical calculations of the half-widths for Raman Q lines of the N(sub 2)-N(sub 2) pair have demonstrated that effects from the line coupling are important. In comparison with values derived from the RB formalism, new calculated values for these lines are significantly reduced. A recent study has shown that in comparison with the measurements and the most accurate close coupling calculations, the RB formalism overestimates the half-widths by a large amount. As a result, the refinement of the RB formalism goes in the right direction and these new calculated half-widths become closer to the "true" values.

spectra↗
DOE OSTI2025

Solar cities: A case study analysis of city-level enablers of expanded solar energy access

Rooftop solar photovoltaic (PV) adoption can benefit households by reducing electricity bills and enhancing energy resiliency. Low and moderate-income (LMI) households have been less likely to adopt PV and experience these benefits in the United States than higher-income households. Adopter income trends are often explored through quantitative analysis with limited explanatory power. Our quantitative analysis only explains around one-third of city-level variation in LMI adoption trends through socioeconomic factors such as median home values and income inequality and PV market factors such as cumulative adoption and incentives. We implement semi-structured interviews in three case studies of cities with relatively high rates of LMI PV adoption to better understand the factors that explain PV adopter income trends. The case studies partly reiterate findings from quantitative analysis, such as the role of PV incentives. The case studies reveal a broader set of LMI adoption drivers that are missed in quantitative analyses. The case studies show how city contexts can affect LMI adoption, such as the role of supportive city governments. The case studies also reveal the importance of partnerships, such as partnerships between city governments and state LMI PV program implementers. Finally, interviewees emphasized the importance of building trust among prospective LMI PV adopters. Interviewees suggested that partnerships, outreach, and consumer protection measures were crucial to building trust in PV installers among LMI households.

Adoption↗
NASA NTRS1993

The case of missing CO-13 in mergers

We present a comparison of the CO-13 and CO-12 emissions of six systems of merging galaxies: NGC 828, NGC 3256, NGC 4194, NGC 6240, Arp 220, and Arp 299. The observations were made in both J=1-0 and J=2-1 transitions with the IRAM 30 m and SEST 15 m telescopes. In all galaxies but NGC 828, the CO-13 is much weaker than in spiral galaxies. The average emissivity ratios measured at the few kiloparsec scale are: CO-12(1-0)/CO-13(1-0) approx. equals 30, CO-12(2-1)/CO-13(2-1) approx. equals 40. These values are significantly larger than those usually measured in normal spiral galaxies, which are always between 5 and 15 for the J=1-0 line. We show that such a peculiar behavior cannot be interpreted as due to the dominant presence of diffuse gas and it cannot be attributed to optically thin CO emission either. The faint CO-13 emission of mergers must result from either an underabundance of CO-13 or an overabundance of CO-12. They may be accounted for by different mechanisms: CO-13 molecules are more easily photodissociated than CO-12 ones - however our observations seem to rule out physical conditions characteristic of photodissociation regions; in the course of the merging, large amounts of unprocessed gas, with a high (CO-12/CO-13) abundance ratio, are driven from the external regions of the progenitor galaxies to the new nucleus - this could increase the (CO-12/CO-13) ratio by a factor of 2; the interstellar medium can also be quickly enriched in C-12 (thus in CO-12) by a factor of 2, due to selective nucleosynthesis of this isotope (vs. C-13) in the massive stars born during the starburst. The last two processes can thus significantly contribute to the weakness of the CO-13 lines. The high CO-13-to-CO-12 line ratios that we have measured are due to the deep transformations that take place in the interstellar medium during the merging and the starburst: indeed, the only object with nearly normal line ratios, NGC 828, appears less perturbed and active in star formation than the other sources in the sample. In most mergers, the faintness of the CO-13 emission may indicate that the conversion factor from CO-12 emissivities to H2 column densities could differ substantially from the standard galactic value.

Casoli, Fabienne↗
NASA NTRS1990

Experimental and theoretical aerodynamic characteristics of a high-lift semispan wing model

Experimental and theoretical aerodynamic characteristics were compared for a high-lift, semispan wing configuration that incorporated a slightly modified version of the NASA Advanced Laminar Flow Control airfoil section. The experimental investigation was conducted in the Langley 14- by 22-Foot Subsonic Tunnel at chord Reynolds numbers of 2.36 and 3.33 million. A two-dimensional airfoil code and a three-dimensional panel code were used to obtain aerodynamic predictions. Two-dimensional data were corrected for three-dimensional effects. Comparisons between predicted and measured values were made for the cruise configuration and for various high-lift configurations. Both codes predicted lift and pitching moment coefficients that agreed well with experiment for the cruise configuration. These parameters were overpredicted for all high-lift configurations. Drag coefficient was underpredicted for all cases. Corrected two-dimensional pressure distributions typically agreed well with experiment, while the panel code overpredicted the leading-edge suction peak on the wing. One important feature missing from both of these codes was a capability for separated flow analysis. The major cause of disparity between the measured data and predictions presented herein was attributed to separated flow conditions.

Applin, Zachary T.↗
NASA NTRS2004

Spatially Complete Global Surface Albedos Derived from Terra/MODIS Data

Spectral land surface albedo is an important parameter for describing the radiative properties of the Earth. Accordingly it reflects the consequences of natural and human interactions, such as anthropogenic, meteorological, and phenological effects, on global and local climatological trends. Consequently, albedos are integral parts in a variety of research areas, such as general circulation models (GCMs), energy balance studies, modeling of land use and land use change, and biophysical, oceanographic, and meteorological studies. Recent production of land surface anisotropy, diffuse bihemispherical (white-sky) albedo and direct beam directional hemispherical (black-sky) albedo from observations acquired by the MODIS instruments aboard NASA s Terra and Aqua satellite platforms have provided researchers with unprecedented spatial, spectral, and temporal information on the land surface's radiative characteristics. Cloud cover, which cutails retrievals, and the presence of ephemeral and seasonal snow limit the snow-free data to approximately half the global land surfaces on an annual equal-angle basis. This precludes the MOD43B3 albedo products from being used in some remote sensing and ground-based applications, climate models, and global change research projects. An ecosystem-dependent temporal interpolation technique is described that has been developed to fill missing or seasonally snow-covered data in the official MOD43B3 albedo product. The method imposes pixel-level and local regional ecosystem-dependent phenological behavior onto retrieved pixel temporal data in such a way as to maintain pixel-level spatial and spectral detail and integrity. The phenological curves are derived from statistics based on the MODIS MOD12Q1 IGBP land cover classification product geolocated with the MOD43B3 data. The resulting snow-free value-added products provide the scientific community with spatially and temporally complete global white- and black-sky surface albedo maps and statistics. These products are stored on 1'(approximately 10 km) and coarser resolution equal-angle grids, and are computed for the first seven MODIS wavelengths, ranging from 0.47 through 2.1 microns, and for three broadband wavelengths, 0.3-0.7,0.3-5.0 and 0.7-5.0 microns.

King, Michael D.↗
DOE OSTI2026

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
NASA NTRS2012

Detecting Falling Snow from Space

There is an increased interest in detecting and estimating the amount of falling snow reaching the Earth's surface in order to fully capture the atmospheric water cycle. An initial step toward global spaceborne falling snow algorithms includes determining the thresholds of detection for various active and passive sensor channel configurations, snow event cloud structures and microphysics, snowflake particle electromagnetic properties, and surface types. In this work, cloud resolving model simulations of a lake effect and synoptic snow event were used to determine the minimum amount of snow (threshold) that could be detected by the following instruments: the W -band radar of CloudSat, Global Precipitation Measurement (GPM) Dual-frequency Precipitation Radar (DPR) Ku and Ka band, and the GPM Microwave Imager (GMI) channels from 10 to 183 plus or minus 7 GHz. Eleven different snowflake shapes were used to compute radar reflectivities and passive brightness temperatures. Notable results include: (1) the W-Band radar has detection thresholds more than an order of magnitude lower than the future GPM sensors, (2) the cloud structure macrophysics influences the thresholds of detection for passive channels, (3) the snowflake microphysics plays a large role in the detection threshold for active and passive instruments, (4) with reasonable assumptions, "the passive 166 GHz channel has detection threshold values comparable to the GPM DPR Ku and Ka band radars with approximately 0.05 g per cubic meter detected at the surface, or an approximately 0.5-1 millimeter per hr. melted snow rate (equivalent to 0.5-2 centimeters per hr. solid fluffy snowflake rate). With detection levels of falling snow known, we can focus current and future retrieval efforts on detectable storms and concentrate advances on achievable results. We will also have an understanding of the light snowfall events missed by the sensors and not captured in the global estimates.

Gail Skofronick-Jackson↗
NASA NTRS2005

Completed Beltrami-Michell Formulation in Polar Coordinates

A set of conditions had not been formulated on the boundary of an elastic continuum since the time of Saint-Venant. This limitation prevented the formulation of a direct stress calculation method in elasticity for a continuum with a displacement boundary condition. The missed condition, referred to as the boundary compatibility condition, is now formulated in polar coordinates. The augmentation of the new condition completes the Beltrami-Michell formulation in polar coordinates. The completed formulation that includes equilibrium equations and a compatibility condition in the field as well as the traction and boundary compatibility condition is derived from the stationary condition of the variational functional of the integrated force method. The new method is illustrated by solving an example of a mixed boundary value problem for mechanical as well as thermal loads.

Patnaik, Surya N.↗