SEARCH · Engineering Papers

Results for “Smooth approximation”

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 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS2014

A Revised Calibration Function and Results for the Phoenix Mission TECP Relative Humidity Sensor

The original calibration function of the R(sub H) sensor on the Phoenix Thermal and Electrical Conductivity Sensor (TECP) has been revised in order to extend the range of the valid calibration, and to improve accuracy. The original function returned non-physical R(sub H) values at the lowest temperatures. To resolve this, and because the original calibration was performed against a pair of hygrometers that measured frost point (T(sub f)), the revised calibration equation is also cast in terms of frost point. Because of the complexity of maintaining very low temperatures and high R(sub H) in the laboratory, no calibration data exists at T is greater than 203K. However, sensor response during the mission was smooth and continuous down to 181 K. Therefore we have opted to include flight data in the calibration data set; selection was limited to data acquired during periods when the atmosphere is known to have been saturated. T(sub f) remained below 210 K throughout the mission(P is greater than 0.75 Pa). R(sub H), conversely, ranged from 1 to well under 0.01 diurnally, due to approximately 50 K temperature variations. To first order, both vapor pressure and its variance are greater during daylight hours. Variance in overnight humidity is almost entirely explained by temperature, while atmospheric turbulence contributes substantial variance to daytime humidity. Likewise, data gathered with the TECP aloft reflect higher H2O abundances than at the surface, as well as greater variance. There is evidence for saturation of the atmosphere overnight throughout much of the mission. In virtually every overnight observation, once the atmosphere cooled to T(sub f), water vapor begins to decrease, and tracks air temperature. There is no evidence for substantial decreases in water vapor prior to saturation, as expected for adsorptive exchange. Likewise, there is no evidence of local control of vapor by phases such as perchlorate hydrates hydrated minerals. The daytime average H2O pressure does not change substantially over the course of the mission, although the H2O column abundance varies by a factor of 2. Column abundances calculated from TECP data are lower than orbital measurements if one assumes that H2O is uniformly mixed through a single scale height. These results argue that the vertical distribution of H2O begins to change well in advance of surface concentrations as northern autumn approaches.

Phoenix Thermal and Electrical Conductivity Sensor↗

NASA NTRS2000

Hybrid Active/Passive Control of Sound Radiation from Panels with Constrained Layer Damping and Model Predictive Feedback Control

There has been considerable interest over the past several years in applying feedback control methods to problems of structural acoustics. One problem of particular interest is the control of sound radiation from aircraft panels excited on one side by a turbulent boundary layer (TBL). TBL excitation appears as many uncorrelated sources acting on the panel, which makes it difficult to find a single reference signal that is coherent with the excitation. Feedback methods have no need for a reference signal, and are thus suited to this problem. Some important considerations for the structural acoustics problem include the fact that the required controller bandwidth can easily extend to several hundred Hertz, so a digital controller would have to operate at a few kilohertz. In addition, aircraft panel structures have a reasonably high modal density over this frequency range. A model based controller must therefore handle the modally dense system, or have some way to reduce the bandwidth of the problem. Further complicating the problem is the fact that the stiffness and dynamic properties of an aircraft panel can vary considerably during flight due to altitude changes resulting in significant resonant frequency shifts. These considerations make the tradeoff between robustness to changes in the system being controlled and controller performance especially important. Recent papers concerning the design and implementation of robust controllers for structural acoustic problems highlight the need to consider both performance and robustness when designing the controller. While robust control methods such as H1 can be used to balance performance and robustness, their implementation is not easy and requires assumptions about the types of uncertainties in the plant being controlled. Achieving a useful controller design may require many tradeoff studies of different types of parametric uncertainties in the system. Another approach to achieving robustness to plant changes is to make the controller adaptive. For example, a mathematical model of the plant could be periodically updated as the plant changes, and the feedback gains recomputed from the updated model. To be practical, this approach requires a simple plant model that can be updated quickly with reasonable computational requirements. A recent paper by the authors discussed one way to simplify a feedback controller, by reducing the number of actuators and sensors needed for good performance. The work was done on a tensioned aircraft-style panel excited on one side by TBL flow in a low speed wind tunnel. Actuation was provided by a piezoelectric (PZT) actuator mounted on the center of the panel. For sensing, the responses of four accelerometers, positioned to approximate the response of the first radiation mode of the panel, were summed and fed back through the controller. This single input-single output topology was found to have nearly the same noise reduction performance as a controller with fifteen accelerometers and three PZT patches. This paper extends the previous results by looking at how constrained layer damping (CLD) on a panel can be used to enhance the performance of the feedback controller thus providing a more robust and efficient hybrid active/passive system. The eventual goal is to use the CLD to reduce sound radiation at high frequencies, then implement a very simple, reduced order, low sample rate adaptive controller to attenuate sound radiation at low frequencies. Additionally this added damping smoothes phase transitions over the bandwidth which promotes robustness to natural frequency shifts. Experiments were conducted in a transmission loss facility on a clamped-clamped aluminum panel driven on one side by a loudspeaker. A generalized predictive control (GPC) algorithm, which is suited to online adaptation of its parameters, was used in single input-single output and multiple input-single output configurations. Because this was a preliminary look at the potential constrained layer damping for adaptive control, static feedback control with no online adaptation was used. Two configurations of CLD in addition to a bare panel configuration were studied. For each CLD configuration, two sensor arrangements for the feedback controller were compared. The first arrangement used fifteen accelerometers on the panel to estimate the responses of the first six radiation modes of the panel. The second sensor arrangement was simpler, using the summed responses of only four accelerometers to approximate the response of the first radiation mode of the panel. In all cases a PZT patch was mounted at the center of the panel for control input. The performance of the controller was quantified using the responses of the fifteen accelerometers on the panel to estimate radiated sound power. The paper begins with a brief discussion of the GPC algorithm and the experimental setup. The experimental results are discussed next, comparing the CLD and sensor configurations, followed by discussion and conclusions.

Cabell, Randolph H.↗

DOE OSTI2021

Schrödinger evolution of two-dimensional black holes

This paper systematically treats the evolving quantum state for two-dimensional black holes, with particular focus on the CGHS model, but also elucidating features generalizing to higher dimensions. This is done in Schrödinger picture(s), to exhibit the dynamic evolution of the state at intermediate times. After a review of classical solutions, also connecting to descriptions of higher-dimensional black holes, it overviews the canonical quantum treatment of the full evolution, including gravitational dynamics. Derived in an approximation to this, following conversion to “perturbation picture”, is the evolution of the quantum matter on the background geometry. Features of the evolving matter state are described, based on choice of a time slicing to put the evolution into ADM form. The choices of slicing as well as coordinates on the slices result in different quantum “pictures” for treating the evolution. If such a description is based on smooth trans-horizon slices, that avoids explicit reference to ultra-planckian modes familiar from traditional treatments, and exhibits the Hawking excitations as emerging from a “quantum atmosphere” with thickness comparable to the inverse temperature. Detailed study of the state exhibits the entanglement structure between Hawking quanta and the partner excitations inside the black hole, and the corresponding “missing information”. This explicit description also allows direct study of the evolution and features, e.g. as seen by infalling observers, of these partner excitations, helping to address various puzzles with them. Explicit treatment of the evolving state, and its extension to higher dimensions, provides further connections to information theory and a starting point for study of corrections that can unitarize evolution, arising from new quantum gravity effects — whether wormholes or something entirely different.

2D Gravity↗

NASA NTRS1999

Detection of Crystalline Hematite Mineralization on Mars by the Thermal Emission Spectrometer: Evidence for Near-surface Water

The Thermal Emission Spectrometer (TES) instrument on the Mars Global Surveyor (MGS) mission has discovered a remarkable accumulation of crystalline hematite ((alpha-Fe2O3) that covers an area with very sharp boundaries approximately 350 by 350-750 km in size centered near 2 S latitude between 0 and 5 W longitude (Sinus Meridiani). Crystalline hematite is uniquely identified by the presence of fundamental vibrational absorption features centered near 300, 450, and >525/cm, and by the absence of silicate fundamentals in the 1000/cm region. Spectral features resulting from atmospheric CO2, dust, and water ice were removed using a radiative transfer model. The spectral properties unique to Sinus Meridiani were emphasized by removing the average spectrum of the surrounding region. The depth and shape of the hematite fundamental bands show that the hematite is crystalline and relatively coarse grained (>5-10 micron). Diameters up to and greater than 100s of micrometers are permitted within the instrumental noise and natural variability of hematite spectra. Hematite particles <5-10 micron in diameter (either as an unpacked or hard-packed powders) fail to match the TES spectra. The spectrally-derived areal abundance of hematite varies with particle size from approximately 10% for particles >30 micron in diameter to 40-60% for unpacked 10 micron powders. The hematite in Sinus Meridiani is thus distinct from the fine-grained (diameter <5-10 micron), red, crystalline hematite considered, on the basis of visible, near-IR data, to be a minor spectral component in Martian bright regions like Olympus-Amazonis. Sinus Meridiani hematite is closely associated with a smooth, layered, friable surface that is interpreted to be sedimentary in origin. This material may be the uppermost surface in the region, indicating that it could be a late-stage sedimentary unit, or it could be a layered portion of the heavily cratered plains units. We consider five possible mechanisms for the formation of coarse-grained, crystalline hematite. These processes fall into two classes depending on whether they require a significant amount of near-surface water: (1) chemical precipitation that includes origin by (a) precipitation from oxygenated, Fe-rich water (iron formations), (b) hydrothermal extraction and crystal growth.

Christensen, P. R.↗

DOE OSTI2024

The vertical-velocity skewness in the inertial sublayer of turbulent wall flows

Empirical evidence is provided that within the inertial sublayer (i.e. logarithmic region) of adiabatic turbulent flows over smooth walls, the skewness of the vertical-velocity component 𝑆 𝑤 displays universal behaviour, being a positive constant and constrained within the range 𝑆 𝑤 ≈ 0.1–0.16, regardless of flow configuration and Reynolds number. A theoretical model is then proposed to explain this behaviour, including the observed range of variations of 𝑆 𝑤 . The proposed model clarifies why 𝑆 𝑤 cannot be predicted from down-gradient closure approximations routinely employed in large-scale meteorological and climate models. The proposed model also offers an alternative and implementable approach for such large-scale models.

turbulent boundary layers↗

NASA NTRS1974

Aerodynamic influence coefficient method using singularity splines

A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.

Mercer, J. E.↗

NASA NTRS1987

Fixed lag smoothers for carrier phase and frequency tracking

The application of fixed lag smoothing algorithms are presented for the problem of estimation of the phase and frequency of a sinusoidal carrier received in the presence of process noise and additive observation noise. A suboptimal structure consists of a phase-locked loop (PLL) followed by a post loop correction to the phase and frequency estimates. When the PLL is operating under a high signal-to-noise ratio, the phase detector is approximately linear, and the smoother equations then correspond to the optimal linear equations for an equivalent linear signal model. The performance of such a smoother can be predicted by the linear filtering theory. However, if the PLL is operating near the threshold region of the signal to noise ratio, the phase detector cannot be assumed to be linear. Then the actual performance of the smoother can significantly differ from that predicted by linear filtering theory. Both the theoretical and simulated performance of such smoothers derived on the basis of various models for the phase of frequency processes are presented.

Kumar, R.↗

NASA NTRS1994

Development of an adaptive hp-version finite element method for computational optimal control

In this research effort, the usefulness of hp-version finite elements and adaptive solution-refinement techniques in generating numerical solutions to optimal control problems has been investigated. Under NAG-939, a general FORTRAN code was developed which approximated solutions to optimal control problems with control constraints and state constraints. Within that methodology, to get high-order accuracy in solutions, the finite element mesh would have to be refined repeatedly through bisection of the entire mesh in a given phase. In the current research effort, the order of the shape functions in each element has been made a variable, giving more flexibility in error reduction and smoothing. Similarly, individual elements can each be subdivided into many pieces, depending on the local error indicator, while other parts of the mesh remain coarsely discretized. The problem remains to reduce and smooth the error while still keeping computational effort reasonable enough to calculate time histories in a short enough time for on-board applications.

Hodges, Dewey H.↗

NASA NTRS2009

Gender-related Changes in Dorsal Hand and Foot Vein Function Following 60 Days of Head Down Bed Rest

It is well known that female astronauts are more likely to experience post-flight orthostatic hypotension and presyncope compared to male astronauts. It has been suggested that the disproportionally higher incidence of presyncope (83% of female vs. 20% male crewmembers) may be due to sex-related differences in vascular function between the upper and lower limbs. However, much of this evidence is specific to changes in resistance vessels. Given that more than 70% of the circulating blood volume resides in compliance vessels, it is conceivable that even small changes in venous function may contribute to post-flight orthostatic hypotension. In spite of this, little is currently known regarding the influence of microgravity exposure on venous function between males and females. PURPOSE: To determine the influence of 60 days of HDBR on dorsal foot and hand vein function between healthy males (M) and females (F). METHODS: Using 2-D ultrasound, dorsal hand and foot vein diameter responses to intravenous infusions phenylephrine (PE), acetylcholine (ACh), and nitroglycerine (NTG) were determined in 26 adults; 10 females (age:37 +/- 2 yr ) and 16 males (age:34 +/- 2 yr ). Changes in venous function were calculated as the difference between diameter at baseline and following each venoactive drug. Differences in venous function between limb and sexes across HDBR were determined using mixed-effects linear regression. RESULTS: In response to 60 days of HDBR, the change in venousconstrictor response to PE in the dorsal hand veins was not significantly different between M and F. Interestingly, the change in constrictor response in the dorsal foot veins (compared to pre HDBR) was approximately 30% greater in the F, whereas the constrictor response was approximately 45% less in the M (p=0.026). HDBR had no influence on the change in dilator response to ACh, or NTG between M and F and between vascular beds. CONCLUSION: These results demonstrate that 60 days of HDBR contributes to sex-related differences in venous function. Importantly, the changes are specific to dorsal foot veins and do not appear to be related to changes in endothelial or smooth muscle function. Previous data from our laboratory indicate that females have an impaired sympathetic response following HDBR. However, the results from the current investigation imply that F maintain the ability to respond to sympathetic stimulation but may not be able to adequately activate a proper response following HDBR.

Westby, Christian M.↗

NASA NTRS2017

Growth Angle: A Microscopic View

A microscopic continuum mechanical model of the growth angle is proposed. It is based on the van der Waals type framework that is used for surface force phenomena. The obtained augmented Laplace type integro-differential equations are, in general, difficult to analyze. Here we focused primarily on the particular case of equal melt and crystal surface energies. We derived an approximate equation for the meniscus shape, and obtained an analytical relationship between the contact and the growth angle. Interestingly, the same result can be obtained using the macroscopic model of Herring. The case of a macroscopically sharp corner is also considered. For this case, the macroscopic angle is not defined and it can be any angle between the contact angles of both flat surfaces. The microscopic model yields the smooth shape for the meniscus that also is not unique, but depends on the initial position of the meniscus.

Mazuruk, Konstantin↗

DOE OSTI2022

Application of Recursive Estimation to Heat Tracing for Groundwater/Surface‐Water Exchange

Abstract We present and demonstrate a recursive‐estimation framework to infer groundwater/surface‐water exchange based on temperature time series collected at different vertical depths below the sediment/water interface. We formulate the heat‐transport problem as a state‐space model (SSM), in which the spatial derivatives in the convection/conduction equation are approximated using finite differences. The SSM is calibrated to estimate time‐varying specific discharge using the Extended Kalman Filter (EKF) and Extended Rauch‐Tung‐Striebel Smoother (ERTSS). Whereas the EKF is suited to real‐time (“online”) applications and uses only the past and current measurements for estimation (filtering), the ERTSS is intended for near‐real time or batch‐processing (“offline”) applications and uses a window of data for batch estimation (smoothing). The two algorithms are demonstrated with synthetic and field‐experimental data and are shown to be efficient and rapid for the estimation of time‐varying flux over seasonal periods; further, the recursive approaches are effective in the presence of rapidly changing flux and (or) nonperiodic thermal boundary conditions, both of which are problematic for existing approaches to heat tracing of time‐varying groundwater/surface‐water exchange.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1977

The stress corrosion resistance and the cryogenic temperature mechanical properties of annealed Nitronic 60 bar material

Ambient and cryogenic temperature mechanical properties and the ambient temperature stress corrosion properties of annealed, straightened, and centerless ground Nitronic 60 stainless steel alloy bar material are presented. The mechanical properties of longitudinal specimens were evaluated at test temperatures from ambient to liquid hydrogen. The tensile test data indicated increasing strength with decreasing temperature to -196 C. Below liquid nitrogen temperature the smooth tensile and notched tensile strengths decreased slightly while the elongation and reduction of area decreased drastically. The Charpy V-notched impact energy decreased steadily with decreasing test temperature. Stress corrosion tests were performed on longitudinal tensile specimens and transverse C-ring specimens exposed to: alternate immersion in a 3.5% NaCl bath; humidity cabinet; and a 5% salt spray atmosphere. The longitudinal tensile specimens experienced no corrosive attack. Approximately 3/4 of the transverse C-rings exposed to alternate immersion and to salt spray experienced a pitting attack on the top and bottom ends. Additional stress corrosion tests were performed on transverse tensile specimens. No failures occurred in the 90% stressed specimens exposed for 90 days in the alternate immersion and salt spray environments

Montano, J. W. L.↗

NASA NTRS1975

Long-term differential energy spectrum for solar-flare iron-group particles

A long-term solar-flare differential energy spectrum for iron-group nuclei from approximately 0.1 to approximately 600 MeV/amu is derived from track density profile measurements in sample 64455 and sample 68815. Measurements from uneroded surfaces were obtained from quench crystals of plagioclase in 64455, and a Kr-81/Kr method indicates that the exposure age of this sample is 2,010,000 yrs. The power laws which best fit the normalized track density data are reported; the energy spectrum consists of two power law curves smoothly joined together which in turn are smoothly connected to a modulated galactic cosmic-ray spectrum. Standard track production versus depth profiles can be used to determine solar-flare track exposure ages and erosion rates for lunar samples.

Blanford, G. E.↗

DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗

DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗

NASA NTRS1991

High order filtering methods for approximating hyperbolic systems of conservation laws

The essentially nonoscillatory (ENO) schemes, while potentially useful in the computation of discontinuous solutions of hyperbolic conservation-law systems, are computationally costly relative to simple central-difference methods. A filtering technique is presented which employs central differencing of arbitrarily high-order accuracy except where a local test detects the presence of spurious oscillations and calls upon the full ENO apparatus to remove them. A factor-of-three speedup is thus obtained over the full-ENO method for a wide range of problems, with high-order accuracy in regions of smooth flow.

Lafon, F.↗

NASA NTRS1972

Applicability of a diffusion model to lateral transport in the terrestrial and lunar exospheres.

Kinetic theory is used to determine a series expansion of the vertical flux of particles in an exosphere in terms of time and space derivatives of particle concentration, exobase velocity, and temperature. For sufficiently large scale variations of these parameters in time and space, the series can be truncated to a form that is similar to a diffusion equation. Owing to this analogy, it is possible to unite the mathematical description of molecular diffusion, which governs thermospheric flow, and the corresponding exospheric equation by using effective transport coefficients which change smoothly with altitude through the transition from thermosphere to exosphere. A new definition of the exobase for lateral flow emerges from the analogy of exospheric and thermospheric diffusion, as the altitude where the horizontal mean free path length equals the mean horizontal extent of ballistic trajectories of the transported gas, as opposed to the scale height of the dominant gas which determines the exobase for escape. It is shown that the approximation of exospheric lateral flow as a diffusion process is applicable to global scale problems concerning terrestrial helium and heavier gases, and lunar gases heavier than helium.

Hodges, R. R., Jr.↗

NASA NTRS1986

A comparison of the bounded derivative and the normal-mode initialization methods using real data

Application of the bounded-derivative and normal-mode methods to a simple linear barotropic model at a typical middle latitude shows that the two methods lead to identical constraints up to a certain degree of approximation. Beyond this accuracy the two methods may differ from each other. When applied to a global nonlinear barotropic model using real data, again the two methods lead to similar balanced initial states. The gravity oscillations in the unbalanced height field, which have amplitudes of up to 60 m with a dominant periodicity of about 5 to 6 h, are practically eliminated by both initialization methods. The rotational wind component is smooth even for the unbalanced initial state. The small-scale spatial features of the irrotational wind component are drastically reduced by initialization. Both the nonlinear normal-mode and the bounded-derivative initialization methods yield similar divergence fields centered around the areas of highest orography. The comparison shows that there is no significant loss of information in the mass and momentum fields, despite the fact that the bounded-derivative method employs only the original, rotational wind component to construct a balanced initial state compared to the normal-mode method, which, in addition, makes use of the unbalanced divergent wind and height fields.

Semazzi, F. H. M.↗