SEARCH · Engineering Papers

Results for “structured model reduction”

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 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Cathodoluminescence as an effective probe of carrier transport and deep level defects in droop-mitigating InGaN/GaN quantum well heterostructures

Commercial InGaN/GaN light emitting diodes continue to suffer from efficiency droop at high current densities, and urgently require enhanced structural-optical toolsets for active region characterization. We measure delayed (tens of seconds) cathodoluminescence which is influenced by carrier transport and deep level defects. Further, we observe that drops in efficiency are not correlated with quantum well (QW) width fluctuations. To explain the rise dynamics, we propose a model involving filling of deep level defects and simultaneous reduction of built-in field within the multi-QW active region. These measurements yield insights into carrier transport, efficiency-reducing defects, and QW band structure.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Isotropic in-plane quenched disorder and dilution induce a robust nematic state in electron-doped pnictides

The phase diagram of electron-doped pnictides is studied varying the temperature, electronic density, and isotropic in-plane quenched disorder strength and dilution by means of computational techniques applied to a three-orbital (xz,yz,xy) spin-fermion model with lattice degrees of freedom. In experiments, chemical doping introduces disorder but in theoretical studies the relationship between electronic doping and the randomly located dopants, with their associated quenched disorder, is difficult to address. Moreover, in this publication, the use of computational techniques allows us to study independently the effects of electronic doping, regulated by a global chemical potential, and impurity disorder at randomly selected sites. Surprisingly, our Monte Carlo simulations reveal that the fast reduction with doping of the N eel T N and the structural T S transition temperatures, and the concomitant stabilization of a robust nematic state, is primarily controlled in our model by the magnetic dilution associated with the in-plane isotropic disorder introduced by Fe substitution. In the doping range studied, changes in the Fermi surface produced by electron doping affect only slightly both critical temperatures. Our results also suggest that the specific material-dependent phase diagrams experimentally observed could be explained as a consequence of the variation in disorder profiles introduced by the different dopants. Finally, our findings are also compatible with neutron scattering and scanning tunneling microscopy, unveiling a patchy network of locally magnetically ordered clusters with anisotropic shapes, even though the quenched disorder is locally isotropic. Our study reveals a remarkable and unexpected degree of complexity in pnictides: the fragile tendency to nematicity intrinsic of translational invariant electronic systems needs to be supplemented by quenched disorder and dilution to stabilize the robust nematic phase experimentally found in electron-doped 122 compounds.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Rotary-wing aeroelasticity with application to VTOL vehicles

This survey presents an assessment of the state of the art in rotary-wing aeroelasticity as applied to conventional helicopters as well as other VTOL vehicles such as tilting prop-rotors, the X-wing and a hybrid heavy lift vehicle. The objective is to enable the reader to develop an awareness of what has been accomplished, what remains to be done, and where to find more comprehensive treatments of the various topics discussed. The main topics discussed are: (1) structural modeling; (2) unsteady aerodynamic modeling; (3) formulation of the equations of motion and their solutions; (4) illustrative results for isolated blades in hover and forward flight; (5) illustrative results for coupled rotor/fuselage problems; (6) active control of aeromechanical and aeroelastic problems; (7) active controls for vibration reduction; (8) structural optimization with aeroelastic constraints; (9) gust response analysis of rotors; and (10) aeroelastic problems in special VTOL vehicles. These topics are reviewed with different levels of detail and some useful observation on potentially rewarding areas of future research are made.

Friedmann, Peretz P.↗

NASA NTRS1999

Experiments on Diffusion Flame Structure of a Laminar Vortex Ring

The study of flame-vortex interactions provides one of the means to better understand turbulent combustion, and allows for canonical configurations that contain the fundamental elements found in turbulent flames, These include concentrated vorticity, entrainment and mixing, strain and nonequilibrium phenomena, diffusion and differential diffusion, partial premixing and diluent effects, and heat release effects. In flame- vortex configurations, these fundamental elements can be studied under more controlled conditions than is possible in direct investigations of turbulent flames. Since the paper of Marble, the problem of the flame-vortex interaction has received considerable attention theoretically, numerically and experimentally. Several configurations exist for study of the premixed flame/vortex ring interaction but more limited results have been obtained to date for the diffusion flame/vortex ring case. The setup of Chen and Dahm, which is conceptually similar to that of Karagozian and Manda and Karagozian, Suganuma and Strom where the ring is composed of fuel and air and combustion begins during the ring formation process, is used in the current study. However, it is essential to conduct the experiments in microgravity to remove the asymmetries caused by buoyancy and thus obtain highly symmetric and repeatable interactions. In previous studies it was found that the flame structure of the vortex ring was similar to that obtained analytically by Karagozian and Manda. Dilution of propane with nitrogen led mainly to a reduction in flame luminosities, flame burnout times were affected by both fuel volumes and amount of dilution, and a simple model of the burnout times was developed. In this paper, a discussion on reacting ring displacement and flame burnout time will be given, and the flame structures of vortex rings containing ethane and air will be compared to those of propane reacting in air.

Chen, Shin-Juh↗

DOE OSTI2018

Reductive Cleavage Method for Quantitation of Monolignols and Low‐Abundance Monolignol Conjugates

Abstract As interest in biomass utilization has grown, the manipulation of lignin biosynthesis has received significant attention, such that recent work has demanded more robust lignin analytical methods. As the derivatization followed by reductive cleavage (DFRC) method is particularly effective for structurally characterizing natively acylated lignins, we used an array of synthetic β‐ether γ‐acylated model compounds to determine theoretical yields for all monolignol conjugates currently known to exist in lignin, and we synthesized a new set of deuterated analogs as internal standards for quantification using GC–MS/MS. Yields of the saturated ester conjugates ranged from 40 to 90 %, and NMR analysis revealed the presence of residual unsaturated conjugates in yields of 20 to 35 %. In contrast to traditional selected‐ion‐monitoring, we demonstrated the superior sensitivity and accuracy of multiple‐reaction‐monitoring detection methods, and further highlighted the inadequacy of traditional standards relative to isotopically labeled analogs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1980

Energy Efficient High-Pressure Turbine Leakage Technology Report

The leakage test program was one of such supporting technology programs structured to provide guidance to the Energy Efficient Engine High Pressure Turbine Component Design Effort. Leakage reduction techniques were identified and evaluated. Test models were used to simulate component leak paths and to evaluate leakage reduction techniques. These models simulated the blade/disk attachment, the vane inner platform attachment, and the vane outer platform attachment combined with the blade outer airseal. Disk blade attachment testing indicated that leakage in this area could be reduced to very low levels by paying careful attention to the tolerances along the contact surface between the blade vibration damper and the blade platform contact surface. The aim of feather seal testing was to achieve a goal for an effective leakage gap of one mil (.001 inch) per inch of feather seal length. Results indicated that effective gaps even below the goal level were achievable by (1) maintaining close tolerances between feather seals and their slots to minimize end gaps and limit seal rotation, (2) avoiding feather seal overlap, and (3) minimizing feather seal intersections. W seals were shown to be effective leakage control devices. Wire rope, in its present state of development, was shown not to be an effective sealing concept for application to the component design.

Gardner, W. B.↗

NASA NTRS1985

Fundamental studies of structure borne noise for advanced turboprop applications

The transmission of sound generated by wing-mounted, advanced turboprop engines into the cabin interior via structural paths is considered. The structural model employed is a beam representation of the wing box carried into the fuselage via a representative frame type of carry through structure. The structure for the cabin cavity is a stiffened shell of rectangular or cylindrical geometry. The structure is modelled using a finite element formulation and the acoustic cavity is modelled using an analytical representation appropriate for the geometry. The structural and acoustic models are coupled by the use of hard wall cavity modes for the interior and vacuum structural modes for the shell. The coupling is accomplished using a combination of analytical and finite element models. The advantage is the substantial reduction in dimensionality achieved by modelling the interior analytically. The mathematical model for the interior noise problem is demonstrated with a simple plate/cavity system which has all of the features of the fuselage interior noise problem.

Eversman, W.↗

NASA NTRS1989

A continuum model for dynamic analysis of the Space Station

Dynamic analysis of the International Space Station using MSC/NASTRAN had 1312 rod elements, 62 beam elements, 489 nodes and 1473 dynamic degrees of freedom. A realtime, man-in-the-loop simulation of such a model is impractical. This paper discusses the mathematical model for realtime dynamic simulation of the Space Station. Several key questions in structures and structural dynamics are addressed. First, to achieve a significant reduction in the number of dynamic degrees of freedom, a continuum equivalent representation of the Space Station truss structure which accounted for the unsymmetry of the basic configuration and resulted in the coupling of extensional and transverse deformation, is developed. Next, dynamic equations for the continuum equivalent of the Space Station truss structure are formulated using a matrix version of Kane's dynamical equations. Flexibility is accounted for by using a theory that accommodates extension, bending in two principal planes and shear displacement. Finally, constraint equations suitable for dynamic analysis of flexible bodies with closed loop configuration are developed and solution of the resulting system of equations is based on the zero eigenvalue theorem.

Thomas, Segun↗

NASA NTRS2016

Probabilistic Prognosis of Non-Planar Fatigue Crack Growth

Quantifying the uncertainty in model parameters for the purpose of damage prognosis can be accomplished utilizing Bayesian inference and damage diagnosis data from sources such as non-destructive evaluation or structural health monitoring. The number of samples required to solve the Bayesian inverse problem through common sampling techniques (e.g., Markov chain Monte Carlo) renders high-fidelity finite element-based damage growth models unusable due to prohibitive computation times. However, these types of models are often the only option when attempting to model complex damage growth in real-world structures. Here, a recently developed high-fidelity crack growth model is used which, when compared to finite element-based modeling, has demonstrated reductions in computation times of three orders of magnitude through the use of surrogate models and machine learning. The model is flexible in that only the expensive computation of the crack driving forces is replaced by the surrogate models, leaving the remaining parameters accessible for uncertainty quantification. A probabilistic prognosis framework incorporating this model is developed and demonstrated for non-planar crack growth in a modified, edge-notched, aluminum tensile specimen. Predictions of remaining useful life are made over time for five updates of the damage diagnosis data, and prognostic metrics are utilized to evaluate the performance of the prognostic framework. Challenges specific to the probabilistic prognosis of non-planar fatigue crack growth are highlighted and discussed in the context of the experimental results.

Leser, Patrick E.↗

NASA NTRS2015

A Vehicle Management End-to-End Testing and Analysis Platform for Validation of Mission and Fault Management Algorithms to Reduce Risk for NASA's Space Launch System

The development of the Space Launch System (SLS) launch vehicle requires cross discipline teams with extensive knowledge of launch vehicle subsystems, information theory, and autonomous algorithms dealing with all operations from pre-launch through on orbit operations. The characteristics of these systems must be matched with the autonomous algorithm monitoring and mitigation capabilities for accurate control and response to abnormal conditions throughout all vehicle mission flight phases, including precipitating safing actions and crew aborts. This presents a large complex systems engineering challenge being addressed in part by focusing on the specific subsystems handling of off-nominal mission and fault tolerance. Using traditional model based system and software engineering design principles from the Unified Modeling Language (UML), the Mission and Fault Management (M&FM) algorithms are crafted and vetted in specialized Integrated Development Teams composed of multiple development disciplines. NASA also has formed an M&FM team for addressing fault management early in the development lifecycle. This team has developed a dedicated Vehicle Management End-to-End Testbed (VMET) that integrates specific M&FM algorithms, specialized nominal and off-nominal test cases, and vendor-supplied physics-based launch vehicle subsystem models. The flexibility of VMET enables thorough testing of the M&FM algorithms by providing configurable suites of both nominal and off-nominal test cases to validate the algorithms utilizing actual subsystem models. The intent is to validate the algorithms and substantiate them with performance baselines for each of the vehicle subsystems in an independent platform exterior to flight software test processes. In any software development process there is inherent risk in the interpretation and implementation of concepts into software through requirements and test processes. Risk reduction is addressed by working with other organizations such as S&MA, Structures and Environments, GNC, Orion, the Crew Office, Flight Operations, and Ground Operations by assessing performance of the M&FM algorithms in terms of their ability to reduce Loss of Mission and Loss of Crew probabilities. In addition, through state machine and diagnostic modeling, analysis efforts investigate a broader suite of failure effects and detection and responses that can be tested in VMET and confirm that responses do not create additional risks or cause undesired states through interactive dynamic effects with other algorithms and systems. VMET further contributes to risk reduction by prototyping and exercising the M&FM algorithms early in their implementation and without any inherent hindrances such as meeting FSW processor scheduling constraints due to their target platform - ARINC 653 partitioned OS, resource limitations, and other factors related to integration with other subsystems not directly involved with M&FM. The plan for VMET encompasses testing the original M&FM algorithms coded in the same C++ language and state machine architectural concepts as that used by Flight Software. This enables the development of performance standards and test cases to characterize the M&FM algorithms and sets a benchmark from which to measure the effectiveness of M&FM algorithms performance in the FSW development and test processes. This paper is outlined in a systematic fashion analogous to a lifecycle process flow for engineering development of algorithms into software and testing. Section I describes the NASA SLS M&FM context, presenting the current infrastructure, leading principles, methods, and participants. Section II defines the testing philosophy of the M&FM algorithms as related to VMET followed by section III, which presents the modeling methods of the algorithms to be tested and validated in VMET. Its details are then further presented in section IV followed by Section V presenting integration, test status, and state analysis. Finally, section VI addresses the summary and forward directions followed by the appendices presenting relevant information on terminology and documentation.

Trevino, Luis↗

NASA NTRS1983

Acoustic responses of coupled fluid-structure system by acoustic-structural analogy

The use of an analogy between structural mechanics and acoustics makes it possible to solve fluid-structural interaction (FSI) problems using an existing structural analysis computer program. This method was implemented in MSC/NASTRAN program and the FSI analysis was performed using two dimensional coupled fluid beam model to assess and evaluate the adequacy of this approach. The coupled modal analysis of 3-D model is also briefly discussed. The normal mode, modal frequency response and transient response analysis of 2-D coupled fluid beam system is presented. The significant reduction of the acoustic pressure response at the fluid structure interface is observed as a result of fluid structure interaction.

Shin, Y. S.↗

NASA NTRS2019

Simulated Weightlessness Alters Cardiomyocyte Structure and Transcriptional Regulation of Mediators Related to Immunity and Cardiovascular Disease

Spaceflight and the ensuing fluid shifts, together with an overall reduction in physical activity, lead to acute and latent effects on the cardiovascular system. This current study makes use of the rodent hindlimb unloading (HU) model to determine how factors such as sex, age, and duration of exposure impact cardiac responses to weightlessness. We hypothesize that extended exposure to simulated weightlessness and the ensuing recovery alters cardiac structure and expression of select genes, including those involved in redox signaling which together, negatively impact long-term cardiac tissue health. To begin to test this hypothesis, male and female rats underwent HU at various durations up to 90 days, with a subset reambulated after 90 days of HU. Physiological stress or contractility changes lead to alterations in ventricular cardiomyocyte size and ventricular wall thickness to adapt to greater functional demand and mitigate mechanical stress to ventricular tissue; under certain conditions, these changes also may mark progression to cardiac failure. Hence, left ventricular cardiomyocyte size (cardiomyocyte cross sectional area, CSA) was quantified to determine if HU leads to structural adaptation responses in cardiac tissue and if age and sex had any impact on this outcome. Cardiomyocyte CSA of older males (9 months) were altered by HU in a time-dependent manner, where HU led to decreases in CSA at 14 days and increases at 90 days. In contrast, younger males (3 months) did not show any changes at day 14 of HU. CSA of females (3 months) was increased in response to short-term HU (14 days) suggesting sex-dependence of structural changes. In older HU males, cardiomyocyte CSA was comparable to controls after 90 days of re-ambulation. Levels of the DNA oxidative damage marker, 8-hydroxydeoxyguanosine (8-OHdG) were greater in left ventricular tissue of females that underwent HU compared to sex-matched controls, while there were no such differences in older or younger males. To gain insight into the signals that drive cardiac adaptations to HU, global transcriptomic analysis (RNAseq) was performed on left ventricular tissue of older males that underwent 14 days of HU. Short-term simulated weightlessness led to differential expression of genes involved in immune and pro-inflammatory signaling. A subset of these genes play a role in autoimmune and cardiovascular disease and are targets of current drugs used to treat bradycardia, hypertension, atherosclerosis and rheumatoid arthritis, amongst others. Oxidative damage/redox signaling pathways were not enriched at the timepoint tested in older males. Since young females displayed greater oxidative damage to DNA, activation of oxidative stress responses at earlier or later time points cannot be ruled out. In summary, simulated weightlessness in adult rats caused changes in cardiomyocyte structure in a sex and age-dependent manner, and the transcriptional regulation of key mediators of immunity and cardiovascular disease, meriting further study to define cardiac risks for interplanetary travel of human crew. Our findings also confirm the value of the rat HU model for cardiac health and countermeasure research.

Tahimic, Candice↗

DOE OSTI2026

Glovebox CFD Analysis: Low-Flow, Laminar CFD Analysis of Glovebox Ventilation and Humidity Transport in Ansys Fluent 2025 R1

There are gloveboxes that frequently exceed the 30% relative humidity level required for operations. These exceedances result in alternative operations needing to be utilized. This report documents a computational fluid dynamics (CFD) analysis performed to characterize airflow and humidity transport within a representative glovebox under nominal, low-flow conditions. The analysis focuses on establishing the flow regime, identifying transport-limited regions, and evaluating the implications of laminar flow on humidity removal. A single-phase, laminar flow model with a species transport model was used to represent air and water vapor mixing. Results demonstrate that, at the nominal flow rates, airflow is laminar with weak mixing, which will result in long residence times, tens of minutes long or more, in low-velocity regions. These conclusions suggests that humidity reduction is transport-limited. This analysis provides qualitative insight into airflow structure and humidity transport trends. In addition to characterization of glovebox airflow and humidity transport, results are intended to inform experimental validation and future transient modeling efforts.

42 ENGINEERING↗

DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗

NASA NTRS1996

Differences in Femoral Geometry and Structure Due to Immobilization

Reduction in bone mass of the lower extremity is well documented in individuals with paralysis resulting from spinal cord injury (SCI). The consequent osteopenia leads to elevated fracture risk with fractures occurring more commonly in the femoral shaft and supracondylar regions than the hip. A model has recently been described to estimate geometry and structure of the femoral midshaft from whole body scans by dual X-ray absorptiometry (DXA). Increases in femoral geometric and structural properties during growth were primarily related to mechanical loading as reflected by body mass. In this study, we investigate the relationship between body mass and femoral geometry and structure in adults with normal habitual mechanical loading patterns and those with severely reduced loading. The subjects were 78 ambulatory men (aged 20-72 yrs) and 113 men with complete paralysis from SCI of more than 4 years duration (aged 21 73 yrs). Subregional analysis was performed on DXA whole body scans to obtain bone mineral content (BMC, g), cortical thickness (cm), crosssectional moment of inertia (CSMI, cm4), and section modulus (cm3) of the femoral midshaft. All measured bone variables were significantly lower in SCI compared with ambulatory subjects: -29% (BMC), -33% (cortical thickness), -23% (CSMI), and -22% (section modulus) while body mass was not significantly different. However, the associations between body mass and bone properties were notably different; r2 values were higher for ambulatory than SCI subjects in regressions of body mass on BMC (0.48 vs 0.20), CSMI (0.59 vs 0.32), and section modulus (0.59 vs 0.31). No association was seen between body mass and cortical thickness for either group. The greatest difference between groups is in the femoral cortex, consistent with reduced bone mass via endosteal expansion. The relatively lesser difference in geometric and structural properties implies that there is less effect on mechanical integrity than would be expected from bone mass results alone. The reduced association in SCI subjects between body mass and bone properties is not unexpected. Although mean body mass differs little between ambulatory and SCI individuals, the association between body mass and in vivo skeletal loading is no longer present, as mechanical influences are removed except for transfer activities. The residual association is probably attributable to the strength of this influence during growth. These results highlight the importance of examining geometry and structure in conjunction with bone mass.

Kiratli, Beatrice Jenny↗

NASA NTRS2010

Effect of Turbulence Modeling on an Excited Jet

The flow dynamics in a high-speed jet are dominated by unsteady turbulent flow structures in the plume. Jet excitation seeks to control these flow structures through the natural instabilities present in the initial shear layer of the jet. Understanding and optimizing the excitation input, for jet noise reduction or plume mixing enhancement, requires many trials that may be done experimentally or computationally at a significant cost savings. Numerical simulations, which model various parts of the unsteady dynamics to reduce the computational expense of the simulation, must adequately capture the unsteady flow dynamics in the excited jet for the results are to be used. Four CFD methods are considered for use in an excited jet problem, including two turbulence models with an Unsteady Reynolds Averaged Navier-Stokes (URANS) solver, one Large Eddy Simulation (LES) solver, and one URANS/LES hybrid method. Each method is used to simulate a simplified excited jet and the results are evaluated based on the flow data, computation time, and numerical stability. The knowledge gained about the effect of turbulence modeling and CFD methods from these basic simulations will guide and assist future three-dimensional (3-D) simulations that will be used to understand and optimize a realistic excited jet for a particular application.

Brown, Clifford A.↗

NASA NTRS1992

Eigenvalue error analysis of viscously damped structures using a Ritz reduction method

The efficient solution of the eigenvalue problem that results from inserting passive dampers with variable stiffness and damping coefficients into a structure is addressed. Eigenanalysis of reduced models obtained by retaining a number of normal modes augmented with Ritz vectors corresponding to the static solutions resulting from the load patterns introduced by the dampers has been empirically shown to yield excellent approximations to the full eigenvalue problem. An analysis of this technique in the case of a single damper is presented. A priori and a posteriori error estimates are generated and tested on numerical examples. Comparison theorems with modally truncated models and a Markov parameter matching reduced-order model are derived. These theorems corroborate the heuristic that residual flexibility methods improve low-frequency approximation of the system. The analysis leads to other techniques for eigenvalue approximation. Approximate closed-form solutions are derived that include a refinement to eigenvalue derivative methods for approximation. An efficient Newton scheme is also developed. A numerical example is presented demonstrating the effectiveness of each of these methods.

Chu, Cheng-Chih↗

DOE OSTI2021

A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models

In the present study, the multiphase volume distribution problem, where there can be an arbitrary number of phases, is addressed using a consistent and conservative volume distribution algorithm. The proposed algorithm satisfies the summation constraint, the conservation constraint, and the consistency of reduction. The first application of the volume distribution algorithm is to determine the Lagrange multipliers in multiphase Phase-Field models that enforce the mass conservation, and a multiphase conservative Allen-Cahn model that satisfies the consistency of reduction is developed. A corresponding consistent and conservative numerical scheme is developed for the model. The multiphase conservative Allen-Cahn model has a better ability than the multiphase Cahn-Hilliard model to preserve under-resolved structures. The second application is to develop a numerical procedure, called the boundedness mapping, to map the order parameters, obtained numerically from a multiphase model, into their physical interval, and at the same time to preserve the physical properties of the order parameters. Along with the consistent and conservative schemes for the multiphase Phase-Field models, the numerical solutions of the order parameters are reduction consistent, conservative, and bounded, which are theoretically analyzed and numerically validated. Then, the multiphase Phase-Field models are coupled with the momentum equation by satisfying the consistency of mass conservation and the consistency of mass and momentum transport, thanks to the consistent formulation. Finally, it is demonstrated that the proposed model and scheme converge to the sharp-interface solution and are capable of capturing the complicated multiphase dynamics even when there is a large density and/or viscosity ratio.

42 ENGINEERING↗