Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational challenges”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Discrete element modeling of deformable pinewood chips in cyclic loading test

The design of efficient lignocellulosic biomass feedstock systems is challenging, as current laboratory characterization and design methods were developed primarily for fine powders with relatively low compressibility. Here, the discrete element method (DEM) is gaining prominence as an alternative method for modeling the bulk flow and transport of particulate materials in hoppers and feeders. However, prior DEM simulations investigating the flow of wood chips modeled the particles as simple rigid geometries such as spheres, rods or blocks, and neglected the effects of particle deformability and irregular shapes. As a consequence, those simplified DEM models may not provide enough key diagnosis to help improve the design of biomass feeding and handling equipment. This work presents a bonded-sphere DEM approach for characterizing the mechanical behavior of bulk flexible, deformable pinewood chips in a cyclic stress loading test. Clustered spheres that can bend and twist via elastic bonds are used to model irregular-shape particles in real pinewood chip samples. An axial compression tester, which contains 0.06 million bonded-sphere particles (1.35 million spheres) in a quarter cylinder, is simulated with a domain size similar to the experiments. With careful calibration, the simulations have delivered the bulk densities and the bulk moduli of elasticity that are in good agreement with those measured in the corresponding experiments. However, it is also been found challenging for the present DEM model to accurately predict the overall stress-strain behavior of bulk pinewood chips, especially the large sustained plastic deformation during unloading. Additional numerical tests have shown that the adjustment in certain contact parameters (e.g. bond stiffness) can lead to profound solution improvement, but meanwhile will induce extra challenges such as increased computing time. Future work will include an elasto-plastic particle bond model to enhance the simulation fidelity.

42 ENGINEERING↗
DOE OSTI2021

Multi-fidelity kinetic theory-based approach for the prediction of particle attrition: Application to jet cup attrition system

The timescale difference between the “fast” flow dynamics of fluidized bed reactors and the relatively “slow” rate of particle degradation makes the direct computational prediction of attrition challenging. An approach to this challenge is a multi-fidelity strategy where a high fidelity model for the flow dynamics is coupled with a lower fidelity model for the long-time resolution of the bulk attrition of the reactor inventory. We implement this approach using high-fidelity kinetic theory simulations to calculate the flow dynamics which are post-processed to calculate the frequency and intensity of the particle-particle and particle-wall collisions (e.g. collision energy spectra). This is combined with the particle breakage properties to construct the coefficients for a low-fidelity model [e.g. Monazam et al., 2018, Powder Technology 340, p. 528-536]. Simulations are performed of a jet cup attrition system containing Canadian hematite (Monazam et al. 2018). Furthermore, these are first analyzed using the collision energy spectra. Quantitative predictions of the mass loss are made using a low-fidelity model derived from the collision-spectra and a calibrated material breakage coefficient. The results are found to compare favorably with the experimental measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Urban building energy modeling (UBEM) tools: A state-of-the-art review of bottom-up physics-based approaches

Regulations corroborate the importance of retrofitting existing building stocks or constructing new energy-efficient districts. There is, thus, a need for modeling tools to evaluate energy scenarios to better manage and design cities, and numerous methodologies and tools have been developed. Among them, Urban Building Energy Modelling (UBEM) tools allow the energy simulation of buildings at large scales. Choosing an appropriate UBEM tool, balancing the level of complexity, accuracy, usability, and computing needs, remains a challenge for users. The review focuses on the main bottom-up physics-based UBEM tools, comparing them from a user-oriented perspective. Five categories are used: (i) the required inputs, (ii) the reported outputs, (iii) the exploited workflow, (iv) the applicability of each tool, and (v) the potential users. Moreover, a critical discussion is proposed, focusing on interests and trends in research and development. The results highlighted major differences between UBEM tools that must be considered to choose the proper one for an application. Finally, barriers of adoption of UBEM tools include the needs of a standardized ontology, a common three-dimensional city model, a standard procedure to collect data, and a standard set of test cases. This feeds into future development of UBEM tools to support cities’ sustainability goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Predictive CRISPR-mediated gene downregulation for enhanced production of sustainable aviation fuel precursor in Pseudomonas putida

CRISPR interference (CRISPRi) has emerged as a valuable tool for redirecting metabolic flux to enhance bioproduction. However, its application is often constrained by two challenges: (i) rationally identifying effective gene targets for downregulation and (ii) efficiently constructing multiplexed CRISPRi systems. In this study, we address both challenges by integrating a computational prioritization tool with a versatile assembly method for building multiplexed CRISPRi systems. FluxRETAP (Flux-Reaction Target Prioritization) accurately identified gene targets whose knockdown led to substantial increase of isoprenol titers in Pseudomonas putida KT2440, outperforming a conventional non-computational, pathway-guided target selection. The highest isoprenol titer of nearly 1.5 g/L was achieved by knocking down PP_4118 (a gene encoding α-ketoglutarate dehydrogenase). The use of VAMMPIRE (Versatile Assembly Method for MultiPlexing CRISPRi-mediated downREgulation) enabled accurate assembly of CRISPRi constructs containing up to five sgRNA arrays, reducing context dependency and achieving uniform, position-independent gene downregulation. The integration of FluxRETAP and VAMMPIRE has the potential to advance metabolic engineering by rapidly identifying CRISPRi-mediated knockdowns and knockdown combinations that enhance bioproduction titers, with potential applicability to other microbial systems.

CRISPR interference↗
DOE OSTI2019

Search for Catalysts by Inverse Design: Artificial Intelligence, Mountain Climbers, and Alchemists

In silico catalyst design is a grand challenge of chemistry. Traditional computational approaches have been limited by the need to compute properties for an intractably large number of possible catalysts. Recently, inverse design methods have emerged, starting from a desired property and optimizing a corresponding chemical structure. Techniques used for exploring chemical space include gradient-based optimization, alchemical transformations, and machine learning. Though the application of these methods to catalysis is in its early stages, further development will allow for robust computational catalyst design. This review provides an overview of the evolution of inverse design approaches and their relevance to catalysis. In conclusion, the strengths and limitations of existing techniques are highlighted, and suggestions for future research are provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Distinct Mechanisms and Hydricities of Cp*Ir-Based CO 2 Hydrogenation Catalysts in Basic Water

Transition metal-catalyzed reversible hydrogenation of CO 2 to formate in aqueous solutions under ambient conditions is an attractive environmentally friendly strategy for the storage and transportation of H 2 in a liquid chemical form. Here, mechanistic details for the CO 2 hydrogenation by a series of pentamethylcyclopentadienyl iridium(III) (Cp*Ir) complexes with picolinamidate ligands are investigated through density functional theory and kinetic isotope effects (KIEs) calculations. Additionally, adopting a speciation approach at a pH of 8.2, CO 2 hydrogenation mechanisms of 37 distinct protonation states and conformers of Cp*Ir type complexes are investigated. This imposes a challenging test for the computational modeling in terms of providing a consistent correlation with experimental kinetics and KIE studies of multiple catalysts instead of a single complex. Overall, H 2 heterolysis to generate an iridium hydride intermediate was demonstrated to be the rate determining step, which proceeds via distinct pathways depending on the stereoelectronic properties of the ligand. A peculiar iridium dihydride route was also uncovered, which could be optimized to accelerate the H 2 heterolysis. We have further computed the thermodynamic and kinetic hydricities of 74 complexes with [Cp*Ir(L)(H)] q ( q = -1, 0, 1 depending on the ligand charge by deprotonation or protonation events) general formula. Hydricity values do not show any noticeable correlation with the thermodynamics (Δ G ) of [Cp*Ir(L)(HCO 2 – )] formation but exhibit linear correlation with the kinetics (Δ G ‡ ) of electrophilic CO 2 attack to iridium hydride species, especially when the charges and local structural effects of the metal hydrides are considered. Insights gained in this study on (i) the factors determining the preferred pathway for the rate-limiting H 2 heterolysis step, (ii) the correlation between thermodynamic hydricity and the kinetics of CO 2 insertion to iridium hydride species, and (iii) mechanistic analysis via combination of free energy profiles and KIE studies for a series of iridium catalysts will provide guidelines for the design of next-generation catalysts for the reversible H 2 storage through CO 2 utilization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A topological superconductor tuned by electronic correlations

A topological superconductor, characterized by either a chiral order parameter or a topological surface state in proximity to bulk superconductivity, is foundational to topological quantum computing. A key open challenge is whether electron-electron interactions can tune such topological superconducting phases. Here, we provide experimental signatures of a unique topological superconducting phase in competition with electronic correlations in 10-unit-cell thick FeTe x Se 1-x films grown on SrTiO 3 substrates. When the Te content x exceeds 0.7, we observe a topological transition marked by the emergence of a superconducting surface state. Near the FeTe limit, the system undergoes another transition where the surface state disappears, and superconductivity is suppressed. Theory suggests that electron-electron interactions in the odd-parity xy− band drives this second topological transition. The flattening and eventual decoherence of d xy -derived bands track the superconducting dome, linking correlation effects directly to superconducting coherent transport. Our work establishes many-body electronic correlations as a sensitive knob for tuning topology and superconductivity, offering a pathway to engineer new topological phases in correlated materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Large-scale integration of artificial atoms in hybrid photonic circuits

A central challenge in developing quantum computers and long-range quantum networks is the distribution of entanglement across many individually controllable qubits. Colour centres in diamond have emerged as leading solid-state ‘artificial atom’ qubits because they enable on-demand remote entanglement, coherent control of over ten ancillae qubits with minute-long coherence times and memory-enhanced quantum communication. A critical next step is to integrate large numbers of artificial atoms with photonic architectures to enable large-scale quantum information processing systems. So far, these efforts have been stymied by qubit inhomogeneities, low device yield and complex device requirements. Here we introduce a process for the high-yield heterogeneous integration of ‘quantum microchiplets’—diamond waveguide arrays containing highly coherent colour centres—on a photonic integrated circuit (PIC). We use this process to realize a 128-channel, defect-free array of germanium-vacancy and silicon-vacancy colour centres in an aluminium nitride PIC. Photoluminescence spectroscopy reveals long-term, stable and narrow average optical linewidths of 54 megahertz (146 megahertz) for germanium-vacancy (silicon-vacancy) emitters, close to the lifetime-limited linewidth of 32 megahertz (93 megahertz). We show that inhomogeneities of individual colour centre optical transitions can be compensated in situ by integrated tuning over 50 gigahertz without linewidth degradation. The ability to assemble large numbers of nearly indistinguishable and tunable artificial atoms into phase-stable PICs marks a key step towards multiplexed quantum repeaters and general-purpose quantum processors

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct evidence of microstructure dependence of magnetic flux trapping in niobium

Abstract Elemental type-II superconducting niobium is the material of choice for superconducting radiofrequency cavities used in modern particle accelerators, light sources, detectors, sensors, and quantum computing architecture. An essential challenge to increasing energy efficiency in rf applications is the power dissipation due to residual magnetic field that is trapped during the cool down process due to incomplete magnetic field expulsion. New SRF cavity processing recipes that use surface doping techniques have significantly increased their cryogenic efficiency. However, the performance of SRF Nb accelerators still shows vulnerability to a trapped magnetic field. In this manuscript, we report the observation of a direct link between flux trapping and incomplete flux expulsion with spatial variations in microstructure within the niobium. Fine-grain recrystallized microstructure with an average grain size of 10–50 µm leads to flux trapping even with a lack of dislocation structures in grain interiors. Larger grain sizes beyond 100–400 µm do not lead to preferential flux trapping, as observed directly by magneto-optical imaging. While local magnetic flux variations imaged by magneto-optics provide clarity on a microstructure level, bulk variations are also indicated by variations in pinning force curves with sequential heat treatment studies. The key results indicate that complete control of the niobium microstructure will help produce higher performance superconducting resonators with reduced rf losses 1 related to the magnetic flux trapping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Noise reduction in X-ray photon correlation spectroscopy with convolutional neural networks encoder–decoder models

Abstract Like other experimental techniques, X-ray photon correlation spectroscopy is subject to various kinds of noise. Random and correlated fluctuations and heterogeneities can be present in a two-time correlation function and obscure the information about the intrinsic dynamics of a sample. Simultaneously addressing the disparate origins of noise in the experimental data is challenging. We propose a computational approach for improving the signal-to-noise ratio in two-time correlation functions that is based on convolutional neural network encoder–decoder (CNN-ED) models. Such models extract features from an image via convolutional layers, project them to a low dimensional space and then reconstruct a clean image from this reduced representation via transposed convolutional layers. Not only are ED models a general tool for random noise removal, but their application to low signal-to-noise data can enhance the data’s quantitative usage since they are able to learn the functional form of the signal. We demonstrate that the CNN-ED models trained on real-world experimental data help to effectively extract equilibrium dynamics’ parameters from two-time correlation functions, containing statistical noise and dynamic heterogeneities. Strategies for optimizing the models’ performance and their applicability limits are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Visualizing heterogeneous dipole fields by terahertz light coupling in individual nano-junctions

The challenge underlying superconducting quantum computing is to remove materials bottleneck for highly coherent quantum devices. The nonuniformity and complex structural components in the underlying quantum circuits often lead to local electric field concentration, charge scattering, dissipation and ultimately decoherence. Here we visualize interface dipole heterogeneous distribution of individual Al/AlO$_{x}$/Al junctions employed in transmon qubits by broadband terahertz scanning near-field microscopy that enables the non-destructive and contactless identification of defective boundaries in nano-junctions at an extremely precise nanoscale level. Our THz nano-imaging tool reveals an asymmetry across the junction in electromagnetic wave-junction coupling response that manifests as hot (high intensity) vs cold (low intensity) spots in the spatial electrical field structures and correlates with defected boundaries from the multi-angle deposition processes in Josephson junction fabrication inside qubit devices. The demonstrated local electromagnetic scattering method offers high sensitivity, allowing for reliable device defect detection in the pursuit of improved quantum circuit fabrication for ultimately optimizing coherence times.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Citrate stabilized gold nanoparticles interfere with amyloid fibril formation: D76N and ΔN6 β2-microglobulin variants

Here, protein aggregation including the formation of dimers and multimers in solution, underlies an array of human diseases such as systemic amyloidosis which is a fatal disease caused by misfolding of native globular proteins damaging the structure and function of affected organs. Different kind of interactors can interfere with the formation of protein dimers and multimers in solution. A very special class of interactors are nanoparticles thanks to the extremely efficient extension of their interaction surface. In particular citrate-coated gold nanoparticles (cit-AuNPs) were recently investigated with amyloidogenic protein β 2 -microglobulin (β 2 m). Here we present the computational studies on two challenging models known for their enhanced amyloidogenic propensity, namely ΔN6 and D76N β 2 m naturally occurring variants, and disclose the role of cit-AuNPs on their fibrillogenesis. The proposed interaction mechanism lies in the interference of the cit-AuNPs with the protein dimers at the early stages of aggregation, that induces dimer disassembling. As a consequence, natural fibril formation can be inhibited. Relying on the comparison between atomistic simulations at multiple levels (enhanced sampling molecular dynamics and Brownian dynamics) and protein structural characterisation by NMR, we demonstrate that the cit-AuNPs interactors are able to inhibit protein dimer assembling. As a consequence, the natural fibril formation is also inhibited, as found in experiment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Iterative map-making with two-level preconditioning for polarized cosmic microwave background data sets: A worked example for ground-based experiments

Context. An estimation of the sky signal from streams of time ordered data (TOD) acquired by the cosmic microwave background (CMB) experiments is one of the most important steps in the context of CMB data analysis referred to as the map-making problem. The continuously growing CMB data sets render the CMB map-making problem progressively more challenging in terms of computational cost and memory in particular in the context of ground-based experiments with their operational limitations as well as the presence of contaminants. Aims. We study a recently proposed, novel class of the Preconditioned Conjugate Gradient (PCG) solvers which invoke two-level preconditioners in the context of the ground-based CMB experiments. We compare them against the PCG solvers commonly used in the map-making context considering their precision and time-to-solution. Methods. We compare these new methods on realistic, simulated data sets reflecting the characteristics of current and forthcoming CMB ground-based experiments. We develop a divide-and-conquer implementation of the approach where each processor performs a sequential map-making for a subset of the TOD. Results. We find that considering the map level residuals, the new class of solvers permits us to achieve a tolerance that is better than the standard approach by up to three orders of magnitude, where the residual level often saturates before convergence is reached. This often corresponds to an important improvement in the precision of the recovered power spectra in particular on the largest angular scales. The new method also typically requires fewer iterations to reach a required precision and therefore shorter run times are required for a single map-making solution. However, the construction of an appropriate two-level preconditioner can be as costly as a single standard map-making run. Nevertheless, if the same problem needs to be solved multiple times, for example, as in Monte Carlo simulations, this cost is incurred only once, and the method should be competitive, not only as far as its precision is concerned but also its performance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Exploring Future Storage Options for ATLAS at the BNL/SDCC facility

The ATLAS experiment is expected to deliver an unprecedented amount of scientific data in the High Luminosity(HL-LHC) era. As the demand for disk storage capacity in ATLAS continues to rise steadily, the BNL Scientific Data and Computing Center (SDCC) faces challenges in terms of cost implications for maintaining multiple disk copies and adapting to the coming ATLAS storage requirements. To address these challenges, the SDCC Storage team has undertaken a thorough analysis of the ATLAS experiment's requirements, matching them to suitable storage options and strategies, and has explored alternatives to enhance or replace the current storage solution. This paper aims to present the main challenges encountered while supporting big data experiments such as ATLAS. We describe the experiment's specific requirements and priorities, particularly focusing on the critical storage system characteristics of the high-luminosity run and how the key storage components provided by the Storage team work together: the dCache disk storage system; its archival back-end, HPSS; and its OS-level backend Storage. Specifically, we investigate a novel approach to integrate Lustre and XRootD. In this setup, Lustre serves as backend storage and XRootD acts as an access layer frontend, supporting various grid access protocols. Additionally, we also describe the validation and commissioning tests, including the performance comparison between dCache and XRootd. Furthermore, we provide a performance and cost analysis comparing OpenZFS and LINUX MD RAID, evaluate different storage software stacks, and showcase stress tests conducted to validate Third Party Copy (TPC) functionality.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2014

Properties of reactive oxygen species by quantum Monte Carlo

The electronic properties of the oxygen molecule, in its singlet and triplet states, and of many small oxygen-containing radicals and anions have important roles in different fields of chemistry, biology, and atmospheric science. Nevertheless, the electronic structure of such species is a challenge for ab initio computational approaches because of the difficulties to correctly describe the statical and dynamical correlation effects in presence of one or more unpaired electrons. Only the highest-level quantum chemical approaches can yield reliable characterizations of their molecular properties, such as binding energies, equilibrium structures, molecular vibrations, charge distribution, and polarizabilities. In this work we use the variational Monte Carlo (VMC) and the lattice regularized Monte Carlo (LRDMC) methods to investigate the equilibrium geometries and molecular properties of oxygen and oxygen reactive species. Quantum Monte Carlo methods are used in combination with the Jastrow Antisymmetrized Geminal Power (JAGP) wave function ansatz, which has been recently shown to effectively describe the statical and dynamical correlation of different molecular systems. In particular, we have studied the oxygen molecule, the superoxide anion, the nitric oxide radical and anion, the hydroxyl and hydroperoxyl radicals and their corresponding anions, and the hydrotrioxyl radical. Overall, the methodology was able to correctly describe the geometrical and electronic properties of these systems, through compact but fully-optimised basis sets and with a computational cost which scales as N{sup 3} − N{sup 4}, where N is the number of electrons. This work is therefore opening the way to the accurate study of the energetics and of the reactivity of large and complex oxygen species by first principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Tuned range separated hybrid functionals for solvated low bandgap oligomers

The description of charge transfer excitations has long been a challenge to time dependent density functional theory. The recently developed concept of “optimally tuned range separated hybrid (OT-RSH) functionals” has proven to describe charge transfer excitations accurately in many cases. However, describing solvated or embedded systems is yet a challenge. This challenge is not only computational but also conceptual, because the tuning requires identifying a specific orbital, typically the highest occupied one of the molecule under study. For solvated molecules, this orbital may be delocalized over the solvent. We here demonstrate that one way of overcoming this problem is to use a locally projected self-consistent field diagonalization on an absolutely localized molecular orbital expansion. We employ this approach to determine ionization energies and the optical gap of solvated oligothiophenes, i.e., paradigm low gap systems that are of relevance in organic electronics. Dioxane solvent molecules are explicitly represented in our calculations, and the ambiguities of straightforward parameter tuning in solution are elucidated. We show that a consistent estimate of the optimal range separated parameter (ω) at the limit of bulk solvation can be obtained by gradually extending the solvated system. In particular, ω is influenced by the solvent beyond the first coordination sphere. For determining ionization energies, a considerable number of solvent molecules on the first solvation shell must be taken into account. We demonstrate that accurately calculating optical gaps of solvated systems using OT-RSH can be done in three steps: (i) including the chemical environment when determining the range-separation parameter, (ii) taking into account the screening due to the solvent, and (iii) using realistic molecular geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A simulation method for the phase diagram of complex fluid mixtures

We present that the phase behavior of complex fluid mixtures is of continuing interest, but obtaining the phase diagram from computer simulations can be challenging. In the Gibbs ensemble method, for example, each of the coexisting phases is simulated in a different cell, and ensuring the equality of chemical potentials of all components requires the transfer of molecules from one cell to the other. For complex fluids such as polymers, successful insertions are rare. An alternative method is to simulate both coexisting phases in a single simulation cell, with an interface between them. The challenge here is that the interface position moves during the simulation, making it difficult to determine the concentration profile and coexisting concentrations. In this work, we propose a new method for single cell simulations that uses a spatial concentration autocorrelation function to (spatially) align instantaneous concentration profiles from different snapshots. This allows one to obtain average concentration profiles and hence the coexisting concentrations. We test the method by calculating the phase diagrams of two systems: the Widom-Rowlinson model and the symmetric blends of freely jointed polymer molecules for which phase diagrams from conventional methods are available. Excellent agreement is found, except in the neighborhood of the critical point where the interface is broad and finite size effects are important. Finally, the method is easy to implement and readily applied to any mixture of complex fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗