SEARCH · Engineering Papers

Results for “Weight function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Optimal rocket thrust profile shaping using third degree spline function interpolation

Optimal solid-rocket thrust profiles for the parallel-burn, solid-rocket-assisted space shuttle are investigated. Solid-rocket thrust profiles are simulated by using third-degree spline functions, with the values of the thrust ordinates defined as parameters. The profiles are optimized parametrically, using the Davidon-Fletcher-Powell penalty function method, by minimizing propellant weight subject to state and control inequality constraints and to terminal boundary conditions. This study shows that optimizing a control variable parametrically by using third-degree spline function interpolation allows the control to be shaped so that inequality constraints are strictly adhered to and all corners are eliminated. The absence of corners, which is realistic in nature, makes this method attractive from the viewpoint of solid rocket grain design.

Johnson, I. L.↗

NASA NTRS1995

Pros and Cons of Using Water Immersion to Simulate Physiological Responses to Microgravity

Head-out water immersion (HOI) has been employed as a remedial treatment for various ills and ailments for many millennia, and total body immersion even longer as protective encapsulation for the mammalian fetus. Two discrete differences between stimuli induced by true microgravity (10(exp -4) g) and HOI are readily apparent. External water pressure on the skin and accompanying negative pressure breathing cause blood to shift headward. Secondly, the gravitational force is ever present during immersion and microgravity, but its effect is essentially neutralized during Earth orbital flight. Thus, the physiological responses to immersion should not be expected to match those during microgravity. Immersion has been used mainly to study and understand kidney function and associated cardiovascular responses for control of body fluid volume and osmotic content, with some application to and simulation of microgravity responses. There is a plethora of data from human HOI studies, but relatively few controlled data from microgravity studies. In general, it appears that physiological responses occur more quickly with water immersion than in microgravity, but this may be due to less rigorous control (voluntary and involuntary) of the preflight state of crew members. The central venous pressure-vasopressin (Gauer-Henry) reflex control for fluid balance may not be of prime importance in microgravity. Gross functions such as reduced body weight and water, level of hypovolemia, decreased isokinetic strength, and lower nitrogen balance found during immersion are qualitatively similar in microgravity, but the mechanisms controlling these and other functions are, for the most part, unclear. Only acquisition of data from well-controlled microgravity experiments will resolve this discrepancy.

Greenleaf, J. E.↗

NASA NTRS2015

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗

DOE OSTI2022

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗

NASA NTRS1917

Thermodynamic efficiency of present types of internal combustion engines for aircraft

Report presents requirements of internal combustion engines suitable for aircraft. Topics include: (1) service requirements for aeronautic engines - power versus weight, reliability, and adaptability factors, (2) general characteristics of present aero engines, (3) aero engine processes and functions of parts versus power-weight ratio, reliability, and adaptability factors, and (4) general arrangement, form, proportions, and materials of aero parts - power-weight ratio, reliability, and adaptability.

Lucke, Charles E↗

DOE OSTI2020

Spectroscopic signatures of next-nearest-neighbor hopping in the charge and spin dynamics of doped one-dimensional antiferromagnets

In this work, we study the impact of next-nearest-neighbor (nnn) hopping on the low-energy collective excitations of strongly correlated doped antiferromagnetic cuprate spin chains. Specifically, we use exact diagonalization and the density matrix renormalization group method to study the single-particle spectral function, the dynamical spin and charge structure factors, and the Cu L-edge resonant inelastic x-ray scattering (RIXS) intensity of the doped t-t'-J model for a set of t' values. We find evidence for the breakdown of spin-charge separation as |t'| increases and identify its fingerprints in the dynamical response functions. The inclusion of nnn hopping couples the spinon and holon excitations, resulting in the formation of a spin-polaron, where a ferromagnetic spin polarization cloud dresses the doped carrier. The spin-polaron manifests itself as additional spectral weight in the dynamical correlation functions, which appear simultaneously in the spin- and charge-sensitive channels. We also demonstrate that RIXS can provide a unique view of the spin-polaron, due to its sensitivity to both the spin and charge degrees of freedom.

1D spin chains↗

DOE OSTI2020

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1998

Technical Challenges for Next Generation Liquid Rocket Engine Nozzles

Next generation launch and propulsion systems face significant challenges of providing increased performance at lower cost with shorter development cycles. The use of existing materials in new application areas and the development and application of "advanced" material systems are enabling for the achievement of these aggressive goals. An area which offers substantial opportunity to decrease liquid rocket engine systems weight, and thus provide a measure of increased performance, is the nozzle. This presentation will examine the technical issues and showstoppers limiting the application of alternate materials to liquid rocket engine nozzles and identify the key material systems which have the potential for high payoff relative to the forcing functions of cost and performance (weight). Existing nozzle material development projects will be cited and opportunities for future joint partnerships to address these materials and processes challenges will be addressed. The presentation will include overviews of activities in low cost ablative nozzles, ceramic matrix composite cooled and uncooled nozzles, and advanced metallics.

Clinton, R. G., Jr.↗

DOE OSTI2025

Filling the Gaps: A Bayesian Mixture Model for Imputing Missing Soil Water Content Data

ABSTRACT Soil water content (SWC) data are central to evaluating how soil moisture varies over time and space and influences critical plant and ecosystem functions, especially in water‐limited drylands. However, sensors that record SWC at high frequencies often malfunction, leading to incomplete timeseries and limiting our understanding of dryland ecosystem dynamics. We developed an analytical approach to impute missing SWC data, which we tested at six eddy flux tower sites along an elevation gradient in the southwestern United States. We impute missing data as a mixture of linearly interpolated SWC between the observed endpoints of a missing data gap and SWC simulated by an ecosystem water balance model (SOILWAT2). Within a Bayesian framework, we allowed the relative utility (mixture weight) of each component (linearly interpolated vs. SOILWAT2) to vary by depth, site and gap characteristics. We explored “fixed” weights versus “dynamic” weights that vary as a function of cumulative precipitation, average temperature, and time since the start of the gap. Both models estimated missing SWC data well ( R 2 = 0.70–0.88 vs. 0.75–0.91 for fixed vs. dynamic weights, respectively), but the utility of linearly interpolated versus SOILWAT2 values depended on site and depth. SOILWAT2 was more useful for more arid sites, shallower depths, longer and warmer gaps and gaps that received greater precipitation. Overall, the mixture model reliably gap‐fills SWC, while lending insight into processes governing SWC dynamics. This approach to impute missing data could be adapted to accommodate more than two mixture components and other types of environmental timeseries.

Ogle, Kiona [School of Informatics, Computing, and↗

NASA NTRS1985

Recovering crater production functions from densely cratered surfaces

Two techniques proposed for the recovery of the crater-production functions on densely cratered surfaces have been appraised. The first technique classifies craters according to their degree of rim preservation, then interprets the statistics of the freshest craters as closely mimicking the full population's production function. Although commonly practiced, in fact, this method is demonstrably misleading for densely cratered surfaces. The second technique uses a weighting scheme to recover the production function. Each crater diameter measured is weighed by the fraction of the crater that remains unoverlapped by subsequent impacts. While this method is theoretically sound even at saturation crater densities, it too fails at high densities due to practical difficulties. For this technique to work effectively, one would have to be able to recognize and accurately measure the diameters of craters having less than 10 percent of their rims remaining, a feat beyond reasonable expectation.

Woronow, A.↗

DOE OSTI2023

Holographic open quantum systems: toy models and analytic properties of thermal correlators

We present a unified picture of open quantum systems, the theory of a system probing a noisy thermal environment, distilling lessons learnt from previous holographic analyses. Our treatment is applicable both when the system is coupled to short-lived (Markovian), and long-lived (non-Markovian) environmental degrees of freedom. The thermal environment is modeled using an asymptotically AdS black hole, and the systems of interest are simple probe field theories. The effective stochastic dynamics of the system is governed by real-time thermal correlators, which we compute using the gravitational Schwinger-Keldysh (grSK) geometry. We describe the structure of arbitrary tree-level contact and exchange Witten diagrams in the grSK geometry. In particular, we argue, that all such diagrams reduce to integrals supported on a single copy of the exterior of the black hole. The integrand is obtained as a multiple discontinuity of a function comprising ingoing boundary-bulk propagators, monodromy functions which appear as radial Boltzmann weights, and vertex factors. These results allow us to deduce the analytic structure of real-time thermal n-point functions in holographic CFTs. We illustrate the general statements by a two-dimensional toy model, dual to fields in the BTZ background, which we argue captures many of the essential features of generic open holographic QFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1991

Frequency weighted system identification and linear quadratic controller design

Application of filters for frequency weighting of Markov parameters (pulse response functions) is described in relation to system/observer identification. The time domain identification approach recovers a model which has a pulse response weighted according to frequency. The identified model is composed of the original system and filters. The augmented system is in a form which can be used directly for frequency weighted linear quadratic controller design. Data from either single or multiple experiments can be used to recover the Markov parameters. Measured acceleration signals from a truss structure are used for system identification and the model obtained is used for frequency weighted controller design. The procedure makes the identification and controler design complementary problems.

Horta, Lucas G.↗

NASA NTRS1985

Space Station man-machine automation trade-off analysis

The man machine automation tradeoff methodology presented is of four research tasks comprising the autonomous spacecraft system technology (ASST) project. ASST was established to identify and study system level design problems for autonomous spacecraft. Using the Space Station as an example spacecraft system requiring a certain level of autonomous control, a system level, man machine automation tradeoff methodology is presented that: (1) optimizes man machine mixes for different ground and on orbit crew functions subject to cost, safety, weight, power, and reliability constraints, and (2) plots the best incorporation plan for new, emerging technologies by weighing cost, relative availability, reliability, safety, importance to out year missions, and ease of retrofit. A fairly straightforward approach is taken by the methodology to valuing human productivity, it is still sensitive to the important subtleties associated with designing a well integrated, man machine system. These subtleties include considerations such as crew preference to retain certain spacecraft control functions; or valuing human integration/decision capabilities over equivalent hardware/software where appropriate.

Zimmerman, W. F.↗

NASA NTRS2004

Real-Time Adaptive Color Segmentation by Neural Networks

Artificial neural networks that would utilize the cascade error projection (CEP) algorithm have been proposed as means of autonomous, real-time, adaptive color segmentation of images that change with time. In the original intended application, such a neural network would be used to analyze digitized color video images of terrain on a remote planet as viewed from an uninhabited spacecraft approaching the planet. During descent toward the surface of the planet, information on the segmentation of the images into differently colored areas would be updated adaptively in real time to capture changes in contrast, brightness, and resolution, all in an effort to identify a safe and scientifically productive landing site and provide control feedback to steer the spacecraft toward that site. Potential terrestrial applications include monitoring images of crops to detect insect invasions and monitoring of buildings and other facilities to detect intruders. The CEP algorithm is reliable and is well suited to implementation in very-large-scale integrated (VLSI) circuitry. It was chosen over other neural-network learning algorithms because it is better suited to realtime learning: It provides a self-evolving neural-network structure, requires fewer iterations to converge and is more tolerant to low resolution (that is, fewer bits) in the quantization of neural-network synaptic weights. Consequently, a CEP neural network learns relatively quickly, and the circuitry needed to implement it is relatively simple. Like other neural networks, a CEP neural network includes an input layer, hidden units, and output units (see figure). As in other neural networks, a CEP network is presented with a succession of input training patterns, giving rise to a set of outputs that are compared with the desired outputs. Also as in other neural networks, the synaptic weights are updated iteratively in an effort to bring the outputs closer to target values. A distinctive feature of the CEP neural network and algorithm is that each update of synaptic weights takes place in conjunction with the addition of another hidden unit, which then remains in place as still other hidden units are added on subsequent iterations. For a given training pattern, the synaptic weight between (1) the inputs and the previously added hidden units and (2) the newly added hidden unit is updated by an amount proportional to the partial derivative of a quadratic error function with respect to the synaptic weight. The synaptic weight between the newly added hidden unit and each output unit is given by a more complex function that involves the errors between the outputs and their target values, the transfer functions (hyperbolic tangents) of the neural units, and the derivatives of the transfer functions.

Duong, Tuan A.↗

DOE OSTI2021

Solution-processable and functionalizable ultra-high molecular weight polymers via topochemical synthesis

Abstract Topochemical polymerization reactions hold the promise of producing ultra-high molecular weight crystalline polymers. However, the totality of topochemical polymerization reactions has failed to produce ultra-high molecular weight polymers that are both soluble and display variable functionality, which are restrained by the crystal-packing and reactivity requirements on their respective monomers in the solid state. Herein, we demonstrate the topochemical polymerization reaction of a family of para -azaquinodimethane compounds that undergo facile visible light and thermally initiated polymerization in the solid state, allowing for the first determination of a topochemical polymer crystal structure resolved via the cryoelectron microscopy technique of microcrystal electron diffraction. The topochemical polymerization reaction also displays excellent functional group tolerance, accommodating both solubilizing side chains and reactive groups that allow for post-polymerization functionalization. The thus-produced soluble ultra-high molecular weight polymers display superior capacitive energy storage properties. This study overcomes several synthetic and characterization challenges amongst topochemical polymerization reactions, representing a critical step toward their broader application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

On the influence of over-parameterization in manifold based surrogates and deep neural operators

Constructing accurate and generalizable approximators (surrogate models) for complex physico-chemical processes exhibiting highly non-smooth dynamics is challenging. The main question is what type of surrogate models we should construct and should these models be under-parameterized or over-parameterized. In this work, we propose new developments and perform comparisons for two promising approaches: manifold-based polynomial chaos expansion (m-PCE) and the deep neural operator (DeepONet), and we examine the effect of over-parameterization on generalization. While m-PCE enables the construction of a mapping by first identifying low-dimensional embeddings of the input functions, parameters, and quantities of interest (QoIs), a neural operator learns the nonlinear mapping via the use of deep neural networks. Here, we demonstrate the performance of these methods in terms of generalization accuracy by solving the 2D time-dependent Brusselator reaction-diffusion system with uncertainty sources, modeling an autocatalytic chemical reaction between two species. We first propose an extension of the m-PCE by constructing a mapping between latent spaces formed by two separate embeddings of the input functions and the output QoIs. To further enhance the accuracy of the DeepONet, we introduce weight self-adaptivity in the loss function. We demonstrate that the performance of m-PCE and DeepONet is comparable for cases of relatively smooth input-output mappings. However, when highly non-smooth dynamics is considered, DeepONet shows higher approximation accuracy. We also find that for m-PCE, modest over-parameterization leads to better generalization, both within and outside of distribution, whereas aggressive over-parameterization leads to over-fitting. In contrast, an even highly over-parameterized DeepONet leads to better generalization for both smooth and non-smooth dynamics. Furthermore, we compare the performance of the above models with another recently proposed operator learning model, the Fourier Neural Operator, and show that its over-parameterization also leads to better generalization. Taken together, our studies show that m-PCE can provide very good accuracy at very low training cost, whereas a highly over-parameterized DeepONet can provide better accuracy and robustness to noise but at higher training cost. In both methods, the inference cost is negligible.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2010

Analysis of Skeletal Muscle Metrics as Predictors of Functional Task Performance

PURPOSE: The ability to predict task performance using physiological performance metrics is vital to ensure that astronauts can execute their jobs safely and effectively. This investigation used a weighted suit to evaluate task performance at various ratios of strength, power, and endurance to body weight. METHODS: Twenty subjects completed muscle performance tests and functional tasks representative of those that would be required of astronauts during planetary exploration (see table for specific tests/tasks). Subjects performed functional tasks while wearing a weighted suit with additional loads ranging from 0-120% of initial body weight. Performance metrics were time to completion for all tasks except hatch opening, which consisted of total work. Task performance metrics were plotted against muscle metrics normalized to "body weight" (subject weight + external load; BW) for each trial. Fractional polynomial regression was used to model the relationship between muscle and task performance. CONCLUSION: LPMIF/BW is the best predictor of performance for predominantly lower-body tasks that are ambulatory and of short duration. LPMIF/BW is a very practical predictor of occupational task performance as it is quick and relatively safe to perform. Accordingly, bench press work best predicts hatch-opening work performance.

Ryder, Jeffrey W.↗

NASA NTRS2017

Robust Acoustic Objective Functions and Sensitivities in Adjoint-Based Design Optimizations

The multidisciplinary design of aircraft typically includes considerations of performance, weight, fuel burn, and noise, among other factors. An objective function is applied to each of these considerations in order to weight the influence of trade-offs between different designs. Higher-order optimization exercises have utilized an adjoint approach to reach an optimal set of objective functions, such as maximum lift and reduced drag. Taking advantage of the adjoint approach significantly reduces the computational time required to find an optimal configuration. Including acoustics in the set of objective functions during an adjoint-based design optimization requires the sensitivity of the acoustic objective function. This document will present an approach for defining the sensitivity of several acoustic metrics and operations that can fill the role of the acoustic objective function. This includes time-integrated metrics such as effective perceived noise level (EPNL) and frequency-integrated metrics such as overall sound pressure level (OASPL). A demonstration case, validation, and details on the implementation in the second generation Aircraft NOise Prediction Program (ANOPP2) are also shown.

Lopes, Leonard V.↗