Search NASA⌕ Search

SEARCH · Search NASA

Results for “APPLES”

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 379 records · Page 21

Adiponectin enhances osteogenic differentiation in human adipose-derived stem cells by activating the APPL1-AMPK signaling pathway

Human adipose-derived stem cells (hASCs) are multipotent progenitor cells with multi-lineage differentiation potential including osteogenesis and adipogenesis. While significant progress has been made in understanding the transcriptional control of hASC fate, little is known about how hASC differentiation is regulated by the autocrine loop. The most abundant adipocytokine secreted by adipocytes, adiponectin (APN) plays a pivotal role in glucose metabolism and energy homeostasis. Growing evidence suggests a positive association between APN and bone formation yet little is known regarding the direct effects of APN on hASC osteogenesis. Therefore, this study was designed to investigate the varied osteogenic effects and regulatory mechanisms of APN in the osteogenic commitment of hASCs. We found that APN enhanced the expression of osteoblast-related genes in hASCs, such as osteocalcin, alkaline phosphatase, and runt-related transcription factor-2 (Runx2, also known as CBFa1), in a dose- and time-dependent manner. This was further confirmed by the higher expression levels of alkaline phosphatase and increased formation of mineralization nodules, along with the absence of inhibition of cell proliferation. Importantly, APN at 1 μg/ml was the optimal concentration, resulting in maximum deposition of calcium nodules, and was significant superior to bone morphogenetic protein 2. Mechanistically, we found for the first time that APN increased nuclear translocation of the leucine zipper motif (APPL)-1 as well as AMP-activated protein kinase (AMPK) phosphorylation, which were reversed by pretreatment with APPL1 siRNA. Our results indicate that APN promotes the osteogenic differentiation of hASCs by activating APPL1-AMPK signaling, suggesting that manipulation of APN is a novel therapeutic target for controlling hASC fate. - Highlights: • Adiponectin enhances osteogenic differentiation in human adipose-derived stem cells. • The knock-down of APPL1 block the enhancement of osteogenic differentiation in hASC. • AMPK signaling cascade is activated by adiponectin through APPL1.

60 APPLIED LIFE SCIENCES↗

Overexpression of a repressor MdMYB15L negatively regulates anthocyanin and cold tolerance in red-fleshed callus

Highlights: • MdMYB15L inhibit the expression of MdCBF2 by binding its promoter. • MdMYB15L interact with MdbHLH33 via its bHLH binding sequence. • MdMYB15L weaken MdbHLH33-induced anthocyanin accumulation by interacting with MdbHLH33. The cold-induced metabolic pathway and anthocyanin biosynthesis play important roles in plant growth. In this study, we identified a bHLH binding motif in the MdMYB15L protein using protein sequence analyses. Yeast two-hybrid and pull-down assays showed that MdMYB15L could interact with MdbHLH33. Overexpressing MdMYB15L in red-fleshed callus inhibited the expression of MdCBF2 and resulted in reduced cold tolerance but did not affect anthocyanin levels. Chip-PCR and EMSA analysis showed that MdMYB15L could bind the type II cis-acting element found in the MdCBF2 promoter. Overexpressing MdMYB15L in red-fleshed callus overexpressing MdbHLH33 also reduced cold tolerance and reduced MdbHLH33-induced anthocyanin biosynthesis. Knocking out the bHLH binding sequence of MdMYB15L (LBSMdMYB15L) prevented LBSMdMYB15L from interacting with MdbHLH33. Overexpressing LBSMdMYB15L in red-fleshed callus overexpressing MdbHLH33 also reduced cold tolerance and reduced MdbHLH33-induced anthocyanin biosynthesis. Together, these results suggested that an apple repressor MdMYB15L might play a key role in the cold signaling and anthocyanin metabolic pathways.

60 APPLIED LIFE SCIENCES↗

Scalable DPG multigrid solver for Helmholtz problems: A study on convergence

This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from $\mathscr{O}$(10 7 ) degrees of freedom (DOFs) to $\mathscr{O}$(10 9 ) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. Finally, the paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast.

97 MATHEMATICS AND COMPUTING↗

Coupled uptake and conversion of C 12 H 26 and NO on Pd/SSZ-13: Experiments and modeling

The uptake of NO in the presence of C 12 H 26 (dodecane) and H 2 O over a Pd/SSZ-13 washcoated Passive NOx Adsorber (PNA) monolith is reported. When a co-feed containing C 12 H 26 and NO is supplied to an unsaturated (with C 12 H 26 ) sample, the NO uptake is unaffected but during the subsequent temperature ramp the release of trapped NO is delayed from 175°C to over 220°C. The release delay is beneficial for PNA performance as the primary NO x aftertreatment technology, Selective Catalytic Reduction (SCR), is not operated below 200°C. However, pre-saturation of C 12 H 26 followed by the same NO and C 12 H 26 co-feed results in a decrease in the NO uptake compared to the NO-only feed. We conjecture that C 12 H 26 pre-adsorbed on the exterior surface of the sample blocks NO access to the pores, decreasing the number of available sites for NO uptake. Oxidation of C 12 H 26 leads to the generation of partial oxidation product CO at lower temperatures (<250°C) and deep oxidation product CO 2 at higher temperatures. Carbon monoxide binds strongly to Pd sites with NO and can delay NO release. A Pd/SSZ-13 washcoated monolith model developed in an earlier study [Ambast et al., Appl. Catal. B. Environmental (2021)] is upgraded to include C 12 H 26 storage, release, and conversion. Through a systematic combination of judicious experiments, model tuning, and validation, we provide evidence for the underlying NO uptake and release mechanisms in the presence of C 12 H 26 .

42 ENGINEERING↗

Self-consistent clustering analysis: an efficient multiscale scheme for inelastic heterogeneous materials

A predictive computational theory is shown for modeling complex, hierarchical materials ranging from metal alloys to polymer nanocomposites. The theory can capture complex mechanisms such as plasticity and failure that span across multiple length scales. This general multiscale material modeling theory relies on sound principles of mathematics and mechanics, and a cutting-edge reduced order modeling method named self-consistent clustering analysis (SCA) [Zeliang Liu, M.A. Bessa, Wing Kam Liu, “Self-consistent clustering analysis: An efficient multi-scale scheme for inelastic heterogeneous materials,” Comput. Methods Appl. Mech. Engrg. 306 (2016) 319–341]. SCA reduces by several orders of magnitude the computational cost of micromechanical and concurrent multiscale simulations, while retaining the microstructure information. This remarkable increase in efficiency is achieved with a data-driven clustering method. Computationally expensive operations are performed in the so-called offline stage, where degrees of freedom (DOFs) are agglomerated into clusters. The interaction tensor of these clusters is computed. In the online or predictive stage, the Lippmann-Schwinger integral equation is solved cluster-wise using a self-consistent scheme to ensure solution accuracy and avoid path dependence. To construct a concurrent multiscale model, this scheme is applied at each material point in a macroscale structure, replacing a conventional constitutive model with the average response computed from the microscale model using just the SCA online stage. A regularized damage theory is incorporated in the microscale that avoids the mesh and RVE size dependence that commonly plagues microscale damage calculations. The SCA method is illustrated with two cases: a carbon fiber reinforced polymer (CFRP) structure with the concurrent multiscale model and an application to fatigue prediction for additively manufactured metals. For the CFRP problem, a speed up estimated to be about 43,000 is achieved by using the SCA method, as opposed to FE2, enabling the solution of an otherwise computationally intractable problem. The second example uses a crystal plasticity constitutive law and computes the fatigue potency of extrinsic microscale features such as voids. This shows that local stress and strain are capture sufficiently well by SCA. This model has been incorporated in a process-structure-properties prediction framework for process design in additive manufacturing.

Liu, Z.↗

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗

An Efficient Algorithm for Automatic Structure Optimization in X-ray Standing-Wave Experiments

We present that x-ray standing-wave photoemission experiments involving multilayered samples are emerging as unique probes of the buried interfaces that are ubiquitous in current device and materials research. Such data require for their analysis a structure optimization process comparing experiment to theory that is not straightforward. In this work, we present a new computer program for optimizing the analysis of standing-wave data, called SWOPT, that automates this trial-and-error optimization process. The program includes an algorithm that has been developed for computationally expensive problems: so-called black-box simulation optimizations. It also includes a more efficient version of the Yang X-ray Optics Program (YXRO) [Yang, S.-H., Gray, A.X., Kaiser, A.M., Mun, B.S., Sell, B.C., Kortright, J.B., Fadley, C.S., J. Appl. Phys. 113, 1 (2013)] which is about an order of magnitude faster than the original version. Human interaction is not required during optimization. We tested our optimization algorithm on real and hypothetical problems and show that it finds better solutions significantly faster than a random search approach. The total optimization time ranges, depending on the sample structure, from minutes to a few hours on a modern laptop computer, and can be up to 100x faster than a corresponding manual optimization. Finally, these speeds make the SWOPT program a valuable tool for realtime analyses of data during synchrotron experiments.

36 MATERIALS SCIENCE↗

Energy efficiency in industrial drying: A hybrid ultrasonic system with a novel dynamic optimization framework

Drying processes are among the most energy-consuming operations in industrial and manufacturing settings, demanding strategic selection, design, and control for enhanced efficiency. Advancing drying technologies is critical for improving sustainability, lowering energy use, reducing carbon emissions, and minimizing waste. This study explores two innovative strategies aimed at transforming drying processes into sustainable, low-carbon systems by reducing energy consumption, minimizing waste, and maintaining a strong emphasis on preserving product quality. The first strategy showcases a sub-pilot scale hybrid ultrasonic-convective dryer for agrifood products. This technology, powered by electricity (process electrification), integrates non-thermal ultrasonic dehydration with convective heating and is presented as a sustainable and energy-efficient solution that enhances eco-friendly practices. The second strategy involves introducing and implementing a novel, multiobjective, mixed integer dynamic optimization technique to determine the optimal time-dependent process parameter values for the drying operation. This optimization technique yields operating conditions that are piecewise constant in time aiming to maximize the energy efficiency of the hybrid ultrasonic-convective dryer while ensuring strict adherence to product quality constraints. By adopting the hybrid ultrasonic-convective dryer, a notable 35% improvement in energy efficiency was achieved compared to conventional hot-air drying systems for drying apple slices. The proposed optimization framework further enhanced energy efficiency by nearly 14% over the most efficient process on the identical testbed, under static operating conditions. The reported enhancements have been experimentally validated. Regarding drying time (thereby improving production yield), the developed hybrid ultrasonic-convective dryer demonstrates as much as a 41% reduction in total processing time, which is further optimized by an additional 10% using our proposed optimization framework. The research outcomes have profound implications for the design and operation of drying systems, encompassing crucial aspects such as process electrification, cost-effectiveness, energy savings, time efficiency, product yield, product quality, and process automation.

Dynamic optimization↗

A nominally second-order cell-centered Lagrangian scheme for simulating elastic–plastic flows on two-dimensional unstructured grids

In this paper, we describe a cell-centered Lagrangian scheme devoted to the numerical simulation of solid dynamics on two-dimensional unstructured grids in planar geometry. This numerical method, utilizes the classical elastic-perfectly plastic material model initially proposed by Wilkins [M.L. Wilkins, Calculation of elastic–plastic flow, Meth. Comput. Phys. (1964)]. In this model, the Cauchy stress tensor is decomposed into the sum of its deviatoric part and the thermodynamic pressure which is defined by means of an equation of state. Regarding the deviatoric stress, its time evolution is governed by a classical constitutive law for isotropic material. The plasticity model employs the von Mises yield criterion and is implemented by means of the radial return algorithm. The numerical scheme relies on a finite volume cell-centered method wherein numerical fluxes are expressed in terms of sub-cell force. The generic form of the sub-cell force is obtained by requiring the scheme to satisfy a semi-discrete dissipation inequality. Sub-cell force and nodal velocity to move the grid are computed consistently with cell volume variation by means of a node-centered solver, which results from total energy conservation. The nominally second-order extension is achieved by developing a two-dimensional extension in the Lagrangian framework of the Generalized Riemann Problem methodology, introduced by Ben-Artzi and Falcovitz [M. Ben-Artzi, J. Falcovitz, Generalized Riemann Problems in Computational Fluid Dynamics, Cambridge Monogr. Appl. Comput. Math. (2003)]. Finally, the robustness and the accuracy of the numerical scheme are assessed through the computation of several test cases.

97 MATHEMATICS AND COMPUTING↗

Theory of magnetic enhancement in strontium hexaferrite through Zn-Sn pair substitution

We study the site occupancy and magnetic properties of Zn-Sn substituted M-type Sr-hexaferrite SrFe12-x(Zn0.5Sn0.5)(x)O-19 with x=1 using first-principles total-energy calculations. We find that in the lowest-energy configuration Zn2+ and Sn4+ ions preferentially occupy the 4f(1) and 4f(2) sites, respectively, in contrast to the model previously suggested by Ghasemi et al. [J. Appl. Phys, 107, 09A734 (2010)], where Zn2+ and Sn4+ ions occupy the 2b and 4f(2) sites. Density-functional theory calculations show that our model has a lower total energy by more than 0.2 eV per unit cell compared to Ghasemi's model. More importantly, the latter does not show an increase in saturation magnetization (M-s) compared to the pure M-type Sr-hexaferrite, in disagreement with the experiment. On the other hand, our model correctly predicts a rapid increase in M-s as well as a decrease in magnetic anisotropy compared to the pure M-type Sr-hexaferrite, consistent with experimental measurements. (c) 2013 Elsevier B.V. All rights reserved.

Liyanage, LSI↗

A three dimensional field formulation, and isogeometric solutions to point and line defects using Toupin's theory of gradient elasticity at finite strains

In this study, we present a field formulation for defects that draws from the classical representation of the cores as force dipoles. We write these dipoles as singular distributions. Exploiting the key insight that the variational setting is the only appropriate one for the theory of distributions, we arrive at universally applicable weak forms for defects in nonlinear elasticity. Remarkably, the standard, Galerkin finite element method yields numerical solutions for the elastic fields of defects that, when parameterized suitably, match very well with classical, linearized elasticity solutions. The true potential of our approach, however, lies in its easy extension to generate solutions to elastic fields of defects in the regime of nonlinear elasticity, and even more notably for Toupin's theory of gradient elasticity at finite strains (Toupin Arch. Ration. Mech. Anal., 11 (1962) 385). In computing these solutions we adopt recent numerical work on an isogeometric analytic framework that enabled the first three-dimensional solutions to general boundary value problems of Toupin's theory (Rudraraju et al. Comput. Methods Appl. Mech. Eng., 278 (2014) 705). We first present exhaustive solutions to point defects, edge and screw dislocations, and a study on the energetics of interacting dislocations. Then, to demonstrate the generality and potential of our treatment, we apply it to other complex dislocation configurations, including loops and low-angle grain boundaries.

36 MATERIALS SCIENCE↗

Third-order elastic constants of diamond determined from experimental data

The pressure derivatives of the second-order elastic constants (SOECs) of diamond were determined by analyzing previous sound velocity measurements under hydrostatic stress [McSkimin and Andreatch, J. Appl. Phys. 43, 294 (1972)]. Furthermore, our analysis corrects an error in the previously reported results.We present a complete and corrected set of third-order elastic constants (TOECs) using the corrected pressure derivatives, together with published data for the nonlinear elastic response of shock compressed diamond [Lang and Gupta, Phys. Rev. Lett. 106, 125502 (2011)] and it differs significantly from TOECs published previously.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward reducing uncertainty quantification costs in DEM models of particulate flow: Testing simple, sensitivity-based, forward uncertainty propagation techniques

The performance of two conceptually-simple uncertainty quantification techniques are tested against the rigorous nested-loop sampling technique of Roy and Oberkampf (Comput Methods Appl Mech Eng, 200: 2131–2144, 2011) (herein called full-sampling) using two very small-scale DEM-based models of particulate flow (one gas-solid flow and one granular flow). The first simplified forward uncertainty propagation technique, reduced-sampling, uses a sensitivity analysis to eliminate uncertain inputs that have little impact on the model output prior to nested-loop sampling. The second technique, boundary-sampling, uses a sensitivity analysis to inform the selection of two bounding cases for each key model output. In conclusion, the uncertainties in the model outputs obtained via the reduced- and boundary-sampling methods agree well with those from full-sampling for both the gas-solid and granular flow models while yielding computational savings of 65–75% (reduced sampling) and 94–97% (boundary sampling).

42 ENGINEERING↗

Metabolism of bilirubin by human cytochrome P450 2A6

The mouse cytochrome P450 (CYP) 2A5 has recently been shown to function as hepatic “Bilirubin Oxidase” (Abu-Bakar, A., et al., 2011. Toxicol. Appl. Pharmacol. 257, 14–22). To date, no information is available on human CYP isoforms involvement in bilirubin metabolism. In this paper we provide novel evidence for human CYP2A6 metabolising the tetrapyrrole bilirubin. Incubation of bilirubin with recombinant yeast microsomes expressing the CYP2A6 showed that bilirubin inhibited CYP2A6-dependent coumarin 7-hydroxylase activity to almost 100% with an estimated K{sub i} of 2.23 μM. Metabolite screening by a high-performance liquid chromatography/electrospray ionisation mass spectrometry indicated that CYP2A6 oxidised bilirubin to biliverdin and to three other smaller products with m/z values of 301, 315 and 333. Molecular docking analyses indicated that bilirubin and its positively charged intermediate interacted with key amino acid residues at the enzyme's active site. They were stabilised at the site in a conformation favouring biliverdin formation. By contrast, the end product, biliverdin was less fitting to the active site with the critical central methylene bridge distanced from the CYP2A6 haem iron facilitating its release. Furthermore, bilirubin treatment of HepG2 cells increased the CYP2A6 protein and activity levels with no effect on the corresponding mRNA. Co-treatment with cycloheximide (CHX), a protein synthesis inhibitor, resulted in increased half-life of the CYP2A6 compared to cells treated only with CHX. Collectively, the observations indicate that the CYP2A6 may function as human “Bilirubin Oxidase” where bilirubin is potentially a substrate and a regulator of the enzyme. -- Highlights: ► Human CYP2A6 interacts with bilirubin with a high affinity. ► Bilirubin docking to the CYP2A6 active site is more stable than biliverdin docking. ► Recombinant CYP2A6 microsomes metabolised bilirubin to biliverdin. ► Bilirubin increased the hepatic CYP2A6 protein and activity levels but not mRNA. ► Co-treatment with a protein synthesis inhibitor prolongs CYP2A6 half-life.

60 APPLIED LIFE SCIENCES↗

Part 2: Forensic attribution profiling of Russian VX in food using liquid chromatography-mass spectrometry

This work is part two of a three-part series in this issue of a Sweden-United States collaborative effort towards the understanding of the chemical attribution signatures of Russian VX (VR) in synthesized samples and complex food matrices. We describe the sourcing of VR present in food based on chemical analysis of attribution signatures by liquid chromatography-tandem mass spectrometry (LC-MS/MS) combined with multivariate data analysis. Analytical data was acquired from seven different foods spiked with VR batches that were synthesized via six different routes in two separate laboratories. The synthesis products were spiked at a lethal dose into seven food matrices: water, orange juice, apple purée, baby food, pea purée, liquid eggs and hot dog. After acetonitrile sample extraction, the samples were analyzed by LC-MS/MS operated in MRM mode. A multivariate statistical calibration model was built on the chemical attribution profiles from 118 VR spiked food samples. Using the model, an external test-set of the six synthesis routes employed for VR production was correctly identified with no observable major impact of the food matrices to the classification. The overall performance of the statistical models was found to be exceptional (94%) for the test set samples retrospectively classified to their synthesis routes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An efficient method for unfolding kinetic pressure driven VISAR data

Velocity Interferometer System for Any Reflector (VISAR) [Barker and Hollenbach, J. Appl. Phys. 43, 4669 (1972)] is a well-known diagnostic that is employed on many shock physics and pulsed-power experiments. With the VISAR diagnostic, the velocity on the surface of any metal flyer can be found. For most experiments employing VISAR, either a kinetic pressure [Grady, Mech. Mater. 29, 181 (1998)] or a magnetic pressure [Lemke et al., Intl J. Impact Eng. 38, 480 (2011)] drives the motion of the flyer. Moreover, reliable prediction of the time-dependent pressure is often a critical component to understanding the physics of these experiments. Although VISAR can provide a precise measurement of a flyer’s surface velocity, the real challenge of this diagnostic implementation is using this velocity to unfold the time-dependent pressure. The purpose of this study is to elucidate a new method for quickly and reliably unfolding VISAR data.

42 ENGINEERING↗

Recent laser upgrades at Sandia’s Z-backlighter facility in order to accommodate new requirements for magnetized liner inertial fusion on the Z-machine

The Z-backlighter laser facility primarily consists of two high energy, high-power laser systems. Z-Beamlet laser (ZBL) (Rambo et al., Appl. Opt. 44, 2421 (2005)) is a multi-kJ-class, nanosecond laser operating at 1054 nm which is frequency doubled to 527 nm in order to provide x-ray backlighting of high energy density events on the Z-machine. Z-Petawatt (ZPW) (Schwarz et al., J. Phys.: Conf. Ser. 112, 032020 (2008)) is a petawatt-class system operating at 1054 nm delivering up to 500 J in 500 fs for backlighting and various short-pulse laser experiments (see also Figure 10 for a facility overview). With the development of the magnetized liner inertial fusion (MagLIF) concept on the Z-machine, the primary backlighting missions of ZBL and ZPW have been adjusted accordingly. As a result, we have focused our recent efforts on increasing the output energy of ZBL from 2 to 4 kJ at 527 nm by modifying the fiber front end to now include extra bandwidth (for stimulated Brillouin scattering suppression). The MagLIF concept requires a well-defined/behaved beam for interaction with the pressurized fuel. Hence we have made great efforts to implement an adaptive optics system on ZBL and have explored the use of phase plates. We are also exploring concepts to use ZPW as a backlighter for ZBL driven MagLIF experiments. Alternatively, ZPW could be used as an additional fusion fuel pre-heater or as a temporally flexible high energy pre-pulse. All of these concepts require the ability to operate the ZPW in a nanosecond long-pulse mode, in which the beam can co-propagate with ZBL. Finally, some of the proposed modifications are complete and most of them are well on their way.

47 OTHER INSTRUMENTATION↗

Differential formulation of the viscous history force on a particle for efficient and accurate computation

It is well known that the computation of the Basset-like history force is very demanding in terms of CPU and memory requirements, since it requires the evaluation of a history integral. We use the recent rational theory of Beylkin & Monzón ( Appl. Comput. Harmon. Anal. , vol. 19, 2005, pp. 17–48) to approximate the history kernel in the form of exponential sums to reformulate the viscous history force in a differential form. This theory allows us to approximate the history kernel in terms of exponential sums to any desired order of accuracy. This removes the need for long-time storage of the acceleration histories of the particle and the fluid. The proposed differential form approximation is applied to compute the history force on a spherical particle in a synthetic turbulent flow and a wall-bounded turbulent channel flow. Particles of various diameters are considered, and results obtained using the present technique are in reasonable agreement with those achieved using the full history integral.

Mechanics↗