Efficient analysis and storage of data is an integral but often challenging task when working with computation fluid dynamics mainly due to the amount of data it can output. Methods centered around the proper orthogonal decomposition were used to analyze, compress, and model various simulation cases. Two different high-fidelity, time-accurate turbomachinery simulations were investigated to show various applications of the analysis techniques. The first turbomachinery example was used to illustrate the extraction of turbulent coherent structures such as traversing shocks, vortex shedding, and wake variation from deswirler and rotor blade passages. Using only the most dominant modes, flow fields were reconstructed and analyzed for error. The reconstructions reproduced the general dynamics within the flow well, but failed to fully resolve shock fronts and smaller vortices. By decomposing the domain into smaller, independent pieces, reconstruction error was reduced by up to 63 percent. A new method of data compression that combined an image compression algorithm and the proper orthogonal decomposition was used to store the reconstructions of the flow field, increasing data compression ratios by a factor of 40.The second turbomachinery simulation studied was a three-stage fan with inlet total pressure distortion. Both the snapshot and repeating geometry methods were used to characterize structures of static pressure fluctuation within the blade passages of the third rotor blade row. Modal coefficients filtered by frequencies relating to the inlet distortion pattern were used to produce reconstructions of the pressure field solely dependent on the inlet boundary condition. A hybrid proper orthogonal decomposition method was proposed to limit burdens on computational resources while providing high temporal resolution analysis.Parametric reduced order models were created from large databases of transient and steady conjugate heat transfer and airfoil simulations. Performance of the models were found to depend heavily on the range of the parameters varied as well as the number of simulations used to traverse that range. The heat transfer models gave excellent predictions for temperature profiles in heated solids for ambitious parameter ranges. Model development for the airfoil case showed that accuracy was highly dependent on modal truncation. The flow fields were predicted very well, especially outside the boundary layer region of the flow.

This thesis describes a novel new method, termed the Tangent-Line-Chord (TLC) method, that can be used to more efficiently model creep deformation dominated by the tertiary regime. Creep deformation is a widespread mechanical mode of failure found in high-stress and temperature mechanical systems. To accurately simulate creep and its effect on structures, researchers utilize finite element analysis (FEA). General purpose FEA packages require extensive amounts of time and computer resources to simulate creep softening in components because of the large deformation rates that continuously evolve. The goal of this research is to employ multi-regime creep models, such as the Kachanov-Rabotnov model, to determine a set of equations that will allow creep to be simulated using as few iterations as possible. The key outcome is the freeing up of computational resources and the saving of time. Because both the number of equations and the value of material constants within the model change depending on the approach used, programming software will be utilized to automate this analytical process. The materials being considered in this research are mainly generic Ni-based superalloys, as they exhibit creep responses that are dominated by secondary and tertiary creep.

O estudo da dinâmica estrutural de risers oceânicos apresenta instigantes desafios aos pesquisadores da área da engenharia de estruturas, uma vez que os meios tradicionais de análises dinâmicas lineares nem sempre se ajustam às suas complexas particularidades. No atual estágio do desenvolvimento científico da área de engenharia de estruturas, a aplicação de técnicas de análise dinâmica não linear, dentro de determinadas hipóteses, mostra-se como uma das alternativas possíveis e viáveis à tradicional análise dinâmica linear. Com vistas a uma nova abordagem do problema, o presente trabalho adota uma metodologia de análise não linear dinâmica de risers oceânicos em configuração de lançamento de catenária, conjugada a uma técnica de processamento de Modelos de Ordem Reduzida para o estudo dos fenômenos dinâmicos manifestados por risers. Trata-se de um método de modelagem local, restrito à região de contato unilateral do riser com o solo, considerado este último um meio viscoelástico. Os resultados da aplicação desta metodologia são demonstrados nos estudos de caso apresentados com comparações com modelos numéricos (Método dos Elementos Finitos) e modelos físicos. / The dynamic study of offshore risers still demands large efforts from structural engineering researchers, since these systems may behave in a way that is not well modeled and understood using simply linear dynamic theories. Nevertheless, the current development stage of non linear dynamic theories gives hope that their use for the analyses of such systems can be of great value, even though, this must be carefully done specially by the analyst. The present work refers to a non linear dynamic methodology application to offshore risers, particularly steel catenary risers, by a technique known as reduced-order modeling, in the study of dynamic phenomena that these structures may present. The model is local, which means that it represents the touch-down zone of the riser-soil system. The soil modeling was presumed to be viscoelastic. The results obtained in case studies are compared with those from numerical (Finite Element Method) and small scale physical models.

Risers são tubos que transportam fluidos do fundo do mar até as plataformas flutuantes e vice-versa. Diversas configurações e materiais são utilizados, porém apenas os steel catenary risers (SCR) são estudados neste trabalho. Os risers são estruturas extremamente esbeltas e, por isso, grande parte de seu trecho suspenso tem comportamento de cabo. Apenas em duas regiões a rigidez flexional é relevante: no hang-off (topo) e na touch-down zone (TDZ), sendo esta última a região mais complexa para análise devido ao contato unilateral com o solo. Em função dos diversos carregamentos dinâmicos a que o riser é submetido, grandes variações na curvatura ocorrem na TDZ, além de impacto e atrito com o solo, que podem reduzir a vida útil da estrutura ou até mesmo por em risco a sua integridade. Por estas razões, este trabalho visa à elaboração de uma metodologia analítica para a construção de um modelo de ordem reduzida (MOR) capaz de analisar o comportamento dinâmico não linear da TDZ de um SCR. Como na TDZ a rigidez flexional predomina sobre a rigidez geométrica, o riser é modelado como uma viga semi-infinita, tendo uma parte suspensa e outra apoiada sobre solo hipoteticamente elástico com contato unilateral. Na extremidade suspensa são impostos deslocamentos verticais e trações dinâmicas que fazem com que a posição do touch-down point (TDP) também varie com o tempo. Trata-se, portanto, de um problema com condições de contorno móveis. A metodologia adotada para a resolução deste problema foi transformá-lo em um problema de condições de contorno fixas por meio de uma transformação de variáveis. Contudo, paga-se um preço por tal transformação, e fortes não linearidades surgem na equação diferencial de movimento, tornando-a extremamente complexa para uma resolução analítica direta. Para o problema de flexão simples, consegue-se obter os modos de vibração não lineares através do método das múltiplas escalas. De posse destes modos, utiliza-se o método de Galerkin não linear para projetar a equação completa em um modo escolhido, transformando o modelo contínuo em um modelo de ordem reduzida com apenas um grau de liberdade, cuja coordenada generalizada modal é o deslocamento horizontal do TDP. Obtida a equação do MOR, nota-se que existem coeficientes que variam com o tempo, como na clássica equação de Mathieu, indicando a possibilidade de ocorrer ressonância paramétrica. Neste tipo de ressonância, entre outras possibilidades, pode ocorrer que a frequência de excitação seja o dobro da frequência natural trata-se da ressonância paramétrica principal. A equação do MOR é integrada numericamente e suas respostas são comparadas com as respostas obtidas por modelos de elementos finitos elaborados em softwares comerciais, como o Abaqus e o Orcaflex. Por fim, discutem-se as potencialidades e limitações do MOR, sendo uma grande vantagem a possibilidade de processar diversos casos facilmente, variando os parâmetros que influem nas respostas. Com este mapeamento das respostas, é possível estimar as amplitudes dos estados estacionários pós-críticos. / Risers are pipes that convey fluids from the seabed up to the floating platforms and vice-versa. Many configurations and materials are used, but only steel catenary risers (SCR) are studied in this work. Risers are extremely slender structures, and for this reason, most of the suspended part has cable behavior. Only in two regions the bending stiffness is important: at the hang-off and at the touch-down zone (TDZ), which is the most complex region for analysis because of the unilateral contact with the seabed. Due to several dynamic loads that the riser is subjected to, great curvature variations occur at the TDZ, apart from impacts and friction with the soil, which can reduce the life time of the structure or even jeopardize its integrity. For these reasons, this work aims at the development of an analytical methodology for the construction of a reduced-order model (ROM) able to analyze the nonlinear dynamic behavior of the TDZ of a SCR. As at the TDZ the bending stiffness prevails over the geometric stiffness, the riser is modeled as a semi-infinite beam, having a suspended part and another one resting on the elastic soil with unilateral contact. At the end of the suspended part, vertical displacements and dynamic tensions are imposed, that cause the TDPs position to vary with time. It is, therefore, a problem with moving boundary conditions. The methodology adopted for solving this problem was to transform it into a problem with fixed boundary conditions via a variable transformation. However, a price is paid for such a transformation, and strong nonlinearities appear in the differential equation of motion, making it extremely complex to solve analytically. For the simple bending problem, nonlinear vibration modes are obtained via the method of multiple scales. In possession of these modes, the nonlinear Galerkin method is used to project the complete equation into a chosen mode, transforming the continuum model into a reduced-order model (ROM) with only one degree of freedom whose modal generalized coordinate is the horizontal displacement of the TDP. After obtaining the ROM, it is noticed that there are coefficients that vary with time, as in the classic Mathieu equation, indicating the possibility of parametric resonance. In this kind of resonance, among other possibilities, the excitation frequency may be twice the natural frequency it is the so-called principal parametric resonance. The ROMs equation is integrated numerically and the responses are compared to those given by finite-element models studied with the help of commercial softwares, like Abaqus and Orcaflex. Finally, the potentialities and limitations of the ROM are discussed. One of the advantages is the possibility of processing several cases easily, changing the parameters that affect the responses. With this response mapping, it is possible to estimate the post-critical steady-state amplitudes that take place.

The use of Computational Fluid Dynamics (CFD) tools is crucial for the development of novel and cost-effective combustion technologies and the minimization of environmental concerns at industrial scale. CFD simulations facilitate scaling-up procedures that otherwise would be complicated by strong interactions between reaction kinetics, turbulence and heat transfer. CFD calculations can be applied directly at the industrial scale of interest, thus avoiding scaling-up from lab-scale experiments. However, this advantage can only be obtained if CFD tools are quantitatively predictive and trusted as so. Despite the improvements in the computational capability, the implementation of detailed physical and chemical models in CFD simulations can still be prohibitive for real combustors, which require large computational grids and therefore significant computational efforts. Advanced simulation approaches like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) guarantee higher fidelity in computational modeling of combustion at, unfortunately, increased computational cost. However, with adequate, reduced, and cost-effective modeling of physical phenomena, such as chemical kinetics and turbulence-chemistry interactions, and state of the art computing, LES will be the tool of choice to describe combustion processes at industrial scale accurately. Therefore, the development of reduced physics and chemistry models with quantified model-form uncertainty is needed to overcome the challenges of performing LES of industrial systems. Reduced-order models must reproduce the main features of the corresponding detailed models. They feature predictivity and capability of bridging scales when validated against a broad range of experiments and targeted by Validation and Uncertainty Quantification (V/UQ) procedures. In this work, V/UQ approaches are applied for reduced-order modeling of pulverized coal devolatilization and subsequent char oxidation, and furthermore for modeling NOx emissions in combustion systems.For coal devolatilization, a benchmark of the Single First-Order Reaction (SFOR) model was performed concerning the accuracy of the prediction of volatile yield. Different SFOR models were implemented and validated against experimental data coming from tests performed in an entrained flow reactor at oxy-conditions, to shed light on their drawbacks and benefits. SFOR models were chosen because of their simplicity: they can be easily included in CFD codes and are very appealing in the perspective of LES of pulverized coal combustion burners. The calibration of kinetic parameters was required to allow the investigated SFOR model to be predictive and reliable for different heating rates, hold temperatures and coal types. A comparison of several calibration approaches was performed to determine if one-step models can be adaptive and able to bridge scales, without losing accuracy, and to select the calibration method to employ for wider ranges of coal rank and operating conditions. The analysis pointed out that the main drawback of the SFOR models is the assumption of a constant ultimate volatile yield, equal to the value from the coal proximate analysis. To overcome this drawback, a yield model, i.e. a simple functional form that relates the ultimate volatile yield to the particle temperature, was proposed. The model depends on two parameters that have a certain degree of uncertainty. The performances of the yield model were assessed using a collaboration of experiments and simulations of a pilot-scale entrained flow reactor. A consistency analysis, based on the Bound-to-Bound Data Collaboration (B2B-DC) approach, and a Bayesian method, based on Gaussian Process Regression (GPR), were employed for the investigation of experiments and simulations. In Bound-to- Bound Data Collaboration the model output, evaluated at specified values of the model parameters, is compared with the experimental data: if the prediction of the model falls within the experimental uncertainty, the corresponding parameter values would be included in the so-called feasible set. The existence of a non-empty feasible set signifies consistency between the experiments and the simulations, i.e. model-data agreement. Consistency was indeed found when a relative error of 19% for all the experimental data was applied. Hence, a feasible set of the two SFOR model parameters was provided. A posterior state of knowledge, indicating potential model forms that could be explored in yield modeling, was obtained by Gaussian Process Regression. The model form evaluated through the consistency analysis is included within the posterior derived from GPR, indicating that it can satisfactorily match the experimental data and provide reliable estimation in almost every range of temperatures. CFD simulations were carried out using the proposed yield model with first-order kinetics, as in the SFOR model. Results showed promising agreement between predicted and experimental conversion for all the investigated cases.Regarding char combustion modeling, the consistency analysis has been applied to validate a reduced-order model and quantify the uncertainty in the prediction of char conversion. The model capability to address heterogeneous reaction between char carbon and O2, CO2 and H2O reagents, mass transport of species in the particle boundary layer, pore diffusion, and internal surface area changes was assessed by comparison with a large number of experiments performed in air and oxy-coal conditions. Different model forms had been considered, with an increasing degree of complexity, until consistency between model outputs and experimental results was reached. Rather than performing forward propagation of the model-form uncertainty on the predictions, the reduction of the parameter uncertainty of a selected model form was pursued and eventually achieved. The resulting 11-dimensional feasible set of model parameters allows the model to predict the experimental data within almost ±10% uncertainty. Due to the high dimensionality of the problem, the employed surrogate models resulted in considerable fitting errors, which led to a spoiled UQ inverse problem. Different strategies were taken to reduce the discrepancy between the surrogate outputs and the corresponding predictions of the simulation model, in the frameworks of constrained optimization and Bayesian inference. Both strategies succeeded in reducing the fitting errors and also resulted in a least-squares estimate for the simulation model. The variety of experimental gas environments ensured the validity of the consistent reduced model for both conventional and oxy-conditions, overcoming the differences in mass transport and kinetics observed in several experimental campaigns.The V/UQ-aided modeling of coal devolatilization and char combustion was done in the framework of the Predictive Science Academic Alliance Program II (PSAAP-II) funded by the US Department of Energy. One of the final goals of PSAAP-II is to develop high-fidelity simulation tools that ensure 5% uncertainty in the incident heat flux predictions inside a 1.2GW Ultra-Super-Critical (USC) coal-fired boiler. The 5% target refers to the expected predictivity of the full-scale simulation without considering the uncertainty in the scenario parameters. The data-driven approaches used in this Thesis helped to improve the predictivity of the investigated models and made them suitable for LES of the 1.2GW USC coal-fired boiler. Moreover, they are suitable for scale-bridging modeling of similar multi-phase processes involved in the conversion of solid renewable sources, such as biomass.In the final part of the Thesis, the sensitivity to finite-rate chemistry combustion models and kinetic mechanisms on the prediction of NO emissions was assessed. Moreover, the forward propagation of the uncertainty in the kinetics of the NNH route (included in the NOx chemistry) on the predictions of NO was investigated to reveal the current state of the art of kinetic modeling of NOx formation. The analysis was carried out on a case where NOx formation comes from various formation routes, both conventional (thermal and prompt) and unconventional ones. To this end, a lab-scale combustion system working in Moderate and Intense Low-oxygen Dilution (MILD) conditions was selected. The results showed considerable sensitivity of the NO emissions to the uncertain kinetic parameters of the rate-limiting reactions of the NNH pathway when a detailed kinetic mechanism is used. The analysis also pointed out that the use of one-step global rate schemes for the NO formation pathways, necessary when a skeletal kinetic mechanism is employed, lacks the required chemical accuracy and dims the importance of the NNH pathway in this combustion regime. An engineering modification of the finite-rate combustion model was proposed to account for the different chemical time scales of the fuel-oxidizer reactions and NOx formation pathways. It showed an equivalent impact on the emissions of NO than the uncertainty in the kinetics of the NNH route. At the cost of introducing a small mass imbalance (of the order of ppm), the adjustment led to improved predictions of NO. The investigation established a possibility for the engineering modeling of NO formation in MILD combustion with a finite-rate chemistry combustion model that can incorporate a detailed mechanism at affordable computational costs. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

Assessing the impact of inlet flow distortion in turbomachinery is desired early in the design cycle. This thesis introduces and validates the use of methods based on the Proper Orthogonal Decomposition (POD) to analyze clean and 1/rev static pressure distortion simulation results at design and near stall operating condition. The value of POD comes in its ability to efficiently extract both quantitative and qualitative information about dominant spatial flow structures as well as information about temporal fluctuations in flow properties. Observation of the modes allowed qualitative identification of shock waves as well as quantification of their location and range of motion. Modal coefficients revealed the location of the passage shock at a given angular location. Distortion amplification and attenuation between rotors was also identified. A relationship was identified between how distortion manifests itself based on downstream conditions. POD provides an efficient means for extracting the most meaningful information from large CFD simulation data. Static pressure and axial velocity were analyzed to explore the flow physics of 3 rotors of a compressor with a distorted inlet. Based on the results of the analysis of static pressure using the POD modes, it was concluded that there was a decreased range of motion in passage shock oscillation. Analysis of axial velocity POD modes revealed the presence of a separated region on the low pressure surface of the blade which was most dynamic in rotor 1. The thickness of this structure decreased in the near stall operating condition. The general conclusion is made that as the fan approaches stall the apparent effects of distortion are lessened which leads to less variation in the operating condition. This is due to the change in operating condition placing the fan at a different position on the speedline such that distortion effects are less pronounced. POD modes of entropy flux were used to identify three distinct levels of entropy flux in the blade row passage. The separated region was the region with the highest entropy due to the irreversibilities associated with separation.

In this paper we describe an approach to generating low-order models of spiral inductors that accurately capture the dependence on both frequency and geometry (width and spacing) parameters. The approach is based on adapting a multiparameter Krylov-subspace based moment matching method to reducing an integral equation for the three dimensional electromagnetic behavior of the spiral inductor. The approach is demonstrated on a typical on-chip rectangular inductor. / Singapore-MIT Alliance (SMA)

Data centers are computing infrastructure facilities used by industries with large data processing needs and the rapid increase in power density of high performance computing equipment has caused many thermal issues in these facilities. Systems-level thermal management requires modeling and analysis of complex fluid flow and heat transfer processes across several decades of length scales. Conventional computational fluid dynamics and heat transfer techniques for such systems are severely limited as a design tool because their large model sizes render parameter sensitivity studies and optimization impractically slow. The traditional proper orthogonal decomposition (POD) methodology has been reformulated to construct physics-based models of turbulent flows and forced convection. Orthogonal complement POD subspaces were developed to parametrize inhomogeneous boundary conditions and greatly extend the use of the existing POD methodology beyond prototypical flows with fixed parameters. A flux matching procedure was devised to overcome the limitations of Galerkin projection methods for the Reynolds-averaged Navier-Stokes equations and greatly improve the computational efficiency of the approximate solutions. An implicit coupling procedure was developed to link the temperature and velocity fields and further extend the low-dimensional modeling methodology to conjugate forced convection heat transfer. The overall reduced-order modeling framework was able to reduce numerical models containing 105 degrees of freedom (DOF) down to less than 20 DOF, while still retaining greater that 90% accuracy over the domain. Rigorous a posteriori error bounds were formulated by using the POD subspace to partition the error contributions and dual residual methods were used to show that the flux matching procedure is a computationally superior approach for low-dimensional modeling of steady turbulent convection. To efficiently model large-scale systems, individual reduced-order models were coupled using flow network modeling as the component interconnection procedure. The development of handshaking procedures between low-dimensional component models lays the foundation to quickly analyze and optimize the modular systems encountered in electronics thermal management. This modularized approach can also serve as skeletal structure to allow the efficient integration of highly-specialized models across disciplines and significantly advance simulation-based design.

In der vorliegenden Arbeit wird eine Prozedur zur Ordnungsreduktion von Finite Elemente Modellen mikromechanischer Struktur mit elektrostatischem Wirkprinzip entwickelt und analysiert. Hintergrund der Ordnungsreduktion ist eine Koordinatentransformation von lokalen Finite Elemente Koordinaten in globale Koordinaten. Die globalen Koordinaten des reduzierten Modells werden durch einige wenige Formfunktionen beschrieben. Damit wird das Makromodell nicht mehr durch lokale Knotenverschiebungen beschrieben, sondern durch globale Formfunktionen, welche die gesamte Deformation der Struktur beeinflussen. Es wird gezeigt, dass Eigenvektoren der linearisierten mechanischen Struktur einfache und effiziente Formfunktionen darstellen. Weiterhin kann diese Methode für bestimmte Nichtlinearitäten und für verschiedene in Mikrosystemen auftretende Lasten angewendet werden. Das Ergebnis sind Makromodelle, die über Klemmen in Systemsimulatoren eingebunden werden können, die Genauigkeiten einer Finite Elemente Analyse erreichen und für Systemsimulationen typische Laufzeitverhalten besitzen.

One of the challenges for Airbus preparing a new freighter development process was the design of a solid freighter barrier, which separates the courier area from the cargo compartment. The major task of such a barrier is to protect the passengers against all risks caused due to cargo impact by a justifiable design. These risks may result from all kind of survivable incident and accident scenarios. Real aircraft crashes were analyzed to get away from a static book-case and come to a more realistic dynamic crash scenario. A reduced-order simulation model was built up to investigate and simulate the dynamic effects during crash. The simulation model considers the highly nonlinear stiffness and damping characteristics of all critical cargo types and also includes their energy absorption potentials. A series of full scale container crash tests have been performed at accredited car crash facilities. The test campaigns were complemented by numerous component tests to study also general crash principles. The critical simulation parameters were identified and implemented into the simulation model. The subsequent validation process showed a close agreement between simulation and test. The simulation environment has turned out to be a reliable basis to simulate all critical barrier loads with respect to the specific aircraft loading distributions. The essence of this investigation is an adequate understanding of the real crash effects. The proposed dynamic crash approach is more realistic than the static condition and results in an optimized safety barrier wall concept. This dynamic approach provides equivalent safety compared to the existing devices and is accepted by FAA and EASA.