Search NASA⌕ Search

SEARCH · Search NASA

Results for “collocation techniques”

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 145 records · Page 8

Use of TEMPO as a Proxy for Hyperspectral Geostationary Ocean Color Measurements from the GeoXO OCX Instrument: Harnessing Machine Learning and Principal Component Techniques for Atmospheric and Glint Correction

Retrievals of ocean color from space are important for better understanding the ocean ecosystem. The launch of atmospheric geostationary hyperspectral sensors such as TEMPO, provides a unique opportunity to examine the diurnal variability in ocean ecology. While TEMPO does not have as high spatial resolution or full spectral coverage as planned coastal ocean sensors such as the Geosynchronous Littoral Imaging and Monitoring Radiometer (GLIMR) or GeoXO Ocean Color instrument (OCX), its hourly measurements provide coverage of regions such as Lake Erie and the Gulf of Mexico at spatial scales of approximately 5 km. These data can be useful for testing new algorithms. We will apply our newly developed machine learning based atmospheric correction approach for ocean color retrievals to TEMPO data. Our approach begins by decomposing measured radiances from hyperspectral sensors into spectral features that describe the scattering and absorption of the atmosphere as well as the underlying surface reflectance. The coefficients of the principal components are then used to train a neural network to predict ocean color properties derived from collocated MODIS/VIIRS physically-based retrievals. This machine learning approach does not rely on radiative transfer modeling, and the use of MODIS/VIIRS data for training accounts for possible calibration b in hyperspectral data. Previously, we applied our approach using blue and UV wavelengths with the Ozone Monitoring Instrument (OMI) and TROPOspheric Monitoring Instrument (TROPOMI) to show that it can estimate ocean color properties in less-than-ideal conditions such as lightly to moderately clouded conditions as well as sun glint and thus improve the spatial coverage of ocean color measurements. TEMPO provides an opportunity to improve on this approach since it will provide collocated measurements at green and red wavelengths that were not available from OMI and TROPOMI and are important particularly for coastal waters. Additionally, our technique can be applied early in the mission and has potential to demonstrate the value of near real time ocean color products that are important for monitoring of harmful algae blooms and other oceanic phenomena.

Zachary Fasnacht↗

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗

Numerical derivative techniques for trajectory optimization

The adoption of robust numerical optimization techniques in trajectory simulation programs has resulted in powerful design and analysis tools. These trajectory simulation/optimization programs are widely used, and a representative list includes the GTS system, the POST program, and newer collocation methods such as OTIS and FONPAC. All of these programs rely on optimization algorithms which require objective function and constraint gradient data during the iteration process. However, most trajectory optimization problems lack simple analytical expressions for these derivatives. In the general case a function evaluation involves integrating aerodynamic, propulsive, and gravity forces over multiple trajectory phases with complex control models. With the newer collocation methods, the integration is replaced by defect constraints and cubic approximations for the state. While analytic gradient expressions can sometimes be derived for trajectory optimization problems, the derivation is cumbersome, time consuming, and prone to mistakes. Fortunately, an alternate method exists for the gradient evaluation, namely finite difference approximations. In this paper some finite difference gradient techniques developed for use with the GTS system are presented. These techniques include methods for computing first and second partial derivatives of single and multiple sets of functions. A key feature of these methods is an error control mechanism which automatically adjusts the perturbation size to obtain accurate derivative values.

Hallman, Wayne P.↗

Supercooled Liquid Water Detection Capabilities from Ka-Band Doppler Profiling Radars: Moment-Based Algorithm Formulation and Assessment

The occurrence of supercooled liquid water in mixed-phase cloud (MPC) affects their cloud microphysical and radiative properties. The prevalence of MPCs in the mid- and high latitudes translates these effects to significant contributions to Earth’s radiative balance and hydrological cycle. The current study develops and assesses a radar-only, moment-based phase partition technique for the demarcation of supercooled liquid water volumes in arctic, MPC conditions. The study utilizes observations from the Ka band profiling radar, the collocated high spectral resolution lidar, and ambient temperature profiles from radio sounding deployments following a statistical analysis of 5.5 years of data (January 2014–May 2019) from the Atmospheric Radiation Measurement observatory at the North Slope of Alaska. The ice/liquid phase partition occurs via a per-pixel, neighborhood-dependent algorithm based on the premise that the partitioning can be deduced by examining the mean values of locally sampled probability distributions of radar-based observables and then compare those against the means of climatologically derived, per-phase probability distributions. Analyzed radar observables include linear depolarization ratio (LDR), spectral width, and vertical gradients of reflectivity factor and radial velocity corrected for vertical air motion. Results highlight that the optimal supercooled liquid water detection skill levels are realized for the radar variable combination of spectral width and reflectivity vertical gradient, suggesting that radar-based polarimetry, in the absence of full LDR spectra, is not as critical as Doppler capabilities. The cloud phase masking technique is proven particularly reliable when applied to cloud tops with an Equitable Threat Score (ETS) of 65%; the detection of embedded supercooled layers remains much more uncertain (ETS = 27%).

54 ENVIRONMENTAL SCIENCES↗

Estimating Ice Water Content for Winter Storms from Millimeter-Wavelength Radar Measurements Using a Synthesis of Polarimetric and Dual-Frequency Radar Observations

The potential of millimeter-wavelength radar-based ice water content (IWC) estimation is demonstrated using a Ka-band Scanning Polarimetric Radar (KASPR) for the U.S. northeast coast winter storms. Two IWC relations for Ka-band polarimetric radar measurements are proposed: one that uses a combination of the radar reflectivity Z and the estimated total number concentration of snow particles N t and the other based on the joint use of Z, specific differential phase K DP , and the degree of riming f rim . A key element of the algorithms is to obtain the “Rayleigh-equivalent” value of Z measured at the Ka band, i.e., the corresponding Z at a longer radar wavelength for which Rayleigh scattering takes place. This is achieved via polarimetric retrieval of the mean volume diameter D m and incorporating the relationship between the dual-wavelength ratio DWR S/Ka and D m . Those techniques allow for retrievals from single millimeter-wavelength radar measurements and do not necessarily require the dual-wavelength ratio (DWR) measurements, if the DWR–D m relation and Rayleigh assumption for Ka-band K DP are valid. Comparison between the quasivertical profile product obtained from KASPR and the columnar vertical profile product generated from the nearby WSR-88D S-band radar measurements demonstrates that the DWR S/Ka can be estimated from the two close radars without the need for collocated radar beams and synchronized antenna scanning and can be used for determining the Rayleigh-equivalent value of Z. As a result, the performance of the suggested techniques is evaluated for seven winter storms using surface disdrometer and snow accumulation measurements.

Cloud retrieval↗

An adaptive identification and control scheme for large space structures

A unified identification and control scheme capable of achieving space at form performance objectives under nominal or failure conditions is described. Preliminary results are also presented, showing that the methodology offers much promise for effective robust control of large space structures. The control method is a multivariable, adaptive, output predictive controller called Model Predictive Control (MPC). MPC uses a state space model and input reference trajectories of set or tracking points to adaptively generate optimum commands. For a fixed model, MPC processes commands with great efficiency, and is also highly robust. A key feature of MPC is its ability to control either nonminimum phase or open loop unstable systems. As an output controller, MPC does not explicitly require full state feedback, as do most multivariable (e.g., Linear Quadratic) methods. Its features are very useful in LSS operations, as they allow non-collocated actuators and sensors. The identification scheme is based on canonical variate analysis (CVA) of input and output data. The CVA technique is particularly suited for the measurement and identification of structural dynamic processes - that is, unsteady transient or dynamically interacting processes such as between aerodynamics and structural deformation - from short, noisy data. CVA is structured so that the identification can be done in real or near real time, using computationally stable algorithms. Modeling LSS dynamics in 1-g laboratories has always been a major impediment not only to understanding their behavior in orbit, but also to controlling it. In cases where the theoretical model is not confirmed, current methods provide few clues concerning additional dynamical relationships that are not included in the theoretical models. CVA needs no a priori model data, or structure; all statistically significant dynamical states are determined using natural, entropy-based methods. Heretofore, a major limitation in applying adaptive controllers to LSS problems has been the requirement for accurate models. We feel that MPC coupled with the CVA identification scheme will resolve this limitation. The results are based on research performed under an SBIR grant, administered by the Air Force Office of Scientific Research.

Carroll, J. V.↗

A Comprehensive Forward Model for Spaceborne Radar Instruments

We present the development and validation of a comprehensive forward model designed to enhance remote sensing capabilities of spaceborne radar instruments. To overcome limitations in existing models, we integrated a Discrete Dipole Approximation (DDA) cloud scattering database into our Radiative Transfer Model (RTM), focusing on microwave frequencies. By simulating the optical properties of non-spherical frozen hydrometeors, the DDA technique effectively reduced discrepancies between simulated and observed values, surpassing traditional Mie tables. The evaluation of DDA lookup tables involved comparisons with a collocated dataset comprising short-term forecasts and satellite microwave data, providing evidence of their superiority. Additionally, we address the challenges of assimilating active radar measurements, which offer vertically resolved insights into clouds and precipitation. We explored the assimilation of spaceborne radar measurements in Numerical Weather Prediction (NWP) models by integrating a forward radar model, along with its adjoint and tangent linear, into the data assimilation system. Evaluation using CloudSat measurements demonstrated promising agreement between simulations and observations, particularly when the input hydrometeor profiles aligned with the measured reflectivity profiles, showcasing the potential of the developed forward radar model. Moreover, we discuss other challenges in radar measurement assimilation within NWP models, including potential observation errors and biases.

Isaac Moradi↗

Coexistence Analysis of Adjacent Long Term Evolution (LTE) Systems

As the licensing and deployment of Long term evolution (LTE) systems are ramping up, the study of coexistence of LTE systems is an essential topic in civil and military applications. In this paper, we present a coexistence study of adjacent LTE systems aiming at evaluating the effect of inter-system interference on system capacity and performance as a function of some of the most common mitigation techniques: frequency guard band, base station (BS) antenna coupling loss, and user equipment (UE) antenna spacing. A system model is constructed for two collocated macro LTE networks. The developed model takes into consideration the RF propagation environment, power control scheme, and adjacent channel interference. Coexistence studies are performed for a different combination of time/frequency division duplex (TDD/FDD) systems under three different guard-bands of 0MHz, 5MHz, and 10MHz. Numerical results are presented to advice the minimum frequency guard band, BS coupling loss, and UE antenna isolation required for a healthy system operation.

adjacent channel interference↗

Creating Seasonal Climatologies of Aerosol Lidar Ratios Over Ocean Using MODIS AOD Constrained Retrievals and GOCART Model Simulations

The current CALIPSO algorithms assign one lidar ratio (i.e., extinction-to-backscatter ratio; LR) value globally for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms through the development of regional and seasonal LR climatologies. In this work, aerosol LRs are inferred through CALIOP backscatter profiles constrained by collocated aerosol optical depth (AOD) from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data. This analysis is subsampled for those profiles that are cloud-free and contain only one CALIOP-classified aerosol type (e.g., marine). The CALIOP profiles are then collocated with aerosol volume fractions obtained through Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulations. In this poster, we show twelve-year (June 2006-August 2018) mean spatial distributions of inferred aerosol LRs for CALIOP-classified marine aerosols and how they correspond inversely with patterns of GOCART sea salt volume fraction (SSVF). Near land masses, smaller SSVFs and larger LRs are found (due to the influence of over-land aerosols). In the remote ocean regions (likely less impacted by over-land aerosols), larger SSVFs and smaller LRs are found. The developed relationship between the GOCART model SSVFs and MODIS AOD constrained LRs is used to create model-assisted seasonal LR maps. Additionally, we show maps of inferred LRs from constrained retrievals using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) product and comparisons with those from the MODIS analyses. The technique demonstrated here benefits CALIPSO in the near-term, but similar methods can also be applied to the next generation space-based elastic backscatter lidars with collocated passive sensors, such as those of the upcoming NASA Atmosphere Observing System (AOS).

Travis Toth↗

Mapping Aerosol Lidar Ratios Over Ocean Using Constrained Retrievals and A Global Aerosol Model

The current NASA Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) algorithms assign one lidar ratio (i.e., extinction-to-backscatter ratio; LR) value globally for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms through the development of regional and seasonal LR climatologies. In this work, aerosol LRs are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth (AOD) from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data. This analysis is subsampled for those profiles that are cloud-free and contain only one CALIOP-classified aerosol type (e.g., marine). The CALIOP profiles are then collocated with aerosol volume fractions obtained through Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulations. In this presentation, we show that the twelve-year (June 2006-August 2018) mean spatial distributions of inferred aerosol LRs for CALIOP-classified marine aerosols correspond inversely with patterns of GOCART sea salt volume fraction (SSVF). For example, smaller SSVFs (< 65%) and larger LRs (> 55 sr), are found near land masses. This is indicative of the influence of over-land aerosols (e.g., pollution and biomass burning smoke). In the remote oceans (i.e., regions likely less impacted by non-sea salt aerosols), larger SSVFs (> 95%) and smaller LRs (< 25 sr) are found. The developed relationship between the GOCART SSVFs and MODIS AOD constrained LRs (polynomial fit intersect values of ~56 sr for SSVF of 0% and ~21 sr for SSVF of 100%) is used to produce model-assisted climatological LR maps on seasonal scales. Additionally, we show maps of inferred LRs from constrained retrievals using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) product and comparisons with those from the MODIS analyses. The technique demonstrated in this study not only benefits CALIPSO in the near-term, but similar methods can be applied to future spaceborne elastic backscatter lidars with collocated passive sensors, such as those associated with the upcoming NASA Atmosphere Observing System (AOS).

Travis Toth↗

Solubility and Diffusivity: Important Metrics in the Search for the Root Cause of Light- and Elevated Temperature-Induced Degradation

Light- and elevated temperature-induced degradation (LeTID) is a detrimental effect observed under operating conditions in p-type multicrystalline silicon (mc-Si) solar cells. In this paper, we employ synchrotron-based techniques to study the dissolution of precipitates due to different firing processes at grain boundaries in LeTID-affected mc-Si. The synchrotron measurements show clear dissolution of collocated metal precipitates during firing. We compare our observations with degradation behavior in the same wafers. The experimental results are complemented with process simulations to provide insight into the change in bulk point defect concentration due to firing. Several studies have proposed that LeTID is caused by metal-rich precipitate dissolution during contact firing, and we find that the solubility and diffusivity are promising screening metrics to identify metals that are compatible with this hypothesis. While slower and less soluble elements (e.g., Fe and Cr) are not compatible according to our simulations, the point defect concentrations of faster and more soluble elements (e.g., Cu and Ni) increase after a high-temperature firing process, primarily due to emitter segregation rather than precipitate dissolution. Finally, these results are a useful complement to lifetime spectroscopy techniques, and can be used to evaluate additional candidates in the search for the root cause of LeTID.

14 SOLAR ENERGY↗

Trajectory Optimization: OTIS 4

The latest release of the Optimal Trajectories by Implicit Simulation (OTIS4) allows users to simulate and optimize aerospace vehicle trajectories. With OTIS4, one can seamlessly generate optimal trajectories and parametric vehicle designs simultaneously. New features also allow OTIS4 to solve non-aerospace continuous time optimal control problems. The inputs and outputs of OTIS4 have been updated extensively from previous versions. Inputs now make use of objectoriented constructs, including one called a metastring. Metastrings use a greatly improved calculator and common nomenclature to reduce the user s workload. They allow for more flexibility in specifying vehicle physical models, boundary conditions, and path constraints. The OTIS4 calculator supports common mathematical functions, Boolean operations, and conditional statements. This allows users to define their own variables for use as outputs, constraints, or objective functions. The user-defined outputs can directly interface with other programs, such as spreadsheets, plotting packages, and visualization programs. Internally, OTIS4 has more explicit and implicit integration procedures, including high-order collocation methods, the pseudo-spectral method, and several variations of multiple shooting. Users may switch easily between the various methods. Several unique numerical techniques such as automated variable scaling and implicit integration grid refinement, support the integration methods. OTIS4 is also significantly more user friendly than previous versions. The installation process is nearly identical on various platforms, including Microsoft Windows, Apple OS X, and Linux operating systems. Cross-platform scripts also help make the execution of OTIS and post-processing of data easier. OTIS4 is supplied free by NASA and is subject to ITAR (International Traffic in Arms Regulations) restrictions. Users must have a Fortran compiler, and a Python interpreter is highly recommended.

Riehl, John P.↗

Biomass Burning Aerosol Absorption Measurements with MODIS Using the Critical Reflectance Method

This research uses the critical reflectance technique, a space-based remote sensing method, to measure the spatial distribution of aerosol absorption properties over land. Choosing two regions dominated by biomass burning aerosols, a series of sensitivity studies were undertaken to analyze the potential limitations of this method for the type of aerosol to be encountered in the selected study areas, and to show that the retrieved results are relatively insensitive to uncertainties in the assumptions used in the retrieval of smoke aerosol. The critical reflectance technique is then applied to Moderate Resolution Imaging Spectrometer (MODIS) data to retrieve the spectral aerosol single scattering albedo (SSA) in South African and South American 35 biomass burning events. The retrieved results were validated with collocated Aerosol Robotic Network (AERONET) retrievals. One standard deviation of mean MODIS retrievals match AERONET products to within 0.03, the magnitude of the AERONET uncertainty. The overlap of the two retrievals increases to 88%, allowing for measurement variance in the MODIS retrievals as well. The ensemble average of MODIS-derived SSA for the Amazon forest station is 0.92 at 670 nm, and 0.84-0.89 for the southern African savanna stations. The critical reflectance technique allows evaluation of the spatial variability of SSA, and shows that SSA in South America exhibits higher spatial variation than in South Africa. The accuracy of the retrieved aerosol SSA from MODIS data indicates that this product can help to better understand 44 how aerosols affect the regional and global climate.

Zhu, Li↗

Vibration attenuation of the NASA Langley evolutionary structure experiment using H(sub infinity) and structured singular value (micron) robust multivariable control techniques

The use is studied of active control to attenuate structural vibrations of the NASA Langley Phase Zero Evolutionary Structure due to external disturbance excitations. H sub infinity and structured singular value (mu) based control techniques are used to analyze and synthesize control laws for the NASA Langley Controls Structures Interaction (CSI) Evolutionary Model (CEM). The CEM structure experiment provides an excellent test bed to address control design issues for large space structures. Specifically, control design for structures with numerous lightly damped, coupled flexible modes, collocated and noncollocated sensors and actuators and stringent performance specifications. The performance objectives are to attenuate the vibration of the structure due to external disturbances, and minimize the actuator control force. The control design problem formulation for the CEM Structure uses a mathematical model developed with finite element techniques. A reduced order state space model for the control design is formulated from the finite element model. It is noted that there are significant variations between the design model and the experimentally derived transfer function data.

Balas, Gary J.↗

Improving Lidar Turbulence Estimates for Wind Energy

Remote sensing devices (e.g., lidars) are quickly becoming a cost-effective and reliable alternative to meteorological towers for wind energy applications. Although lidars can measure mean wind speeds accurately, these devices measure different values of turbulence intensity (TI) than an instrument on a tower. In response to these issues, a lidar TI error reduction model was recently developed for commercially available lidars. The TI error model first applies physics-based corrections to the lidar measurements, then uses machine-learning techniques to further reduce errors in lidar TI estimates. The model was tested at two sites in the Southern Plains where vertically profiling lidars were collocated with meteorological towers. This presentation primarily focuses on the physics-based corrections, which include corrections for instrument noise, volume averaging, and variance contamination. As different factors affect TI under different stability conditions, the combination of physical corrections applied in L-TERRA changes depending on the atmospheric stability during each 10-minute time period. This stability-dependent version of L-TERRA performed well at both sites, reducing TI error and bringing lidar TI estimates closer to estimates from instruments on towers. However, there is still scatter evident in the lidar TI estimates, indicating that there are physics that are not being captured in the current version of L-TERRA. Two options are discussed for modeling the remainder of the TI error physics in L-TERRA: machine learning and lidar simulations. Lidar simulations appear to be a better approach, as they can help improve understanding of atmospheric effects on TI error and do not require a large training data set.

17 WIND ENERGY↗

A harmonic balance method for PDEs

In this report, we describe some approaches to calculate the non-linear system of equations prescribed by the harmonic balance method (HB), a frequency domain analysis technique for modelling a non-linear system of partial differential equations (PDEs). The approach which we ultimately pursue can be seen as a time-collocation approach, except that the harmonic balance equations are obtained weakly (in the sense used in the calculus of variations). This weak formulation allows us to adapt existing transient or stationary PDEs models in the Panzer/Trilinos framework for frequency domain analysis via the harmonic balance method. We begin with a motivatation for the harmonic balance method and outline its mathematical formulation. We then describe some approaches to calculate the harmonic balance formulae, and their means of implementation through the modification of a Panzer tutorial problem - a stationary Helmholtz equation with a constant Dirichlet boundary condition and a non-linear source. For each of these approaches, we outline the necessary adaptations to solve the corresponding (periodically) transient Helmholtz equation with a (temporally) periodic Dirichlet boundary condition and non-linear source.

97 MATHEMATICS AND COMPUTING↗

GEOS-2 refraction program summary document

Data from an extensive array of collocated instrumentation at the Wallops Island test facility were intercompared in order to (1) determine the practical achievable accuracy limitations of various tropospheric and ionospheric correction techniques; (2) examine the theoretical bases and derivation of improved refraction correction techniques; and (3) estimate internal systematic and random error levels of the various tracking stations. The GEOS 2 satellite was used as the target vehicle. Data were obtained regarding the ionospheric and tropospheric propagation errors, the theoretical and data analysis of which was documented in some 30 separate reports over the last 6 years. An overview of project results is presented.

Mallinckrodt, A. J.↗

Incompressible Navier-Stokes and parabolized Navier-Stokes solution procedures and computational techniques

Recent developments with finite-difference techniques are emphasized. The quotation marks reflect the fact that any finite discretization procedure can be included in this category. Many so-called finite element collocation and galerkin methods can be reproduced by appropriate forms of the differential equations and discretization formulas. Many of the difficulties encountered in early Navier-Stokes calculations were inherent not only in the choice of the different equations (accuracy), but also in the method of solution or choice of algorithm (convergence and stability, in the manner in which the dependent variables or discretized equations are related (coupling), in the manner that boundary conditions are applied, in the manner that the coordinate mesh is specified (grid generation), and finally, in recognizing that for many high Reynolds number flows not all contributions to the Navier-Stokes equations are necessarily of equal importance (parabolization, preferred direction, pressure interaction, asymptotic and mathematical character). It is these elements that are reviewed. Several Navier-Stokes and parabolized Navier-Stokes formulations are also presented.

Rubin, S. G.↗