VFUP

VFUP, Action, Opt, OPT2, OPT3, OPT4
Controls view factor updating at the substep level for a coupled-field radiation analysis that includes large-deflection effects.

Valid Products: Pro | Premium | Enterprise | PrepPost | Solver | AS add-on

Action

Action to be performed:

DEFINE

 — 

Define view factor updating parameters (default).

STAT

 — 

Print all specified values for VFUP command arguments.

DEFA

 — 

Resets all Opt options to their default values. OPT2, OPT3, and OPT4 are ignored.

Opt

Option to enable view factor updates at the substep level or specify updating criteria.

ON

 — 

Enable view factor updates according to the last criteria specified by the VFUP command. The default criteria is VFUP,,FRQU,1,0 if you only issue VFUP,,ON without issuing another VFUP command to choose between frequency-based (Opt = FRQU) or kinematic-based criteria (Opt = KNEX, KNIM, or KNI1). OPT2, OPT3, and OPT4 are ignored.

OFF

 — 

Do not update view factors at the substep level (default). OPT2, OPT3, and OPT4 are ignored.

FRQU

 — 

Key to specify a desired value for the substep frequency (OPT2) and the pass frequency (OPT3) at which view factors are updated. For guidelines on setting these values and example command usage, see Controlling the frequency of view factor updates. OPT4 is ignored.

KNEX

 — 

Key to specify kinematic-based view factor updating at the end of a substep. The program tests all radiation facets and updates view factors only if the relative change in facet spacing, area, or angle exceeds a threshold specified by OPT2, OPT3, and OPT4. For details on how these values are calculated, see the table and figures below and Kinematic-based algorithm in the notes. For guidelines on the timing of view factor updates, see Controlling the frequency of view factor updates.

KNIM

 — 

Key to specify kinematic-based view factor updating within the radiation multipass loop, after the Newton-Raphson (RP) loop has produced a converged coupled-field solution. The program tests all radiation facets and updates view factors only if the relative change in facet spacing, area, or angle exceeds a threshold specified by OPT2, OPT3, and OPT4. For details on how these values are calculated, see the table and figures below and Kinematic-based algorithm in the notes. For guidelines on the timing of view factor updates, see Controlling the frequency of view factor updates.

KNI1

 — 

Key to specify kinematic-based view factor updating within the first radiation multipass loop, after the Newton-Raphson (RP) loop has produced a converged coupled-field solution. The program tests all radiation facets and updates view factors only if the relative change in facet spacing, area, or angle exceeds a threshold specified by OPT2, OPT3, and OPT4. For details on how these values are calculated, see the following table and figures and Kinematic-based algorithm in the notes. For guidelines on the timing of view factor updates, see Controlling the frequency of view factor updates.

Table 231: Definitions of OPT2, OPT3, and OPT4 for different Opt labels

Opt labelOPT2OPT3OPT4
FRQU

(frequency-based updating)

Substep frequency for view factor updates. This value must be a positive integer (default = 1). If OPT2 = 3, view factors are updated every third substep.Valid values are -1, 0, or a positive integer (default = 0).
If OPT3 = -1, view factors are updated only in the first pass of the multipass loop in the substeps that view factors are being updated (specified by OPT2).
If OPT3 = 0, view factors are not updated in the radiation multipass loop.
If OPT3 is a positive integer, it specifies the pass frequency for view factor updates in the multipass loop.
Ignored
KNEX, KNIM, or KNI1

(kinematic-based updating)

Relative change in facet spacing threshold (default = 0.05).

Set OPT2 = 0 to ignore this threshold test.

Relative change in facet area threshold (default = 0.05).

Set OPT3 = 0 to ignore this threshold test.

Relative change in facet angle threshold (default = 0.05).

Set OPT4 = 0 to ignore this threshold test.

2-surface test:

view factors update if

OPT2

2-surface test:

view factors update if

OPT3

2-surface test:

view factors update if

OPT4

1-surface test:

view factors update if

OPT2

1-surface test:

view factors update if

OPT3

1-surface test:

view factors update if

OPT4


The following figure depicts the values used in kinematic-based tests to determine if thresholds are surpassed and view factors must be updated. There are two tests: a 2-surface test and a 1-surface test. The formulas for each of these tests are listed in the previous table, and the variables in these formulas are explained in the following figures. For a detailed description of the algorithm used to test all radiation surfaces and determine if thresholds are exceeded, see Kinematic-based algorithm in the notes.

Figure 31: Definition of Values Used in 2-surface (A) and 1-surface (B) Kinematic-based Tests

Definition of Values Used in 2-surface (A) and 1-surface (B) Kinematic-based Tests


RefLength in the denominator of the 1-surface facet spacing threshold test is a length scale representative of the radiation facet spacing for an enclosure. It is automatically calculated by putting a bounding box on the radiation facets and computing the diagonal of the bounding box. RefLength is dynamic and updates when view factors are updated. See the following figure for a depiction of how RefLength is computed.

Figure 32: Example Demonstrating How RefLength is Calculated

Example Demonstrating How RefLength is Calculated


For details on how the kinematic-based thresholds are tested to determine if view factors are updated, see Kinematic-based algorithm in the notes.

Command Default

Although the command is not required for solution, if you do not issue VFUP, view factors are not updated at the substep level. Consequently, results of thermal radiation modeling may be inaccurate because they do not account for changes in the view factors of radiation surfaces as they deform or move.

Notes

Deformation and/or motion of radiating surfaces during an anlysis can cause inaccuracies in radiation calculations if the view factors are not updated to account for these changes. Issue VFUP,DEFINE,ON to enable view factor updates at the substep level and improve simulation accuracy for an analysis that includes all of the following conditions:

Updates occur only if all of these conditions are met:

  • VFUP,DEFINE,ON has been issued.

  • NLGEOM,ON has been issued.

  • You have specified KEYOPT(1) = 1 for any surface elements (SURF251 or SURF252) that have been generated via the RSURF command.

When view factor updating is enabled (VFUP,,ON), radiation facet areas are recalculated along with the view factors according to the specified criteria. To specify the updating criteria, you must re-issue the VFUP command and choose one of the following:

  • frequency-based criteria (Opt = FRQU) to have updates occur at defined substep (OPT2) and pass (OPT3) frequencies, regardless of geometric changes to radiation facets,

    or

  • kinematic-based criteria (Opt = KNEX, KNIM, or KNI1) to have updates occur only if relative changes in facet spacing, areas, or angles exceed thresholds defined by OPT2, OPT3, and OPT4, respectively.

Controlling the frequency of view factor updates

The full solution method (THOPT,FULL) with radiosity has three nested loops, which occur for each load step: the substep (or time step) loop, the multipass loop, and the Newton-Raphson (NR) loop. Displacements and temperatures, in addition to any other DOFs that may be included in your coupled-field analysis, are computed in the NR loop. The view factors / radiation areas can be recalculated at different points during the solution, as depicted in Figure 33: View Factor Update Timing. There are three timing strategies, Explicit, Implicit, and Implicit-1, and the implications of each timing strategy are detailed in the following discussion.

Figure 33: View Factor Update Timing

View Factor Update Timing


As seen in the figure, implicit updates occur within the multipass loop, and updated view factors are used in the current substep. Instead, explicit updates occur at the end of the time step, and updated view factors are used in the next time step. An additional update occurs before entering the substep loop for a restart (ANTYPE,, RESTART) or if you have imposed initial conditions on the displacements (IC) as indicated.

There is a trade-off between accuracy and computational cost as updating the view factors at the substep level improves accuracy but also increases the solution time. To control the timing of view factor updates and whether updates are based on geometric changes to the radiation facets due to kinematic motion or deformation during solution, you must choose one of these options:

  • Opt = FRQU: Specify the substep (OPT2) and pass (OPT3) frequency at which view factor updates occur. Updates happen at the specified frequency, regardless of whether deformation or motion causes significant changes to the radiation facets.

  • Opt = KNEX, KNIM, or KNI1: The view factors are updated only if changes in the radiation facets due to motion/deformation exceed a specified threshold (OPT2, OPT3, or OPT4). As shown in Figure 33: View Factor Update Timing and described in detail below, the difference between these kinematic-based update options is when the testing and updates occur.

There are three timing strategies:

  • Explicit - View factors are updated at the end of a substep after a converged coupled-field solution is achieved, and the updated view factors are used in the next time step.

    For Opt = FRQU  —  You chose the explicit strategy by setting OPT3 = 0, and the value specified for the substep frequency, OPT2, determines which substeps include view factor updating. By default, OPT2 = 1, and view factors are updated at the end of every substep.

    For example, if you issue:View factors are updated
    VFUP,DEFINE,FRQU,3,0at the end of substeps 3,6,9,...

    For Opt = KNEX  —  At the end of every substep, the facets are tested (see Kinematic-based algorithm for details) and view factors updated only if one of the three kinematic-based thresholds (OPT2, OPT3, OPT4) are exceeded.

    For example, if you issue:Radiation facets are tested and view factors are updated
    VFUP,DEFINE,KNEXat the end of every substep. Updates occur only if one of three kinematic-based thresholds (facet spacing, area, or angle) exceeds the default value of 0.05
    VFUP,DEFINE,KNEX,0.1,,0at the end of every substep. Updates occur only if the relative change in facet spacing threshold exceeds 0.1 or the relative change in facet area exceeds the default value of 0.05.

    Explicit is the most stable timing strategy. If you have convergence difficulties, explicit updates are recommended. For an accurate solution, the substep size must be small. Explicit updates may lead to an energy imbalance in the radiosity solution since the view factor facet areas were calculated in the preceding time step. For further details, see the discussion related to Figure 37: View Factor Updates and Reciprocity.

  • Implicit - View factors are updated within the radiation multipass loop, after the Newton-Raphson (RP) loop has produced a converged coupled-field solution, and are used in the current pass and substep.

    For Opt = FRQU  —  The value of the pass frequency, OPT3, determines which passes will update the view factors. To have implicit view factor updating, you must set OPT3 to something other than its default value of zero.

    For example, if you issue:View factors are updated
    VFUP,DEFINE,FRQU,5,3in passes 3,6,9,... for substeps 5,10,15...

    For Opt = KNIM  —  The testing and view factor updating if a kinematic-based threshold is exceeded (see Kinematic-based algorithm for details) occurs within the radiation multipass loop, after the Newton-Raphson (RP) loop has produced a converged coupled-field solution. Updated view factors are used in the current pass and substep.

    For example, if you issue:Radiation facets are tested and view factors are updated
    VFUP,DEFINE,KNIMin every pass for every substep. Updates occur only if one of three kinematic-based thresholds (facet spacing, area, or angle) exceeds the default value of 0.05
    VFUP,DEFINE,KNIM,0.1,0,0.1in every pass for every substep. Updates occur only if the relative change in the facet spacing or angle threshold exceeds 0.1.

    Implicit updates are useful when you have strong nonlinearities, strong structural-thermal coupling, or rapid changes in surface geometry. However, implicit updates can be computationally expensive as view factors are recalculated multiple times within a single substep. Also, Implicit updates of view factors can lead to instabilities and convergence difficulties.

  • Implicit-1 - View factor updating occurs onlly in the first pass of the radiation multipass loop.

    For Opt = FRQU  —  View factors are updated in the substeps prescribed by the value of the substep frequency, OPT2, but they are updated only in the first pass of the radiation multipass loop. To have implicit-1 view factor updating, specify -1 for OPT3.

    For example, if you issue:View factors are updated
    VFUP, DEFINE, FRQU, 5, -1only in pass 1 for substeps 5,10,15...

    For Opt = KNI1  —  The facet testing and view-factor updating if a kinematic-based threshold is exceeded (see Kinematic-based algorithm for details) occurs only within the first pass of the radiation multipass loop, after the Newton-Raphson (RP) loop has produced a converged coupled-field solution.

    For example, if you issue:Radiation facets are tested and view factors are updated
    VFUP,DEFINE,KNI1in the first pass for every substep. Updates occur only if one of three kinematic-based thresholds (facet spacing, area, or angle) exceeds the default value of 0.05
    VFUP,DEFINE,KNI1,0.1,0,0in the first pass for every substep. Updates occur only if the relative change in the facet spacing threshold exceeds 0.1.

    Use this timing strategy when the motion of bodies in your model is completely specified, as in rigid body motion (translation and rotation). It is most efficient and accurate to update the view factors in the first pass only when the NR loop achieves a converged structural solution in the first pass and the structural solution does not change for subsequent passes.

Kinematic-based algorithm

For Opt = KNEX, KNIM, or KNI1, the kinematic-based criteria are used in the algorithm depicted below to test and update view factors.

Figure 34: Kinematic-based Algorithm

Kinematic-based Algorithm


In most practical situations, S2 will be found and the 2-surface test will be done. In some situations, when the facet only radiates to the environment, the 1-surface test is required. Two examples of when a 1-surface test is encountered are shown in the following figure. For a 2D square with four radiating surfaces around it, the 4x4 view factor matrix is zero, and only single surface tests will be done as the body translates (A in the figure). The second example (B) is two radiating surfaces, one fixed and one rotating. When the rotating surface reaches an 180 degree rotation, the view factor matrix becomes zero, and the 1-surface test will be encountered.

Figure 35: Examples Scenarios that Require a 1-Surface Test

Examples Scenarios that Require a 1-Surface Test


Notes on kinematic-based and frequency-based updates

Kinematic-based updates (Opt = KNEX, KNIM, or KNI1) are sensitive to motion while frequency-based updates (Opt = FRQU) are not. For example, consider a 2 plate problem with total loadstep time = 100. Plate 1 is fixed and plate to moves only from time 10-20. Frequency-based updates will occur from time = 0 - 100, but kinematic-based updates will occur from time 10-20.

If an enclosure is undergoing rigid-body translation or rotation, and there is no change to the radiation facets:

  • If you are using kinematic-based updates, view factors will not be updated if the 2-surface test is encountered. However, if a 1-surface surface is encountered, they will be updated.

  • If you are using frequency-based updates, view factors are updated at the prescribed frequency regardless of whether or not the radiation facets are moving.

When you have no knowledge about body deformation prior to the simulation, kinematic-based updates are a good option since view factor updates will only be triggered if the geometry of the radiation surfaces move or deform significantly.

For kinematic-based updating, when bodies are moving away from each other, updates occur less often. Vice-versa, when bodies move toward each other, updates occur more frequently.

Restrictions / limitations

Enclosure status change and penetration must be avoided  —  Simulations using the VFUP to update view factors at the substep level are invalid if an enclosure status changes or the surface of a radiating body penetrates the surface of another radiating body. Figure 36: Examples of Radiation Problems that Undergo an Enclosure Status Change illustrates some examples of enclosure status changes that must be avoided. For example, if two radiating bodies in close proximity penetrate each other due to thermal expansion during the analysis, their view factors become zero, and results using the VFUP command to update view factors are unpredictable and inaccurate. It is better to model bodies within close proximity using thermal contact radiation (contact and target elements: CONTA172, CONTA174, CONTA175, CONTA178,TARGE169, and TARGE170).

You must also be careful to avoid an enclosure status change, for example if movement or deformation in the analysis causes an open enclosure to become closed or vice versa. If an enclosure is undergoing a status change in your analysis, split the load step into multiple load steps such that:

  • There is no status change in the individual load steps.

  • VFSM settings are appropriate for the enclosure status in individual load steps.

Figure 36: Examples of Radiation Problems that Undergo an Enclosure Status Change

Examples of Radiation Problems that Undergo an Enclosure Status Change

View factor Updates and Reciprocity  —  For explicit updates, view factors are recalculated at the end of the substep and applied in the next substep. Consequently, they violate the reciprocity rule in the next substep as the radiating areas change since the view factors were calculated from the areas of the previous substep. This may lead to energy imbalance issues. The same problem exists if you specify a substep (OPT2) or pass (OPT3) frequency other than one. Figure 37: View Factor Updates and Reciprocity illustrates this issue for different updating strategies. Reciprocity / row sum properties are only satisfied in substeps where view factors are updated in every pass (OPT3 = 1) provided you have used the VFSM command appropriately (see View Factors in the Theory Reference).

Figure 37: View Factor Updates and Reciprocity

View Factor Updates and Reciprocity

Supported/restricted commands at the substep level with VFUP,DEFINE,ON

The following commands are not supported at the substep level:

The following commands may be used to query and print view factors or average view factors at the substep level:

  • View factor queries (VFQUERY) are based on the last update.

  • Any *GET command with Entity = RAD is based on the last view factor update. To retrieve the latest average view factor, issue *GET,Par,RAD,0,VFAVG.

  • To print view factors to the output file at the substep level, issue MPRINT before SOLVE.

The following radiation commands are valid at the substep level, and may be used with view factor updating enabled (VFUP,DEFINE,ON):

  • HEMIOPT, RADOPT, V2DOPT, SPCTEMP, SPCNOD, and VFCO

  • Settings on the VFSM command to adjust view factors (view factor smoothing) for different enclosures are valid for all substeps. However, if there is an enclosure status change during a load step, VFSM settings may become inappropriate and results inaccurate. See earlier discussion on inaccurate results if motion or deformation causes an enclosure status change (Figure 36: Examples of Radiation Problems that Undergo an Enclosure Status Change). To model a situation where deformation changes the enclosure status, you must solve the problem in multiple load steps and change VFSM settings in each load step to accurately model the state of the enclosure.

Menu Paths

This command cannot be accessed from a menu.