Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel optimization”

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 757 records · Page 42

An Efficient Objective Analysis System for Parallel Computers

A new atmospheric objective analysis system designed for parallel computers will be described. The system can produce a global analysis (on a 1 X 1 lat-lon grid with 18 levels of heights and winds and 10 levels of moisture) using 120,000 observations in 17 minutes on 32 CPUs (SGI Origin 2000). No special parallel code is needed (e.g. MPI or multitasking) and the 32 CPUs do not have to be on the same platform. The system is totally portable and can run on several different architectures at once. In addition, the system can easily scale up to 100 or more CPUS. This will allow for much higher resolution and significant increases in input data. The system scales linearly as the number of observations and the number of grid points. The cost overhead in going from 1 to 32 CPUs is 18%. In addition, the analysis results are identical regardless of the number of processors used. This system has all the characteristics of optimal interpolation, combining detailed instrument and first guess error statistics to produce the best estimate of the atmospheric state. Static tests with a 2 X 2.5 resolution version of this system showed it's analysis increments are comparable to the latest NASA operational system including maintenance of mass-wind balance. Results from several months of cycling test in the Goddard EOS Data Assimilation System (GEOS DAS) show this new analysis retains the same level of agreement between the first guess and observations (O-F statistics) as the current operational system.

Stobie, J.↗

Artificial Neural Network Modeling for Airline Disruption Management

Since the 1970s, most airlines have incorporated computerized support for managing disruptions during flight schedule execution. However, existing platforms for airline disruption management (ADM) employ monolithic system design methods that rely on the creation of specific rules and requirements through explicit optimization routines, before a system that meets the specifications is designed. Thus, current platforms for ADM are unable to readily accommodate additional system complexities resulting from the introduction of new capabilities, such as the introduction of unmanned aerial systems (UAS), operations and infrastructure, to the system. To this end, we use historical data on airline scheduling and operations recovery to develop a system of artificial neural networks (ANNs), which describe a predictive transfer function model (PTFM) for promptly estimating the recovery impact of disruption resolutions at separate phases of flight schedule execution during ADM. Furthermore, we provide a modular approach for assessing and executing the PTFM by employing a parallel ensemble method to develop generative routines that amalgamate the system of ANNs. Our modular approach ensures that current industry standards for tardiness in flight schedule execution during ADM are satisfied, while accurately estimating appropriate time-based performance metrics for the separate phases of flight schedule execution.

Kolawole Ogunsina↗

A linear decomposition method for large optimization problems. Blueprint for development

A method is proposed for decomposing large optimization problems encountered in the design of engineering systems such as an aircraft into a number of smaller subproblems. The decomposition is achieved by organizing the problem and the subordinated subproblems in a tree hierarchy and optimizing each subsystem separately. Coupling of the subproblems is accounted for by subsequent optimization of the entire system based on sensitivities of the suboptimization problem solutions at each level of the tree to variables of the next higher level. A formalization of the procedure suitable for computer implementation is developed and the state of readiness of the implementation building blocks is reviewed showing that the ingredients for the development are on the shelf. The decomposition method is also shown to be compatible with the natural human organization of the design process of engineering systems. The method is also examined with respect to the trends in computer hardware and software progress to point out that its efficiency can be amplified by network computing using parallel processors.

Sobieszczanski-Sobieski, J.↗

An Efficient Objective Analysis System for Parallel Computers

A new objective analysis system designed for parallel computers will be described. The system can produce a global analysis (on a 2 x 2.5 lat-lon grid with 20 levels of heights and winds and 10 levels of moisture) using 120,000 observations in less than 3 minutes on 32 CPUs (SGI Origin 2000). No special parallel code is needed (e.g. MPI or multitasking) and the 32 CPUs do not have to be on the same platform. The system Ls totally portable and can run on -several different architectures at once. In addition, the system can easily scale up to 100 or more CPUS. This will allow for much higher resolution and significant increases in input data. The system scales linearly as the number of observations and the number of grid points. The cost overhead in going from I to 32 CPus is 18%. in addition, the analysis results are identical regardless of the number of processors used. T'his system has all the characteristics of optimal interpolation, combining detailed instrument and first guess error statistics to produce the best estimate of the atmospheric state. It also includes a new quality control (buddy check) system. Static tests with the system showed it's analysis increments are comparable to the latest NASA operational system including maintenance of mass-wind balance. Results from a 2-month cycling test in the Goddard EOS Data Assimilation System (GEOS DAS) show this new analysis retains the same level of agreement between the first guess and observations (0-F statistics) throughout the entire two months.

Stobie, James G.↗

Developing Extended Reality Projects in Support of Design, Fabrication and Procedure

The goal of this internship was to improve and create virtual reality simulations and visualizations for use in parallel with the design, fabrication, and analysis of flight ready hardware for areas like the Environmental Control and Life Support Systems (ECLSS) and also Space Systems. Specifically, my work was done in the XRSpace lab at Marshall Space Flight Center (MSFC) with assistance directly and indirectly from workers at KSC, JSC and LaRC. Led by David Reynolds, the XRSpace lab develops products for various entities at NASA. The work done in the XRSpace lab focuses on Extended Reality (XR) solutions for both simulations and visualization capabilities. The goal of the lab is to support the larger systems of NASA and to help find ways that XR technologies can streamline and optimize the design process. Extended Reality is an umbrella term that encompasses Mixed Reality, Augmented Reality, and Virtual Reality. In this capacity, I was able to complete several elements in the design, building, testing, and deployment for a variety of immersive experiences, including a VR procedure simulation, visualization aids, and a 360-image capture tool.

Joseph Chambers↗

Five myofibrillar lesion types in eccentrically challenged, unloaded rat adductor longus muscle--a test model

Sarcomere disruptions are observed in the adductor longus (AL) muscles following voluntary reloading of spaceflown and hindlimb suspension unloaded (HSU) rat, which resemble lesions in eccentrically challenged muscle. We devised and tested an eccentric contraction (ECCON) test system for the 14-day HSU rat AL. Six to 7 hours following ECCON, ALs were fixed to allow immunostaining and electron microscopy (EM). Toluidine blue-stained histology semithin sections were screened for lesion density (#/mm2). Serial semithin sections from the ECCON group were characterized for myosin immunointensity of lesions. Five myofibrillar lesion types were identified in histological semithin sections: focal contractions; wide A-bands; opaque areas; missing A-bands; and hyperstretched sarcomeres. Lesion density by type was greater for ECCON than NonECCON ALs (P< or =0.05; focal contractions and opaque regions). Lesion density (#-of-all-five-types/mm2) was significantly different (ECCON: 23.91+/-10.58 vs. NonECCON: 5.48+/-1.28, P< or =0.05; ECCON vs. SHAM: 0.00+/-0.00; P< or = 0.025). PostECCON optimal tension decreased (Poi-drop, 17.84+/-4.22%) and was correlated to lesion density (R2=0.596), but prestretch tension demonstrated the highest correlation with lesion density (R2=0.994). In lesions, the darkly staining A-band lost the normally organized thick filament alignment to differing degrees across the different lesion types. Ranking the five lesion types by a measure of lesion length deformation (hypercontracted to hyperstretched) at the light microscopy level, related to the severity of thick filament registry loss across the lesion types at the electron microscopic level. This ranking suggested that the five lesion types seen in semithin sections at the light level represented a lesion progression sequence and paralleled myosin immunostaining loss as the distorted A-band filaments spread across the hyperlengthening lesion types. Lesion ultrastructure indicated damage involved calcium homeostasis loss (focal contraction lesions) and "thick-filament-centering" failure of titin (wide A-band lesions) in the early stages of lesion development.

NASA Discipline Musculoskeletal↗

Design Evolution and Verification of the A-3 Chemical Steam Generator

Following is an overview of the Chemical Steam Generator system selected to provide vacuum conditions for a new altitude test facility, the A-3 Test Stand at Stennis Space Center (SSC) in Bay St. Louis, MS. A-3 will serve as NASA s primary facility for altitude testing of the J-2X rocket engine, to be used as the primary propulsion device for the upper stages of the Ares launch vehicles. The Chemical Steam Generators (CSGs) will produce vacuum conditions in the test cell through the production and subsequent supersonic ejection of steam into a diffuser downstream of the J-2X engine nozzle exit. The Chemical Steam Generators chosen have a rich heritage of operation at rocket engine altitude test facilities since the days of the Apollo program and are still in use at NASA White Sands Test Facility (WSTF) in New Mexico. The generators at WSTF have been modified to a degree, but are still very close to the heritage design. The intent for the A-3 implementation is to maintain this heritage design as much as possible, making minimal updates only where necessary to substitute for obsolete parts and to increase reliability. Reliability improvements are especially desired because the proposed system will require 27 generators, which is nine times the largest system installed in the 1960s. Improvements were suggested by the original design firm, Reaction Motors, by NASA SSC and NASA WSTF engineers, and by the A-3 test stand design contractor, Jacobs Technology, Inc. (JTI). This paper describes the range of improvements made to the design to date, starting with the heritage generator and the minor modifications made over time at WSTF, to the modernized configuration which will be used at A-3. The paper will discuss NASA s investment in modifications to SSC s E-2 test facility fire a full-scale Chemical Steam Generator in advance of the larger steam system installation at A-3. Risk mitigation testing will be performed in early 2009 at this test facility to verify that the CSGs operate as expected. The generator which will undergo this testing is of the most recent A-3 configuration, and will be instrumented far in excess of what is normally required for operation. The extra data will allow for easier troubleshooting and more complete knowledge of expected generator performance. In addition, the early testing will give SSC personnel experience in operating the CSG systems, which will expedite the process of installation and activation at A-3. Each Chemical Steam Generator is supported by a complement of valves, instruments, and flow control devices, with the entire assembly called a "module." The generators will be installed in groups of three, historically called "units". A module is so called because of its modular ability to be replaced or serviced without disturbing the other two modules installed on the same unit. A module is pictured in Figure 1, shown with its generator secured by white bands in its shipping (vs. installed) configuration. The heritage system at WSTF is composed of a single unit (three generator modules), pictured in Figure 2 as it was installed in 1965. In contrast, A-3 will have nine units operating in parallel to achieve vacuum conditions appropriate for testing the J-2X engine. Each of the combustors operates in two modes and achieves the so-called "full-steam" mode after all three of its stages ignite. Ignition of the first stage is achieved by exciting a spark plug; the second stage and main stage are lit by the flame front of the previous stage. The main stage burns approximately 97% of the total propellant flow and uses the heat energy to vaporize water into superheated steam. While the main stage remains unlit, the combustor is in so-called "idle" mode. In the WSTF system, this idle mode is not optimized for water usage, and does not need to be, as the water is pumped from a large reservoir. The water supply at A-3 will be contained in tanks with finite volume, so water optimization is preferred for the modnized configuration. Multiple solutions for this issue have been proposed, with the leading concept being a change to the operational definition of "idle mode," with the generator running in a lower heat flux condition.

Kirchner, Casey K.↗

Self regulating body bias generator

The back bias voltage on a functional circuit is controlled through a closed loop process. A delay element receives a clock pulse and produces a delay output. The delay element is advantageously constructed of the same materials as the functional circuit so that the aging and degradation of the delay element parallels the degradation of the functional circuit. As the delay element degrades, the transistor switching time increases, increasing the time delay of the delay output. An AND gate compares a clock pulse to an output pulse of the delay element, the AND output forming a control pulse. A duty cycle of the control pulse is determined by the delay time between the clock pulse and the delay element output. The control pulse is received at the input of a charge pump. The charge pump produces a back bias voltage which is then applied to the delay element and to the functional circuit. If the time delay produced by the delay element exceeds the optimal delay, the duty cycle of the control pulse is shortened, and the back bias voltage is lowered, thereby increasing the switching speed of the transistors in the delay element and reducing the time delay. If the throughput of the delay element is too fast, the duty cycle of the control pulse is lengthened, raising the back bias voltage produced by the charge pump. This, in turn, lowers the switching speed of the transistors in both the delay element and the functional circuit. The slower switching speed in the delay element increases time delay. In this manner, the switching speed of the delay element, and of the functional circuit, is maintained at a constant level over the life of the circuit.

Hass, Kenneth↗

Reusable and Extensible High Level Data Distributions

This paper presents a reusable design of a data distribution framework for data parallel high performance applications. We are implementing the design in the context of the Chapel high productivity programming language. Distributions in Chapel are a means to express locality in systems composed of large numbers of processor and memory components connected by a network. Since distributions have a great effect on,the performance of applications, it is important that the distribution strategy can be chosen by a user. At the same time, high productivity concerns require that the user is shielded from error-prone, tedious details such as communication and synchronization. We propose an approach to distributions that enables the user to refine a language-provided distribution type and adjust it to optimize the performance of the application. Additionally, we conceal from the user low-level communication and synchronization details to increase productivity. To emphasize the generality of our distribution machinery, we present its abstract design in the form of a design pattern, which is independent of a concrete implementation. To illustrate the applicability of our distribution framework design, we outline the implementation of data distributions in terms of the Chapel language.

patterns↗

Nonrecursive formulations of multibody dynamics and concurrent multiprocessing

Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Kurdila, Andrew J.↗

Experimental Study of Split-Path Transmission Load Sharing

Split-path transmissions are promising, attractive alternatives to the common planetary transmissions for helicopters. The split-path design offers two parallel paths for transmitting torque from the engine to the rotor. Ideally, the transmitted torque is shared equally between the two load paths; however, because of manufacturing tolerances, the design must be sized to allow for other than equal load sharing. To study the effect of tolerances, experiments were conducted using the NASA split-path test gearbox. Two gearboxes, nominally identical except for manufacturing tolerances, were tested. The clocking angle was considered to be a design parameter and used to adjust the load sharing of an otherwise fixed design. The torque carried in each path was measured for a matrix of input torques and clocking angles. The data were used to determine the optimal value and a tolerance for the clocking angles such that the most heavily loaded split path carried no greater than 53 percent of an input shaft torque of 367 N-m. The range of clocking angles satisfying this condition was -0.0012 +/- 0.0007 rad for box 1 and -0.0023 +/- 0.0009 rad for box 2. This study indicates that split-path gearboxes can be used successfully in rotorcraft and can be manufactured with existing technology.

Krantz, Timothy L.↗

A Generalized Framework for Constrained Design Optimization of General Supersonic Configurations Using Adjoint Based Sensitivity Derivatives

The Aeronautics Research Mission Directorate (ARMD) sent out an NASA Research Announcement (NRA) for proposals soliciting research and technical development. The proposed research program was aimed at addressing the desired milestones and outcomes of ROA (ROA-2006) Subtopic A.4.1.1 Advanced Computational Methods. The second milestone, SUP.1.06.02 Robust, validated mesh adaptation and error quantification for near field Computational Fluid Dynamics (CFD), was addressed by the proposed research. Additional research utilizing the direct links to geometry through a CAD interface enabled by this work will allow for geometric constraints to be applied and address the final milestone, SUP2.07.06 Constrained low-drag supersonic aerodynamic design capability. The original product of the proposed research program was an integrated system of tools that can be used for the mesh mechanics required for rapid high fidelity analysis and for design of supersonic cruise vehicles. These Euler and Navier-Stokes volume grid manipulation tools were proposed to efficiently use parallel processing. The mesh adaptation provides a systematic approach for achieving demonstrated levels of accuracy in the solutions. NASA chose to fund only the mesh generation/adaptation portion of the proposal. So this report describes the completion of the proposed tasks for mesh creation, manipulation and adaptation as it pertains to sonic boom prediction of supersonic configurations.

Karman, Steve L., Jr.↗

Instrumentation and Control Needs for Reliable Operation of Lunar Base Surface Nuclear Power Systems

As one of the near-term goals of the President's Vision for Space Exploration, establishment of a multi-person lunar base will require high-endurance power systems which are independent of the sun, and can operate without replenishment for several years. These requirements may be obtained using nuclear power systems specifically designed for use on the lunar surface. While it is envisioned that such a system will generally be supervised by humans, some of the evolutions required maybe semi or fully autonomous. The entire base complement for near-term missions may be less than 10 individuals, most or all of which may not be qualified nuclear plant operators and may be off-base for extended periods thus, the need for power system autonomous operation. Startup, shutdown, and load following operations will require the application of advanced control and health management strategies with an emphasis on robust, supervisory, coordinated control of, for example, the nuclear heat source, energy conversion plant (e.g., Brayton Energy Conversion units), and power management system. Autonomous operation implies that, in addition to being capable of automatic response to disturbance input or load changes, the system is also capable of assessing the status of the integrated plant, determining the risk associated with the possible actions, and making a decision as to the action that optimizes system performance while minimizing risk to the mission. Adapting the control to deviations from design conditions and degradation due to component failures will be essential to ensure base inhabitant safety and mission success. Intelligent decisions will have to be made to choose the right set of sensors to provide the data needed to do condition monitoring and fault detection and isolation because of liftoff weight and space limitations, it will not be possible to have an extensive set of instruments as used for earth-based systems. Advanced instrumentation and control technologies will be needed to enable this critical functionality of autonomous operation. It will be imperative to consider instrumentation and control requirements in parallel to system configuration development so as to identify control-related, as well as integrated system-related, problem areas early to avoid potentially expensive work-arounds . This paper presents an overview of the enabling technologies necessary for the development of reliable, autonomous lunar base nuclear power systems with an emphasis on system architectures and off-the-shelf algorithms rather than hardware. Autonomy needs are presented in the context of a hypothetical lunar base nuclear power system. The scenarios and applications presented are hypothetical in nature, based on information from open-literature sources, and only intended to provoke thought and provide motivation for the use of autonomous, intelligent control and diagnostics.

Turso, James↗

New laser materials for laser diode pumping

The potential advantages of laser diode pumped solid state lasers are many with high overall efficiency being the most important. In order to realize these advantages, the solid state laser material needs to be optimized for diode laser pumping and for the particular application. In the case of the Nd laser, materials with a longer upper level radiative lifetime are desirable. This is because the laser diode is fundamentally a cw source, and to obtain high energy storage, a long integration time is necessary. Fluoride crystals are investigated as host materials for the Nd laser and also for IR laser transitions in other rare earths, such as the 2 micron Ho laser and the 3 micron Er laser. The approach is to investigate both known crystals, such as BaY2F8, as well as new crystals such as NaYF8. Emphasis is on the growth and spectroscopy of BaY2F8. These two efforts are parallel efforts. The growth effort is aimed at establishing conditions for obtaining large, high quality boules for laser samples. This requires numerous experimental growth runs; however, from these runs, samples suitable for spectroscopy become available.

Jenssen, H. P.↗

Mission and Vehicle-Level Updates for the Parallel Electric-Gas Architecture with Synergistic Utilization Scheme (PEGASUS) Concept Aircraft

NASA created the PEGASUS concept with the goal of lowering mission energy (a surrogate for operating cost) compared to other regional aircraft by leveraging electrified aircraft propulsion (EAP). Since its inception, researchers have explored multiple facets of PEGASUS in varying fidelity but have not completed a rigorous, integrated design. The goal of this memorandum is to provide an updated design using recent studies and improved methods. This memorandum explores the initial vehicle concept and concept of operations, while considering ways to improve both the mission concept of operations and the integrated vehicle-level performance. Additionally, the design and analysis methodologies for EAP-enabled aircraft concepts are improved in several areas. This research incorporates new propulsion-airframe integration and wing weight surrogates to model the impacts of wingtip propulsors on the configuration. Detailed weight and balance calculations enable calculating dynamic stability and flight qualities within the conceptual design environment. Ultimately, the vehicle is optimized to reduce "well-to-wake" equivalent CO2, CO2e, rather than minimizing either fuel (or total energy) consumption or maximum takeoff weight. Using fuel/energy or takeoff weight leads to conflicting optima for hybrid-electric aircraft. Two aircraft are developed to provide points of comparison for PEGASUS: an advanced conventional turboprop vehicle and a hybrid-electric variant. The results show that the PEGASUS concept can reduce CO2e relative to the advanced turboprop or a hybrid-electric propulsion architecture, albeit with an increase in maximum takeoff weight. PEGASUS's maximum takeoff weight is 55% heavier than the advanced conventional turboprop but releases 18% less CO2e for a 400 nmi mission. Over the same mission, PEGASUS's maximum takeoff weight is 47% heavier than the comparator hybrid-electric vehicle but releases 12% less CO2e. This study shows that the PEGASUS configuration reduces CO2e through its use of wingtip propulsors and that its benefit is not solely a result of switching to a hybrid-electric propulsion architecture. PEGASUS achieves this reduction in CO2e while maintaining satisfactory Level 1 or 2 flight qualities for all of its longitudinal- and lateral-directional modes.

Hybrid-electric↗

Optimization of the oxidant supply system for combined cycle MHD power plants

An in-depth study was conducted to determine what, if any, improvements could be made on the oxidant supply system for combined cycle MHD power plants which could be reflected in higher thermal efficiency and a reduction in the cost of electricity, COE. A systematic analysis of air separation process varitions which showed that the specific energy consumption could be minimized when the product stream oxygen concentration is about 70 mole percent was conducted. The use of advanced air compressors, having variable speed and guide vane position control, results in additional power savings. The study also led to the conceptual design of a new air separation process, sized for a 500 MW sub e MHD plant, referred to a internal compression is discussed. In addition to its lower overall energy consumption, potential capital cost savings were identified for air separation plants using this process when constructed in a single large air separation train rather than multiple parallel trains, typical of conventional practice.

Juhasz, A. J.↗

A simple, first step to the optimization of regenerator geometry

This paper presents a simplified set of equations for calculating and optimizing regenerator geometries. A number of competing parameters must be accounted for in the design of a regenerator. To obtain high values of effectiveness, there is need for a large surface area for heat transfer and high heat capacity. The void volume should be small to maintain the pressure ratio in the entire system. Moreover, the pressure drop should be small compared to the absolute pressure. With the assumptions made here, the calculations can be done with a hand calculator. The equations show that the optimum regenerator length, hydraulic radius, and porosity are independent of mass flow rate and the gas cross sectional area is proportional to the mass flow rate. It is shown that using gas gaps between parallel plates produces a significantly better regenerator than is possible with packed spheres or screens, particularly at temperatures below approximately 50%.

Radebaugh, R.↗

The Fundamental Gain Limitation of High-Rate Gaseous Detectors and Optimum Designs for High-Rate Applications

We have performed a systematic study of the breakdown mechanism of high-rate detectors: MSGC, MGC, MICROMEGAS and GEM, recently chosen or considered as candidates for high-luminosity applications, together with newly elaborated designs such as high-rate wire chambers, thin gap parallel-plate chamber and super-high-rate RPC's. It was found for all these that the maximum breakdown-limited dramatically with incident flux and further decreases in the presence of alpha particles which are typical of the backgrounds in high-energy experiments. It was determined that the key parameter that limits a detector's performance is the effective current density in the avalanche, and that above a certain value breakdowns occur. At these limits the subsequent breakdowns take place by a mechanism which does not seem to have been previously reported in the literature. We present the results of this study together with a qualitative theory of this new type of breakdown and suggestions for optimizing high-rate detectors which we have verified experimentally.

Fonte, Paulo↗