Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partitioned methods”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A parallel algorithm for multi-level logic synthesis using the transduction method

The Transduction Method has been shown to be a powerful tool in the optimization of multilevel networks. Many tools such as the SYLON synthesis system (X90), (CM89), (LM90) have been developed based on this method. A parallel implementation is presented of SYLON-XTRANS (XM89) on an eight processor Encore Multimax shared memory multiprocessor. It minimizes multilevel networks consisting of simple gates through parallel pruning, gate substitution, gate merging, generalized gate substitution, and gate input reduction. This implementation, called Parallel TRANSduction (PTRANS), also uses partitioning to break large circuits up and performs inter- and intra-partition dynamic load balancing. With this, good speedups and high processor efficiencies are achievable without sacrificing the resulting circuit quality.

Lim, Chieng-Fai↗
NASA NTRS1978

The development and application of strainrange partitioning as a tool in the treatment of high temperature metal fatigue

The effects of frequency, stress and strain hold periods, stress and strain ramping, creep-fatigue interspersion, and thermomechanical cycling were studied. Of special interest was the establishment of a set of universalized life relations which were normalized by the ductility of the material. These relations together with a set of rules that were devised for treating multiaxiality by the method, were, for example, useful in predicting torsion behavior, and in particular in establishing a new relationship between torsion and tension that had previously not been considered. They also were useful in the development of a framework for treatment of environmental effects.

Manson, S. S.↗
NASA NTRS1998

Nitrogen Species in the Post-Pinatubo Stratosphere: Model Analysis Utilizing UARS Measurements

We present an analysis of the impact of heterogeneous chemistry on the partitioning of nitrogen species measured by the Upper Atmosphere Research Satellite (UARS) instruments. The UARS measurements utilized include: N2O, HNO3 and ClONO2 (Cryogen Limb Array Etalon Spectrometer (CLAES), version 7), temperature, methane, ozone, H2O, HCI, NO and NO2 (HALogen Occultation Experiment (HALOE), version 18). The analysis is carried out for the data from January 1992 to September 1994 in the 100-1 mbar (approximately 17-47 km) altitude range and over 10 degree latitude bins from 70 deg S to 70 deg N. Temporal-spatial evolution of aerosol surface area density (SAD) is adopted according to the Stratospheric Aerosol and Gas Experiment (SAGE) II data. A diurnal steady-state photochemical box model, constrained by the temperature, ozone, H2O, CH4, aerosol SAD and columns of O2 and O3 above the point of interest, has been used as the main tool to analyze these data. Total inorganic nitrogen (NOy) is obtained by three different methods: (1) as a sum of the UARS measured NO, NO2, HNO3, and CIONO2; (2) from the N2O-NOy correlation, (3) from the CH4-NOy correlation. To validate our current understanding of stratospheric heterogeneous chemistry for post-Pinatubo conditions, the model-calculated NOx/NOy ratios and the NO, NO2, and HNO3 profiles are compared to the UARS-derived data. In general, the UARS-constrained box model captures the main features of nitrogen species partitioning in the post-Pinatubo years. However, the model underestimates the NO2 content, particularly, in the 30-7 mbar (approximately 23-32 km) range. Comparisons of the calculated temporal behavior of the partial columns of NO2 and HNO3 and ground based measurements at 45 deg S and 45 deg N are also presented. Our analysis indicates that ground-based and HALOE v.18 measurements of the NO2 vertical columns are consistent within the range of their uncertainties and are systematically higher (up to 50%) than the model results at mid-latitudes in both hemispheres. Reasonable agreement is obtained for HNO3 columns at 45 deg S suggesting some problems with nitrogen species partitioning in the model. Outstanding uncertainties are discussed.

Danilin, M. Y.↗
NASA NTRS1994

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods for the analysis of complete jet engines. We have entitled this program by applying the two dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition, and solution capabilities were successfully tested. We then focused attention on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion that results from these structural displacements. This is treated by a new arbitrary Lagrangian-Eulerian (ALE) technique that models the fluid mesh motion as that of a fictitious mass-spring network. New partitioned analysis procedures to treat this coupled three-component problem are developed. These procedures involved delayed corrections and subcycling. Preliminary results on the stability, accuracy, and MPP computational efficiency are reported.

Felippa, C. A.↗
NASA NTRS2006

Strict Constraint Feasibility in Analysis and Design of Uncertain Systems

This paper proposes a methodology for the analysis and design optimization of models subject to parametric uncertainty, where hard inequality constraints are present. Hard constraints are those that must be satisfied for all parameter realizations prescribed by the uncertainty model. Emphasis is given to uncertainty models prescribed by norm-bounded perturbations from a nominal parameter value, i.e., hyper-spheres, and by sets of independently bounded uncertain variables, i.e., hyper-rectangles. These models make it possible to consider sets of parameters having comparable as well as dissimilar levels of uncertainty. Two alternative formulations for hyper-rectangular sets are proposed, one based on a transformation of variables and another based on an infinity norm approach. The suite of tools developed enable us to determine if the satisfaction of hard constraints is feasible by identifying critical combinations of uncertain parameters. Since this practice is performed without sampling or partitioning the parameter space, the resulting assessments of robustness are analytically verifiable. Strategies that enable the comparison of the robustness of competing design alternatives, the approximation of the robust design space, and the systematic search for designs with improved robustness characteristics are also proposed. Since the problem formulation is generic and the solution methods only require standard optimization algorithms for their implementation, the tools developed are applicable to a broad range of problems in several disciplines.

Crespo, Luis G.↗
NASA NTRS1991

Large constraint length high speed viterbi decoder based on a modular hierarchial decomposition of the deBruijn graph

A method of formulating and packaging decision-making elements into a long constraint length Viterbi decoder which involves formulating the decision-making processors as individual Viterbi butterfly processors that are interconnected in a deBruijn graph configuration. A fully distributed architecture, which achieves high decoding speeds, is made feasible by novel wiring and partitioning of the state diagram. This partitioning defines universal modules, which can be used to build any size decoder, such that a large number of wires is contained inside each module, and a small number of wires is needed to connect modules. The total system is modular and hierarchical, and it implements a large proportion of the required wiring internally within modules and may include some external wiring to fully complete the deBruijn graph. pg,14.

Collins, Oliver↗
NASA NTRS2008

Experimental Constraints on the Partitioning and Valence of V and Cr in Garnet and Coexisting Glass

A series of experiments with garnet and coexisting melt have been carried out across a range of oxygen fugacities (near hematite-magnetite (HM) to below the iron-wustite (IW) buffers) at 1.7 GPa to study the partitioning and valence of Cr and V in both phases. Experiments were carried out in a non end loaded piston cylinder apparatus, and the run products were analyzed with electron microprobe and xray absorption near edge structure (XANES) analysis at beamline 13-ID at the Advanced Photon Source of Argonne National Lab. The valence of vanadium and chromium were determined using the position and intensity of the Ka pre-edge peaks, calibrated on a series of Cr and Vbearing standard glasses. This technique has been applied to V and Cr in glasses and V in spinels previously, and in these isotropic phases there are no orientational effects on the XANES spectra (Righter et al., 2006, Amer. Mineral. 91, 1643-1656). We also now demonstrate this to be true for V and Cr in garnet. Also, previous work has shown that V has a higher valence in the glass (or melt) than in the coexisting spinel. This is also true for V in garnet-glass pairs in this study. Vanadium valence in garnets varies from 2.7 below the IW buffer to 3.7 near HM, and for coexisting glass it varies from 3.2 to 4.3. Vanadium valence measured in some natural garnets from mantle localities indicates V in the more reduced range at 2.5. Comparisons will be made between fO2 estimated from V valence and other methods for garnet-bearing mantle samples. In contrast, Cr valence measured in garnet and coexisting glass for all experimental and natural samples is 2.9- 3.0, suggesting that the valence of Cr does not vary within either phase across a large fO2 range. These results demonstrate that while V varies from 2+ to 3+ to 4+ in garnet-melt systems, Cr does not, and this will ultimately affect the partitioning behavior of these two elements in natural systems. Garnet/melt D(Cr) are between 12 and 17 across this range of fO2, whereas D(V) has the highest partition coefficient approx.3, near the IW buffer where the valence of V is almost entirely 3+.

Righter, K.↗
NASA NTRS1998

Nitrogen Species in the Post-Pinatubo Stratosphere: Model Analysis Utilizing UARS Measurements

We present an analysis of the impact of heterogeneous chemistry on the partitioning of nitrogen species measured by the Upper Atmosphere Research Satellite (UARS) instruments. The UARS measurements utilized include: N2O, HNO3 and ClONO2 (Cryogen Limb Array Etalon Spectrometer (CLAES), version 7), temperature, methane, ozone, H2O, HCl, NO and NO2 (HALogen Occultation Experiment (HALOE), version 18). The analysis is carried out for the data from January 1992 to September 1994 in the 100-1 mbar (approx.17-47 km) altitude range and over 10 degree latitude bins from 70degS to 70degN. Temporal-spatial evolution of aerosol surface area density (SAD) is adopted according to the Stratospheric Aerosol and Gas Experiment (SAGE) 11 data. A diurnal steady-state photochemical box model, constrained by the temperature, ozone, H2O, CH4, aerosol SAD and columns of O2 and O3 above the point of interest, has been used as the main tool to analyze these data. Total inorganic nitrogen (NO(y)) is obtained by three different methods: (1) as a sum of the UARS measured NO, NO2, HNO3, and ClONO2; (2) from the N2O-NO(y) correlation, and (3) from the CH4-NO(y) correlation. To validate our current understanding of stratospheric heterogeneous chemistry for post-Pinatubo conditions, the model-calculated NO(x)/NO(y) ratios and the NO, NO2, and HNO3 profiles are compared to the UARS-derived data. In general, the UARS-constrained box model captures the main features of nitrogen species partitioning in the post-Pinatubo years. However, the model underestimates the NO2 content, particularly, in the 30-7 mbar (approx. 23-32 km) range. Comparisons of the calculated temporal behavior of the partial columns of NO2 and HNO3 and ground based measurements at 45degS and 45degN are also presented. Our analysis indicates that ground-based and HALOE v. 18 measurements of the NO2 vertical columns are consistent within the range of their uncertainties and are systematically higher (up to 50%) than the model results at mid-latitudes in both hemispheres. Reasonable agreement is obtained for HNO3 columns at 45degS suggesting some problems with nitrogen species partitioning in the model. Outstanding uncertainties are discussed.

Danilin, M. Y.↗
NASA NTRS1994

Are models of catalytic removal of O3 by HO(x) accurate? Constraints from in situ measurements of the OH to HO2 ratio

Measurements of the ratio OH/HO2, NO, O3, ClO, and BrO were obtained at altitudes from 15-20 km and latitudes from 15-60 deg N. A method is presented for interpreting the rates of chemical transformations that (1) are responsible for over half the ozone removal rate in the lower stratosphere via reactions of HO2; and (2) control the abundance of HO2 through coupling to nitrogen and halogen radicals. The results show our understanding of the chemical reactions controlling the partitioning of OH and HO2 is complete and accurate and that the potential effects of 'missing chemistry' are strickly constrained in the region of the atmosphere encompassed by the observations. The analysis demonstrates that the sensitivity of the ratio OH/HO2 to changes in NO is described to within 12% by current models. This reduces by more than a factor of 2 the effect of uncertainty in the coupling of hydrogen and nitrogen radicals on the analysis of the potential effects of perturbations to odd notrogen in the lower statosphere.

Cohen, R. C.↗
NASA NTRS2022

Atmospheric Energy Spectra in Global Kilometre-Scale Models

Eleven 40-day long integrations of five different global models with horizontal resolutions of less than 9 km are compared in terms of their global energy spectra. The method of normal-mode function decomposition is used to distinguish between balanced (Rossby wave; RW) and unbalanced (inertia-gravity wave; IGW) circulation. The simulations produce the expected canonical shape of the spectra, but their spectral slopes at mesoscales, and the zonal scale at which RW and IGW spectra intersect differ significantly. The partitioning of total wave energies into RWs an IGWs is most sensitive to the turbulence closure scheme and this partitioning is what determines the spectral crossing scale in the simulations, which differs by a factor of up to two. It implies that care must be taken when using simple spatial filtering to compare gravity wave phenomena in storm-resolving simulations, even when the model horizontal resolutions are similar. In contrast to the energy partitioning between the RWs and IGWs, changes in turbulence closure schemes do not seem to strongly affect spectral slopes, which only exhibit major differences at mesoscales. Despite their minor contribution to the global (horizontal kinetic plus potential available) energy, small scales are important for driving the global mean circulation. Our results support the conclusions of previous studies that the strength of convection is a relevant factor for explaining discrepancies in the energies at small scales. The models studied here produce the major large-scale features of tropical precipitation patterns. However, particularly at large horizontal wavenumbers, the spectra of upper tropospheric vertical velocity, which is a good indicator for the strength of deep convection, differ by factors of three or more in energy. High vertical kinetic energies at small scales are mostly found in those models that do not use any convective parameterisation.

DYAMOND↗
NASA NTRS1999

Nitrogen Species in the Post-Pinatubo Stratosphere: Model Analysis Utilizing UARS Measurements

We present an analysis of the impact of heterogeneous chemistry on the partitioning of nitrogen species measured by the Upper Atmosphere Research Satellite (UARS) instruments. The UARS measurements utilized include N2O, HNO3, and ClONO2 from the cryogenic limb array etalon spectrometer (CLAES), version 7 (v.7), and temperature, methane, ozone, H2O, HCl, NO and NO2 from the halogen occultation experiment (HALOE), version 18. The analysis is carried out for the UARS data obtained between January 1992 and September 1994 in the 100-to 1-mbar (approx. 17-47 km) altitude range and over 10 degrees latitude bins from 70 S to 70 N. The spatiotemporal evolution of aerosol surface area density (SAD) is adopted from analysis of the Stratospheric Aerosol and Gas Experiment (SAGE) II data. A diurnal steady state photochemical box model, constrained by the temperature, ozone, H2O, CH4, aerosol SAD, and columns of O2 and O3 above the point of interest, has been used as the main tool to analyze these data. Total inorganic nitrogen (NOy) is obtained by three different methods: (1) as a sum of the UARS-measured NO, NO2, HNO3, and ClONO2; (2) from the N2O-NOy correlation, and (3) from the CH4-NOy correlation. To validate our current understanding of stratospheric heterogeneous chemistry for post-Pinatubo conditions, the model-calculated monthly averaged NOx/NOy ratios and the NO, NO2, and HNO3 profiles are compared with the UARS-derived data. In general, the UARS-constrained box model captures the main features of nitrogen species partitioning in the post-Pinatubo years, such as recovery of NOx after the eruption, their seasonal variability and vertical profiles. However, the model underestimates the NO2 content, particularly in the 30- to 7-mbar (approx.23-32 km) range. Comparisons of the calculated temporal behavior of the partial columns of NO2 and HNO3 and ground-based measurements at 45 S and 45 N are also presented. Our analysis indicates that ground-based and HALOE v.18 measurements of the NO2 vertical columns are consistent within the range of their uncertainties and are systematically higher (up to 50%) than the model results at midlatitudes in both hemispheres. Reasonable agreement is obtained for HNO3 columns at 45 S, suggesting some problems with nitrogen species partitioning in the model. Outstanding uncertainties are discussed.

Danilin, Michael Y.↗
NASA NTRS1989

Surface energy and radiation balance systems - General description and improvements

Surface evaluation of sensible and latent heat flux densities and the components of the radiation balance were desired for various vegetative surfaces during the ASCOT84 experiment to compare with modeled results and to relate these values to drainage winds. Five battery operated data systems equipped with sensors to determine the above values were operated for 105 station days during the ASCOT84 experiment. The Bowen ratio energy balance technique was used to partition the available energy into the sensible and latent heat flux densities. A description of the sensors and battery operated equipment used to collect and process the data is presented. In addition, improvements and modifications made since the 1984 experiment are given. Details of calculations of soil heat flow at the surface and an alternate method to calculate sensible and latent heat flux densities are provided.

Fritschen, Leo J.↗
NASA NTRS1996

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗
NASA NTRS1997

Progress in Parallel Schur Complement Preconditioning for Computational Fluid Dynamics

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a non-overlapping Schur complement procedure for arbitrary triangulated domains. The ultimate goal of this research is to develop scalable preconditioning algorithms for fluid flow discretizations on parallel computing architectures. In our implementation of the Schur complement preconditioning technique, the triangulation is first partitioned into a number of subdomains using the METIS multi-level k-way partitioning code. This partitioning induces a natural 2X2 partitioning of the p.d.e. discretization matrix. By considering various inverse approximations of the 2X2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 and the SGI Power Challenge array using MPI message passing protocol. A number of example CFD calculations will be presented to illustrate and assess various Schur complement approximations.

Barth, Timothy J.↗
NASA NTRS2021

Spectral Derivatives of Optical Depth for Partitioning Aerosol Type and Loading

Quantifying aerosol compositions (e.g., type, loading) from remotely sensed measurements by spaceborne, suborbital and ground-based platforms is a challenging task. In this study, the first and second-order spectral derivatives of aerosol optical depth (AOD) with respect to wavelength are explored to determine the partitions of the major components of aerosols based on the spectral dependence of their particle optical size and complex refractive index. With theoretical simulations from the Second Simulation of a Satellite Signal in the Solar Spectrum (6S) model, AOD spectral derivatives are characterized for collective models of aerosol types, such as mineral dust (DS) particles, biomass-burning (BB) aerosols and anthropogenic pollutants (AP), as well as stretching out to the mixtures among them. Based on the intrinsic values from normalized spectral derivatives, referenced as the Normalized Derivative Aerosol Index (NDAI), a unique pattern is clearly exhibited for bounding the major aerosol components; in turn, fractions of the total AOD (fAOD) for major aerosol components can be extracted. The subtlety of this NDAI method is examined by using measurements of typical aerosol cases identified carefully by the ground-based Aerosol Robotic Network (AERONET) sun–sky spectroradiometer. The results may be highly practicable for quantifying fAOD among mixed-type aerosols by means of the normalized AOD spectral derivatives.

aerosol partition↗
NASA NTRS1999

Nitrogen Species in the Post-Pinatubo Stratosphere: Model Analysis Utilizing UARS Measurements

We present an analysis of the impact of heterogeneous chemistry on the partitioning of nitrogen species measured by the Upper Atmosphere Research Satellite (UARS) instruments. The UARS measurements utilized include N2O, HNO3, and ClONO2 from the cryogenic limb array etalon spectrometer (CLAES), version 7 (v.7), and temperature, methane, ozone, H2O, HCl, NO and NO2 from the halogen occultation experiment (HALOE), version 18. The analysis is carried out for the UARS data obtained between January 1992 and September 1994 in the 100- to 1-mbar (approx. 17-47 km) altitude range and over 10 deg latitude bins from 70 deg S to 70 deg N. The spatiotemporal evolution of aerosol surface area density (SAD) is adopted from analysis of the Stratospheric Aerosol and Gas Experiment (SAGE) II data. A diurnal steady state photochemical box model, constrained by the temperature, ozone, H2O, CH4, aerosol SAD, and columns of O2 and O3 above the point of interest, has been used as the main tool to analyze these data. Total inorganic nitrogen (NOY) is obtained by three different methods: (1) as a sum of the UARS-measured NO, NO2, HNO3, and ClONO2; (2) from the N2O-NOY correlation; and (3) from the CH4-NOY correlation. To validate our current understanding of stratospheric heterogeneous chemistry for post-Pinatubo conditions, the model-calculated monthly averaged NO(x)/NO(y) ratios and the NO, NO2, and HNO3 profiles are compared with the UARS-derived data. In general, the UARS-constrained box model captures the main features of nitrogen species partitioning in the post-Pinatubo years, such as recovery of NO(x) after the eruption, their seasonal variability and vertical profiles. However, the model underestimates the NO2 content, particularly in the 30- to 7-mbar (approx. 23-32 km) range. Comparisons of the calculated temporal behavior of the partial columns of NO2 and HNO3 and ground-based measurements at 45 deg S and 45 deg N are also presented. Our analysis indicates that ground-based and HALOE v.18 measurements of the NO2 vertical columns are consistent within the range of their uncertainties and are systematically higher (up to 50%) than the model results at midlatitudes in both hemispheres. Reasonable agreement is obtained for HNO3 columns at 45 deg S, suggesting some problems with nitrogen species partitioning in the model. Outstanding uncertainties are discussed.

Danilin, Michael Y.↗
NASA NTRS1983

Determination of aerodynamic parameters of a fighter airplane from flight data at high angles of attack

A procedure for the estimation of airplane model structure and parameters is applied to data from a modern fighter airplane operating within an angle of attack range of 5 to 60 deg. The paper briefly describes the airplane, flight and wind tunnel data available, and the estimation method. The results presented contain basic longitudinal characteristics of the airplane and the estimates of aerodynamic parameters in the yawing-moment equations. These estimates are obtained from small and large amplitude maneuvers. Because the latter set of data was not suitable for airplane identification, some of the large amplitude maneuvers were joined together and then partitioned into subsets according to the values of angle of attack. Each subset was then analyzed as a separate data set. Most of the estimated parameters and functions are in good agreement with the wind tunnel measurements. The estimated lateral parameters in the model equations also demonstrate good prediction capabilities.

Klein, V.↗
NASA NTRS1995

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a by-pass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by an ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled 3-component problem were developed in 1994. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 has been developed. It is planned to use the steady-state global solution provided by ENG10 as input to a localized three-dimensional FSI analysis for engine regions where aeroelastic effects may be important.

Felippa, C. A.↗